Protein Denaturation

Thermodynamics of modulation of interaction of α-helix inducer 2, 2, 2-trifluoroethanol with lysozyme in presence of cationic, anionic and non-ionic surfactants

The interactions of anionic sodium dodecyl sulphate (SDS), cationic cetyltrimethylammonium bromide (CTAB) and nonionic triton X-100 (TX-100) surfactants with lysozyme at pH = 2.4 have been studied individually as well as in combination with 2,2,2-trifluoroetanol (TFE). Urea has also been used in combination with surfactants. By using these combinations, efforts have been made to obtain partially folded conformations of the protein in the presence of surfactants and effect of α-helix inducer 2,2,2-trifluoroethanol on these intermediate states. Thermodynamic analysis of all these interactions has been done employing a combination of UV-visible, fluorescence and circular dichroism spectroscopies. The results have been correlated with each other and characterized qualitatively as well as quantitatively. At lower concentration of surfactant, the thermodynamic parameters indicated the destabilizing effect of SDS, stabilizing effect of CTAB and unappreciable destabilizing impact of TX-100 on lysozyme. The enhancement in destabilization effect or reduction in stabilization effect of surfactants on lysozyme in the presence of TFE and urea has also been indicated.Communicated by Ramaswamy H. Sarma.

Solvent accessible surface area-assessed molecular basis of osmolyte-induced protein stability

In solvent-modulated protein folding, under certain physiological conditions, an equilibrium exists between the unfolded and folded states of the protein without any need to break or make a covalent bond. In this process, interactions between various protein groups (peptides) and solvent molecules are known to play a major role in determining the directionality of the chemical reaction. However, an understanding of the mechanism of action of the co(solvent) by a generic theoretical underpinning is lacking. In this study, a generic solvation model is developed based on statistical mechanics and the thermodynamic transfer free energy model by considering the microenvironment polarity of the interacting co(solvent)-protein system. According to this model, polarity and the fractional solvent-accessible surface areas contribute to the interaction energies. The present model includes various orientations of participating interactant solvent surfaces of suitable areas. As model systems, besides the backbone we consider naturally occurring amino acid residues solvated in ten different osmolytes, small organic compounds known to modulate protein stability. The present model is able to predict the correct trend of the osmolyte-peptide interactions ranging from stabilizing to destabilizing not only for the backbone but also for side chains. Our model predicts Asn, Gln, Asp, Glu, Arg and Pro to be highly stable in most of the protecting osmolytes while Ala, Val, Ile, Leu, Thr, Met, Lys, Phe, Trp and Tyr are predicted to be moderately stable, and Ser, Cys and Histidine are predicted to be least stable. However, in denaturing solvents, both backbone and side chain models show similar stabilities in urea and guanidine. One of the important aspects of this model is that it is essentially parameter-free and consistent with the electrostatics of the interaction partners that make this model suitable for estimating any solute-solvent interaction energies.

Disordered proteins interact with the chemical environment to tune their protective function during drying

The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins combined with the exposure of their residues accounts for this sensitivity. One context in which IDPs play important roles that is concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family, synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet, the mechanisms underlying this synergy differ between IDP families.

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Probing protein aggregation through spectroscopic insights and multimodal approaches: A comprehensive review for counteracting neurodegenerative disorders

Aberrant accumulation of protein misfolding can cause aggregation and fibrillation and is one of the primary characteristic features of neurodegenerative diseases. Because they are disordered, misfolded, and aggregated proteins pose a significant setback in drug designing. The structural study of intermediate steps in these kinds of aggregated proteins will allow us to determine the conformational changes as well as the probable pathways encompassing various neurodegenerative disorders. The analysis of protein aggregates involved in neurodegenerative diseases relies on a diverse toolkit of biophysical techniques, encompassing both morphological and non-morphological methods. Additionally, Thioflavin T (ThT) assays and Circular Dichroism (CD) spectroscopy facilitate investigations into aggregation kinetics and secondary structure alterations. The collective application of these biophysical techniques empowers researchers to comprehensively unravel the intricate nature of protein aggregates associated with neurodegeneration. Furthermore, the topics covered in this review have summed up a handful of well-established techniques used for the structural analysis of protein aggregation. This multifaceted approach advances our fundamental understanding of the underlying mechanisms driving neurodegenerative diseases and informs potential therapeutic strategies.

Do Ionic Liquids Exhibit the Required Characteristics to Dissolve, Extract, Stabilize, and Purify Proteins? Past-Present-Future Assessment

Proteins are highly labile molecules, thus requiring the presence of appropriate solvents and excipients in their liquid milieu to keep their stability and biological activity. In this field, ionic liquids (ILs) have gained momentum in the past years, with a relevant number of works reporting their successful use to dissolve, stabilize, extract, and purify proteins. Different approaches in protein-IL systems have been reported, namely, proteins dissolved in (i) neat ILs, (ii) ILs as co-solvents, (iii) ILs as adjuvants, (iv) ILs as surfactants, (v) ILs as phase-forming components of aqueous biphasic systems, and (vi) IL-polymer-protein/peptide conjugates. Herein, we critically analyze the works published to date and provide a comprehensive understanding of the IL-protein interactions affecting the stability, conformational alteration, unfolding, misfolding, and refolding of proteins while providing directions for future studies in view of imminent applications. Overall, it has been found that the stability or purification of proteins by ILs is bispecific and depends on the structure of both the IL and the protein. The most promising IL-protein systems are identified, which is valuable when foreseeing market applications of ILs, e.g., in "protein packaging" and "detergent applications". Future directions and other possibilities of IL-protein systems in light-harvesting and biotechnology/biomedical applications are discussed.

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

