Search

Markov State Models (MSMs) have been widely applied to understand protein folding mechanisms by predicting long time scale dynamics from ensembles of short molecular simulations. Most MSM estimators enforce detailed balance, assuming that trajectory data are sampled at an equilibrium. This is rarely the case for ab initio folding studies, however, and as a result, MSMs can severely underestimate protein folding stabilities from such data. To remedy this problem, we have developed an enhanced-sampling protocol in which (1) unbiased folding simulations are performed and sparse tICA is used to obtain features that best capture the slowest events in folding, (2) umbrella sampling along this reaction coordinate is performed to observe folding and unfolding transitions, and (3) the thermodynamics and kinetics of folding are estimated using multiensemble Markov models (MEMMs). Using this protocol, folding pathways, rates, and stabilities of a designed α-helical hairpin, Z34C, can be predicted in good agreement with experimental measurements. These results indicate that accurate simulation-based estimates of absolute folding stabilities are within reach, with implications for the computational design of folded miniproteins and peptidomimetics.

ConspectusInsulin has remained indispensable in the treatment of diabetes since it was first discovered in 1921. Unlike small molecular drugs, insulin and other protein drugs are prone to degradation when exposed to elevated temperatures, mechanical agitation during transportation, and prolonged storage periods. Therefore, strict cold-chain management is crucial for the insulin supply, requiring significant resources, which can limit the access to insulin, particularly in low-income areas. Moreover, although insulin formulations have advanced tremendously in the last century, insulin treatment still imposes a challenging regimen and provides suboptimal outcomes for the majority of patients. There is an increasing focus on pursuing improved pharmacology, specifically on safer, more user-friendly insulin therapies that minimize the self-management burden. These challenges underscore the need for developing novel insulin formulations with improved stability that are compatible with advanced insulin therapy.Insulin stabilization can be achieved through either chemical modification of insulin or formulation component design. Inspired by insulin-like peptides from invertebrates, we have developed novel stable insulin analogs based on a fundamental understanding of the insulin receptor engagement for insulin bioactivity. We created a novel four-disulfide insulin analog with high aggregation stability and potency by introducing a fourth disulfide bond between a C-terminal extended insulin A-chain and residues near the C-terminus of the B-chain. In an effort to stabilize insulin in its monomeric state to develop ultrafast-acting insulin with rapid absorption upon injection, we have developed a series of structurally miniaturized yet fully active insulin analogs that do not form dimers due to the lack of the canonical B-chain C-terminal octapeptide. Additionally, our study provided strategies for expanding the scope of cucurbit[7]uril (CB[7])-assisted insulin stabilization by engineering safe and biodegradable CB[7]-zwitterionic polypeptide excipients. We also explored insulin N-terminal substitution methods to achieve pH-dependent insulin stabilization without prolonging the duration of action.This Account describes our exploration of engineering stable insulin analogs and formulation design strategies for stabilizing insulin in aqueous solutions. Beyond conventional stabilization strategies for insulin injections, the unmet challenges and recent innovations in insulin stabilization are discussed, addressing the growing demand for alternative, less invasive routes of insulin administration. Additionally, we aim to provide a thorough overview of insulin stabilization from the perspective of commercially available insulin drugs and common pharmaceutical engineering practices in the industry. We also highlight unresolved insulin stabilization challenges and ongoing research strategies. We anticipate that further emphasis on collective efforts of protein engineering, pharmaceutical formulation design, and drug delivery will inform the development of stable and advanced insulin therapy.

Ammonium acetate is widely used in native mass spectrometry to provide adequate ionic strength without adducting to protein ions, but different ions can preferentially stabilize or destabilize the native form of proteins in solution. The stability of bovine serum albumin (BSA) was investigated in 50 mM solutions of a variety of salts using electrospray emitters with submicron tips to desalt protein ions. The charge-state distribution of BSA is narrow (+14 to +18) in ammonium acetate (AmmAc), whereas it is much broader (+13 to +42) in solutions containing sodium acetate (NaAc), ammonium chloride (AmmCl), potassium chloride (KCl), and sodium chloride (NaCl). The average charge state and percent of unfolded protein increase in these respective solutions, indicating greater extents of protein destabilization and conformational changes. In contrast, no high charge states of either bovine carbonic anhydrase II or IgG1 were formed in AmmAc or NaCl despite their similar melting temperatures to BSA, indicating that the presence of unfolded BSA in some of these solutions is not an artifact of the electrospray ionization process. The charge states formed from the nonvolatile salt solutions do not change significantly for up to 7 min of electrospray, but higher charging occurs after 10 min, consistent with solution acidification. Formation of unfolded BSA in NaAc but not in AmmAc indicates that the cation identity, not acidification, is responsible for structural differences in these two solutions. Temperature-dependent measurements show both increased charging and aggregation at lower temperatures in NaCl:Tris than in AmmAc, consistent with lower protein stability in the former solution. These results are consistent with the order of these ions in the Hofmeister series and indicate that data on protein stability in AmmAc may not be representative of solutions containing nonvolatile salts that are directly relevant to biology.

A subunit of factor XIII (FXIII-A) contains a unique activation peptide (AP) that protects the catalytic triad and prevents degradation. In plasma, FXIII is activated proteolytically (FXIII-A*) by thrombin and Ca 2+ cleaving AP, while in cytoplasm, it is activated nonproteolytically (FXIII-A°) with increased Ca 2+ concentrations. This study aimed to elucidate the role of individual parts of the FXIII-A AP in protein stability, thrombin activation, and transglutaminase activity. Recombinant FXIII-A AP variants were expressed, and SDS-PAGE was used to monitor thrombin hydrolysis at the AP cleavage sites R37-G38. Transglutaminase activities were assessed by cross-linking lysine mimics to Fbg αC (233-425, glutamine-substrate) and monitoring reactions by mass spectrometry and in-gel fluorescence assays. FXIII-A AP variants, S19P, E23K, and D24V, degraded during purification, indicating their vital role in FXIII-A 2 stability. Mutation of P36 to L36/F36 abolished the proteolytic cleavage of AP and thus prevented activation. FXIII-A N20S and P27L exhibited slower thrombin activation, likely due to the loss of key interdomain H-bonding interactions. Except N20S and P15L/P16L, all activatable FXIII-A* variants (P15L, P16L, S19A, and P27L) showed similar cross-linking activity to WT. By contrast, FXIII-A° P15L, P16L, and P15L/P16L had significantly lower cross-linking activity than FXIII-A° WT, suggesting that loss of these prolines had a greater structural impact. In conclusion, FXIII-A AP residues that play crucial roles in FXIII-A stability, activation, and activity were identified. The interactions between these AP amino acid residues and other domains control the stability and activity of FXIII.

