Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

Bio-Chemoinformatics-Driven Analysis of nsp7 and nsp8 Mutations and Their Effects on Viral Replication Protein Complex Stability

The nonstructural proteins 7 and 8 (nsp7 and nsp8) of SARS-CoV-2 are highly important proteins involved in the RNA-dependent polymerase (RdRp) protein replication complex. In this study, we analyzed the global mutation of nsp7 and nsp8 in 2022 and 2023 and analyzed the effects of mutation on the viral replication protein complex using bio-chemoinformatics. Frequently occurring variants are found to be single amino acid mutations for both nsp7 and nsp8. The most frequently occurring mutations for nsp7 which include L56F, L71F, S25L, M3I, D77N, V33I and T83I are predicted to cause destabilizing effects, whereas those in nsp8 are predicted to cause stabilizing effects, with the threonine to isoleucine mutation (T89I, T145I, T123I, T148I, T187I) being a frequent mutation. A conserved domain database analysis generated critical interaction residues for nsp7 (Lys-7, His-36 and Asn-37) and nsp8 (Lys-58, Pro-183 and Arg-190), which, according to thermodynamic calculations, are prone to destabilization. Trp-29, Phe-49 of nsp7 and Trp-154, Tyr-135 and Phe-15 of nsp8 cause greater destabilizing effects to the protein complex based on a computational alanine scan suggesting them as possible new target sites. This study provides an intensive analysis of the mutations of nsp7 and nsp8 and their possible implications for viral complex stability.

Distinct chikungunya virus polymerase palm subdomains contribute to virus replication and virion assembly

Alphaviruses encode an error-prone RNA-dependent RNA polymerase (RdRp), nsP4, required for genome synthesis, yet how the RdRp functions in the complete alphavirus life cycle is not well-defined. Previous work using chikungunya virus (CHIKV) has established the importance of the nsP4 residue cysteine 483 in maintaining viral genetic fidelity. Given the location of residue C483 in the nsP4 palm domain, we hypothesized that other residues within this domain and surrounding subdomains would also contribute to polymerase function. To test this hypothesis, we designed a panel of nsP4 variants via homology modeling based on the Coxsackievirus B3 3 polymerase. We rescued each variant in both mammalian and mosquito cells and discovered that the palm domain and ring finger subdomain contribute to polymerase host-specific replication and genetic stability. Surprisingly, in mosquito cells, these variants in the ring finger and palm domain were replication competent and produced viral structural proteins, but they were unable to produce infectious progeny, indicating a yet uncharacterized role for the polymerase in viral assembly. Finally, we have identified additional residues in the nsP4 palm domain that influence the genetic diversity of the viral progeny, potentially via an alteration in NTP binding and/or discrimination by the polymerase. Taken together, these studies highlight that distinct nsP4 subdomains regulate multiple processes of the alphavirus life cycle, placing nsP4 in a central role during the switch from RNA synthesis to packaging and assembly.

Tyrosinase inhibitory peptides: Structure-activity relationship study on peptide chemical properties, terminal preferences and intracellular regulation of melanogenesis signaling pathways

BACKGROUND

Bioactive peptides have gained attention as potential alternatives to chemical-based skin lightening agents. Based on literature search, the reported articles focused mainly on the sources and preparation methods of tyrosinase inhibitory peptides and there is lacking information regarding the structure-activity relationship (SAR) between peptide property and tyrosinase inhibition. It was hypothesized that peptide properties such as hydrophobicity/hydrophilicity and the amino acid type and position/arrangement at the terminal positions could affect peptide mode of binding hence result in various degrees of tyrosinase inhibition.

METHODS

In this study, the sequences of 128 tyrosinase inhibitory peptides were collected from peer-reviewed articles. The hydrophobicity/hydrophilicity property and the amino acid profile of peptides at the N- and C-terminals were analyzed using bioinformatics tools. Molecular docking analysis was employed to further elucidate the roles of reactive amino acids in tyrosinase-peptide binding interactions. The peptide-regulated intracellular melanogenesis pathways were also compiled and discussed.

RESULTS

It was found that hydrophobic and/or polar neutral properties may facilitate or stabilize peptide binding with tyrosinase. Moreover, short peptides featuring a cysteine and tyrosine at the N- and C- terminal ultimate positions tend to bind to the active site of tyrosinase whereas positively charged amino acid such as arginine at the N-terminal does not favor peptide binding to tyrosinase.

CONCLUSIONS

These findings provide detailed explanation on how peptide/amino acid structures are related and what function they play in tyrosinase inhibition. It could also inspire researchers to account for tyrosinase-peptide SAR and the underlying anti-melanogenesis mechanisms in formulating peptide-based treatments or strategies against skin hyperpigmentation.

Nitrogen regulates the synthesis of hydrophobic amino acids to improve protein structural and gel properties in common buckwheat