N-acetyl amino acid amide solubility in aqueous 1,6-hexanediol solutions: Insights into the protein droplet deformation mechanism

Proteinaceous liquid droplets, generated by liquid-liquid phase separation, function as membraneless compartments that are essential for diverse biological functions. Studies addressing droplet generation have used 1,6-hexanediol (1,6-HD) as a droplet-discerning agent owing to its capacity to induce droplet deformation. Despite the empirical utility of 1,6-HD, the mechanism underlying 1,6-HD-induced droplet deformation remains unknown. In this study, the solubilities of N-acetyl amino acid amides, which correspond to proteinogenic amino acid residues, were examined in the presence of 1,6-HD at 25 °C. Other solvents included ethanol, 1-propanol, and amides. Remarkably, 1,6-HD effectively solubilized hydrophobic species (particularly aromatic species) and exhibited reduced efficacy in solubilizing hydrophilic species and peptide bond moieties. These solubilizing effects are reflected in changes in protein solubility and structure. Specifically, 1,6-HD primarily targets the hydrophobic regions of a protein, increasing protein solubility without causing substantial structural changes. This solubilization mechanism is essential for elucidating the role of 1,6-HD as a droplet-discerning agent and recognizing its potential limitations.

pH-dependent interactions of biologically important metal ions with hen egg white lysozyme based on its hydration properties: Thermodynamic and mechanistic insights

Importance of metal ion selectivity in biomolecules and their key role in proteins are widely explored. However, understanding the thermodynamics of how hydrated metal ions alter the protein hydration and their conformation is also important. In this study, the interaction of some biologically important Ca2+, Mn2+, Co2+, Cu2+, and Zn2+ ions with hen egg white lysozyme at pH 2.1, 3.0, 4.5 and 7.4 has been investigated. Intrinsic fluorescence studies have been employed for metal ion-induced protein conformational changes analysis. Thermostability based on protein hydration has been investigated using differential scanning calorimetry (DSC). Thermodynamic parameters emphasizing on metal ion-protein binding mechanistic insights have been well discussed using isothermal titration calorimetry (ITC). Overall, these experiments have reported that their interactions are pH-dependent and entropically driven. This research also reports the strongly hydrated metal ions as water structure breaker unlike osmolytes based on DSC studies. These experimental results have highlighted higher concentrations of different metal ions effect on the protein hydration and thermostability which might be helpful in understanding their interactions in aqueous solutions.

Actual Amount Adsorbed as Estimated from the Surface Excess Isotherm

The amount of adsorption at equilibrium is commonly used for reporting solid/solution isotherms, despite the admonishment by the International Union of Pure and Applied Chemistry (IUPAC) against equating the surface excess (i.e., the measurable quantity for sorption, signifying the competitive sorption of adsorbate and solvent) with the actual amount adsorbed. The consensus, more generally stated, is that the surface excess cannot be divided into individual isotherms for sorbate and solvent unless simplifying model assumptions are introduced. Here we show, contrary to the IUPAC report, that there exists a simple method for assigning the total isotherm to the sorbate's actual amount adsorbed and to the individual solute isotherm. This requires a combination of isotherm and volumetric measurements. For dilute sorbates, we establish criteria to show if the total isotherm is dominated by the amount of sorption at the interface, in agreement with the common assumption in the practical literature. In the absence of the volume data, we propose an approximate yet more versatile method based on the specific surface area to carry out order-of-magnitude analysis to examine whether the actual amount adsorbed dominates surface excess. Application of our methods to the adsorption of sodium decyl sulfate on polystyrene latex, malachite green on activated carbons, and thiophenes on a metal-organic framework all demonstrated the dominance of the actual amount adsorbed, significantly simplifying isotherm analysis in terms of the underlying interactions (i.e., surface-sorbate and net self-interactions at the interface), eliminating the need for excess surface quantities. Analysis of fully miscible solvent-sorbate isotherms (e.g., the mixtures of organic solvents adsorbed on mesoporous silica and carbonaceous adsorbents) indicates the contributions from both sorbate and solvent isotherms.

K-Pro: Kinetics Data on Proteins and Mutants

The study of protein folding plays a crucial role in improving our understanding of protein function and of the relationship between genetics and phenotypes. In particular, understanding the thermodynamics and kinetics of the folding process is important for uncovering the mechanisms behind human disorders caused by protein misfolding. To address this issue, it is essential to collect and curate experimental kinetic and thermodynamic data on protein folding. K-Pro is a new database designed for collecting and storing experimental kinetic data on monomeric proteins, with a two-state folding mechanism. With 1,529 records from 62 proteins corresponding to 65 structures, K-Pro contains various kinetic parameters such as the logarithm of the folding and unfolding rates, Tanford's β and the ϕ values. When available, the database also includes thermodynamic parameters associated with the kinetic data. K-Pro features a user-friendly interface that allows browsing and downloading kinetic data of interest. The graphical interface provides a visual representation of the protein and mutants, and it is cross-linked to key databases such as PDB, UniProt, and PubMed. K-Pro is open and freely accessible through https://folding.biofold.org/k-pro and supports the latest versions of popular browsers.