The design of stable formulations remains a major challenge for protein therapeutics, particularly the need to minimize aggregation. Experimental formulation screens are typically based on thermal transition midpoints ( T m ), and forced degradation studies at elevated temperatures. Both approaches give limited predictions of long-term storage stability, particularly at low temperatures. Better understanding of the mechanisms of action for formulation of excipients and buffers could lead to improved strategies for formulation design. Here, we identified a complex impact of glycine concentration on the experimentally determined stability of an antibody Fab fragment and then used molecular dynamics simulations to reveal mechanisms that underpin these complex behaviors. T m values increased monotonically with glycine concentration, but associated Δ S vh measurements revealed more complex changes in the native ensemble dynamics, which reached a maximum at 30 mg/mL. The aggregation kinetics at 65 °C were similar at 0 and 20 mg/mL glycine, but then significantly slower at 50 mg/mL. These complex behaviors indicated changes in the dominant stabilizing mechanisms as the glycine concentration was increased. MD revealed a complex balance of glycine self-interaction, and differentially preferred interactions of glycine with the Fab as it displaced hydration-shell water, and surface-bound water and citrate buffer molecules. As a result, glycine binding to the Fab surface had different effects at different concentrations, and led from preferential interactions at low concentrations to preferential exclusion at higher concentrations. During preferential interaction, glycine displaced water from the Fab hydration shell, and a small number of water and citrate molecules from the Fab surface, which reduced the protein dynamics as measured by root-mean-square fluctuation (RMSF) on the short time scales of MD. By contrast, the native ensemble dynamics increased according to Δ S vh , suggesting increased conformational changes on longer time scales. The aggregation kinetics did not change at low glycine concentrations, and so the opposing dynamics effects either canceled out or were not directly relevant to aggregation. During preferential exclusion at higher glycine concentrations, glycine could only bind to the Fab surface through the displacement of citrate buffer molecules already favorably bound on the Fab surface. Displacement of citrate increased the flexibility (RMSF) of the Fab, as glycine formed fewer bridging hydrogen bonds to the Fab surface. Overall, the slowing of aggregation kinetics coincided with reduced flexibility in the Fab ensemble at the very highest glycine concentrations, as determined by both RMSF and Δ S vh , and occurred at a point where glycine binding displaced neither water nor citrate. These final interactions with the Fab surface were driven by mass action and were the least favorable, leading to a macromolecular crowding effect under the regime of preferential exclusion that stabilized the dynamics of Fab.

Computational methods including machine learning and molecular dynamics simulations have strong potential to characterize, understand, and ultimately predict the properties of proteins relevant to their stability and function as therapeutics. Such methods would streamline the development pathway by minimizing the current experimental testing required for many protein variants and formulations. The molecular understanding of thermostability and aggregation propensity has advanced significantly along with predictive algorithms based on the sequence-level or structural-level information on a protein. However, these approaches focus largely on a comparison of protein sequence variations to correlate the properties of proteins to their stability, solubility, and aggregation propensity. For therapeutic protein development, it is of equal importance to take into account the impact of the formulation conditions to elucidate and predict the stability of the antibody drugs. At the macroscopic level, changing temperature, pH, ionic strength, and the addition of excipients can significantly alter the kinetics of protein aggregation. The mechanisms controlling aggregation kinetics have been traced back to a combination of molecular features, including conformational stability, partial unfolding to aggregation-prone states, and the colloidal stability governed by surface charges and hydrophobicity. However, very little has been done to evaluate these features in the context of protein dynamics in different formulations. In this work, we have combined a range of molecular features calculated from the Fab A33 protein sequence and molecular dynamics simulations. Using the power of advanced, yet interpretable, statistical tools, it has been possible to uncover greater insights into the mechanisms behind protein stability, validating previous findings, and also develop models that can predict the aggregation kinetics within a range of 49 different solution conditions.

Thermostability, which is essential for the functional performance of enzymes, is largely determined by intramolecular physical interactions. Although many tools have been developed, existing computational methods have struggled to find the universal principles of protein thermostability. Recent advancements in structural proteomics have been driven by the introduction of deep neural networks such as AlphaFold2 and ESMFold. These innovations have enabled the characterization of protein structures with unprecedented speed and accuracy. Here, we introduce qProtein, a Python-implemented workflow designed for the quantitative analysis of physical interactions on the scale of structural proteomics. This platform accepts protein sequences as input and produces four structural features, including hydrophobic clusters, hydrogen bonds, electrostatic interactions, and disulfide bonds. To demonstrate the use of qProtein, we investigate the structural features related to protein thermostability in six glycoside hydrolase (GH) families, comprising a total of 3,811 protein structures. Our results indicate that in five enzyme families (GH11, GH12, GH5_2, GH10, and GH48), the thermophilic enzymes have a larger average area of hydrophobic clusters compared to the nonthermophilic enzymes within each family. Furthermore, our analysis of the local-structure regions reveals that the hydrophobic clusters are predominantly distributed in the distal regions of the GH11 enzymes. In addition, the average hydrophobic cluster area of the thermophilic enzymes is significantly higher than that of the nonthermophilic enzymes in the distal regions of the GH11 enzymes. Therefore, qProtein is a well-suited platform for analyzing the structural features of thermal stability at the level of structural proteomics. We provide the source code for qProtein at https://github.com/bj600800/qProtein, and the web server is available at http://qProtein.sdu.edu.cn:8888.

We study the influence of urea on the stability of monoclonal antibodies (mAbs) using molecular dynamics (MD) simulations in combination with differential scanning fluorimetry (DSF). We show that a denaturing cosolute such as urea binds strongly to the protein, which can lead to denaturation and enhanced aggregation behavior at high temperatures. The interaction between protein and urea crucially depends on the surface properties of the individual mAb domains and therefore affects the general binding to the protein differently. The study of these mechanisms for proteins with multiple domains, such as mAbs, encounters significant limitations in experimental analysis methods due to their complexity. Using computational and experimental methods, we are able to separate the protein-urea interaction by domain and show that Lennard-Jones interactions are mainly responsible for significant binding effects. Our results emphasize the potential of MD simulations in combination with Kirkwood-Buff theory to study the interactions between proteins with multiple domains and cosolutes as formulation excipients for drug discovery and development.

Since the proposition of the Hofmeister series, guanidinium (Gdm) salts hold a special mention in protein science owing to their contrasting effect on protein(s) depending on the counteranion(s). For example, while GdmCl is known to act as a potential protein denaturant, Gdm 2 SO 4 offers minimal effect on protein structure. Despite the fact that theoretical studies reckon the formation of ion-pairing to be responsible for such behavior, experimental validation of this hypothesis is still in sparse. In this study, we combine electrochemical impedance spectroscopy (EIS) and THz spectroscopy to underline the effect of GdmCl and Gdm 2 SO 4 on a model amide molecule N -methylacetamide (NMA). Molecular dynamics (MD) simulation studies predict that Gdm 2 SO 4 forms heteroion pairing in water, which inhibits Gdm + ions to approach NMA molecules, while in case of GdmCl, Gdm + ions directly interact with NMA. The experimental findings on ion hydration, specifically the detailed analysis of the ion-water rattling mode, which appears in the THz frequency domain, unambiguously endorse this hypothesis. Our study establishes the fact that the propensity of ion-pairing in Gdm salts dictates their (de)stabilization effect on proteins.

The aberrant accumulation of the transactive response deoxyribonucleic acid (DNA)-binding protein of 43 kDa (TDP-43) aggregates in the cytoplasm of motor neurons is the main pathological hallmark of amyotrophic lateral sclerosis (ALS). Previous experiments reported that adenosine triphosphate (ATP), the universal energy currency for all living cells, could induce aggregation and enhance the folding of TDP-43 fibrillar aggregates. However, the significance of ATP on TDP-43 fibrillation and the mechanism behind it remain elusive. In this work, we conducted multiple atomistic molecular dynamics (MD) simulations totaling 20 μs to search the critical nucleus size of TDP-43 282-360 and investigate the impact of ATP molecules on preformed protofibrils. The results reveal that the trimer is the critical nucleus for TDP-43 282-360 fibril formation and the tetramer is the minimal stable nucleus. When ATP molecules bind to the TDP-43 282-360 trimer and tetramer, they can consolidate the TDP-43 282-360 protofibrils by increasing the content of the β-sheet structure and promoting the formation of hydrogen bonds (H-bonds). Binding site analyses show that the N-terminus of TDP-43 282-360 protofibrils is the main binding site of ATP, and R293 dominates the direct binding of ATP. Further analyses reveal that the π-π, cation-π, salt bridge, and H-bonding interactions together contribute to the binding of ATP to TDP-43 282-360 protofibrils. This study decoded the detailed stabilization mechanism of protofibrillar TDP-43 282-360 oligomers by ATP, and may provide new avenues for the development of drug design against ALS.