Nitrogen (N) fertilizer impacts the grain quality of common buckwheat, but the effects and regulatory mechanisms of N on various protein parameters of buckwheat are not fully understood. The purpose of this study was to investigate the particle morphology, structural and gel properties, and regulation mechanism of buckwheat protein under four N levels. The bulk density, surface hydrophobicity, particle size, and thermal properties of the buckwheat protein were maximized through the optimal N application (180 kg N/ha), further enhancing the thermal stability of the protein. N application increased the β-sheet content and reduced the random coil content. Appropriate N fertilizer input enhanced the tertiary structure stability and gel elasticity of buckwheat protein by promoting hydrophobic interactions, disulfide bonds, ionic bonds, storage modulus and loss modulus. The differentially expressed proteins induced by N are primarily enriched in small ribosomal subunit and ribosome, improving protein quality mainly by promoting the synthesis of hydrophobic amino acids. Future agriculture should pay attention to the hydrophobic amino acid content of buckwheat to effectively improve protein quality. This study further advances the application of buckwheat protein in the field of food processing and provides a theoretical basis for the extensive development and utilization of buckwheat protein.

Design of Plant-Based Food: Influences of Macronutrients and Amino Acid Composition on the Techno-Functional Properties of Legume Proteins

Imitating animal-based products using vegetable proteins is a technological challenge that can be mastered based on their techno-functional properties. These properties of legume proteins can be influenced by multiple factors, among which the macronutrients and amino acid contents play an important role. Therefore, the question arises as to what extent the techno-functional properties are related to these factors. The water- and oil-holding capacities and the emulsion and foaming properties of commercially available legume protein powders were analyzed. Correlations between macronutrient, amino acid content, steric structure, and techno-functional properties were conducted. However, the protein concentration is the focus of techno-functional properties, as well as the type of protein and the interaction with the non-protein ingredients. The type of protein is not always quantified by the quantity of amino acids or by their spatial arrangement. In this study, the effects of the three-dimensional structure were observed by the used purification method, which overshadow the influencing factors of the macronutrients and amino acid content. In summary, both the macronutrient and amino acid contents of legume proteins provide a rough indication but not a comprehensive statement about their techno-functional properties and classification in an adequate product context.

Exploring the Binding Mechanism of NRG1-ERBB3 Complex and Discovery of Potent Natural Products to Reduce Diabetes-Assisted Breast Cancer Progression

Diabetes mellitus significantly contributes to breast cancer progression, where hyperglycemia upregulates specific genes, leading to more aggressive tumor growth. In patients with BC that develop diabetes, neuregulin 1 (NRG1) and epidermal growth factor receptor 3 (ERBB3) overexpression exacerbate tumor growth and progression. Since the interaction between NRG1 and ERBB3 is critical for tumor growth, understanding the molecular mechanisms underlying NRG1-ERBB3 complex formation is essential for elucidating diabetes-assisted breast cancer progression. However, the key residues forming the NRG1-ERBB3 complex remain unknown. Here, we substituted specific residues in NRG1 with alanine and studied its interactions with ERBB3 using computational structural biology tools. We further screened the South African natural compounds database to target the complex's interface residues to discover potential inhibitors. The conformational stability and dynamic features of NRG1-WT, -H2A, -L3A, and -K35A complexed with ERBB3 were subjected to 400 ns molecular dynamics simulations. The free binding energies of all NRG1-ERBB3 complexes were calculated using the molecular mechanics-generalized Born surface area (MM/GBSA). The H2 and L3 alanine substitutions caused a loss of interaction with ERBB3 residue D73, weakening the interaction with ERBB3. Screening 1300 natural compounds identified four (SANC00643, SANC00824, SANC00975, and SANC00335) with the best potential to inhibit ERRB3-NRG1 coupling. The binding free energies for each complex were - 48.55 kcal/mol for SANC00643, - 47.68 kcal/mol for SANC00824, - 46.04 kcal/mol for SANC00975, and - 45.29 kcal/mol for SANC00335, showing their overall stronger binding with ERBB3 than NRG1 and their potential to act as ERBB3-NRG1 complex inhibitors. In conclusion, this complex may represent a residue-specific drug target to inhibit BC progression.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Piperidine-based natural products targeting Type IV pili antivirulence: A computational approach

Type IV (T4) pilus is among the virulence factors with a key role in serious bacterial diseases. Specifically, in Neisseria meningitidis and Pseudomonas aeruginosa, it determines pathogenicity and causes infection. Here, a computational approach has been pursued to find piperidine-based inhibitor molecules against the elongation ATPase of T4 pili in these two selected pathogens. Using the modeled structures of the PilF and PilB ATPases of N. meningitidis and P. aeruginosa, virtual library screening via molecular docking has returned inhibitor molecule candidates. The dynamics of the best three binders have further been investigated in detail via molecular dynamic simulations. Among these, ligands with COCONUT IDs CNP0030078 and CNP0051517 were found to have higher potential in the inhibition of ATPases based on molecular dynamic simulation analysis and biological activity information. The obtained results will guide future efforts in antivirulence drug development against T4 pili of N. meningitidis and P. aeruginosa.

Structure of the AlgKX modification and secretion complex required for alginate production and biofilm attachment in Pseudomonas aeruginosa