Molecular mechanism of the common and opposing cosolvent effects of fluorinated alcohol and urea on a coiled coil protein

Alcohols and urea are widely used as effective protein denaturants. Among monohydric alcohols, 2,2,2-trifluoroethanol (TFE) has large cosolvent effects as a helix stabilizer in proteins. In contrast, urea efficiently denatures ordered native structures, including helices, into coils. These opposing cosolvent effects of TFE and urea are well known, even though both preferentially bind to proteins; however, the underlying molecular mechanism remains controversial. Cosolvent-dependent relative stability between native and denatured states is rigorously related to the difference in preferential binding parameters (PBPs) between these states. In this study, GCN4-p1 with two-stranded coiled coil helices was employed as a model protein, and molecular dynamics simulations for the helix dimer and isolated coil were conducted in aqueous solutions with 2 M TFE and urea. As 2 M cosolvent aqueous solutions did not exhibit clustering of cosolvent molecules, we were able to directly investigate the molecular origin of the excess PBP without considering the enhancement effect of PBPs arising from the concentration fluctuations. The calculated excess PBPs of TFE for the helices and those of urea for the coils were consistent with experimentally observed stabilization of helix by TFE and that of coil by urea. The former was caused by electrostatic interactions between TFE and side chains of the helices, while the latter was attributed to both electrostatic and dispersion interactions between urea and the main chains. Unexpectedly, reverse-micelle-like orientations of TFE molecules strengthened the electrostatic interactions between TFE and the side chains, resulting in strengthening of TFE solvation.

JAK3 Y841 Autophosphorylation Is Critical for STAT5B Activation, Kinase Domain Stability and Dimer Formation

Janus tyrosine kinase 3 (JAK3) is primarily expressed in immune cells and is needed for signaling by the common gamma chain (γc) family of cytokines. Abnormal JAK3 signal transduction can manifest as hematological disorders, e.g., leukemia, severe combined immunodeficiency (SCID) and autoimmune disease states. While regulatory JAK3 phosphosites have been well studied, here a functional proteomics approach coupling a JAK3 autokinase assay to mass spectrometry revealed ten previously unreported autophosphorylation sites (Y105, Y190, Y238, Y399, Y633, Y637, Y738, Y762, Y824, and Y841). Of interest, Y841 was determined to be evolutionarily conserved across multiple species and JAK family members, suggesting a broader role for this residue. Phospho-substitution mutants confirmed that Y841 is also required for STAT5 tyrosine phosphorylation. The homologous JAK1 residue Y894 elicited a similar response to mutagenesis, indicating the shared importance for this site in JAK family members. Phospho-specific Y841-JAK3 antibodies recognized activated kinase from various T-cell lines and transforming JAK3 mutants. Computational biophysics analysis linked Y841 phosphorylation to enhanced JAK3 JH1 domain stability across pH environments, as well as to facilitated complementary electrostatic JH1 dimer formation. Interestingly, Y841 is not limited to tyrosine kinases, suggesting it represents a conserved ubiquitous enzymatic function that may hold therapeutic potential across multiple kinase families.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape - all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs - questioning the differences and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the overall consensus protein is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted region is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order couplings using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Analysis of the Shear Stresses in a Filling Line of Parenteral Products: The Role of Tubing

Parenteral products appear to be sensitive to process conditions in bioprocessing steps, such as interfacial stress and shear stress. The combination of these elements is widely believed and proven to influence product stability, but the defined roles of these players in the product damage process have not yet been identified. The present work addresses a current industrial problem, by focusing on the analysis of shear stress on protein-based therapeutics flowing in tubing by means of Computational Fluid Dynamics simulations. The purpose of this article is not to pinpoint the mechanism triggering the damage of the product, but it represents the first step towards wider experimental investigations and introduces a new strategy to quantify the average shear stress. The field of scale-down approaches, used to scale the commercial process down to the laboratory level, is also explored. Since quality control is critical in the pharmaceutical realm, it is essential that the scale-down approach preserves the same stress exposure as the commercial scale, which in the present work is considered to be that resulting from shear effects. Therefore, a new approach for scaling down the commercial process is proposed, which has been compared with traditional approaches and shown to provide greater representativeness between the two scales.

The B domain of protein A retains residual structures in 6 M guanidinium chloride as revealed by hydrogen/deuterium-exchange NMR spectroscopy

The characterization of residual structures persistent in unfolded proteins is an important issue in studies of protein folding, because the residual structures present, if any, may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the residual structures of the isolated B domain (BDPA) of staphylococcal protein A in 6 M guanidinium chloride. BDPA is a small three-helix-bundle protein, and until recently its folding/unfolding reaction has been treated as a simple two-state process between the native and the fully unfolded states. We employed a dimethylsulfoxide (DMSO)-quenched hydrogen/deuterium (H/D)-exchange 2D NMR techniques with the use of spin desalting columns, which allowed us to investigate the H/D-exchange behavior of individually identified peptide amide (NH) protons. We obtained H/D-exchange protection factors of the 21 NH protons that form an α-helical hydrogen bond in the native structure, and the majority of these NH protons were significantly protected with a protection factor of 2.0-5.2 in 6 M guanidinium chloride, strongly suggesting that these weakly protected NH protons form much stronger hydrogen bonds under native folding conditions. The results can be used to deduce the structure of an early folding intermediate, when such an intermediate is shown by other methods. Among three native helical regions, the third helix in the C-terminal side was highly protected and stabilized by side-chain salt bridges, probably acting as the folding initiation site of BDPA. The present results are discussed in relation to previous experimental and computational findings on the folding mechanisms of BDPA.

Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome

Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG^0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H ⁺ activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG^0.3 is between 3 and 4 k _B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG^0.3 ∼ 10k _B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG^0.3 = 2.2k _B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG^0.3 = 5.2k _B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