The stability of proteins from rates of oxidation (SPROX), thermal proteome profiling (TPP), and limited proteolysis (LiP) techniques were used to profile the stability of ∼2500 proteins in hippocampus tissue cell lysates from 2- and 8-months-old wild-type (C57BL/6J; n = 7) and transgenic (5XFAD; n = 7) mice with five Alzheimer's disease (AD)-linked mutations. Approximately 200-500 protein hits with AD-related stability changes were detected by each technique at each age point. The hit overlap from technique to technique was low, and all of the techniques generated protein hits that were more numerous and largely different from those identified in protein expression level analyses, which were also performed here. The hit proteins identified by each technique were enriched in a number of the same pathways and biological processes, many with known connections to AD. The protein stability hits included 25 high-value conformation biomarkers with AD-related stability changes detected using at least 2 techniques at both age points. Also discovered were subunit- and age-specific AD-related stability changes in the proteasome, which had reduced function at both age points. The different folding stability profiles of the proteasome at the two age points are consistent with a different mechanism for proteasome dysfunction at the early and late stages of AD.

α-N-Methylation (Nα-methylation), catalyzed by protein N-terminal methyltransferases (NTMTs), constitutes a crucial post-translational modification involving the transfer of a methyl group from S -adenosyl-l-methionine (SAM) to the Nα-terminal amino group of substrate proteins. NTMT1/2 are known to methylate canonical Nα sequences, such as X-P-K/R. With over 300 potential human protein substrates, only a small fraction has been validated, and even less is known about the functions of Nα-methylation. This study delves into the characterizations of protein arginine deiminase 1 (PAD1) as a substrate of NTMT1. By employing biochemical and cellular assays, we demonstrated NTMT1-mediated Nα-methylation of PAD1, leading to an increase in protein half-life and the modulation of protein-protein interactions in HEK293T cells. The methylation of PAD1 appears nonessential to its enzymatic activity or cellular localization. Proteomic studies revealed differential protein interactions between unmethylated and Nα-methylated PAD1, suggesting a regulatory role for Nα-methylation in modulating PAD1's protein-protein interactions. These findings shed light on the intricate molecular mechanisms governing PAD1 function and expand our knowledge of Nα-methylation in regulating protein function.

Exenatide (Ex4), a GLP-1 incretin mimetic polypeptide, is an effective therapeutic agent against diabetes and obesity. We highlight the indirect role of Ex4's structure-stabilizing Trp-cage (Tc) motif in governing GLP-1 receptor (GLP-1R) signal transduction. We use various Ex4 derivatives to explore how Tc compactness influences thermal stability, aggregation, enhancement of insulin secretion, and GLP-1R binding. We found that Ex4 variants decorated with fortified Tc motifs exhibit increased resistance to unfolding and aggregation but show an inverse relationship between the bioactivity and stability. Molecular dynamics simulations coupled with a rigid-body segmentation protocol to analyze dynamic interconnectedness revealed that the constrained Tc motifs remain intact within the receptor-ligand complexes but interfere with one of the major stabilizing contacts and recognition loci on the extracellular side of GLP-1R, dislodging the N-terminal activating region of the hormone mimetics, and restrict the free movement of TM6, the main signal transduction device of GLP-1R.

Macromolecular crowding bridges in vivo and in vitro studies by simulating cellular complexities such as high viscosity and limited space while maintaining the experimental feasibility. Over the last two decades, the impact of macromolecular crowding on protein stability and activity has been a significant topic of study and discussion, though still lacking a thorough mechanistic understanding. This article investigates the role of associated water dynamics on protein stability and activity within crowded environments, using bromelain and Ficoll-70 as the model systems. Traditional crowding theory primarily attributes protein stability to entropic effects (excluded volume) and enthalpic interactions. However, our recent findings suggest that water structure modulation plays a crucial role in a crowded environment. In this report, we strengthen the conclusion of our previous study, i.e., rigid-associated water stabilizes proteins via entropy and destabilizes them via enthalpy, while flexible water has the opposite effect. In the process, we addressed previous shortcomings with a systematic concentration-dependent study using a single-domain protein and component analysis of solvation dynamics. More importantly, we analyze bromelain's hydrolytic activity using the Michaelis-Menten model to understand kinetic parameters like maximum velocity ( V max ) achieved by the system and the Michaelis-Menten coefficient ( K ) complex formation, where an increase in the microviscosity makes the ES complex formation less favorable. On the other hand, flexible associated water dynamics were found to favor the rate of product formation significantly from the ES complex, while rigid associated water hinders it. This study improves our understanding of protein stability and activity in crowded environments, highlighting the critical role of associated water dynamics.M ). Results indicate that microviscosity (not the bulk viscosity) controls the enzyme-substrate ( ES ) complex formation, where an increase in the microviscosity makes the ES complex formation less favorable. On the other hand, flexible associated water dynamics were found to favor the rate of product formation significantly from the ES complex, while rigid associated water hinders it. This study improves our understanding of protein stability and activity in crowded environments, highlighting the critical role of associated water dynamics.

Virus-like particles (VLPs) from bacteriophage MS2 provide a platform to study protein self-assembly and create engineered systems for drug delivery. Here, we aim to understand the impact of intersubunit interface mutations on the local and global structure and function of MS2-based VLPs. In previous work, our lab identified locally supercharged double mutants [T71K/G73R] that concentrate positive charge at capsid pores, enhancing uptake into mammalian cells. To study the effects of particle size on cellular internalization, we combined these double mutants with a single point mutation [S37P] that was previously reported to switch particle geometry from T = 3 to T = 1 icosahedral symmetry. These new variants retained their enhanced cellular uptake activity and could deliver small-molecule drugs with efficacy levels similar to our first-generation capsids. Surprisingly, these engineered triple mutants exhibit increased thermostability and unexpected geometry, producing T = 3 particles instead of the anticipated T = 1 assemblies. Transmission electron microscopy revealed various capsid assembly states, including wild-type ( T = 3), T = 1, and rod-like particles, that could be accessed using different combinations of these point mutations. Molecular dynamics experiments recapitulated the structural rationale in silico for the single point mutation [S37P] forming a T = 1 virus-like particle and showed that this assembly state was not favored when combined with mutations that favor rod-like architectures. Through this work, we investigated how interdimer interface dynamics influence VLP size and morphology and how these properties affect particle function in applications such as drug delivery.

The use of protein biomarkers in blood for clinical settings is limited by the cost and accessibility of traditional venipuncture sampling. The dried blood spot (DBS) technique offers a less invasive and more accessible alternative. However, protein stability in DBS has not been well evaluated. Herein, we deployed a quantitative LC-MS/MS system to construct proteomic atlases of whole blood, DBSs, plasma, and blood cells. Approximately 4% of detected proteins' abundance was significantly altered during blood drying into blood spots, with overwhelming disturbances in cytoplasmic fraction. We also reported a novel finding suggesting a decrease in the level of membrane/cytoskeletal proteins (SLC4A1, RHAG, DSC1, DSP, and JUP) and an increase in the level of proteins (ATG3, SEC14L4, and NRBP1) related to intracellular trafficking. Furthermore, we identified 19 temporally dynamic proteins in DBS samples stored at room temperature for up to 6 months. There were three declined cytoskeleton-related proteins (RDX, SH3BGRL3, and MYH9) and four elevated proteins (XPO7, RAN, SLC2A1, and SLC29A1) involved in cytoplasmic transport as representatives. The instability was governed predominantly by hydrophilic proteins and enhanced significantly with an increasing storage time. Our analyses provide comprehensive knowledge of both short- and long-term storage stability of DBS proteins, forming the foundation for the widespread use of DBS in clinical proteomics and other analytical applications.

The indole ring of tryptophan can form NH/π hydrogen bonds, acting both as a hydrogen donor at the NH group in the pyrrole subring and as a hydrogen acceptor at the benzene subring. In the structural core of the trimeric stable protein Pholiota squarrosa lectin (PhoSL), three indoles are symmetrically arranged and form NH/π hydrogen bonds among each other. Here, we conducted quantum chemical calculations on this indole triad by using various methods and basis sets. The analyses revealed cooperativity in triad formation, with the many-body effect contributing approximately -2 kcal mol -1 , which significantly stabilizes this protein. Symmetry-adapted perturbation theory ascribed this effect to the induced polarization. The electrostatic potential and atomic charges indeed revealed a charge redistribution through the NH/π hydrogen bond, which was favorable for triad formation.