Synthase-dependent secretion systems are a conserved mechanism for producing exopolysaccharides in Gram-negative bacteria. Although widely studied, it is not well understood how these systems are organized to coordinate polymer biosynthesis, modification, and export across both membranes and the peptidoglycan. To investigate how synthase-dependent secretion systems produce polymer at a molecular level, we determined the crystal structure of the AlgK-AlgX (AlgKX) complex involved in Pseudomonas aeruginosa alginate exopolysaccharide acetylation and export. We demonstrate that AlgKX directly binds alginate oligosaccharides and that formation of the complex is vital for polymer production and biofilm attachment. Finally, we propose a structural model for the AlgEKX outer membrane modification and secretion complex. Together, our study provides insight into how alginate biosynthesis proteins coordinate production of a key exopolysaccharide involved in establishing persistent Pseudomonas lung infections.

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

Comparative Analysis of the Mitochondrial Genome of Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae) From a Depth of 9,462 m at the Mariana Trench

The hadal zone, which represents the deepest marine habitat on Earth (6,000–11,000 m), is a harsh environment mainly characterized by extremely high hydrostatic pressure, and this habitat is believed to have a high degree of endemism. The deep-sea anemone family Galatheanthemidae comprises two valid species exclusively from the hadal; however, no other information about this family is currently available. In the present study, a sea anemone was collected from a depth of 9,462 m at the Mariana Trench and was defined as Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae). The mitochondrial genome of Galatheanthemum sp. MT-2020 was circular, was 16,633 bp in length, and contained two ribosomal RNA genes, 13 protein-coding genes and two transfer RNA genes. The order of the genes of Galatheanthemum sp. MT-2020 was identical to that of the majority of the species of the order Actiniaria. The value of the AT-skew was the lowest in the whole mitochondrial genome, with a positive GC skew value for the atp8 gene, while other species, except Antholoba achates, had the negative values of the GC skew. Galatheanthemum sp. MT-2020 was clustered with another abyssal species, Paraphelliactis xishaensis, in the phylogenetic tree, and these species diverged in the early Jurassic approximately 200 Mya from the shallow-sea species. The usage ratio of valine, which is one of the five amino acids with the strongest barophilic properties, in the mitochondrial genomes of the two abyssal species was significantly higher than that in other species with habitats above the depth of 3,000 m. The ω (dN/dS) ratio of the genomes was 2.45-fold higher than that of the shallow-sea species, indicating a slower evolutionary rate. Overall, the present study is the first to provide a complete mitogenome of sea anemones from the hadal and reveal some characteristics that may be associated with adaptation to an extreme environment.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔG_u°'(H₂O), of UBA(1) is 2.4 kcal mol^-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, k_f, in the absence of a denaturant of 13,000 s^-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone ¹⁵NH, ¹³CO, and ¹³Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Heterologous expression and biological characteristics of UGPases from Lactobacillus acidophilus

Herein, two genes (LBA0625 and LBA1719) encoding UGPases (UDP-glucose pyrophosphorylase) in Lactobacillus acidophilus (L. acidophilus) were successfully transformed into Escherichia coli BL21 (DE3) to construct recombinant overexpressing strains (E-0625, E-1719) to investigate the biological characteristics of UGPase-0625 and UGPase-1719. The active sites, polysaccharide yield, and anti-freeze-drying stress of L. acidophilus ATCC4356 were also detected. UGPase-0625 and UGPase-1719 belong to the nucleotidyltransferase of stable hydrophilic proteins; contain 300 and 294 amino acids, respectively; and have 20 conserved active sites by prediction. Αlpha-helixes and random coils were the main secondary structures, which constituted the main skeleton of UGPases. The optimal mixture for the high catalytic activity of the two UGPases included 0.5 mM UDP-Glu (uridine diphosphate glucose) and Mg²⁺ at 37 °C, pH 10.0. By comparing the UGPase activities of the mutant strains with the original recombinant strains, A10, L130, and L263 were determined as the active sites of UGPase-0625 (P < 0.01) and A11, L130, and L263 were determined as the active sites of UGPase-1719 (P < 0.01). In addition, UGPase overexpression could increase the production of polysaccharides and the survival rates of recombinant bacteria after freeze-drying. This is the first study to determine the enzymatic properties, active sites, and structural simulation of UGPases from L. acidophilus, providing in-depth understanding of the biological characteristics of UGPases in lactic acid bacteria.Key points• We detected the biological characteristics of UGPases encoded by LBA0625 and LBA1719.• We identified UGPase-0625 and UGPase-1719 active sites.• UGPase overexpression elevates polysaccharide levels and post-freeze-drying survival.

Biochemical and structural studies of two tetrameric nucleoside 2'-deoxyribosyltransferases from psychrophilic and mesophilic bacteria: Insights into cold-adaptation

Nucleoside 2'-deoxyribosyltransferases (NDTs) catalyze the cleavage of glycosidic bonds of 2'-deoxynucleosides and the following transfer of the 2'-deoxyribose moiety to acceptor nucleobases. Here, we report the crystal structures and biochemical properties of the first tetrameric NDTs: the type I NDT from the mesophilic bacterium Enterococcus faecalis V583 (EfPDT) and the type II NDT from the bacterium Desulfotalea psychrophila (DpNDT), the first psychrophilic NDT. This novel structural and biochemical data permitted an exhaustive comparative analysis aimed to shed light into the basis of the high global stability of the psychrophilic DpNDT, which has a higher melting temperature than EfPDT (58.5 °C versus 54.4 °C) or other mesophilic NDTs. DpNDT possesses a combination of unusual structural motifs not present neither in EfPDT nor any other NDT that most probably contribute to its global stability, in particular, a large aliphatic isoleucine-leucine-valine (ILV) bundle accompanied by a vicinal disulfide bridge and also an intersubunit disulfide bridge, the first described for an NDT. The functional and structural features of DpNDT do not fit the standard features of psychrophilic enzymes, which lead us to consider the implication of (sub)cellular levels together with the protein level in the adaptation of enzymatic activity to low temperatures.

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