Dry Molten Globule-Like Intermediates in Protein Folding, Function, and Disease

The performance of a protein depends on its correct folding to the final functional native form. Hence, understanding the process of protein folding has remained an important field of research for the scientific community for the past five decades. Two important intermediate states, namely, wet molten globule (WMG) and dry molten globule (DMG), have emerged as critical milestones during protein folding-unfolding reactions. While much has been discussed about WMGs as a common unfolding intermediate, the evidence for DMGs has remained elusive owing to their near-native features, which makes them difficult to probe using global structural probes. This Review puts together the available literature and new evidence on DMGs to give a broader perspective on the universality of DMGs and discuss their significance in protein folding, function, and disease.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

The pH Robustness of Bacterial Sensing

Bacteria have evolved many different signal transduction systems to sense and respond to changing environmental conditions. Signal integration is mainly achieved by signal recognition at extracytosolic ligand-binding domains (LBDs) of receptors. Hundreds of different LBDs have been reported, and our understanding of their sensing properties is growing. Receptors must function over a range of environmental pH values, but there is little information available on the robustness of sensing as a function of pH. Here, we have used isothermal titration calorimetry to determine the pH dependence of ligand recognition by nine LBDs that cover all major LBD superfamilies, of periplasmic solute-binding proteins, and cytosolic LBDs. We show that periplasmic LBDs recognize ligands over a very broad pH range, frequently stretching over eight pH units. This wide pH range contrasts with a much narrower pH response range of the cytosolic LBDs analyzed. Many LBDs must be dimeric to bind ligands, and analytical ultracentrifugation studies showed that the LBD of the Tar chemoreceptor forms dimers over the entire pH range tested. The pH dependences of Pseudomonas aeruginosa motility and chemotaxis were bell-shaped and centered at pH 7.0. Evidence for pH robustness of signaling in vivo was obtained by Förster Resonance Energy Transfer (FRET) measurements of the chemotaxis pathway responses in Escherichia coli. Bacteria have evolved several strategies to cope with extreme pH, such as periplasmic chaperones for protein refolding. The intrinsic pH resistance of periplasmic LBDs appears to be another strategy that permits bacteria to survive under adverse conditions.

The influence of cations on α-lactalbumin amyloid aggregation

There is limited knowledge regarding α-lactalbumin amyloid aggregation and its mechanism. We examined the formation of α-lactalbumin amyloid fibrils (α-LAF) in the presence of cations (Mg²⁺, Ca²⁺, Na⁺, K⁺, NH₄⁺, and Cs⁺) in the form of chloride salts at two concentrations. We have shown that studied cations affect the conformation of α-lactalbumin, the kinetics of its amyloid formation, morphology, and secondary structure of α-LAF in a different manner. The higher salts concentration significantly accelerated the aggregation process. Both salt concentrations stabilized α-lactalbumin's secondary structure. However, the presence of divalent cations resulted in shorter fibrils with less β-sheet content. Moreover, strongly hydrated Mg²⁺ significantly altered α-lactalbumin's tertiary structure, followed by Na⁺, NH₄⁺, K⁺, and weakly hydrated Cs⁺. On the other hand, Ca²⁺, despite being also strongly hydrated, stabilized the tertiary structure, supposedly due to its high affinity towards α-lactalbumin. Yet, Ca²⁺ was not able to inhibit α-lactalbumin amyloid aggregation.

Protein stability [determination] problems

Human health depends on the correct folding of proteins, for misfolding and aggregation lead to diseases. An unfolded (denatured) protein can refold to its original folded state. How does this occur is known as the protein folding problem. One of several related questions to this problem is that how much more stable is the folded state than the unfolded state. There are several measures of protein stability. In this article, protein stability is given a thermodynamic definition and is measured by Gibbs free energy change (ΔGD0) associated with the equilibrium, native (N) conformation ↔ denatured (D) conformation under the physiological condition usually taken as dilute buffer (or water) at 25 °C. We show that this thermodynamic quantity (ΔGD0), where subscript D represents transition between N and D states, and superscript 0 (zero) represents the fact that the transition occurs in the absence of denaturant, can be neither measured nor predicted under physiological conditions. However, ΔGD can be measured in the presence of strong chemical denaturants such as guanidinium chloride and urea which are shown to destroy all noncovalent interactions responsible for maintaining the folded structure. A problem with this measurement is that the estimate of ΔGD0 comes from the analysis of the plot of ΔGDversus denaturant concentration, which requires a long extrapolation of values of ΔGD, and all the three methods of extrapolation give three different values of ΔGD0 for a protein. Thus, our confidence in the authentic value of ΔGD0 is eroded. Another problem with this in vitro measurement of ΔGD0 is that it is done on the pure protein sample in dilute buffer which is a very large extrapolation of the in vivo conditions, for the crowding effect on protein stability is ignored.

A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study