Thermal inactivation is a major bottleneck to the scalable production, storage, and transportation of protein-based reagents and therapies. Failures in temperature control both compromise protein bioactivity and increase the risk of microorganismal contamination. Herein, we report the rational design of fluorochemical additives that promiscuously bind to and coat the surfaces of proteins to enable their stable dispersion within fluorous solvents. By replacing traditional aqueous liquids with fluorinated media, this strategy conformationally rigidifies proteins to preserve their structure and function at extreme temperatures (≥90 °C). We show that fluorous protein formulations resist contamination by bacterial, fungal, and viral pathogens, which require aqueous environments for survival, and display equivalent serum bioavailability to standard saline samples in animal models. Importantly, by designing dispersants that decouple from the protein surface in physiologic solutions, we deliver a fluorochemical formulation that does not alter the pharmacologic function or safety profile of the functionalized protein in vivo . As a result, this nonaqueous protein storage paradigm is poised to open technological opportunities in the design of shelf-stable protein reagents and biopharmaceuticals.

The objective of this work was to investigate the effect of succinylation treatment on the physicochemical properties of black bean proteins (BBPI), and the relationship mechanism between BBPI structure and gel properties was further analyzed. The results demonstrated that the covalent formation of higher-molecular-weight complexes with BBPI could be achieved by succinic anhydride (SA). With the addition of SA at 10% (v/v), the acylation of proteins amounted to 92.53 ± 1.10%, at which point there was a minimized particle size of the system (300.90 ± 9.57 nm). Meanwhile, the protein structure was stretched with an irregular curl content of 34.30% and the greatest processable flexibility (0.381 ± 0.004). The dense three-dimensional mesh structure of the hydrogel as revealed by scanning electron microscopy was the fundamental prerequisite for the ability to resist external extrusion. The thermally induced hydrogels of acylated proteins with 10% (v/v) addition of SA showed excellent gel elastic behavior (1.44 ± 0.002 nm) and support capacity. Correlation analysis showed that the hydrogel strength and stability of hydrogels were closely related to the changes in protein conformation. This study provides theoretical guidance for the discovery of flexible proteins and their application in hydrogels.

The HMGB1 protein typically serves as a DNA chaperone that assists DNA-repair enzymes and transcription factors but can translocate from the nucleus to the cytoplasm or even to extracellular space upon some cellular stimuli. One of the factors that triggers the translocation of HMGB1 is its phosphorylation near a nuclear localization sequence by protein kinase C (PKC), although the exact modification sites on HMGB1 remain ambiguous. In this study, using spectroscopic methods, we investigated the HMGB1 phosphorylation and its impact on the molecular properties of the HMGB1 protein. Our nuclear magnetic resonance (NMR) data on the full-length HMGB1 protein showed that PKC specifically phosphorylates the A-box domain, one of the DNA binding domains of HMGB1. Phosphorylation of S46 and S53 was particularly efficient. Over a longer reaction time, PKC phosphorylated some additional residues within the HMGB1 A-box domain. Our fluorescence-based binding assays showed that the phosphorylation significantly reduces the binding affinity of HMGB1 for DNA. Based on the crystal structures of HMGB1-DNA complexes, this effect can be ascribed to electrostatic repulsion between the negatively charged phosphate groups at the S46 side chain and DNA backbone. Our data also showed that the phosphorylation destabilizes the folding of the A-box domain. Thus, phosphorylation by PKC weakens the DNA-binding affinity and folding stability of HMGB1.

Fc-fusion proteins are an emerging class of protein therapeutics that combine the properties of biological ligands with the unique properties of the fragment crystallizable (Fc) domain of an immunoglobulin G (IgG). Due to their diverse higher-order structures (HOSs), Fc-fusion proteins remain challenging characterization targets within biopharmaceutical pipelines. While high-resolution biophysical tools are available for HOS characterization, they frequently demand extended time frames and substantial quantities of purified samples, rendering them impractical for swiftly screening candidate molecules. Herein, we describe the development of ion mobility-mass spectrometry (IM-MS) and collision-induced unfolding (CIU) workflows that aim to fill this technology gap, where we focus on probing the HOS of a model Fc-Interleukin-10 (Fc-IL-10) fusion protein engineered using flexible glycine-serine linkers. We evaluate the ability of these techniques to probe the flexibility of Fc-IL-10 in the absence of bulk solvent relative to other proteins of similar size, as well as localize structural changes of low charge state Fc-IL-10 ions to specific Fc and IL-10 unfolding events during CIU. We subsequently apply these tools to probe the local effects of glycine-serine linkers on the HOS and stability of IL-10 homodimer, which is the biologically active form of IL-10. Our data reveals that Fc-IL-10 produces significantly more structural transitions during CIU and broader IM profiles when compared to a wide range of model proteins, indicative of its exceptional structural dynamism. Furthermore, we use a combination of enzymatic approaches to annotate these intricate CIU data and localize specific transitions to the unfolding of domains within Fc-IL-10. Finally, we detect a strong positive, quadratic relationship between average linker mass and fusion protein stability, suggesting a cooperative influence between glycine-serine linkers and overall fusion protein stability. This is the first reported study on the use of IM-MS and CIU to characterize HOS of Fc-fusion proteins, illustrating the practical applicability of this approach.

Accurately predicting protein-ligand interactions is essential in computational molecular biochemistry and in silico drug development. Monitoring changes in molecular dipole moments through molecular dynamics simulations provides valuable insights into dipole-dipole interactions, which are critical for understanding protein structure stability and predicting protein-ligand binding affinity. In this study, we propose a novel method to monitor changes in the interangle between dipole vectors of residue molecules within proteins and ligand molecules, aiming to evaluate the strength and consistency of interactions within the complex. Additionally, we extend the concept of positional root-mean-square fluctuation (RMSF), commonly used for protein structure stability analysis, to dipole moments, thus defining dipole moment RMSF. This enables us to analyze the stability of dipole moments for each residue within the protein and compare them across residues and between binding and non-binding complexes. Using the CRBP1-retinoic acid complex as our model system, we observed a significant difference in the interangle change of dipole moments for the key residue at the residue-level between the non-binding and binding complexes. Furthermore, we found that the dipole moment RMSF value of the non-binding complex was substantially larger than that of the binding complex, indicating greater dipole moment instability in the non-binding complex. Leveraging the concept of scalability inherent in the calculation of dipole moment vectors, we systematically expanded the residues within the protein's primary secondary structure. Our dipole moment analysis approach can provide valuable predictive insights into complex candidates, especially in situations where experimental comparisons are challenging.

PGM1 -linked congenital disorder of glycosylation (PGM1-CDG) is an autosomal recessive disease characterized by several phenotypes, some of which are life-threatening. Research focusing on the disease-related variants of the α-D-phosphoglucomutase 1 (PGM1) protein has shown that several are insoluble in vitro and expressed at low levels in patient fibroblasts. Due to these observations, we hypothesized that some disease-linked PGM1 protein variants are structurally destabilized and subject to protein quality control (PQC) and rapid intracellular degradation. Employing yeast-based assays, we show that a disease-associated human variant, PGM1 L516P, is insoluble, inactive, and highly susceptible to ubiquitylation and rapid degradation by the proteasome. In addition, we show that PGM1 L516P forms aggregates in S. cerevisiae and that both the aggregation pattern and the abundance of PGM1 L516P are chaperone-dependent. Finally, using computational methods, we perform saturation mutagenesis to assess the impact of all possible single residue substitutions in the PGM1 protein. These analyses identify numerous missense variants with predicted detrimental effects on protein function and stability. We suggest that many disease-linked PGM1 variants are subject to PQC-linked degradation and that our in silico site-saturated data set may assist in the mechanistic interpretation of PGM1 variants.