Transcriptomics and Metabolomics Analyses Provide Novel Insights into Glucose-Induced Trophic Transition of the Marine Diatom Nitzschia laevis

Diatoms have important ecological roles and are natural sources of bioactive compounds. Nitzschia laevis is a member of marine diatoms that accumulates high-value products including fucoxanthin and eicosapentaenoic acid (EPA). In this study, physiological data showed that comparing to autotrophic growth, mixotrophic cultivation with glucose supplementation led to a decrease of chlorophyll and fucoxanthin content in N. laevis, and an increase of biomass density and EPA yield. To further examine the metabolic barriers for fucoxanthin and EPA biosynthesis, comparative transcriptomic and metabolome analyses were conducted, with a focus on the genes related to carotenoids biosynthesis and fatty acid metabolism. The results indicated that phytoene desaturase (PDS) and zeta-carotene isomerase (ZISO) could be the rate-limiting enzymes in carotenoid biosynthesis. The transcription regulation of 3-ketoacyl-CoA synthase (KCS) and elongation of very long chain fatty acids protein (EVOVL) are important contributors associated with polyunsaturated fatty acids (PUFAs) accumulation. Furthermore, we also investigated the glucose-associated regulatory genes using weighted gene co-expression network analysis, and identified potential hub genes linked with cell cycle, carbohydrate metabolism, purine biosynthesis, and lipid metabolism. This study offers a high-quality transcriptome resource for N. laevis and provides a molecular framework for further metabolic engineering studies on fucoxanthin and EPA production.

ProteinTools: a toolkit to analyze protein structures

The experimental characterization and computational prediction of protein structures has become increasingly rapid and precise. However, the analysis of protein structures often requires researchers to use several software packages or web servers, which complicates matters. To provide long-established structural analyses in a modern, easy-to-use interface, we implemented ProteinTools, a web server toolkit for protein structure analysis. ProteinTools gathers four applications so far, namely the identification of hydrophobic clusters, hydrogen bond networks, salt bridges, and contact maps. In all cases, the input data is a PDB identifier or an uploaded structure, whereas the output is an interactive dynamic web interface. Thanks to the modular nature of ProteinTools, the addition of new applications will become an easy task. Given the current need to have these tools in a single, fast, and interpretable interface, we believe that ProteinTools will become an essential toolkit for the wider protein research community. The web server is available at https://proteintools.uni-bayreuth.de.

A conserved folding nucleus sculpts the free energy landscape of bacterial and archaeal orthologs from a divergent TIM barrel family

Orthologous proteins from the three superkingdoms have conserved their structures and functions over evolutionary time. We ask whether their folding mechanisms and the structures of their partially folded states are similarly conserved, using bacterial and archaeal representatives of the IGPS TIM barrel enzyme. Comparison of circular dichroism and fluorescence spectroscopic studies reveal a highly conserved mechanism, and hydrogen–deuterium exchange mass spectrometry analyses highlight similar cores of stability in regions dominated by clusters of branched aliphatic side chains. A bioinformatics analysis of hundreds of IGPS sequences from each superkingdom shows a very highly conserved sequence, V/ILLI, that nucleates the formation of a misfolded, microsecond intermediate and has existed since the last universal common ancestor of the IGPS family of proteins.

The amino acid sequences of proteins have evolved over billions of years, preserving their structures and functions while responding to evolutionary forces. Are there conserved sequence and structural elements that preserve the protein folding mechanisms? The functionally diverse and ancient (βα)_1–8 TIM barrel motif may answer this question. We mapped the complex six-state folding free energy surface of a ∼3.6 billion y old, bacterial indole-3-glycerol phosphate synthase (IGPS) TIM barrel enzyme by equilibrium and kinetic hydrogen–deuterium exchange mass spectrometry (HDX-MS). HDX-MS on the intact protein reported exchange in the native basin and the presence of two thermodynamically distinct on- and off-pathway intermediates in slow but dynamic equilibrium with each other. Proteolysis revealed protection in a small (α1β2) and a large cluster (β5α5β6α6β7) and that these clusters form cores of stability in I_a and I_bp. The strongest protection in both states resides in β4α4 with the highest density of branched aliphatic side chain contacts in the folded structure. Similar correlations were observed previously for an evolutionarily distinct archaeal IGPS, emphasizing a key role for hydrophobicity in stabilizing common high-energy folding intermediates. A bioinformatics analysis of IGPS sequences from the three superkingdoms revealed an exceedingly high hydrophobicity and surprising α-helix propensity for β4, preceded by a highly conserved βα-hairpin clamp that links β3 and β4. The conservation of the folding mechanisms for archaeal and bacterial IGPS proteins reflects the conservation of key elements of sequence and structure that first appeared in the last universal common ancestor of these ancient proteins.