We developed a method for implicit solvent molecular dynamics simulations of proteins in solvent mixtures (model with implicit solvation thermodynamics, MIST). The MIST method introduces experimental group transfer free energies to the generalized Born formulation for generating molecular trajectories without the need for developing rigorous explicit-solvent force fields for multicomponent solutions. As a test case, we studied the urea-induced denaturation of the Trp-cage miniprotein in water. We demonstrate that our method allows efficient exploration of the conformational space of the protein in only a few hundreds of nanoseconds of all-atom unbiased simulations. Furthermore, selective implementation of the transfer free energies of specific peptide groups, backbone, and side chains enables us to decouple their specific energetic contributions to the conformational changes of the protein. The approach herein developed can readily be extended to the investigation of complex matrices as well as to the characterization of protein aggregation. The MIST method is implemented in Plumed (ver. 2.8) as a separate module called SASA.

Quantification of Entropic Excluded Volume Effects Driving Crowding-Induced Collapse and Folding of a Disordered Protein

We investigate the conformational properties of the intrinsically disordered DNA-binding domain of CytR in the presence of the polymeric crowder polyethylene glycol (PEG). Integrating circular dichroism, nuclear magnetic resonance, and single-molecule Förster resonance energy transfer measurements, we demonstrate that disordered CytR populates a well-folded minor conformation in its native ensemble, while the unfolded ensemble collapses and folds with an increase in crowder density independent of the crowder size. Employing a statistical-mechanical model, the effective reduction in the accessible conformational space of a residue in the unfolded state is estimated to be 10% at 300 mg/mL PEG8000, relative to dilute conditions. The experimentally consistent PEG-temperature phase diagram thus constructed reveals that entropic effects can stabilize disordered CytR by 10 kJ mol^-1, driving the equilibrium toward folded conformations under physiological conditions. Our work highlights the malleable conformational landscape of CytR, the presence of a folded conformation in the disordered ensemble, and proposes a scaling relation for quantifying excluded volume effects on protein stability.

Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a relaxed state at pH 4-8 (with the typical alignment of the crystallographic structure) and a compact state at pH 3-4 (with a tighter association and rotated alignment). This work shows that an approach based on constant-pH MD simulations can produce rich detailed pictures of pH-dependent protein associations, as illustrated for BLG dimerization.

The pH Effects on SARS-CoV and SARS-CoV-2 Spike Proteins in the Process of Binding to hACE2

COVID-19 has been threatening human health since the late 2019, and has a significant impact on human health and economy. Understanding SARS-CoV-2 and other coronaviruses is important to develop effective treatments for COVID-19 and other coronavirus-caused diseases. In this work, we applied multi-scale computational approaches to study the electrostatic features of spike (S) proteins for SARS-CoV and SARS-CoV-2. From our results, we found that SARS-CoV and SARS-CoV-2 have similar charge distributions and electrostatic features when binding with the human angiotensin-converting enzyme 2 (hACE2). Energy pH-dependence calculations revealed that the complex structures of hACE2 and the S proteins of SARS-CoV/SARS-CoV-2 are stable at pH values ranging from 7.5 to 9. Three independent 100 ns molecular dynamics (MD) simulations were performed using NAMD to investigate the hydrogen bonds between S proteins RBD and hACE2 RBD. From MD simulations, we found that SARS-CoV-2 forms 19 pairs (average of three simulations) of hydrogen bonds with high occupancy (>50%) to hACE2, compared to 16 pairs between SARS-CoV and hACE2. Additionally, SARS-CoV viruses prefer sticking to the same hydrogen bond pairs, while SARS-CoV-2 tends to have a larger range of selections on hydrogen bonds acceptors. We also labelled key residues involved in forming the top five hydrogen bonds that were found in all three independent 100 ns simulations. This identification is important to potential drug designs for COVID-19 treatments. Our work will shed the light on current and future coronavirus-caused diseases.

Classification of protein solubilizing additives by fluorescence assay

Aromatic interaction plays a crucial role in controlling protein interaction by additives. Here we investigated the interaction of protein salting-in (solubilizing) additives with tryptophan (Trp), tyrosine (Tyr), indole, and proteins based on their fluorescence spectra. Five salting-in additives, i.e., arginine (Arg), urea, guanidine (Gdn), ethylene glycol (EG), and magnesium chloride (MgCl₂), showed different effects on the fluorescence properties of Trp and Tyr. Arg significantly reduced fluorescence intensity of Trp and Tyr, as was the case for glycine to a lesser extent. MgCl₂ and calcium chloride (CaCl₂) showed little effect on the aromatic fluorescence spectra. Gdn also showed little effect on the aromatic fluorescence spectra of Trp and Tyr even at high concentrations. EG increased the aromatic fluorescence intensity of Trp and Tyr with blue-shifted emission wavelength. Urea enhanced fluorescence of Trp and Tyr without altering emission wavelength. These results indicate that the protein solubilizing additives interact with aromatic groups differently.

Immobilization of horseradish peroxidase on lysine-functionalized gum Arabic-coated Fe3O4 nanoparticles for cholesterol determination

BACKGROUND

Horseradish Peroxidase (HRP) is ranked as one of the most important industrial enzymes that is extensively used in industry. Cholesterol is routinely detected indirectly by cholesterol oxidase in the presence of O₂, liberating H₂O₂ as a by-product. The H₂O₂ content is determined through the HRP activity in the presence of a redox dye, producing a red colored quinoneimine which can be measured quantitatively. Herein, we have designed a magnetic nanoparticle for reusing and easily separating HRP as the most expensive compartment for the low-cost cholesterol assay.