Human γD-crystallin belongs to a crucial family of proteins known as crystallins located in the fiber cells of the human lens. Since crystallins do not undergo any turnover after birth, they need to possess remarkable thermodynamic stability. However, their sporadic misfolding and aggregation, triggered by environmental perturbations or genetic mutations, constitute the molecular basis of cataracts, which is the primary cause of blindness in the globe according to the World Health Organization. Here, we investigate the impact of high pressure on the conformational landscape of wild-type HγD-crystallin using replica exchange molecular dynamics simulations augmented with principal component analysis. We find pressure to have a modest impact on global measures of protein stability, such as root-mean-square displacement and radius of gyration. Upon projecting our trajectories along the first two principal components from principal component analysis, however, we observe the emergence of distinct free energy basins at high pressures. By screening local order parameters previously shown or hypothesized as markers of HγD-crystallin stability, we establish correlations between a tyrosine-tyrosine aromatic contact within the N-terminal domain and the protein's end-to-end distance with projections along the first and second principal components, respectively. Furthermore, we observe the simultaneous contraction of the hydrophobic core and its intrusion by water molecules. This exploration sheds light on the intricate responses of HγD-crystallin to elevated pressures, offering insights into potential mechanisms underlying its stability and susceptibility to environmental perturbations, crucial for understanding cataract formation.

The disaccharide trehalose is generally acknowledged as a superior stabilizer of proteins and other biomolecules in aqueous environments. Despite many theories aiming to explain this, the stabilization mechanism is still far from being fully understood. This study compares the stabilizing properties of trehalose with those of the structurally similar disaccharide sucrose. The stability has been evaluated for the two proteins, lysozyme and myoglobin, at both low and high temperatures by determining the glass transition temperature, T g , and the denaturation temperature, T den . The results show that the sucrose-containing samples exhibit higher T den than the corresponding trehalose-containing samples, particularly at low water contents. The better stabilizing effect of sucrose at high temperatures may be explained by the fact that sucrose, to a greater extent, binds directly to the protein surface compared to trehalose. Both sugars show T den elevation with an increasing sugar-to-protein ratio, which allows for a more complete sugar shell around the protein molecules. Finally, no synergistic effects were found by combining trehalose and sucrose. Conclusively, the exact mechanism of protein stabilization may vary with the temperature, as influenced by temperature-dependent interactions between the protein, sugar, and water. This variability can make trehalose to a superior stabilizer under some conditions and sucrose under others.

The T-cell receptor (TCR) is a crucial molecule in cellular immunity. The single-chain T-cell receptor (scTCR) is a potential format in TCR therapeutics because it eliminates the possibility of αβ-TCR mispairing. However, its poor stability and solubility impede the in vitro study and manufacturing of therapeutic applications. In this study, some conserved structural motifs are identified in variable domains regardless of germlines and species. Theoretical analysis helps to identify those unfavored factors and leads to a general strategy for stabilizing scTCRs by substituting residues at exact IMGT positions with beneficial propensities on the consensus sequence of germlines. Several representative scTCRs are displayed to achieve stability optimization and retain comparable binding affinities with the corresponding αβ-TCRs in the range of μM to pM. These results demonstrate that our strategies for scTCR engineering are capable of providing the affinity-enhanced and specificity-retained format, which are of great value in facilitating the development of TCR-related therapeutics.

Acyl capping groups stabilize α-helices relative to free N-termini by providing one additional C═O i ···H i +4 -N hydrogen bond. The electronic properties of acyl capping groups might also directly modulate α-helix stability: electron-rich N-terminal acyl groups could stabilize the α-helix by strengthening both i / i + 4 hydrogen bonds and i / i + 1 n → π* interactions. This hypothesis was tested in peptides X-AKAAAAKAAAAKAAGY-NH 2 , where X = different acyl groups. Surprisingly, the most electron-rich acyl groups (pivaloyl and iso -butyryl) strongly destabilized the α-helix. Moreover, the formyl group induced nearly identical α-helicity to that of the acetyl group, despite being a weaker electron donor for hydrogen bonds and for n → π* interactions. Other acyl groups exhibited intermediate α-helicity. These results indicate that the electronic properties of the acyl carbonyl do not directly determine the α-helicity in peptides in water. In order to understand these effects, DFT calculations were conducted on α-helical peptides. Using implicit solvation, α-helix stability correlated with acyl group electronics, with the pivaloyl group exhibiting closer hydrogen bonds and n → π* interactions, in contrast to the experimental results. However, DFT and MD calculations with explicit water solvation revealed that hydrogen bonding to water was impacted by the sterics of the acyl capping group. Formyl capping groups exhibited the closest water-amide hydrogen bonds, while pivaloyl groups exhibited the longest. In α-helices in the PDB, the highest frequency of close amide-water hydrogen bonds is observed when the N-cap residue is Gly. The combination of experimental and computational results indicates that solvation (hydrogen bonding of water) to the N-terminal amide groups is a central determinant of α-helix stability.

HdeA and HdeB are dimeric ATP-independent acid-stress chaperones, which protect the periplasmic proteins of enteric bacteria at pH 2.0 and 4.0, respectively, during their passage through the acidic environment of the mammalian stomach. Despite being structurally similar, they exhibit distinct functional pH optima and conformational prerequisite for their chaperone action. HdeA undergoes a dimer-to-monomer transition at pH 2.0, whereas HdeB remains dimeric at pH 4.0. The monomerization of HdeA exposes its hydrophobic motifs, which facilitates its interaction with the partially folded substrates. How HdeB functions despite maintaining its dimeric conformation has been poorly elucidated in the literature. Herein, we characterized the conformational states and stability of HdeB at its physiologically relevant pH and compared the data with those of HdeA. At pH 4.0, HdeB exhibited distinct spectroscopic signatures and higher stability against heat and guanidine-HCl-induced denaturation than at pH 7.5. We affirm that the pH 4.0 conformer of HdeB was distinct from that at pH 7.5 and that these two conformational states were hierarchically unrelated. Salt-bridge mutations that perturbed HdeB's intersubunit interactions resulted in the loss of its stability and function at pH 4.0. In contrast, mutations affecting intrasubunit interactions enhanced its function, albeit with a reduction in stability. These findings suggest that, unlike HdeA, HdeB acts as a noncanonical chaperone, where pH-dependent stability and conformational rearrangement at pH 4.0 play a core role in its chaperone function rather than its surface hydrophobicity. Such rearrangement establishes a stability-function trade-off that allows HdeB to function while maintaining its stable dimeric state.

NANOG protein levels correlate with stem cell pluripotency. NANOG concentrations fluctuate constantly with low NANOG levels leading to spontaneous cell differentiation. Previous literature implicated Pin1, a phosphorylation-dependent prolyl isomerase, as a key player in NANOG stabilization. Here, using NMR spectroscopy, we investigate the molecular interactions of Pin1 with the NANOG unstructured N-terminal domain that contains a PEST sequence with two phosphorylation sites. Phosphorylation of NANOG PEST peptides increases affinity to Pin1. By systematically increasing the amount of cis PEST conformers, we show that the peptides bind tighter to the prolyl isomerase domain (PPIase) of Pin1. Phosphorylation and cis Pro enhancement at both PEST sites lead to a 5-10-fold increase in NANOG binding to the Pin1 WW domain and PPIase domain, respectively. The cis- populated NANOG PEST peptides can be potential inhibitors for disrupting Pin1-dependent NANOG stabilization in cancer stem cells.