Control over the fibrillization yield by varying the oligomeric nucleation propensities of self-assembling peptides

Self-assembling peptides are an exemplary class of supramolecular biomaterials of broad biomedical utility. Mechanistic studies on the peptide self-assembly demonstrated the importance of the oligomeric intermediates towards the properties of the supramolecular biomaterials being formed. In this study, we demonstrate how the overall yield of the supramolecular assemblies are moderated through subtle molecular changes in the peptide monomers. This strategy is exemplified with a set of surfactant-like peptides (SLPs) with different β-sheet propensities and charged residues flanking the aggregation domains. By integrating different techniques, we show that these molecular changes can alter both the nucleation propensity of the oligomeric intermediates and the thermodynamic stability of the fibril structures. We demonstrate that the amount of assembled nanofibers are critically defined by the oligomeric nucleation propensities. Our findings offer guidance on designing self-assembling peptides for different biomedical applications, as well as insights into the role of protein gatekeeper sequences in preventing amyloidosis.

The role of gut microbiota (GM) and GM-related metabolites in diabetes and obesity. A review of analytical methods used to measure GM-related metabolites in fecal samples with a focus on metabolites' derivatization step

Disruption of gut microbiota (GM) composition is increasingly related to the pathogenesis of various metabolic diseases. Additionally, GM is responsible for the production and transformation of metabolites involved in the development of metabolic disorders, such as obesity and type 2 diabetes mellitus (T2DM). The current state of knowledge regarding the composition of GM and GM-related metabolites in relation to the progress and development of obesity and T2DM is presented in this review. To understand the relationships between GM-related metabolites and the development of metabolic disorders, their accurate qualitative and quantitative measurement in biological samples is needed. Feces represent a valuable biological matrix which composition may reflect the health status of the lower gastrointestinal tract and the whole organism. Mass spectrometry (MS), mainly in combination with gas chromatography (GC) or liquid chromatography (LC), is commonly used to measure fecal metabolites. However, profiling metabolites in such a complex matrix as feces is challenging from both analytical chemistry and biochemistry standpoints. Chemical derivatization is one of the most effective methods used to overcome these problems. In this review, we provide a comprehensive summary of the derivatization methods of GM-related metabolites prior to GC-MS or LC-MS analysis, which have been published in the last five years (2015-2020). Additionally, analytical methods used for the analysis of GM-related metabolites without the derivatization step are also presented.

CYP20-3 deglutathionylates 2-CysPRX A and suppresses peroxide detoxification during heat stress

In plants, growth-defense trade-offs occur because of limited resources, which demand prioritization towards either of them depending on various external and internal factors. However, very little is known about molecular mechanisms underlying their occurrence. Here, we describe that cyclophilin 20-3 (CYP20-3), a 12-oxo-phytodienoic acid (OPDA)-binding protein, crisscrosses stress responses with light-dependent electron reactions, which fine-tunes activities of key enzymes in plastid sulfur assimilations and photosynthesis. Under stressed states, OPDA, accumulates in the chloroplasts, binds and stimulates CYP20-3 to convey electrons towards serine acetyltransferase 1 (SAT1) and 2-Cys peroxiredoxin A (2CPA). The latter is a thiol-based peroxidase, protecting and optimizing photosynthesis by reducing its toxic byproducts (e.g., H2O2). Reduction of 2CPA then inactivates its peroxidase activity, suppressing the peroxide detoxification machinery, whereas the activation of SAT1 promotes thiol synthesis and builds up reduction capacity, which in turn triggers the retrograde regulation of defense gene expressions against abiotic stress. Thus, we conclude that CYP20-3 is a unique metabolic hub conveying resource allocations between plant growth and defense responses (trade-offs), ultimately balancing optimal growth phonotype.

The hot sites of α-synuclein in amyloid fibril formation

The role of alpha-synuclein (αS) amyloid fibrillation has been recognized in various neurological diseases including Parkinson's Disease (PD). In early stages, fibrillation occurs by the structural transition from helix to extended states in monomeric αS followed by the formation of beta-sheets. This alpha-helix to beta-sheet transition (αβT) speeds up the formation of amyloid fibrils through the formation of unstable and temporary configurations of the αS. In this study, the most important regions that act as initiating nuclei and make unstable the initial configuration were identified based on sequence and structural information. In this regard, a Targeted Molecular Dynamics (TMD) simulation was employed using explicit solvent models under physiological conditions. Identified regions are those that are in the early steps of structural opening. The trajectory was clustered the structures characterized the intermediate states. The findings of this study would help us to better understanding of the mechanism of amyloid fibril formation.

The acidic tail of HMGB1 regulates its secondary structure and conformational flexibility: A circular dichroism and molecular dynamics simulation study