METHODS

The gum Arabic coated magnetic nanoparticles were functionalized with L-lysine linker for maintaining protein flexibility on nanoparticle. Enzyme-loaded nanoparticles were characterized by TEM, FTIR, DLS, VSM and XRD analysis.

RESULTS

The immobilization efficiency was ∼65% and the immobilized HRP retained 60% of its activity after 8 times reuse. The optimum pH and thermal stability shifted from 7.0 to 8.0 and 60 to 70 °C after immobilization, respectively. Storage stability of HRP was improved by 10%, at 4 °C for 60 days. Immobilized HRP showed more catalytic activity in presence of Fe²⁺, Ca²⁺ and Na⁺. The designed system has cholesterol detection linearity range from 0.2 to 5.0 mM and detection limit of 0.08 mM and acceptable correlation coefficient of 0.9973 and 0.9982 on sample serum using both chromogens.

CONCLUSION

The HRP-loaded magnetic nanoparticles are capable of being used as a cost-effective system for cholesterol determination in laboratory due to its reusability and stability benefits.

Novel variations in spermatogenic transcription regulators RFX2 and TAF7 increase risk of azoospermia

PURPOSE

Genetic etiology of idiopathic male infertility is enigmatic owing to involvement of multiple gene regulatory networks in spermatogenesis process. Any change in optimal function of the transcription factors involved in this process owing to polymorphisms/mutations may increase the risk of infertility. We investigated polymorphisms/mutations of spermatogenic transcription regulators TAF7 and RFX2 and analysed their association with incidence of azoospermia among the men from West Bengal, India.

METHODS

Genotyping was carried by Sanger's dideoxy sequencing of 130 azoospermic men who were detected negative in Y chromosome microdeletion screening and 140 healthy controls. Association study was done by suitable statistical methods. In silico analysis was performed to infer the intuitive damaging effects of detected variants at transcripts and protein level.

RESULTS

We found significant association of TAF7 C16T (MW827584 G > A), RFX2 562delT (MZ560629delA), rs11547633 A > C, rs17606721 A > G, MW827583 C > T, and MZ379836 C > T variants with the incidence of azoospermia. In silico analysis predicted that the variants either alter the natural splice junctions of the transcript or cause probable damage in the structure of proteins of respective genes.

CONCLUSION

Polymorphisms/mutations of TAF7 and RFX2 genes increase risk of male infertility in Bengali population. The novel variants may be used as markers for male infertility screening in ART practise.

Machine learning implicit solvation for molecular dynamics

Accurate modeling of the solvent environment for biological molecules is crucial for computational biology and drug design. A popular approach to achieve long simulation time scales for large system sizes is to incorporate the effect of the solvent in a mean-field fashion with implicit solvent models. However, a challenge with existing implicit solvent models is that they often lack accuracy or certain physical properties compared to explicit solvent models as the many-body effects of the neglected solvent molecules are difficult to model as a mean field. Here, we leverage machine learning (ML) and multi-scale coarse graining (CG) in order to learn implicit solvent models that can approximate the energetic and thermodynamic properties of a given explicit solvent model with arbitrary accuracy, given enough training data. Following the previous ML-CG models CGnet and CGSchnet, we introduce ISSNet, a graph neural network, to model the implicit solvent potential of mean force. ISSNet can learn from explicit solvent simulation data and be readily applied to molecular dynamics simulations. We compare the solute conformational distributions under different solvation treatments for two peptide systems. The results indicate that ISSNet models can outperform widely used generalized Born and surface area models in reproducing the thermodynamics of small protein systems with respect to explicit solvent. The success of this novel method demonstrates the potential benefit of applying machine learning methods in accurate modeling of solvent effects for in silico research and biomedical applications.

Atomic view of cosolute-induced protein denaturation probed by NMR solvent paramagnetic relaxation enhancement

The cosolvent effect arises from the interaction of cosolute molecules with a protein and alters the equilibrium between native and unfolded states. Denaturants shift the equilibrium toward the latter, while osmolytes stabilize the former. The molecular mechanism whereby cosolutes perturb protein stability is still the subject of considerable debate. Probing the molecular details of the cosolvent effect is experimentally challenging as the interactions are very weak and transient, rendering them invisible to most conventional biophysical techniques. Here, we probe cosolute-protein interactions by means of NMR solvent paramagnetic relaxation enhancement together with a formalism we recently developed to quantitatively describe, at atomic resolution, the energetics and dynamics of cosolute-protein interactions in terms of a concentration normalized equilibrium average of the interspin distance, [Formula: see text], and an effective correlation time, τ_c The system studied is the metastable drkN SH3 domain, which exists in dynamic equilibrium between native and unfolded states, thereby permitting us to probe the interactions of cosolutes with both states simultaneously under the same conditions. Two paramagnetic cosolute denaturants were investigated, one neutral and the other negatively charged, differing in the presence of a carboxyamide group versus a carboxylate. Our results demonstrate that attractive cosolute-protein backbone interactions occur largely in the unfolded state and some loop regions in the native state, electrostatic interactions reduce the [Formula: see text] values, and temperature predominantly impacts interactions with the unfolded state. Thus, destabilization of the native state in this instance arises predominantly as a consequence of interactions of the cosolutes with the unfolded state.