The fibrillation of therapeutic peptides can present significant quality concerns and poses challenges for manufacturing and storage. A fundamental understanding of the mechanisms of fibrillation is critical for the rational design of fibrillation-resistant peptide drugs and can accelerate product development by guiding the selection of solution-stable candidates and formulations. The studies reported here investigated the effects of structural modifications on the fibrillation of a 29-residue peptide (PepA) and two sequence modified variants (PepB, PepC). The C-terminus of PepA was amidated, whereas both PepB and PepC retained the carboxylate, and Ser16 in PepA and PepB was substituted with a helix-stabilizing residue, α-aminoisobutyric acid (Aib), in PepC. In thermal denaturation studies by far-UV CD spectroscopy and fibrillation kinetic studies by fluorescence and turbidity measurements, PepA and PepB showed heat-induced conformational changes and were found to form fibrils, whereas PepC did not fibrillate and showed only minor changes in the CD signal. Pulsed hydrogen-deuterium exchange mass spectrometry (HDX-MS) showed a high degree of protection from HD exchange in mature PepA fibrils and its proteolytic fragments, indicating that most of the sequence had been incorporated into the fibril structure and occurred nearly simultaneously throughout the sequence. The effects of the net peptide charge and formulation pH on fibrillation kinetics were investigated. In real-time stability studies of two formulations of PepA at pH's 7.4 and 8.0, analytical methods detected significant changes in the stability of the formulations at different time points during the study, which were not observed during accelerated studies. Additionally, PepA samples were withdrawn from real-time stability and subjected to additional stress (40 °C, continuous shaking) to induce fibrillation; an approach that successfully amplified oligomers or prefibrillar species previously undetected in a thioflavin T assay. Taken together, these studies present an approach to differentiate and characterize fibrillation risk in structurally related peptides under accelerated and real-time conditions, providing a model for rapid, iterative structural design to optimize the stability of therapeutic peptides.

Lactoferrin (LTF) has diverse biological activities and is widely used in functional foods and active additives. Nevertheless, evaluating the proteoform heterogeneity, conformational stability, and activity of LTF remains challenging during its production and storage processes. In this study, we describe the implementation of native mass spectrometry (nMS), glycoproteomics, and an antimicrobial activity assay to assess the quality of LTF. We systematically characterize the purity, glycosylation heterogeneity, conformation, and thermal stability of LTF samples from different sources and transient high-temperature treatments by using nMS and glycoproteomics. Meanwhile, the nMS peak intensity and antimicrobial activity of LTF samples after heat treatment decreased significantly, and the two values were positively correlated. The nMS results provide essential molecular insights into the conformational stability and glycosylation heterogeneity of different LTF samples. Our results underscore the great potential of nMS for LTF quality control and activity evaluation in industrial production.

While synthetic biology has advanced complex capabilities such as sensing and molecular synthesis in aqueous solutions, important applications may also be pursued for biological systems in solid materials. Harsh processing conditions used to produce many synthetic materials such as plastics make the incorporation of biological functionality challenging. One technology that shows promise in circumventing these issues is cell-free protein synthesis (CFPS), where core cellular functionality is reconstituted outside the cell. CFPS enables genetic functions to be implemented without the complications of membrane transport or concerns over the cellular viability or release of genetically modified organisms. Here, we demonstrate that dried CFPS reactions have remarkable tolerance to heat and organic solvent exposure during the casting processes for polymer materials. We demonstrate the utility of this observation by creating plastics that have spatially patterned genetic functionality, produce antimicrobials in situ, and perform sensing reactions. The resulting materials unlock the potential to deliver DNA-programmable biofunctionality in a ubiquitous class of synthetic materials.

Tyrosinase is capable of oxidizing tyrosine residues in proteins, leading to intermolecular protein cross-linking, which could modify the protein network of food and improve the texture of food. To obtain the recombinant tyrosinase with microbial cell factory instead of isolation tyrosinase from the mushroom Agaricus bisporus , a TYR expression cassette was constructed in this study. The expression cassette was electroporated into Trichoderma reesei Rut-C30 and integrated into its genome, resulting in a recombinant strain C30-TYR. After induction with microcrystalline cellulose for 7 days, recombinant tyrosinase could be successfully expressed and secreted by C30-TYR, corresponding to approximately 2.16 g/L tyrosinase in shake-flask cultures. The recombinant TYR was purified by ammonium sulfate precipitation and gel filtration, and the biological activity of purified TYR was 45.6 U/mL. The purified TYR could catalyze the cross-linking of glycinin, and the emulsion stability index of TYR-treated glycinin emulsion was increased by 30.6% compared with the untreated one. The cross-linking of soy glycinin by TYR resulted in altered properties of oil-in-water emulsions compared to emulsions stabilized by native glycinin. Therefore, cross-linking with this recombinant tyrosinase is a feasible approach to improve the properties of protein-stabilized emulsions and gels.

How mutations impact protein stability and structure dynamics is crucial for understanding the pathological process and rational drug design. Herein, we establish a time-resolved native mass spectrometry (TR-nMS) platform via a rapid-mixing capillary apparatus for monitoring the acid-initiated protein unfolding process. The molecular details in protein structure unfolding are further profiled by a 193 nm ultraviolet photodissociation (UVPD) analysis of the structure-informative photofragments. Compared with the wild-type dihydrofolate reductase (WT-DHFR), the M42T/H114R mutant (MT-DHFR) exhibits a significant stability decrease in TR-nMS characterization. UVPD comparisons of the unfolding intermediates and original DHFR forms indicate the special stabilization effect of cofactor NADPH on DHFR structure, and the M42T/H114R mutations lead to a significant decrease in NADPH-DHFR interactions, thus promoting the structure unfolding. Our study paves the way for probing the mutation-induced subtle changes in the stability and structure dynamics of drug targets.

Though protein stability and aggregation have been well characterized in dilute solutions, the influence of a confining environment that exists (e.g., in intercellular and tissue spaces and therapeutic formulations) on the protein structure is largely unknown. Herein, the effects of confinement on stability and aggregation were explored for proteins of different sizes, stability, and hydrophobicity when encapsulated in hydrophilic poly(ethylene glycol) hydrogels. Denaturation curves show linear correlations between confinement size (mesh size) and thermodynamic stability, i.e., unfolding free energy and surface area accessible for solvation (represented by m -value). Two clusters of protein types are identifiable from these correlations; the clusters are defined by differences in protein stability, surface area, and aggregation propensity. Proteins with higher stability, larger surface area, and lower aggregation propensity (e.g., lysozyme and myoglobin) are less affected by the confinement imposed by mesh size than proteins with lower stability, smaller surface area, and higher aggregation propensity (e.g., growth hormone and aldehyde dehydrogenase). According to aggregation kinetics measured by thioflavin T fluorescence, confinement in smaller mesh sizes resulted in slower aggregation rates than that in larger mesh sizes. Compared to that in buffer solution, the confinement of a hydrophobic protein (e.g., human insulin) in the hydrogels accelerates protein aggregation. Conversely, the confinement of a hydrophilic protein (e.g., human amylin) in the hydrogels decelerates or prevents aggregation, with the rates of aggregation inversely proportional to mesh size. These findings provide new insights into protein conformational stability in confined microenvironments relevant to various cellular, tissue, and therapeutics scenarios.