High mobility group box 1 (HMGB1) is a damage-associated molecular pattern (DAMP) molecule that triggers the progression of several pro-inflammatory diseases such as diabetes, Alzheimer's disease and cancer, by inducing signals upon interaction with the receptors such as the receptor for advanced glycation end-products (RAGE) and toll-like receptors (TLRs). The acidic C-terminal tail of HMGB1 is an intrinsically disordered region of the protein which is known to determine the interaction of HMGB1 to DNA and histones. This study characterizes its structural properties using a combination of circular dichroism (CD) and molecular dynamics (MD) simulations. The full-length and tail-less forms of HMGB1 were compared to rationalise the role of the acidic tail in maintaining the stability of the entire structure of HMGB1 in atomistic detail. Consistent with experimental data, the acidic tail was predicted to adopt an extended conformation that allows it to make a range of hydrogen-bonding and electrostatic interactions with the box-like domains that stabilize the overall structure of HMGB1. Absence of the acidic tail was predicted to increase structural fluctuations of all amino acids, leading to changes in secondary structure from α-helical to more hydrophilic turns along with increased exposure of multiple amino acids to the surrounding solvent. These structural changes reveal the intrinsic conformational dynamics of HMGB1 that are likely to affect the accessibility of its receptors.

Revealing the Structural Contributions to Thermal Adaptation of the TATA-Box Binding Protein: Molecular Dynamics and QSPR Analyses

The TATA-box binding protein (TBP) is an important element of the transcription machinery in archaea and eukaryotic organisms. TBP is expressed in organisms adapted to different temperatures, indicating a robust structure, and experimental studies have shown that the mid-unfolding temperature (Tm) of TBP is directly correlated with the optimal growth temperature (OGT) of the organism. To understand which are the relevant structural requirements for its stability, we present the first structural and dynamic computational study of TBPs, combining molecular dynamics (MD) simulations and a quantitative structure-property relationship (QSPR) over a set of TBPs of organisms adapted to different temperatures. We found that the main structural properties of TBP used to adapt to high temperatures are an increase in the ease of desolvation of charged residues at the surface, an increase in the local resiliency, the presence of Leu clusters in the protein core, and an increase in the loss of hydrophobic packing in the N-terminal subdomain. In view of our results, we consider that TBP is a good model to study thermal adaptation, and our analysis opens the possibility of performing protein engineering on TBPs to study transcription at high or low temperatures.

Frustration and folding of a TIM barrel protein

Triosephosphate isomerase (TIM) barrel proteins have not only a conserved architecture that supports a myriad of enzymatic functions, but also a conserved folding mechanism that involves on- and off-pathway intermediates. Although experiments have proven to be invaluable in defining the folding free-energy surface, they provide only a limited understanding of the structures of the partially folded states that appear during folding. Coarse-grained simulations employing native centric models are capable of sampling the entire energy landscape of TIM barrels and offer the possibility of a molecular-level understanding of the readout from sequence to structure. We have combined sequence-sensitive native centric simulations with small-angle X-ray scattering and time-resolved Förster resonance energy transfer to monitor the formation of structure in an intermediate in the Sulfolobus solfataricus indole-3-glycerol phosphate synthase TIM barrel that appears within 50 μs and must at least partially unfold to achieve productive folding. Simulations reveal the presence of a major and 2 minor folding channels not detected in experiments. Frustration in folding, i.e., backtracking in native contacts, is observed in the major channel at the initial stage of folding, as well as late in folding in a minor channel before the appearance of the native conformation. Similarities in global and pairwise dimensions of the early intermediate, the formation of structure in the central region that spreads progressively toward each terminus, and a similar rate-limiting step in the closing of the β-barrel underscore the value of combining simulation and experiment to unravel complex folding mechanisms at the molecular level.

Unfolded states under folding conditions accommodate sequence-specific conformational preferences with random coil-like dimensions

Proteins are marginally stable molecules that fluctuate between folded and unfolded states. Here, we provide a high-resolution description of unfolded states under refolding conditions for the N-terminal domain of the L9 protein (NTL9). We use a combination of time-resolved Förster resonance energy transfer (FRET) based on multiple pairs of minimally perturbing labels, time-resolved small-angle X-ray scattering (SAXS), all-atom simulations, and polymer theory. Upon dilution from high denaturant, the unfolded state undergoes rapid contraction. Although this contraction occurs before the folding transition, the unfolded state remains considerably more expanded than the folded state and accommodates a range of local and nonlocal contacts, including secondary structures and native and nonnative interactions. Paradoxically, despite discernible sequence-specific conformational preferences, the ensemble-averaged properties of unfolded states are consistent with those of canonical random coils, namely polymers in indifferent (theta) solvents. These findings are concordant with theoretical predictions based on coarse-grained models and inferences drawn from single-molecule experiments regarding the sequence-specific scaling behavior of unfolded proteins under folding conditions.

Characterization of TDP-43 RRM2 Partially Folded States and Their Significance to ALS Pathogenesis

The human protein TDP-43 is a major component of the cellular aggregates found in amyotrophic lateral sclerosis and other neurodegenerative diseases. Insoluble cytoplasmic aggregates isolated from the brain of amyotrophic lateral sclerosis and frontotemporal lobar degeneration patients contain ubiquitinated, hyperphosphorylated, and N-terminally truncated TDP-43. Truncated fragments of TDP-43 identified from patient tissues contain part of the second RNA recognition motif (RRM2) and the disordered C-terminus, indicating that both domains can be involved in aggregation and toxicity. Here, we focus on RRM2. Using all-atom replica-averaged metadynamics simulations with NMR chemical shift restraints, we characterized the atomic structure of non-native states of RRM2, sparsely populated under native conditions. These structures reveal the exposure to the solvent of aggregation-prone peptide regions, normally buried in the native state, supporting a role in aggregation for the partially folded states of RRM2.