The Right-Handed Parallel β-Helix Topology of Erwinia chrysanthemi Pectin Methylesterase Is Intimately Associated with Both Sequential Folding and Resistance to High Pressure

The complex topologies of large multi-domain globular proteins make the study of their folding and assembly particularly demanding. It is often characterized by complex kinetics and undesired side reactions, such as aggregation. The structural simplicity of tandem-repeat proteins, which are characterized by the repetition of a basic structural motif and are stabilized exclusively by sequentially localized contacts, has provided opportunities for dissecting their folding landscapes. In this study, we focus on the Erwinia chrysanthemi pectin methylesterase (342 residues), an all-β pectinolytic enzyme with a right-handed parallel β-helix structure. Chemicals and pressure were chosen as denaturants and a variety of optical techniques were used in conjunction with stopped-flow equipment to investigate the folding mechanism of the enzyme at 25 °C. Under equilibrium conditions, both chemical- and pressure-induced unfolding show two-state transitions, with average conformational stability (ΔG° = 35 ± 5 kJ·mol⁻¹) but exceptionally high resistance to pressure (P_m = 800 ± 7 MPa). Stopped-flow kinetic experiments revealed a very rapid (τ < 1 ms) hydrophobic collapse accompanied by the formation of an extended secondary structure but did not reveal stable tertiary contacts. This is followed by three distinct cooperative phases and the significant population of two intermediate species. The kinetics followed by intrinsic fluorescence shows a lag phase, strongly indicating that these intermediates are productive species on a sequential folding pathway, for which we propose a plausible model. These combined data demonstrate that even a large repeat protein can fold in a highly cooperative manner.

Multiple-State Monitoring of SOD1 Amyloid Formation at Single-Residue Resolution by Rheo-NMR Spectroscopy

Formation of protein aggregates or fibrils entails the conversion of soluble native protein monomers via multiple molecular states. No spectroscopic techniques have succeeded in capturing the transient molecular-scale events of fibrillation in situ. Here we report residue- and state-specific real-time monitoring of the fibrillation of amyotrophic lateral sclerosis-related SOD1 by rheology NMR (Rheo-NMR) spectroscopy. Under moderately denaturing conditions, where NMR signals of folded and unfolded monomeric SOD1 are simultaneously observable, the cross-peak intensities of folded monomeric SOD1 decreased faster than those of the unfolded species, and a 3₁₀-helix in folded SOD1 was deformed prior to global unfolding. Furthermore, real-time protein dynamics analysis identified residues involved in the core structure formation of SOD1 oligomers. Our findings provide insight into local and global unfolding events in SOD1 and fibril formation. This Rheo-NMR analysis will be applicable not only to atomic-level monitoring of other amyloidogenic proteins but also to quantification of shear-induced structural changes of non-amyloidogenic proteins and elucidation of shear-enhanced chemical phenomena such as viscosity increase and crystallization of various solution-state compounds.

A Glimpse into the Structural Properties of the Intermediate and Transition State in the Folding of Bromodomain 2 Domain 2 by Φ Value Analysis

Bromodomains (BRDs) are small protein interaction modules of about 110 amino acids that selectively recognize acetylated lysine in histones and other proteins. These domains have been identified in a variety of multi-domain proteins involved in transcriptional regulation or chromatin remodeling in eukaryotic cells. BRD inhibition is considered an attractive therapeutic approach in epigenetic disorders, particularly in oncology. Here, we present a Φ value analysis to investigate the folding pathway of the second domain of BRD2 (BRD2(2)). Using an extensive mutational analysis based on 25 site-directed mutants, we provide structural information on both the intermediate and late transition state of BRD2(2). The data reveal that the C-terminal region represents part of the initial folding nucleus, while the N-terminal region of the domain consolidates its structure only later in the folding process. Furthermore, only a small number of native-like interactions have been identified, suggesting the presence of a non-compact, partially folded state with scarce native-like characteristics. Taken together, these results indicate that, in BRD2(2), a hierarchical mechanism of protein folding can be described with non-native interactions that play a significant role in folding.

A pH-dependent protein stability switch coupled to the perturbed pKa of a single ionizable residue

The contribution of electrostatic interactions in protein stability has not been fully understood. Burial of an ionizable amino acid inside the hydrophobic protein core can affect its ionization equilibrium and shift its pKa differentially in the native (N) and unfolded (U) states of a protein and this coupling between the folding/unfolding cycle and the ionization equilibria of the ionizable residue can substantially influence the protein stability. Here, we studied the coupling of the folding/unfolding cycle with the ionization of a buried ionizable residue in a multi-domain protein, Human Serum Albumin (HSA) using fluorescence spectroscopy. A pH-dependent change in the stability of HSA was observed in the near native pH range (pH 6.0-9.0). The protonation-deprotonation equilibrium of a single thiol residue that is buried in the protein structure was identified to give rise to the pH-dependent protein stability. We quantified the pKa of the thiol residue in the N and the U states. The mean pKa of the thiol in the N state was upshifted by 0.5 units to 8.7 due to the burial of the thiol in the protein structure. Surprisingly, the mean pKa of the thiol in the U state was observed to be downshifted by 1.3 units to 6.9. These results indicate that some charged residues are spatially proximal to the thiol group in the U state. Our results suggest that, in addition to the N state, electrostatic interactions in the U state are important determinants of protein stability.