Single-point mutations in kinase proteins can affect their stability and fitness, and computational analysis of these effects can provide insights into the relationships among protein sequence, structure, and function for this enzyme family. To assess the impact of mutations on protein stability, we used a sequence-based Potts Hamiltonian model trained on a kinase family multiple-sequence alignment (MSA) to calculate the statistical energy (fitness) effects of mutations and compared these against relative folding free energies (ΔΔ G s) calculated from all-atom molecular dynamics free energy perturbation (FEP) simulations in explicit solvent. The fitness effects of mutations in the Potts model (Δ E s) showed good agreement with experimental thermostability data (Pearson r = 0.68), similar to the correlation we observed with ΔΔ G s predicted from structure-based relative FEP simulations. Recognizing the possible advantages of using Potts models to rapidly estimate protein stability effects of kinase mutations seen in cancer genomics data, we used the Potts statistical energy model to estimate the stability effects of 65 conservative and nonconservative mutations across three distinct kinases (Wee1, Abl1, and Cdc7) with somatic mutations reported in the Genomic Data Commons (GDC) database. The Δ E s of these mutations calculated from the Potts model are consistent with the corresponding ΔΔ G s from FEP simulations (Pearson ratio of 0.72). The agreement between these methods suggests that the Potts model may be used as a sequence-based tool for high-throughput screening of mutational effects as part of a computational pipeline for predicting the stability effects of mutations. We also demonstrate how the scalability of the fitness-based Potts model calculations permits analyses that are not easily accessed using FEP simulations. To this end, we employed site-saturation mutagenesis in the Potts model in order to investigate the relative stability effects of mutations seen in different cancer evolutionary scenarios. We used this approach to analyze the effects of drug pressure in Abl kinase by contrasting the relative fitness penalties of somatic mutations seen in miscellaneous cancer types with those calculated for mutations associated with cancer drug resistance. We observed that, in contrast to somatic mutations of Abl seen in various tumors that appear to have evolved neutrally, cancer mutations that evolved under drug pressure in Abl-targeted therapies tend to preserve enzyme stability.

The impact of macromolecular crowding on biological macromolecules has been elucidated through the excluded volume phenomenon and soft interactions. However, it has often been difficult to provide a clear demarcation between the two regions. Here, using temperature-dependent dynamics (local and global) of the multidomain protein human serum albumin (HSA) in the presence of commonly used synthetic crowders (Dextran 40, PEG 8, Ficoll 70, and Dextran 70), we have shown the presence of a transition that serves as a bridge between the soft and hard regimes. The bridging region is independent of the crowder identity and displays no apparent correlation with the critical overlap concentration of the polymeric crowding agents. Moreover, the dynamics of domains I and II and the protein gating motion respond differently, thereby bringing to the fore the asymmetry underlying the crowder influence on HSA. In addition, solvent-coupled and decoupled protein motions indicate the heterogeneity of the dynamic landscape in the crowded milieu. We also propose an intriguing correlation between protein stability and dynamics, with increased global stability being accompanied by eased local domain motion.

Protein thermodynamic stability is essential to clarify the relationships among structure, function, and interaction. Therefore, developing a faster and more accurate method to predict the impact of the mutations on protein stability is helpful for protein design and understanding the phenotypic variation. Recent studies have shown that protein embedding will be particularly powerful at modeling sequence information with context dependence, such as subcellular localization, variant effect, and secondary structure prediction. Herein, we introduce a novel method, ProSTAGE, which is a deep learning method that fuses structure and sequence embedding to predict protein stability changes upon single point mutations. Our model combines graph-based techniques and language models to predict stability changes. Moreover, ProSTAGE is trained on a larger data set, which is almost twice as large as the most used S2648 data set. It consistently outperforms all existing state-of-the-art methods on mutation-affected problems as benchmarked on several independent data sets. The protein embedding as the prediction input achieves better results than the previous results, which shows the potential of protein language models in predicting the effect of mutations on proteins. ProSTAGE is implemented as a user-friendly web server.

Determining changes in the protein's thermal stability following mutations is critical in protein engineering and understanding pathogenic missense mutations. Despite the development of various computational methods to predict the effects of single-point mutations, their accuracy remains limited. In this study, we propose a new computational method, OmeDDG, that more accurately predicts mutation-induced Gibbs free energy changes in protein folding (ΔΔ G ). OmeDDG takes the sequences of wild-type and mutant proteins as input, utilizes OmegaFold to obtain the 3D structure, employs a convolutional neural network to extract structural features, and combines them with protein mutation features and pretraining features to predict the stability of single-point mutations in proteins. We performed a comprehensive comparison between OmeDDG and other available prediction methods on four blind test datasets, confirming that OmeDDG can effectively enhance protein mutation prediction performance. Notably, on the antisymmetric dataset Ssym, OmeDDG achieves the best performance, demonstrating favorable antisymmetry with PCC = 0.79 and RMSE = 0.96 for forward mutations and PCC = 0.77 and RMSE = 0.97 for reverse mutant types.

The oral delivery of protein therapeutics offers numerous advantages for patients but also presents significant challenges in terms of development. Currently, there is limited knowledge available regarding the stability and shelf life of orally delivered protein therapeutics. In this study, a comprehensive assessment of the stability of an orally delivered solid dosage variable domain of heavy-chain antibody (VHH antibody) drug product was conducted. Four stability related quality attributes that undergo change as a result of thermal and humidity stress were identified. Subsequently, these attributes were modeled using an accelerated stability approach facilitated by ASAP prime software. To the best of our knowledge, this is the first time that this approach has been reported for an antibody drug product. We observed overall good model quality and accurate predictions regarding the protein stability during storage. Notably, we discovered that protein aggregation, formed through a degradation pathway, requires additional adjustments to the modeling method. In summary, the ASAP approach demonstrated promising results in predicting the stability of this complex solid-state protein formulation. This study sheds light on the stability and shelf life of orally delivered protein therapeutics, addressing an important knowledge gap in the field.

Hydrogen exchange (HX) has become an important tool to monitor protein structure and dynamics. The interpretation of HX data with respect to protein structure requires understanding of the factors that influence exchange. Simulated protein structures can be validated by comparing experimental deuteration profiles with the profiles derived from the modeled protein structure. To do this, we propose here a new method, POPPeT, for protection factor prediction based on protein motions that enable HX. By comparing POPPeT with two existing methods, the phenomenological approximation and COREX, we show enhanced predictability measured at both protection factor and deuteration level. This method can be subsequently used by modeling strategies for protein structure prediction.ᅟ [ABSTRACT FROM AUTHOR]

Gas-phase hydrogen/deuterium exchange measured by mass spectrometry (gas-phase HDX-MS) is a fast method to probe the conformation of protein ions. The use of gas-phase HDX-MS to investigate the structure and interactions of protein complexes is however mostly unharnessed. Ionizing proteins under conditions that maximize preservation of their native structure (native MS) enables the study of solution-like conformation for milliseconds after electrospray ionization (ESI), which enables the use of ND3-gas inside the mass spectrometer to rapidly deuterate heteroatom-bound non-amide hydrogens. Here, we explored the utility of gas-phase HDX-MS to examine protein-protein complexes and inform on their binding surface and the structural consequences of gas-phase dissociation. Protein complexes ranging from 24 kDa dimers to 395 kDa 24mers were analyzed by gas-phase HDX-MS with subsequent collision-induced dissociation (CID). The number of exchangeable sites involved in complex formation could, therefore, be estimated. For instance, dimers of cytochrome c or α-lactalbumin incorporated less deuterium/subunit than their unbound monomer counterparts, providing a measure of the number of heteroatom-bound side-chain hydrogens involved in complex formation. We furthermore studied if asymmetric charge-partitioning upon dissociation of protein complexes caused intermolecular H/D migration. In larger multimeric protein complexes, the dissociated monomer showed a significant increase in deuterium. This indicates that intermolecular H/D migration occurs as part of the asymmetric partitioning of charge during CID. We discuss several models that may explain this increase deuterium content and find that a model where only deuterium involved in migrating charge can account for most of the deuterium enrichment observed on the ejected monomer. In summary, the deuterium content of the ejected subunit can be used to estimate that of the intact complex with deviations observed for large complexes accounted for by charge migration.ᅟ [ABSTRACT FROM AUTHOR]

The ability to design detergents that are suitable for protein analysis by mass spectrometry (MS) represents an on-going challenge in the field of native MS. Desirable detergent characteristics include charge-reducing properties and low gas-phase stabilities of complexes formed with proteins. In this work, the gas-phase properties of oligoglycerol detergents (OGDs) are optimized by fine tuning of their molecular structure. Furthermore, a tandem mass spectrometry (MS/MS) approach is presented that estimates the gas-phase properties of detergents simply by studying the dissociation behaviour of protein-detergent complexes (PDCs) formed with the soluble protein β-lactoglobulin (BLG).ᅟ [ABSTRACT FROM AUTHOR]