Structural and functional impact of troponin C-mediated Ca2+ sensitization on myofilament lattice spacing and cross-bridge mechanics in mouse cardiac muscle

Acto-myosin cross-bridge kinetics are important for beat-to-beat regulation of cardiac contractility; however, physiological and pathophysiological mechanisms for regulation of contractile kinetics are incompletely understood. Here we explored whether thin filament-mediated Ca²⁺ sensitization influences cross-bridge kinetics in permeabilized, osmotically compressed cardiac muscle preparations. We used a murine model of hypertrophic cardiomyopathy (HCM) harboring a cardiac troponin C (cTnC) Ca²⁺-sensitizing mutation, Ala8Val in the regulatory N-domain. We also treated wild-type murine muscle with bepridil, a cTnC-targeting Ca²⁺ sensitizer. Our findings suggest that both methods of increasing myofilament Ca²⁺ sensitivity increase cross-bridge cycling rate measured by the rate of tension redevelopment (k_TR); force per cross-bridge was also enhanced as measured by sinusoidal stiffness and I_1,1/I_1,0 ratio from X-ray diffraction. Computational modeling suggests that Ca²⁺ sensitization through this cTnC mutation or bepridil accelerates k_TR primarily by promoting faster cross-bridge detachment. To elucidate if myofilament structural rearrangements are associated with changes in k_TR, we used small angle X-ray diffraction to simultaneously measure myofilament lattice spacing and isometric force during steady-state Ca²⁺ activations. Within in vivo lattice dimensions, lattice spacing and steady-state isometric force increased significantly at submaximal activation. We conclude that the cTnC N-domain controls force by modulating both the number and rate of cycling cross-bridges, and that the both methods of Ca²⁺ sensitization may act through stabilization of cTnC's D-helix. Furthermore, we propose that the transient expansion of the myofilament lattice during Ca²⁺ activation may be an additional factor that could increase the rate of cross-bridge cycling in cardiac muscle. These findings may have implications for the pathophysiology of HCM.

A Hyperthermophilic Phage Decoration Protein Suggests Common Evolutionary Origin with Herpesvirus Triplex Proteins and an Anti-CRISPR Protein

Virus capsids are protein shells that protect the viral genome from environmental assaults, while maintaining the high internal pressure of the tightly packaged genome. To elucidate how capsids maintain stability under harsh conditions, we investigated the capsid components of the hyperthermophilic phage P74-26. We determined the structure of capsid protein gp87 and show that it has the same fold as decoration proteins in many other phages, despite lacking significant sequence homology. We also find that gp87 is significantly more stable than mesophilic homologs. Our analysis of the gp87 structure reveals that the core "β tulip" domain is conserved in trimeric capsid components across numerous double-stranded DNA viruses, including Herpesviruses. Moreover, this β barrel domain is found in anti-CRISPR protein AcrIIC1, suggesting a mechanism for the evolution of this Cas9 inhibitor. Our work illustrates the principles for increased stability of gp87, and extends the evolutionary reach of the β tulip domain.

Spontaneous refolding of the large multidomain protein malate synthase G proceeds through misfolding traps

Most protein folding studies until now focus on single domain or truncated proteins. Although great insights in the folding of such systems has been accumulated, very little is known regarding the proteins containing multiple domains. It has been shown that the high stability of domains, in conjunction with inter-domain interactions, manifests as a frustrated energy landscape, causing complexity in the global folding pathway. However, multidomain proteins despite containing independently foldable, loosely cooperative sections can fold into native states with amazing speed and accuracy. To understand the complexity in mechanism, studies were conducted previously on the multidomain protein malate synthase G (MSG), an enzyme of the glyoxylate pathway with four distinct and adjacent domains. It was shown that the protein refolds to a functionally active intermediate state at a fast rate, which slowly produces the native state. Although experiments decoded the nature of the intermediate, a full description of the folding pathway was not elucidated. In this study, we use a battery of biophysical techniques to examine the protein's folding pathway. By using multiprobe kinetics studies and comparison with the equilibrium behavior of protein against urea, we demonstrate that the unfolded polypeptide undergoes conformational compaction to a misfolded intermediate within milliseconds of refolding. The misfolded product appears to be stabilized under moderate denaturant concentrations. Further folding of the protein produces a stable intermediate, which undergoes partial unfolding-assisted large segmental rearrangements to achieve the native state. This study reveals an evolved folding pathway of the multidomain protein MSG, which involves surpassing the multiple misfolding traps during refolding.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Modulations of Keap1-Nrf2 signaling axis by TIIA ameliorated the oxidative stress-induced myocardial apoptosis