Quantitative Analysis of Protein Unfolded State Energetics: Experimental and Computational Studies Demonstrate That Non-Native Side-Chain Interactions Stabilize Local Native Backbone Structure

Proteins fold on relatively smooth free energy landscapes which are biased toward the native state, but even simple topologies which fold rapidly can experience roughness on their free energy landscape. The details of these interactions are difficult to elucidate experimentally. Closely related to the problem of deciphering the details of the free energy landscape is the problem of defining the interactions in the denatured state ensemble (DSE) which is populated under native conditions, that is, under conditions where the native state is stable. The DSE of many proteins deviates from random coil models, but quantifying and defining the energetics of the transiently populated interactions in this ensemble is extremely challenging. Characterization of the DSE of proteins which fold to compact structures is also relevant to studies of intrinsically disordered proteins (IDPs) since interactions in the dynamic ensemble populated by IDPs can modulate their behavior. Here we show how experimental thermodynamic and pK_a measurements can be combined with computational thermodynamic integration to quantify interactions in the DSE. We show that non-native side chain interactions can stabilize native backbone structure in the DSE and demonstrate that that even rapidly folding proteins can form energetically significant non-native interactions in their DSE. As an example, we characterize a non-native salt bridge that stabilizes local native backbone structure in the DSE of a widely studied fast-folding protein, the villin headpiece helical domain. The combined computational experimental approach is applicable to other protein unfolded states and provides insight that is impossible to obtain with either method alone.

Protein-complex stability in cells and in vitro under crowded conditions

Biology is beginning to appreciate the effects of the crowded and complex intracellular environment on the equilibrium thermodynamics and kinetics of protein folding. The next logical step involves the interactions between proteins. We review quantitative, wet-experiment based efforts aimed at understanding how and why high concentrations of small molecules, synthetic polymers, biologically relevant cosolutes and the interior of living cells affect the energetics of protein-protein interactions. We then address popular theories used to explain the effects and suggest expeditious paths for a more methodical integration of experiment and simulation.

MIRRAGGE - Minimum Information Required for Reproducible AGGregation Experiments

Reports on phase separation and amyloid formation for multiple proteins and aggregation-prone peptides are recurrently used to explore the molecular mechanisms associated with several human diseases. The information conveyed by these reports can be used directly in translational investigation, e.g., for the design of better drug screening strategies, or be compiled in databases for benchmarking novel aggregation-predicting algorithms. Given that minute protocol variations determine different outcomes of protein aggregation assays, there is a strong urge for standardized descriptions of the different types of aggregates and the detailed methods used in their production. In an attempt to address this need, we assembled the Minimum Information Required for Reproducible Aggregation Experiments (MIRRAGGE) guidelines, considering first-principles and the established literature on protein self-assembly and aggregation. This consensus information aims to cover the major and subtle determinants of experimental reproducibility while avoiding excessive technical details that are of limited practical interest for non-specialized users. The MIRRAGGE table (template available in Supplementary Information) is useful as a guide for the design of new studies and as a checklist during submission of experimental reports for publication. Full disclosure of relevant information also enables other researchers to reproduce results correctly and facilitates systematic data deposition into curated databases.

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23-34) and an N-terminal β-hairpin (residues 2-16) were significantly protected with a protection factor of 2.1-4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.

Effects of natural mutations (L94I and L94V) on the stability and mechanism of folding of horse cytochrome c: A combined in vitro and molecular dynamics simulations approach

Known crystal structures of 10 cytochromes (cyts) c from different sources led to the conclusion that natural mutations in these proteins does not affect their 3D structure, hence evolution preserved structure for function. A sequence alignment of horse cyt c with all other 284 cyts c led to two important conclusions: (i) Leu at position 94 is conserved in all 30 mammalian known sequences, and (ii) there are 14 other species which have either Val or Ile at 94th position. We asked a question: Is the avoidance of substitution by Val or Ile at position 94 in the mammalian cyts c by design or by chance? To answer this question, we introduced natural substitutes of Leu94 by Val and Ile in horse cyt c using site-directed mutagenesis. Here, from our in vitro and molecular dynamic simulation studies on L94V and L94I mutants, we concluded that (i) although the natural mutations destabilize the wild type cyt c, it does not significantly affect the mechanism of folding of the protein, (ii) urea-induced denaturation of WT cyt c and its mutants is a two-state process, and (iii) denaturation of WT cyt c and its mutants by guanidinium chloride is not a two-state process.

Duplex-tetraplex equilibria in guanine- and cytosine-rich DNA

Noncanonical four-stranded DNA structures, including G-quadruplexes and i-motifs, have been discovered in the cell and are implicated in a variety of genomic regulatory functions. The tendency of a specific guanine- and cytosine-rich region of genomic DNA to adopt a four-stranded conformation depends on its ability to overcome the constraints of duplex base-pairing by undergoing consecutive duplex-to-coil and coil-to-tetraplex transitions. The latter ability is determined by the balance between the free energies of participating ordered and disordered structures. In this review, we present an overview of the literature on the stability of G-quadruplex and i-motif structures and discuss the extent of duplex-tetraplex competition as a function of the sequence context of the DNA and environmental conditions including temperature, pH, salt, molecular crowding, and the presence of G-quadruplex-binding ligands. We outline how the results of in vitro studies can be expanded to understanding duplex-tetraplex equilibria in vivo.

Global analysis of protein stability by temperature and chemical denaturation

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 μl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 μg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.