Native top-down mass spectrometry (MS) and ion mobility spectrometry (IMS) were applied to characterize the interaction of a molecular tweezer assembly modulator, CLR01, with tau, a protein believed to be involved in a number of neurodegenerative disorders, including Alzheimer's disease. The tweezer CLR01 has been shown to inhibit aggregation of amyloidogenic polypeptides without toxic side effects. ESI-MS spectra for different forms of tau protein (full-length, fragments, phosphorylated, etc.) in the presence of CLR01 indicate a primary binding stoichiometry of 1:1. The relatively high charging of the protein measured from non-denaturing solutions is typical of intrinsically disordered proteins, such as tau. Top-down mass spectrometry using electron capture dissociation (ECD) is a tool used to determine not only the sites of post-translational modifications but also the binding site(s) of non-covalent interacting ligands to biomolecules. The intact protein and the protein-modulator complex were subjected to ECD-MS to obtain sequence information, map phosphorylation sites, and pinpoint the sites of inhibitor binding. The ESI-MS study of intact tau proteins indicates that top-down MS is amenable to the study of various tau isoforms and their post-translational modifications (PTMs). The ECD-MS data point to a CLR01 binding site in the microtubule-binding region of tau, spanning residues K294-K331, which includes a six-residue nucleating segment PHF6 (VQIVYK) implicated in aggregation. Furthermore, ion mobility experiments on the tau fragment in the presence of CLR01 and phosphorylated tau reveal a shift towards a more compact structure. The mass spectrometry study suggests a picture for the molecular mechanism of the modulation of protein-protein interactions in tau by CLR01.ᅟ [ABSTRACT FROM AUTHOR]

The initial stages of protein unfolding may reflect the stability of the entire fold and can also reveal which parts of a protein can be perturbed, without restructuring the rest. In this work, we couple UVPD with activated ion mobility mass spectrometry to measure how three model proteins start to unfold. Ubiquitin, cytochrome c and myoglobin ions produced via nESI from salty solutions are subjected to UV irradiation pre-mobility separation; experiments are conducted with a range of source conditions which alter the conformation of the precursor ion as shown by the drift time profiles. For all three proteins, the compact structures result in less fragmentation than more extended structures which emerge following progressive in-source activation. Cleavage sites are found to differ between conformational ensembles, for example, for the dominant charge state of cytochrome c [M + 7H]7+, cleavage at Phe10, Thr19 and Val20 was only observed in activating conditions whilst cleavage at Ala43 is dramatically enhanced. Mapping the photo-cleaved fragments onto crystallographic structures provides insight into the local structural changes that occur as protein unfolding progresses, which is coupled to global restructuring observed in the drift time profiles.Graphical Abstract [ABSTRACT FROM AUTHOR]

Proteins are commonly encapsulated in alginate gels for drug delivery and tissue-engineering applications. However, there is limited knowledge of how encapsulation impacts intrinsic protein properties such as folding stability or unfolding kinetics. Here, we use fast relaxation imaging (FReI) to image protein unfolding in situ in alginate hydrogels after applying a temperature jump. Based on changes in the Förster resonance energy transfer (FRET) response of FRET-labeled phosphoglycerate kinase (PGK), we report the quantitative impact of multiple alginate hydrogel concentrations on protein stability and folding dynamics. The gels stabilize PGK by increasing its melting temperature up to 18.4 °C, and the stabilization follows a nonmonotonic dependence on the alginate density. In situ kinetic measurements also reveal that PGK deviates more from two-state folding behavior in denser gels and that the gel decreases the unfolding rate and accelerates the folding rate of PGK, compared to buffer. Phi-value analysis suggests that the folding transition state of an encapsulated protein is structurally similar to that of folded protein. This work reveals both beneficial and negative impacts of gel encapsulation on protein folding, as well as potential mechanisms contributing to altered stability.

Accurate prediction of protein stability upon a point mutation has important applications in drug discovery and personalized medicine. It remains a challenging issue in computational biology. Existing computational prediction methods, which range from mechanistic to supervised learning approaches, have experienced limited progress over the last few decades. This stagnation is largely due to their heavy reliance on both the quantity and quality of the training data. This is evident in recent state-of-the-art methods that continue to yield substantial errors on two challenging blind test sets: frataxin and p53, with average root-mean-square errors exceeding 3 and 1.5 kcal/mol, respectively, which is still above the theoretical 1 kcal/mol prediction barrier. Rigorous approaches, on the other hand, offer greater potential for accuracy without relying on training data but are computationally demanding and require both wild-type and mutant structure information. Although they showed high accuracy for conserving mutations, their performance is still limited for charge-changing mutation cases. This might be due to the lack of an available mutant structure, often represented by a simplified capped peptide. The recent advances in protein structure prediction methods now make it possible to obtain structures comparable to experimental ones, including complete mutant structure information. In this work, we compare the performance of supervised learning-based methods and rigorous approaches for predicting protein stability on point mutations in difficult targets: frataxin and p53. The rigorous alchemical method significantly surpasses state-of-the-art techniques in terms of both the root-mean-squared error and Pearson correlation coefficient in these two challenging blind test sets. Additionally, we propose an improved alchemical method that employs the pmx double-system/single-box approach to accurately predict the folding free energy change upon both conserving and charge-changing mutations. The enhanced protocol can accurately predict both types of mutations, thereby outperforming existing state-of-the-art methods in overall performance.

Colloidal stability is an important consideration when developing high concentration mAb formulations. PEG-induced protein precipitation is a commonly used assay to assess the colloidal stability of protein solutions. However, the practical usefulness and the current theoretical model for this assay have yet to be verified over a large formulation space across multiple mAbs and mAb-based modalities. In the present study, we used PEG-induced protein precipitation assays to evaluate colloidal stability of 3 mAbs in 24 common formulation buffers at 20 and 5 °C. These prediction assays were conducted at low protein concentration (1 mg/mL). We also directly characterized high concentration (100 mg/mL) formulations for cold-induced phase separation, turbidity, and concentratibility by ultrafiltration. This systematic study allowed analysis of the correlation between the results of low concentration assays and the high concentration attributes. The key findings of this study include the following: (1) verification of the usefulness of three different parameters ( C mid , μ B , and T cloud ) from PEG-induced protein precipitation assays for ranking colloidal stability of high concentration mAb formulations; (2) a new method to implement PEG-induced protein precipitation assay suitable for high throughput screening with low sample consumption; (3) improvement in the theoretical model for calculating robust thermodynamic parameters of colloidal stability (μ B and ε B ) that are independent of specific experimental settings; (4) systematic evaluation of the effects of pH and buffer salts on colloidal stability of mAbs in common formulation buffers. These findings provide improved theoretical and practical tools for assessing the colloidal stability of mAbs and mAb-based modalities during formulation development.

Salts reduce the p K a of weak acids by a mechanism sensitive to ion identity and concentration via charge screening of the deprotonated state. In this study, we utilize constant pH molecular dynamics simulations to understand the molecular mechanism behind the salt-dependent dissociation of aspartic acid (Asp). We calculate the p K a nonmonotonically, interpretable in the context of the law of matching water affinity. The net effect of salt on Asp acidity is governed by an interplay of solvation and competing ion interactions. The proposed mechanism is rather general and can be applicable to several problems in Hofmeister ion chemistry, such as pH effects on protein stability and soft matter interfaces.K a . Conversely, the cation size affects the p K a nonmonotonically, interpretable in the context of the law of matching water affinity. The net effect of salt on Asp acidity is governed by an interplay of solvation and competing ion interactions. The proposed mechanism is rather general and can be applicable to several problems in Hofmeister ion chemistry, such as pH effects on protein stability and soft matter interfaces.