Mounting evidence has strongly implicated oxidative stress in the development of cardiac dysfunction, and myocardial apoptosis contributes to the pathogenesis of heart failure. Quantitative cardiac proteomics data revealed that pressure load by TAC resulted in a significant decline in mitochondrial metabolic activity, where TIIA (Tanshinone IIA sulfonate) treatment reversed it in vivo, which might be mediated by Nrf2. In NRVMs, TIIA treatment ameliorated H₂O₂-induced caspase-3/9 activations through the suppression of p38 and mTOR signaling pathways, where caspase-mediated cleavage of YY1 and PARP resulted in the defects in mitochondrial biogenesis and DNA repair, and this event finally led to cardiomyocyte apoptosis. Mass spectrometry analysis showed that TIIA hydrophobically interacted with Keap1 (the cytoplasmic repressor of Nrf2) and induced its degradation in vitro. Site-directed mutagenesis of Keap1 identified V122/V123/I125 to be the critical residues for the TIIA-induced de-dimerization and degradation of Keap1. Besides, TIIA treatment also epigenetically up-regulated Nrf2 gene transcription, where it hypomethylated the first 5 CpGs of Nrf2 promoter. Furthermore, cardiac-specific Nrf2 knockout mice exhibited the significantly dampened anti-apoptotic effects of TIIA.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Structural Rearrangement upon Fragmentation of the Stability Core of the ALS-Linked Protein TDP-43

Amyotrophic lateral sclerosis (ALS) is the most common adult degenerative motor neuron disease. Experimental evidence indicates a direct role of transactive-response DNA-binding protein 43 (TDP-43) in the pathology of ALS and other neurodegenerative diseases. TDP-43 has been identified as a major component of cytoplasmic inclusions in patients with sporadic ALS; however, the molecular basis of the disease mechanism is not yet fully understood. Fragmentation within the second RNA recognition motif (RRM2) of TDP-43 has been observed in patient tissues and may play a role in the formation of aggregates in disease. To determine the structural and dynamical changes resulting from the truncation that could lead to aggregation and toxicity, we performed molecular dynamics simulations of the full-length RRM2 domain (the stability core of TDP-43) and of a truncated variant (where residues 189-207 are deleted to mimic a site of cleavage within RRM2 found in ALS patients). Our simulations show heterogeneous structural reorganization and decreased stability of the truncated RRM2 domain compared to the full-length domain, consistent with previous experimental results. The decreased stability and structural reorganization in the truncated RRM2 result in a higher probability of protein-protein interactions through altered electrostatic surface charges and increased accessibility of hydrophobic residues (including the nuclear export sequence), providing a rationale for the increased cytoplasmic aggregation of RRM2 fragments seen in sporadic ALS patients.

The Hsp70 interdomain linker is a dynamic switch that enables allosteric communication between two structured domains

Hsp70 molecular chaperones play key roles in cellular protein homeostasis by binding to exposed hydrophobic regions of incompletely folded or aggregated proteins. This crucial Hsp70 function relies on allosteric communication between two well-structured domains: an N-terminal nucleotide-binding domain (NBD) and a C-terminal substrate-binding domain (SBD), which are tethered by an interdomain linker. ATP or ADP binding to the NBD alters the substrate-binding affinity of the SBD, triggering functionally essential cycles of substrate binding and release. The interdomain linker is a well-structured participant in the interdomain interface in ATP-bound Hsp70s. By contrast, in the ADP-bound state, exemplified by the Escherichia coli Hsp70 DnaK, the interdomain linker is flexible. Hsp70 interdomain linker sequences are highly conserved; moreover, mutations in this region compromise interdomain allostery. To better understand the role of this region in Hsp70 allostery, we used molecular dynamics simulations to explore the conformational landscape of the interdomain linker in ADP-bound DnaK and supported our simulations by strategic experimental data. We found that while the interdomain linker samples many conformations, it behaves as three relatively ordered segments connected by hinges. As a consequence, the distances and orientations between the NBD and SBD are limited. Additionally, the C-terminal region of the linker forms previously unreported, transient interactions with the SBD, and the predominant linker-docking site is available in only one allosteric state, that with high affinity for substrate. This preferential binding implicates the interdomain linker as a dynamic allosteric switch. The linker-binding site on the SBD is a potential target for small molecule modulators of the Hsp70 allosteric cycle.

Correlation of fitness landscapes from three orthologous TIM barrels originates from sequence and structure constraints

Sequence divergence of orthologous proteins enables adaptation to environmental stresses and promotes evolution of novel functions. Limits on evolution imposed by constraints on sequence and structure were explored using a model TIM barrel protein, indole-3-glycerol phosphate synthase (IGPS). Fitness effects of point mutations in three phylogenetically divergent IGPS proteins during adaptation to temperature stress were probed by auxotrophic complementation of yeast with prokaryotic, thermophilic IGPS. Analysis of beneficial mutations pointed to an unexpected, long-range allosteric pathway towards the active site of the protein. Significant correlations between the fitness landscapes of distant orthologues implicate both sequence and structure as primary forces in defining the TIM barrel fitness landscape and suggest that fitness landscapes can be translocated in sequence space. Exploration of fitness landscapes in the context of a protein fold provides a strategy for elucidating the sequence-structure-fitness relationships in other common motifs.

The TIM barrel fold is an evolutionarily conserved motif found in proteins with a variety of enzymatic functions. Here the authors explore the fitness landscape of the TIM barrel protein IGPS and uncover evolutionary constraints on both sequence and structure, accompanied by long range allosteric interactions.

Insights from coarse-grained Gō models for protein folding and dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.