Unfolding free energy changes determined by the linear extrapolation method. 1. Unfolding of phenylmethanesulfonyl alpha-chymotrypsin using different denaturants

Survey of chemical unfolding complexity as a unique stability assessment assay for monoclonal antibodies

Seventy-two intentionally sequence-diverse antibody variable regions were selected, expressed as IgG1 antibodies, and evaluated by chemical unfolding to survey the complexities of denaturant induced unfolding behavior. A two-transition fit well described the curves and uncovered a wide range of sensitivities to denaturant. Four general types of unfolding curves were observed: balanced traces (each transition responsible for half of the total unfolding curve), low-unfolding traces (first transition is a majority of the unfolding curve), high-unfolding traces (the second transition is the majority of the unfolding curve), and coincident traces (the two transitions are found close to each other, approximating a single transition). The complexity of the data from this survey indicates that focusing on the first inflection point or fitting a single transition model is likely an over-simplistic method for measuring stability by the chemical unfolding assay. Additionally, other conformational assays, such as thermal and low pH unfolding, showed no correlation with the chemical unfolding results, indicating that each of these assays provide alternate information on the different pathways of antibody conformational stability. These results provide a basis for beginning to better connect unfolding behavior to other physical phenotypic behaviors and production process behaviors.

MIA40 circumvents the folding constraints imposed by TRIAP1 function

The MIA40 relay system mediates the import of small cysteine-rich proteins into the intermembrane mitochondrial space (IMS). MIA40 substrates are synthesized in the cytosol and assumed to be disordered in their reduced state in this compartment. As they cross the outer mitochondrial membrane, MIA40 promotes the oxidation of critical native disulfides to facilitate folding, trapping functional species in the IMS. Here, we study the redox-controlled folding of TRIAP1, a small cysteine-rich protein with moonlighting function: regulating phospholipid trafficking between mitochondrial membranes in the IMS and preventing apoptosis in the cytosol. TRIAP1 dysregulation is connected to oncogenesis. Although TRIAP1 contains a canonical twin CX9C motif, its sequence characteristics and folding pathway deviate from typical MIA40 substrates. In its reduced state, TRIAP1 rapidly populates a hydrophobic collapsed, alpha-helical, and marginally stable molten globule. This intermediate, biases oxidative folding towards a non-native Cys37-Cys47 kinetic trap, slowing the reaction. MIA40 accelerates TRIAP1 folding rate by 30-fold, bypassing the formation of this folding trap. MIA40 drives the oxidation of the inner disulfide bond Cys18-Cys37, and subsequently, it can catalyze the formation of the outer disulfide bond Cys8-Cys47 to attain the native two-disulfide-bridged structure. We demonstrate that, unlike most MIA40 substrates, TRIAP1's folding pathway is strongly constrained by the structural requirements for its function in phospholipid traffic at the IMS. The obligatory population of a reduced, alpha-helical, metastable molten globule in the cytoplasm may explain TRIAP1's connection to the p53-dependent cell survival pathway, constituting a remarkable example of a functional molten globule state.

Intermolecular Misfolding Captured in Parallelly Organized Titin

The giant muscle protein titin is largely responsible for the passive elasticity of the muscles. The I-band part of titin is elastic, and its constitutive immunoglobulin (Ig) domains undergo force-induced unfolding and refolding when the muscle is stretched toward or beyond the end of the physiological range of sarcomere length. Correct folding of the titin Ig domains is essential to the structure and functions of titin. Although our knowledge of titin elasticity at the molecular level has been largely obtained from single molecule experiments, titin does not exist as an isolated molecule. Instead, six titins are parallelly organized in the muscle sarcomeres. It remains unknown what impact such a parallel organization brings on the folding of titin Ig domains and titin elasticity. Using the two-molecule force spectroscopy technique, here, we report the direct observation of the intermolecular misfolding of titin Ig domains that are arranged in parallel. Our results reveal that when parallelly arranged, two I94 domains can misfold into an intermolecular domain-swapped state that is thermally and mechanically stable. Such intermolecular misfolding may play important structural and functional roles in titin organization and elasticity.

Thermodynamic versus kinetic basis for the high conformational stability of nanobodies for therapeutic applications

Nanobodies (NB) are powerful tools for biotechnological and therapeutic applications. They strongly bind to their targets and are very stable. Early studies showed that NB unfolding is reversible and can be analyzed by equilibrium thermodynamics, whereas more recent studies focused on their kinetic stability in very harsh conditions that are far from storage or physiological temperatures (4-37 °C). Here, we show that the thermodynamic view of NB stability holds in a wide range of temperatures (18-100 °C). The thermodynamic stability of three different NBs did not correlate with binding affinity for their target. Alpha-Fold 2 analyses of these NBs showed structural differences in the binding site and hydrogen bond networks. We expect that our approach will be helpful to improve our capacity to enhance structure-function-stability relationships of NB.

Interactions of human serum albumin with phosphate and Tris buffers: impact on paclitaxel binding and nanoparticles self-assembly

AIM

To investigate the conformational changes in human serum albumin (HSA) caused by chemical (CD) and thermal denaturation (TD) at pH 7.4 and 9.9, crucial for designing controlled drug delivery systems with paclitaxel (PTX).

METHODS

Experimental and computational methods, including differential scanning calorimetry (DSC), UV-Vis and intrinsic fluorescence spectroscopy, mean diameter, polydispersity index (PDI), ζ-potential, encapsulation efficiency (EE), in vitro release and protein docking studies were conducted to study the HSA denaturation and nanoparticles (NPs) preparation.

RESULTS

TD at pH 7.4 produced smaller NPs (287.1 ± 12.9 nm) than CD at pH 7.4 with NPs (584.2 ± 47.7 nm). TD at pH 9.9 exhibited high EE (97.3 ± 0.2%w/w) with rapid PTX release (50% within 1h), whereas at pH 7.4 (96.4 ± 2.1%w/w), release only 40%. ζ-potentials were around -30 mV.

CONCLUSION

Buffer type and pH significantly influence NP properties. TD in PBS at pH 7.4, provided optimal conditions for a stable and efficient drug delivery system.

Cooperative Substructure and Energetics of Allosteric Regulation of the Catalytic Core of the E3 Ubiquitin Ligase Parkin by Phosphorylated Ubiquitin

Mutations in the parkin gene product Parkin give rise to autosomal recessive juvenile parkinsonism. Parkin is an E3 ubiquitin ligase that is a critical participant in the process of mitophagy. Parkin has a complex structure that integrates several allosteric signals to maintain precise control of its catalytic activity. Though its allosterically controlled structural reorganization has been extensively characterized by crystallography, the energetics and mechanisms of allosteric regulation of Parkin are much less well understood. Allostery is fundamentally linked to the energetics of the cooperative (sub)structure of the protein. Herein, we examine the mechanism of allosteric activation by phosphorylated ubiquitin binding to the enzymatic core of Parkin, which lacks the antagonistic Ubl domain. In this way, the allosteric effects of the agonist phosphorylated ubiquitin can be isolated. Using native-state hydrogen exchange monitored by mass spectrometry, we find that the five structural domains of the core of Parkin are energetically distinct. Nevertheless, association of phosphorylated ubiquitin destabilizes structural elements that bind the ubiquitin-like domain antagonist while promoting the dissociation of the catalytic domain and energetically poises the protein for transition to the fully activated structure.

Factor defining the effects of tetraalkylammonium chloride on stability, folding, and dynamics of horse cytochrome c

This article describes the molecular mechanism by which tetraalkylammonium chloride ([R4N]Cl: R- = methyl (Me), ethyl (Et), propyl (Pr),butyl (Bu)) modulates the stability, folding, and dynamics of cytochrome c (Cyt c). Analysis of [R4N]Cl effects on thermal/chemical denaturations, millisecond refolding/unfolding kinetics, and slow CO-association kinetics of Cyt c without and with denaturant providing some significant results: (i) [R4N]Cl decreasing the unfolding free energy estimated by thermodynamic and kinetic analysis of thermal/chemical denaturation curves and kinetic chevrons (Log kobs-[GdmCl]) of Cyt c, respectively (ii) hydrophobicity of R-group of [R4N]Cl, preferential inclusion of [R4N]Cl at the protein surface, and destabilizing enthalpic attractive interactions of [Me4N]Cl and steric entropic interactions of [Et4N]Cl,[Pr4N]Cl and [Bu4N]Cl with protein contribute to [R4N]Cl-induced decrease thermodynamic stability of Cyt c (iii) [R4N]Cl exhibits an additive effect with denaturant to decrease thermodynamic stability and refolding rates of Cyt c (iv) low concentrations of [R4N]Cl (≤ 0.5 M) constrain the motional dynamics while the higher concentrations (>0.75 M [R4N]Cl) enhance the structural-fluctuations that denture protein (v) hydrophobicity of R-group of [R4N]Cl alters the [denaturant]-dependent conformational stability, refolding-unfolding kinetics, and CO-association kinetics of Cyt c. Furthermore, the MD simulations depicted that the addition of 1.0 M of [R4N]Cl increased the conformational fluctuations in Cyt c leading to decreased structural stability in the order [Me4N]Cl < [Et4N]Cl < [Pr4N]Cl < [Bu4N]Cl consistent with the experimental results.

Backbone Modification in a Protein Hydrophobic Core

Targeted protein backbone modification can recreate tertiary structures reminiscent of folds found in nature on artificial scaffolds with improved biostability. Incorporation of altered monomers in such entities is typically limited to sites distant from the hydrophobic core to avoid potential disruptions to folding. This is limiting, as it is advantageous in some applications to incorporate artificial connectivity at buried sites. Here, we report an examination of protein backbone modification targeted specifically to hydrophobic core positions and its impacts on tertiary folded structure and fold stability. Different artificial monomer types are placed at core, core-flanking, or solvent-exposed positions in a compact three-helix protein. Effects on structure and folding energetics are assessed by NMR spectroscopy and biophysical methods. Results show that artificial residues can be well accommodated in the hydrophobic core of a defined tertiary fold, with effects on stability only modestly larger than identical changes at solvent-exposed sites. Collectively, these results provide new insights into folding behavior of protein-like artificial chains as well as strategies for the design of such molecules.

Unraveling the impact of the p.R107L mutation on the structure and function of human αB-Crystallin: Implications for cataract formation

To date, several pathogenic mutations have been identified in the primary structure of human α-Crystallin, frequently involving the substitution of arginine with a different amino acid. These mutations can lead to the incidence of cataracts and myopathy. Recently, an important cataract-associated mutation has been reported in the functional α-Crystallin domain (ACD) of human αB-Crystallin protein, where arginine 107 (R107) is replaced by a leucine. In this study, we investigated the structure, chaperone function, stability, oligomerization, and amyloidogenic properties of the p.R107L human αB-Crystallin using a number of different techniques. Our results suggest that the p.R107L mutation can cause significant changes in the secondary, tertiary, and quaternary structures of αB-Crystallin. This cataractogenic mutation led to the formation of protein oligomers with larger sizes than the wild-type protein and reduced the chemical and thermal stability of the mutant chaperone. Both fluorescence and microscopic assessments indicated that this mutation significantly altered the amyloidogenic properties of human αB-Crystallin. Furthermore, the mutant protein indicated an attenuated in vitro chaperone activity. The molecular dynamics (MD) simulation confirmed the experimental results and indicated that p.R107L mutation could alter the proper conformation of human αB-Crystallin dimers. In summary, our results indicated that the p.R107L mutation could promote the formation of larger oligomers, diminish the stability and chaperone activity of human αB-Crystallin, and these changes, in turn, can play a crucial role in the development of cataract disorder.

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

Solvent concentration at 50% protein unfolding may reform enzyme stability ranking and process window identification

As water miscible organic co-solvents are often required for enzyme reactions to improve e.g., the solubility of the substrate in the aqueous medium, an enzyme is required which displays high stability in the presence of this co-solvent. Consequently, it is of utmost importance to identify the most suitable enzyme or the appropriate reaction conditions. Until now, the melting temperature is used in general as a measure for stability of enzymes. The experiments here show, that the melting temperature does not correlate to the activity observed in the presence of the solvent. As an alternative parameter, the concentration of the co-solvent at the point of 50% protein unfolding at a specific temperature T in shortc U 50 T is introduced. Analyzing a set of ene reductases,c U 50 T is shown to indicate the concentration of the co-solvent where also the activity of the enzyme drops fastest. Comparing possible rankings of enzymes according to melting temperature andc U 50 T reveals a clearly diverging outcome also depending on the specific solvent used. Additionally, plots ofc U 50 versus temperature enable a fast identification of possible reaction windows to deduce tolerated solvent concentrations and temperature.

Encapsulation of Transketolase into In Vitro-Assembled Protein Nanocompartments Improves Thermal Stability

Protein compartments offer definitive structures with a large potential design space that are of particular interest for green chemistry and therapeutic applications. One family of protein compartments, encapsulins, are simple prokaryotic nanocompartments that self-assemble from a single monomer into selectively permeable cages of between 18 and 42 nm. Over the past decade, encapsulins have been developed for a diverse application portfolio utilizing their defined cargo loading mechanisms and repetitive surface display. Although it has been demonstrated that encapsulation of non-native cargo proteins provides protection from protease activity, the thermal effects arising from enclosing cargo within encapsulins remain poorly understood. This study aimed to establish a methodology for loading a reporter protein into thermostable encapsulins to determine the resulting stability change of the cargo. Building on previous in vitro reassembly studies, we first investigated the effectiveness of in vitro reassembly and cargo-loading of two size classes of encapsulins Thermotoga maritima T = 1 and Myxococcus xanthus T = 3, using superfolder Green Fluorescent Protein. We show that the empty T. maritima capsid reassembles with higher yield than the M. xanthus capsid and that in vitro loading promotes the formation of the M. xanthus T = 3 capsid form over the T = 1 form, while overloading with cargo results in malformed T. maritima T = 1 encapsulins. For the stability study, a Förster resonance energy transfer (FRET)-probed industrially relevant enzyme cargo, transketolase, was then loaded into the T. maritima encapsulin. Our results show that site-specific orthogonal FRET labels can reveal changes in thermal unfolding of encapsulated cargo, suggesting that in vitro loading of transketolase into the T. maritima T = 1 encapsulin shell increases the thermal stability of the enzyme. This work supports the move toward fully harnessing structural, spatial, and functional control of in vitro assembled encapsulins with applications in cargo stabilization.

Electrostatics introduce a trade-off between mesophilic stability and adaptation in halophilic proteins

Extremophile organisms have adapted to extreme physicochemical conditions. Halophilic organisms, in particular, survive at very high salt concentrations. To achieve this, they have engineered the surface of their proteins to increase the number of short, polar and acidic amino acids, while decreasing large, hydrophobic and basic residues. While these adaptations initially decrease protein stability in the absence of salt, they grant halophilic proteins remarkable stability in environments with extremely high salt concentrations, where non-adapted proteins unfold and aggregate. The molecular mechanisms by which halophilic proteins achieve this, however, are not yet clear. Here, we test the hypothesis that the halophilic amino acid composition destabilizes the surface of the protein, but in exchange improves the stability in the presence of salts. To do that, we have measured the folding thermodynamics of various protein variants with different degrees of halophilicity in the absence and presence of different salts, and at different pH values to tune the ionization state of the acidic amino acids. Our results show that halophilic amino acids decrease the stability of halophilic proteins under mesophilic conditions, but in exchange improve salt-induced stabilization and solubility. We also find that, in contrast to traditional assumptions, contributions arising from hydrophobic effect and preferential ion exclusion are more relevant for haloadaptation than electrostatics. Overall, our findings suggest a trade-off between folding thermodynamics and halophilic adaptation to optimize proteins for hypersaline environments.

Nascent chains derived from a foldable protein sequence interact with specific ribosomal surface sites near the exit tunnel

In order to become bioactive, proteins must be translated and protected from aggregation during biosynthesis. The ribosome and molecular chaperones play a key role in this process. Ribosome-bound nascent chains (RNCs) of intrinsically disordered proteins and RNCs bearing a signal/arrest sequence are known to interact with ribosomal proteins. However, in the case of RNCs bearing foldable protein sequences, not much information is available on these interactions. Here, via a combination of chemical crosslinking and time-resolved fluorescence-anisotropy, we find that nascent chains of the foldable globin apoHmp1-140 interact with ribosomal protein L23 and have a freely-tumbling non-interacting N-terminal compact region comprising 63-94 residues. Longer RNCs (apoHmp1-189) also interact with an additional yet unidentified ribosomal protein, as well as with chaperones. Surprisingly, the apparent strength of RNC/r-protein interactions does not depend on nascent-chain sequence. Overall, foldable nascent chains establish and expand interactions with selected ribosomal proteins and chaperones, as they get longer. These data are significant because they reveal the interplay between independent conformational sampling and nascent-protein interactions with the ribosomal surface.

Quantum Defect Sensitization via Phase-Changing Supercharged Antibody Fragments

Quantum defects in single-walled carbon nanotubes promote exciton localization, which enables potential applications in biodevices and quantum light sources. However, the effects of local electric fields on the emissive energy states of quantum defects and how they can be controlled are unexplored. Here, we investigate quantum defect sensitization by engineering an intrinsically disordered protein to undergo a phase change at a quantum defect site. We designed a supercharged single-chain antibody fragment (scFv) to enable a full ligand-induced folding transition from an intrinsically disordered state to a compact folded state in the presence of a cytokine. The supercharged scFv was conjugated to a quantum defect to induce a substantial local electric change upon ligand binding. Employing the detection of a proinflammatory biomarker, interleukin-6, as a representative model system, supercharged scFv-coupled quantum defects exhibited robust fluorescence wavelength shifts concomitant with the protein folding transition. Quantum chemical simulations suggest that the quantum defects amplify the optical response to the localization of charges produced upon the antigen-induced folding of the proteins, which is difficult to achieve in unmodified nanotubes. These findings portend new approaches to modulate quantum defect emission for biomarker sensing and protein biophysics and to engineer proteins to modulate binding signal transduction.

Mimicking kidney flow shear efficiently induces aggregation of LECT2, a protein involved in renal amyloidosis

Aggregation of leukocyte cell-derived chemotaxin 2 (LECT2) causes ALECT2, a systemic amyloidosis that affects the kidney and liver. Previous studies established that LECT2 fibrillogenesis is accelerated by the loss of its bound zinc ion and stirring/shaking. These forms of agitation create heterogeneous shear conditions, including air-liquid interfaces that denature proteins, that are not present in the body. Here, we determined the extent to which a more physiological form of mechanical stress-shear generated by fluid flow through a network of narrow channels-drives LECT2 fibrillogenesis. To mimic blood flow through the kidney, where LECT2 and other proteins form amyloid deposits, we developed a microfluidic device consisting of progressively branched channels narrowing from 5 mm to 20 μm in width. Shear was particularly pronounced at the branch points and in the smallest capillaries. Aggregation was induced within 24 h by shear levels that were in the physiological range and well below those required to unfold globular proteins such as LECT2. EM images suggested the resulting fibril ultrastructures were different when generated by laminar flow shear versus shaking/stirring. Importantly, results from the microfluidic device showed the first evidence that the I40V mutation accelerated fibril formation and increased both the size and the density of the aggregates. These findings suggest that kidney-like flow shear, in combination with zinc loss, acts in combination with the I40V mutation to trigger LECT2 amyloidogenesis. These microfluidic devices may be of general use for uncovering mechanisms by which blood flow induces misfolding and amyloidosis of circulating proteins.

Bacterial Chaperone Domain Insertions Convert Human FKBP12 into an Excellent Protein-Folding Catalyst-A Structural and Functional Analysis

Many folding enzymes use separate domains for the binding of substrate proteins and for the catalysis of slow folding reactions such as prolyl isomerization. FKBP12 is a small prolyl isomerase without a chaperone domain. Its folding activity is low, but it could be increased by inserting the chaperone domain from the homolog SlyD of E. coli near the prolyl isomerase active site. We inserted two other chaperone domains into human FKBP12: the chaperone domain of SlpA from E. coli, and the chaperone domain of SlyD from Thermococcus sp. Both stabilized FKBP12 and greatly increased its folding activity. The insertion of these chaperone domains had no influence on the FKBP12 and the chaperone domain structure, as revealed by two crystal structures of the chimeric proteins. The relative domain orientations differ in the two crystal structures, presumably representing snapshots of a more open and a more closed conformation. Together with crystal structures from SlyD-like proteins, they suggest a path for how substrate proteins might be transferred from the chaperone domain to the prolyl isomerase domain.

Degradation versus fibrillogenesis, two alternative pathways modulated by seeds and glycosaminoglycans

The mechanism that converts native human transthyretin into amyloid fibrils in vivo is still a debated and controversial issue. Commonly, non-physiological conditions of pH, temperature, or organic solvents are used in in vitro models of fibrillogenesis of globular proteins. Transthyretin amyloid formation can be achieved under physiological conditions through a mechano-enzymatic mechanism involving specific serine proteases such as trypsin or plasmin. Here, we investigate S52P and L111M transthyretin variants, both causing a severe form of systemic amyloidosis mostly targeting the heart at a relatively young age with heterogeneous phenotype among patients. Our studies on thermodynamics show that both proteins are significantly less stable than other amyloidogenic variants. However, despite a similar thermodynamic stability, L111M variant seems to have enhanced susceptibility to cleavage and a lower tendency to form fibrils than S52P in the presence of specific proteases and biomechanical forces. Heparin strongly enhances the fibrillogenic capacity of L111M transthyretin, but has no effect on the S52P variant. Fibrillar seeds similarly affect the fibrillogenesis of both proteins, with a stronger effect on the L111M variant. According to our model of mechano-enzymatic fibrillogenesis, both full-length and truncated monomers, released after the first cleavage, can enter into fibrillogenesis or degradation pathways. Our findings show that the kinetics of the two processes can be affected by several factors, such as intrinsic amyloidogenicity due to the specific mutations, environmental factors including heparin and fibrillar seeds that significantly accelerate the fibrillogenic pathway.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

Thermodynamic characterization of amyloid polymorphism by microfluidic transient incomplete separation

Amyloid fibrils of proteins such as α-synuclein are a hallmark of neurodegenerative diseases and much research has focused on their kinetics and mechanisms of formation. The question as to the thermodynamic stability of such structures has received much less attention. Here, we newly utilize the principle of transient incomplete separation of species in laminar flow in combination with chemical depolymerization for the quantification of amyloid fibril stability. The relative concentrations of fibrils and monomer at equilibrium are determined through an in situ separation of these species based on their different diffusivity inside a microfluidic capillary. The method is highly sample economical, using much less than a microliter of sample per data point and its only requirement is the presence of aromatic residues (W, Y) because of its label-free nature, which makes it widely applicable. Using this method, we investigate the differences in thermodynamic stability between different fibril polymorphs of α-synuclein and quantify these differences for the first time. Importantly, we show that fibril formation can be under kinetic or thermodynamic control and that a change in solution conditions can both stabilise and destabilise amyloid fibrils. Taken together, our results establish the thermodynamic stability as a well-defined and key parameter that can contribute towards a better understanding of the physiological roles of amyloid fibril polymorphism.

Sequence identification and in silico characterization of novel thermophilic lipases from Geobacillus species

Microbial lipases are utilized in various biotechnological areas, including pharmaceuticals, food, biodiesel, and detergents. In this study, we cloned and sequenced Lip21 and Lip33 genes from Geobacillus sp. GS21 and Geobacillus sp. GS33, then we in silico and experimentally analyzed the encoded lipases. For this purpose, Lip21 and Lip33 were cloned, sequenced, and their amino acid sequences were investigated for determination of biophysicochemical characteristics, evolutionary relationships, and sequence similarities. 3D models were built and computationally affirmed by various bioinformatics tools, and enzyme-ligand interactions were investigated by docking analysis using six ligands. Biophysicochemical property of Lip21 and Lip33 was also determined experimentally and the results demonstrated that they had similar isoelectric point (pI) (6.21) and Tm (75.5°C) values as Tm was revealed by denatured protein analysis of the circular dichroism spectrum and pI was obtained by isoelectric focusing. Phylogeny analysis indicated that Lip21 and Lip33 were the closest to lipases from Geobacillus sp. SBS-4S and Geobacillus thermoleovorans, respectively. Alignment analysis demonstrated that S144-D348-H389 was catalytic triad residues in Lip21 and Lip33, and enzymes possessed a conserved Gly-X-Ser-X-Gly motif containing catalytic serine. 3D structure analysis indicated that Lip21 and Lip33 highly resembled each other and they were α/β hydrolase-fold enzymes with large lid domains. BANΔIT analysis results showed that Lip21 and Lip33 had higher thermal stability, compared to other thermostable Geobacillus lipases. Docking results revealed that Lip21- and Lip33-docked complexes possessed common residues (H112, K115, Q162, E163, and S141) that interacted with the substrates, except paranitrophenyl (pNP)-C10 and pNP-C12, indicating that these residues might have a significant action on medium and short-chain fatty acid esters. Thus, Lip21 and Lip33 can be potential candidates for different industrial applications.

Design of complicated all-α protein structures

A wide range of de novo protein structure designs have been achieved, but the complexity of naturally occurring protein structures is still far beyond these designs. Here, to expand the diversity and complexity of de novo designed protein structures, we sought to develop a method for designing 'difficult-to-describe' α-helical protein structures composed of irregularly aligned α-helices like globins. Backbone structure libraries consisting of a myriad of α-helical structures with five or six helices were generated by combining 18 helix-loop-helix motifs and canonical α-helices, and five distinct topologies were selected for de novo design. The designs were found to be monomeric with high thermal stability in solution and fold into the target topologies with atomic accuracy. This study demonstrated that complicated α-helical proteins are created using typical building blocks. The method we developed will enable us to explore the universe of protein structures for designing novel functional proteins.

Differential effects of disulfide bond formation in TEM-1 versus CTX-M-9 β-lactamase

To investigate how disulfide bonds can impact protein energy landscapes, we surveyed the effects of adding or removing a disulfide in two β-lactamase enzymes, TEM-1 and CTX-M-9. The homologs share a structure and 38% sequence identity, but only TEM-1 contains a native disulfide bond. They also differ in thermodynamic stability and in the number of states populated at equilibrium: CTX-M-9 is two-state whereas TEM-1 has an additional intermediate state. We hypothesized that the disulfide bond is the major underlying determinant for these observed differences in their energy landscapes. To test this, we removed the disulfide bridge from TEM-1 and introduced a disulfide bridge at the same location in CTX-M-9. This modest change to sequence modulates the stabilities-and therefore populations-of TEM-1's equilibrium states and, more surprisingly, creates a novel third state in CTX-M-9. Unlike TEM-1's partially folded intermediate, this third state is a higher-order oligomer with reduced cysteines that retains the native fold and is fully active. Sub-denaturing concentrations of urea shifts the equilibrium to the monomeric form, allowing the disulfide bond to form. Interestingly, comparing the stability of the oxidized monomer with a variant lacking cysteines reveals the disulfide is neither stabilizing nor destabilizing in CTX-M-9, in contrast with the observed stabilization in TEM-1. Thus, we can conclude that engineering disulfide bonds is not always an effective stabilization strategy even when analogous disulfides exist in more stable structural homologs. This study also illustrates how homo-oligomerization can result from a small number of mutations, suggesting complex formation might be easily accessed during a protein family's evolution.

Single-molecule analysis of osmolyte-mediated nanomechanical unfolding behavior of a protein domain

The small organic molecules, known as osmolytes being ubiquitously present in different cell types, affect protein folding, stability and aggregation. However, it is unknown how the osmolytes affect the nanomechanical unfolding behavior of protein domain. Here, we show the osmolyte-dependent mechanical unfolding properties of protein titin immunoglobulin-27 (I27) domain using an atomic force microscopy (AFM)-based single-molecule force spectroscopy. We found that amines and methylamines improved the mechanical stability of I27 domain, whereas polyols had no effect. Interestingly, glycine betaine (GB) or trimethylamine-N-oxide (TMAO) increased the average unfolding force of the protein domain. The kinetic parameters analyzed at single-molecule level reveal that stabilizing effect of osmolytes is due to a decrease in the unfolding rate constant of I27, which was confirmed by molecular dynamics simulations. Our study reveals different effects that diverse osmolytes have on the mechanical properties of the protein, and suggests the potential use of osmolytes in modulating the mechanical stability of proteins required for various nano-biotechnological applications.

Engineering neuroglobin nitrite reductase activity based on myoglobin models

Neuroglobin is a hemoprotein expressed in several nervous system cell lineages with yet unknown physiological functions. Neuroglobin presents a very similar structure to that of the related globins hemoglobin and myoglobin, but shows an hexacoordinate heme as compared to the pentacoordinated heme of myoglobin and hemoglobin. While several reactions of neuroglobin have been characterized in vitro, the relative importance of most of those reactions in vivo is yet undefined. Neuroglobin, like other heme proteins, can reduce nitrite to nitric oxide, providing a possible route to generate nitric oxide in vivo in low oxygen conditions. The reaction kinetics are highly dependent on the nature of the distal residue, and replacement of the distal histidine His64(E7) can increase the reaction rate constants by several orders of magnitude. However, mutation of other distal pocket positions such as Phe28(B10) or Val68(E11) has more limited impact on the rates. Computational analysis using myoglobin as template, guided by the structure of dedicated nitrite reductases like cytochrome cd1 nitrite reductase, has pointed out that combined mutations of the residues B10 and CD1 could increase the nitrite reductase activity of myoglobin, by mimicking the environment of the distal heme pocket in cytochrome cd1 nitrite reductase. As neuroglobin shows high sequence and structural homology with myoglobin, we hypothesized that such mutations (F28H and F42Y in neuroglobin) could also modify the nitrite reductase activity of neuroglobin. Here we study the effect of these mutations. Unfortunately, we do not observe in any case an increase in the nitrite reduction rates. Our results provide some further indications of nitrite reductase regulation in neuroglobin and highlight the minor but critical differences between the structure of penta- and hexacoordinate globins.

Amyloid Aggregation Is Potently Slowed Down by Osmolytes Due to Compaction of Partially Folded State

Amyloid aggregation is a key process in amyloidoses and neurodegenerative diseases. Hydrophobicity is one of the major driving forces for this type of aggregation, as an increase in hydrophobicity generally correlates with aggregation susceptibility and rate. However, most experimental systems in vitro and prediction tools in silico neglect the contribution of protective osmolytes present in the cellular environment. Here, we assessed the role of hydrophobic mutations in amyloid aggregation in the presence of osmolytes. To achieve this goal, we used the model protein human muscle acylphosphatase (mAcP) and mutations to leucine that increased its hydrophobicity without affecting its thermodynamic stability. Osmolytes significantly slowed down the aggregation kinetics of the hydrophobic mutants, with an effect larger than that observed on the wild-type protein. The effect increased as the mutation site was closer to the middle of the protein sequence. We propose that the preferential exclusion of osmolytes from mutation-introduced hydrophobic side-chains quenches the aggregation potential of the ensemble of partially unfolded states of the protein by inducing its compaction and inhibiting its self-assembly with other proteins. Our results suggest that including the effect of the cellular environment in experimental setups and predictive softwares, for both mechanistic studies and drug design, is essential in order to obtain a more complete combination of the driving forces of amyloid aggregation.

The variable domain from dynamin-related protein 1 promotes liquid-liquid phase separation that enhances its interaction with cardiolipin-containing membranes

Dynamins are an essential superfamily of mechanoenzymes that remodel membranes and often contain a "variable domain" important for regulation. For the mitochondrial fission dynamin, dynamin-related protein 1, a regulatory role for the variable domain (VD) is demonstrated by gain- and loss-of-function mutations, yet the basis for this is unclear. Here, the isolated VD is shown to be intrinsically disordered and undergo a cooperative transition in the stabilizing osmolyte trimethylamine N-oxide. However, the osmolyte-induced state is not folded and surprisingly appears as a condensed state. Other co-solutes including known molecular crowder Ficoll PM 70, also induce a condensed state. Fluorescence recovery after photobleaching experiments reveal this state to be liquid-like indicating the VD undergoes a liquid-liquid phase separation under crowding conditions. These crowding conditions also enhance binding to cardiolipin, a mitochondrial lipid, which appears to promote phase separation. Since dynamin-related protein 1 is found assembled into discrete punctate structures on the mitochondrial surface, the inference from the present work is that these structures might arise from a condensed state involving the VD that may enable rapid tuning of mechanoenzyme assembly necessary for fission.

Domain swap facilitates structural transitions of spider silk protein C-terminal domains

Domain swap is a mechanism of protein dimerization where the two interacting domains exchange parts of their structure. Web spiders make use of the process in the connection of C-terminal domains (CTDs) of spidroins, the soluble protein building blocks that form tough silk fibers. Besides providing connectivity and solubility, spidroin CTDs are responsible for inducing structural transitions during passage through an acidified assembly zone within spinning ducts. The underlying molecular mechanisms are elusive. Here, we studied the folding of five homologous spidroin CTDs from different spider species or glands. Four of these are domain-swapped dimers formed by five-helix bundles from spidroins of major and minor ampullate glands. The fifth is a dimer that lacks domain swap, formed by four-helix bundles from a spidroin of a flagelliform gland. Spidroins from this gland do not undergo structural transitions whereas the others do. We found a three-state mechanism of folding and dimerization that was conserved across homologues. In chemical denaturation experiments the native CTD dimer unfolded to a dimeric, partially structured intermediate, followed by full unfolding to denatured monomers. The energetics of the individual folding steps varied between homologues. Contrary to the common belief that domain swap stabilizes protein assemblies, the non-swapped homologue was most stable and folded four orders of magnitude faster than a swapped variant. Domain swap of spidroin CTDs induces an entropic penalty to the folding of peripheral helices, thus unfastening them for acid-induced unfolding within a spinning duct, which primes them for refolding into alternative structures during silk formation.

The complex impact of cancer-related missense mutations on the stability and on the biophysical and biochemical properties of MAPK1 and MAPK3 somatic variants

Mitogen-activated protein kinases 1 and 3 (MAPK1 and MAPK3), also called extracellular regulated kinases (ERK2 and ERK1), are serine/threonine kinase activated downstream by the Ras/Raf/MEK/ERK signal transduction cascade that regulates a variety of cellular processes. A dysregulation of MAPK cascade is frequently associated to missense mutations on its protein components and may be related to many pathologies, including cancer. In this study we selected from COSMIC database a set of MAPK1 and MAPK3 somatic variants found in cancer tissues carrying missense mutations distributed all over the MAPK1 and MAPK3 sequences. The proteins were expressed as pure recombinant proteins, and their biochemical and biophysical properties have been studied in comparison with the wild type. The missense mutations lead to changes in the tertiary arrangements of all the variants. The thermodynamic stability of the wild type and variants has been investigated in the non-phosphorylated and in the phosphorylated form. Significant differences in the thermal stabilities of most of the variants have been observed, as well as changes in the catalytic efficiencies. The energetics of the catalytic reaction is affected for all the variants for both the MAPK proteins. The stability changes and the variation in the enzyme catalysis observed for most of MAPK1/3 variants suggest that a local change in a residue, distant from the catalytic site, may have long-distance effects that reflect globally on enzyme stability and functions.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

An automated single-molecule FRET platform for high-content, multiwell plate screening of biomolecular conformations and dynamics

Single-molecule FRET (smFRET) has become a versatile tool for probing the structure and functional dynamics of biomolecular systems, and is extensively used to address questions ranging from biomolecular folding to drug discovery. Confocal smFRET measurements are amongst the widely used smFRET assays and are typically performed in a single-well format. Thus, sampling of many experimental parameters is laborious and time consuming. To address this challenge, we extend here the capabilities of confocal smFRET beyond single-well measurements by integrating a multiwell plate functionality to allow for continuous and automated smFRET measurements. We demonstrate the broad applicability of the multiwell plate assay towards DNA hairpin dynamics, protein folding, competitive and cooperative protein-DNA interactions, and drug-discovery, revealing insights that would be very difficult to achieve with conventional single-well format measurements. For the adaptation into existing instrumentations, we provide a detailed guide and open-source acquisition and analysis software.

Chemical Protein Unfolding - A Simple Cooperative Model

Chemical unfolding with guanidineHCl or urea is a common method to study the conformational stability of proteins. The analysis of unfolding isotherms is usually performed with an empirical linear extrapolation method (LEM). A large positive free energy is assigned to the native protein, which is usually considered to be a minimum of the free energy. The method thus contradicts common expectations. Here, we present a multistate cooperative model that addresses specifically the binding of the denaturant to the protein and the cooperativity of the protein unfolding equilibrium. The model is based on a molecular statistical-mechanical partition function of the ensemble, but simple solutions for the calculation of the binding isotherm and the associated free energy are presented. The model is applied to 23 published unfolding isotherms of small and large proteins. For a given denaturant, the binding constant depends on temperature and pH but shows little protein specificity. Chemical unfolding is less cooperative than thermal unfolding. The cooperativity parameter σ is at least 2 orders of magnitude larger than that of thermal unfolding. The multistate cooperative model predicts zero free energy for the native protein, which becomes strongly negative beyond the midpoint concentration of unfolding. The free energy to unfold a cooperative unit corresponds exactly to the diffusive energy of the denaturant concentration gradient necessary for unfolding. The temperature dependence of unfolding isotherms yields the denaturant-induced unfolding entropy and, in turn, the unfolding enthalpy. The multistate cooperative model provides molecular insight and is as simple to apply as the LEM but avoids the conceptual difficulties of the latter.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Mimicking Kidney Flow Shear Efficiently Induces Aggregation of LECT2, a Protein Involved in Renal Amyloidosis

Aggregation of leukocyte cell-derived chemotaxin 2 (LECT2) causes ALECT2, a systemic amyloidosis that affects the kidney and liver. Homozygosity of the I40V LECT2 mutation is believed to be necessary but not sufficient for the disease. Previous studies established that LECT2 fibrillogenesis is greatly accelerated by loss of its single bound zinc ion and stirring or shaking. These forms of agitation are often used to facilitate protein aggregation, but they create heterogeneous shear conditions, including air-liquid interfaces that denature proteins, that are not present in the body. Here, we determined the extent to which a more physiological form of mechanical stress-shear generated by fluid flow through a network of artery and capillary-sized channels-drives LECT2 fibrillogenesis. To mimic blood flow through the human kidney, where LECT2 and other proteins form amyloid deposits, we developed a microfluidic device consisting of progressively branched channels narrowing from 5 mm to 20 μm in width. Flow shear was particularly pronounced at the branch points and in the smallest capillaries, and this induced LECT2 aggregation much more efficiently than conventional shaking methods. EM images suggested the resulting fibril structures were different in the two conditions. Importantly, results from the microfluidic device showed the first evidence that the I40V mutation accelerated fibril formation and increased both size and density of the aggregates. These findings suggest that kidney-like flow shear, in combination with zinc loss, acts in combination with the I40V mutation to trigger LECT2 amyloidogenesis. These microfluidic devices may be of general use for uncovering the mechanisms by which blood flow induces misfolding and amyloidosis of circulating proteins.

Binding to the conserved and stably folded guide RNA pseudoknot induces Cas12a conformational changes during ribonucleoprotein assembly

Ribonucleoproteins (RNPs) comprise one or more RNA and protein molecules that interact to form a stable complex, which commonly involves conformational changes in the more flexible RNA components. Here, we propose that Cas12a RNP assembly with its cognate CRISPR RNA (crRNA) guide instead proceeds primarily through Cas12a conformational changes during binding to more stable, prefolded crRNA 5' pseudoknot handles. Phylogenetic reconstructions and sequence and structure alignments revealed that the Cas12a proteins are divergent in sequence and structure while the crRNA 5' repeat region, which folds into a pseudoknot and anchors binding to Cas12a, is highly conserved. Molecular dynamics simulations of three Cas12a proteins and their cognate guides revealed substantial flexibility for unbound apo-Cas12a. In contrast, crRNA 5' pseudoknots were predicted to be stable and independently folded. Limited trypsin hydrolysis, differential scanning fluorimetry, thermal denaturation, and CD analyses supported conformational changes of Cas12a during RNP assembly and an independently folded crRNA 5' pseudoknot. This RNP assembly mechanism may be rationalized by evolutionary pressure to conserve CRISPR loci repeat sequence, and therefore guide RNA structure, to maintain function across all phases of the CRISPR defense mechanism.

Enhanced Thermal Stability and Reduced Aggregation in an Antibody Fab Fragment at Elevated Concentrations

The aggregation of protein therapeutics such as antibodies remains a major challenge in the biopharmaceutical industry. The present study aimed to characterize the impact of the protein concentration on the mechanisms and potential pathways for aggregation, using the antibody Fab fragment A33 as the model protein. Aggregation kinetics were determined for 0.05 to 100 mg/mL Fab A33, at 65 °C. A surprising trend was observed whereby increasing the concentration decreased the relative aggregation rate, ln(v) (% day^-1), from 8.5 at 0.05 mg/mL to 4.4 at 100 mg/mL. The absolute aggregation rate (mol L^-1 h^-1) increased with the concentration following a rate order of approximately 1 up to a concentration of 25 mg/mL. Above this concentration, there was a transition to an apparently negative rate order of -1.1 up to 100 mg/mL. Several potential mechanisms were examined as possible explanations. A greater apparent conformational stability at 100 mg/mL was observed from an increase in the thermal transition midpoint (T_m) by 7-9 °C, relative to those at 1-4 mg/mL. The associated change in unfolding entropy (△S_vh) also increased by 14-18% at 25-100 mg/mL, relative to those at 1-4 mg/mL, indicating reduced conformational flexibility in the native ensemble. Addition of Tween or the crowding agents Ficoll and dextran, showed that neither surface adsorption, diffusion limitations nor simple volume crowding affected the aggregation rate. Fitting of kinetic data to a wide range of mechanistic models implied a reversible two-state conformational switch mechanism from aggregation-prone monomers (N*) into non-aggregating native forms (N) at higher concentrations. k_D measurements from DLS data also suggested a weak self-attraction while remaining colloidally stable, consistent with macromolecular self-crowding within weakly associated reversible oligomers. Such a model is also consistent with compaction of the native ensemble observed through changes in T_m and △S_vh.

Collagen-like Motifs of SasG: A Novel Fold for Protein Mechanical Strength

The Staphylococcus aureus surface protein G (SasG) is associated with host colonisation and biofilm formation. As colonisation occurs at the liquid-substrate interface bacteria are subject to a myriad of external forces and, presumably as a consequence, SasG displays extreme mechanical strength. This mechanical phenotype arises from the B-domain; a repetitive region composed of alternating E and G5 subdomains. These subdomains have an unusual structure comprising collagen-like regions capped by triple-stranded β-sheets. To identify the determinants of SasG mechanical strength, we characterised the mechanical phenotype and thermodynamic stability of 18 single substitution variants of a pseudo-wildtype protein. Visualising the mechanically-induced transition state at a residue-level by ϕ-value analysis reveals that the main force-bearing regions are the N- and C-terminal 'Mechanical Clamps' and their side-chain interactions. This is tailored by contacts at the pseudo-hydrophobic core interface. We also describe a novel mechanical motif - the collagen-like region and show that glycine to alanine substitutions, analogous to those found in Osteogenesis Imperfecta (brittle bone disease), result in a significantly reduced mechanical strength.

Enhancing response of a protein conformational switch by using two disordered ligand binding domains

Introduction: Protein conformational switches are often constructed by fusing an input domain, which recognizes a target ligand, to an output domain that establishes a biological response. Prior designs have employed binding-induced folding of the input domain to drive a conformational change in the output domain. Adding a second input domain can in principle harvest additional binding energy for performing useful work. It is not obvious, however, how to fuse two binding domains to a single output domain such that folding of both binding domains combine to effect conformational change in the output domain. Methods: Here, we converted the ribonuclease barnase (Bn) to a switchable enzyme by duplicating a C-terminal portion of its sequence and appending it to its N-terminus, thereby establishing a native fold (OFF state) and a circularly permuted fold (ON state) that competed for the shared core in a mutually exclusive fashion. Two copies of FK506 binding protein (FKBP), both made unstable by the V24A mutation and one that had been circularly permuted, were inserted into the engineered barnase at the junctions between the shared and duplicated sequences. Results: Rapamycin-induced folding of FK506 binding protein stretched and unfolded the native fold of barnase via the mutually exclusive folding effect, and rapamycin-induced folding of permuted FK506 binding protein stabilized the permuted fold of barnase by the loop-closure entropy principle. These folding events complemented each other to turn on RNase function. The cytotoxic switching mechanism was validated in yeast and human cells, and in vitro with purified protein. Discussion: Thermodynamic modeling and experimental results revealed that the dual action of loop-closure entropy and mutually exclusive folding is analogous to an engine transmission in which loop-closure entropy acts as the low gear, providing efficient switching at low ligand concentrations, and mutually exclusive folding acts as the high gear to allow the switch to reach its maximum response at high ligand concentrations.

Assessing and Engineering Antibody Stability Using Experimental and Computational Methods

Engineering increased stability into antibodies can improve their developability. While a range of properties need to be optimized, thermal stability and aggregation are two key factors that affect the antibody yield, purity, and specificity throughout the development and manufacturing pipeline. Therefore, an ideal goal would be to apply protein engineering methods early-on, such as in parallel to affinity maturation, to screen out potential drug molecules with the desired conformational and colloidal stability. This chapter introduces our methods to computationally characterize an antibody Fab fragment, propose stabilizing variants, and then experimentally verify these predictions.

Influence of the Interdomain Interface on Structural and Redox Properties of Multiheme Proteins

Multiheme proteins are important in energy conversion and biogeochemical cycles of nitrogen and sulfur. A diheme cytochrome c4 (c4) was used as a model to elucidate roles of the interdomain interface on properties of iron centers in its hemes A and B. Isolated monoheme domains c4-A and c4-B, together with the full-length diheme c4 and its Met-to-His ligand variants, were characterized by a variety of spectroscopic and stability measurements. In both isolated domains, the heme iron is Met/His-ligated at pH 5.0, as in the full-length c4, but becomes His/His-ligated in c4-B at higher pH. Intradomain contacts in c4-A are minimally affected by the separation of c4-A and c4-B domains, and isolated c4-A is folded. In contrast, the isolated c4-B is partially unfolded, and the interface with c4-A guides folding of this domain. The c4-A and c4-B domains have the propensity to interact even without the polypeptide linker. Thermodynamic cycles have revealed properties of monomeric folded isolated domains, suggesting that ferrous (FeII), but not ferric (FeIII) c4-A and c4-B, is stabilized by the interface. This study illustrates the effects of the interface on tuning structural and redox properties of multiheme proteins and enriches our understanding of redox-dependent complexation.

Stability profile of vimentin rod domain

Intermediate filaments (IFs) form an essential part of the metazoan cytoskeleton. Despite a long history of research, a proper understanding of their molecular architecture and assembly process is still lacking. IFs self-assemble from elongated dimers, which are defined by their central "rod" domain. This domain forms an α-helical coiled coil consisting of three segments called coil1A, coil1B, and coil2. It has been hypothesized that the structural plasticity of the dimer, including the unraveling of some coiled-coil regions, is essential for the assembly process. To systematically explore this possibility, we have studied six 50-residue fragments covering the entire rod domain of human vimentin, a model IF protein. After creating in silico models of these fragments, their evaluation using molecular dynamics was performed. Large differences were seen across the six fragments with respect to their structural variability during a 100 ns simulation. Next, the fragments were prepared recombinantly, whereby their correct dimerization was promoted by adding short N- or C-terminal capping motifs. The capped fragments were subjected to circular dichroism measurements at varying temperatures. The obtained melting temperatures reveal the relative stabilities of individual fragments, which correlate well with in silico results. We show that the least stable regions of vimentin rod are coil1A and the first third of coil2, while the structures of coil1B and the rest of coil2 are significantly more robust. These observations are in line with the data obtained using other experimental approaches, and contribute to a better understanding of the molecular mechanisms driving IF assembly.

Hereditary spastic paraplegia SPG13 mutation increases structural stability and ATPase activity of human mitochondrial chaperonin

Human mitochondrial chaperonin mHsp60 is broadly associated with various human health conditions and the V72I mutation in mHsp60 causes a form of hereditary spastic paraplegia, a neurodegenerative disease. The main function of mHsp60 is to assist folding of mitochondrial proteins in an ATP-dependent manner. In this study, we unexpectedly found that mutant mHsp60^V72I was more stable structurally and more active in the ATPase activity than the wildtype. Analysis of our recently solved cryo-EM structure of mHsp60 revealed allosteric roles of V72I in structural stability and ATPase activity, which were supported by studies including those using the V72A mutation. Despite with the increases in structural stability and ATPase activity, mHsp60^V72I was less efficient in folding malate dehydrogenase, a putative mHsp60 substrate protein in mitochondria and also commonly used in chaperonin studies. In addition, although mHsp60^V72I along with its cochaperonin mHsp10 was able to substitute the E. coli chaperonin system in supporting cell growth under normal temperature of 37 °C, it was unable under heat shock temperature of 42 °C. Our results support the importance of structural dynamics and an optimal ATP turnover that mHsp60 has evolved for its function and physiology. We propose that unproductive energy utilization, or hyperactive ATPase activity and compromised folding function, not mutually exclusive, are responsible for the V72I pathology in neurodegenerative disease.

Folding Mechanism and Aggregation Propensity of the KH0 Domain of FMRP and Its R138Q Pathological Variant

The K-homology (KH) domains are small, structurally conserved domains found in proteins of different origins characterized by a central conserved βααβ “core” and a GxxG motif in the loop between the two helices of the KH core. In the eukaryotic KHI type, additional αβ elements decorate the “core” at the C-terminus. Proteins containing KH domains perform different functions and several diseases have been associated with mutations in these domains, including those in the fragile X mental retardation protein (FMRP). FMRP is an RNA-binding protein crucial for the control of RNA metabolism whose lack or mutations lead to fragile X syndrome (FXS). Among missense mutations, the R138Q substitution is in the KH0 degenerated domain lacking the classical GxxG motif. By combining equilibrium and kinetic experiments, we present a characterization of the folding mechanism of the KH0 domain from the FMRP wild-type and of the R138Q variant showing that in both cases the folding mechanism implies the accumulation of an on-pathway transient intermediate. Moreover, by exploiting a battery of biophysical techniques, we show that the KH0 domain has the propensity to form amyloid-like aggregates in mild conditions in vitro and that the R138Q mutation leads to a general destabilization of the protein and to an increased fibrillogenesis propensity.

Analysis of the effect of 1-Allyl-3-Methylimidazolium chloride on thermodynamic stability, folding kinetics, and motional dynamics of horse cytochrome c

1-allyl-3-methylimidazolium chloride (AMIMCl) acts as a potential green solvent for proteins. The present work provides a possible pathway by which the structural, kinetic, thermodynamic, and folding properties of horse cytochrome c (cyt c) are affected in green aqueous-AMIMCl systems. Analysis of the effect of AMIMCl on thermodynamic stability, refolding/unfolding kinetics, and motional dynamics of cyt c provided important information, (i) AMIMCl decreases the thermodynamic stability of reduced cyt c and also strengthens the guanidinium chloride (GdmCl)-mediated decrease in thermodynamic stability of protein, (ii) AMIMCl reduces the thermal-fluctuation of Met80-containing omega-loop of natively-folded compact state of carbonmonoxycytochrome c (MCO-state) due to polyfunctional interactions between the AMIM⁺ and different groups of protein, (iii) AMIMCl shifts the kinetic chevron plot, ln k_obs[GdmCl] to the lower concentration of GdmCl, (iv) AMIMCl shifts the refolding and unfolding limps to vertically downwards and upwards, respectively, and (v) AMIMCl reducing the unfolding free energy estimated by both thermodynamic and kinetic analysis.

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 ( Δ G D 0 ) 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 ( Δ G D 0 ), 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, Δ G D 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 Δ G D 0 comes from the analysis of the plot of Δ G D versus denaturant concentration, which requires a long extrapolation of values of Δ G D , and all the three methods of extrapolation give three different values of Δ G D 0 for a protein. Thus, our confidence in the authentic value of Δ G D 0 is eroded. Another problem with this in vitro measurement of Δ G D 0 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.

Dynamics of integrin α5β1, fibronectin, and their complex reveal sites of interaction and conformational change

Integrin α5β1 mediates cell adhesion to the extracellular matrix by binding fibronectin (Fn). Selectivity for Fn by α5β1 is achieved through recognition of an RGD motif in the 10th type III Fn domain (Fn10) and the synergy site in the ninth type III Fn domain (Fn9). However, details of the interaction dynamics are unknown. Here, we compared synergy-site and Fn-truncation mutations for their α5β1-binding affinities and stabilities. We also interrogated binding of the α5β1 ectodomain headpiece fragment to Fn using hydrogen-deuterium exchange (HDX) mass spectrometry to probe binding sites and sites of integrin conformational change. Our results suggest the synergistic effect of Fn9 requires both specific residues and a folded domain. We found some residues considered important for synergy are required for stability. Additionally, we show decreases in fibronectin HDX are localized to a synergy peptide containing contacting residues in two β-strands, an intervening loop in Fn9, and the RGD-containing loop in Fn10, indicative of binding sites. We also identified binding sites in the α5-subunit β-propeller domain for the Fn9 synergy site and in the β1-subunit βI domain for Fn10 based on decreases in α5β1 HDX. Interestingly, the dominant effect of Fn binding was an increase in α5β1 deuterium exchange distributed over multiple sites that undergo changes in conformation or solvent accessibility and appear to be sites where energy is stored in the higher-energy, open-integrin conformation. Together, our results highlight regions important for α5β1 binding to Fn and dynamics associated with this interaction.

Insight into the Conformational Transitions of Serine Acetyl Transferase Isoforms in E. histolytica: Implications for Structural and Functional Balance

Serine acetyl transferase (SAT) is one of the crucial enzymes in the cysteine biosynthetic pathway and an essential enzyme for the survival of Entamoeba histolytica, the causative agent of amoebiasis. E. histolytica expresses three isoforms of SAT, where SAT1 and SAT2 are inhibited by the final product cysteine, while SAT3 is not inhibited. SAT3 has a slightly elongated C-terminus compared to SAT1. To understand the stability and conformational transition between two secondary structures of proteins, we measured the effect of urea, a chemical denaturant, on two isoforms of SAT (SAT1 and SAT3) of E. histolytica. The effect of urea on the structure and stability of SAT1 and SAT3 was determined by measuring changes in their far-UV circular dichroism (CD), Trp fluorescence, and near-UV absorption spectra. The urea-induced normal transition curves suggested that the structural transition is reversible and follows a two-state process. Analysis of the urea-induced transition of all optical properties for the stability parameters ΔG _D° (Gibbs free energy change (ΔG _D) in the absence of urea), m (dependence of ΔG _D on urea concentration), and C _m (midpoint of urea transition) suggested that SAT1 is more stable than SAT3. Characterization of the end product of the urea-induced transition of both proteins by the far-UV CD and Trp-fluorescence and near-UV absorbance suggested that urea causes α-helix to β-sheet transition and burial of Trp residues, respectively. To support the in vitro findings, 100 ns molecular dynamics simulations (in silico study) were performed. Both the spectroscopic and molecular dynamics approaches clearly indicated that SAT1 is more stable than SAT3. SAT3 has evolved to escape the feedback inhibition to keep producing cysteine, but in the process, it compromises its structural stability relative to SAT1.

Biophysical Fitness Landscape of the SARS-CoV-2 Delta Variant Receptor Binding Domain

Among the five known SARS-CoV-2 variants of concern, Delta is the most virulent leading to severe symptoms and increased mortality among infected people. Our study seeks to examine how the biophysical parameters of the Delta variant correlate to the clinical observations. Receptor binding domain (RBD) is the first point of contact with the human host cells and is the immunodominant form of the spike protein. Delta variant RBD contains two novel mutations L452R and T478K. We examined the effect of single as well as the double mutations on RBD expression in human Expi293 cells, RBD stability using urea and thermal denaturation, and RBD binding to angiotensin converting enzyme 2 (ACE2) receptor and to neutralizing antibodies using isothermal titration calorimetry. Delta variant RBD showed significantly higher expression compared to the wild-type RBD, and the increased expression is due to L452R mutation. Despite their non-conservative nature, none of the mutations significantly affected RBD structure and stability. All mutants showed similar binding affinity to ACE2 and to Class 1 antibodies (CC12.1 and LY-CoV016) as that of the wild-type. Delta double mutant L452R/T478K showed no binding to Class 2 antibodies (P2B-2F6 and LY-CoV555) and a hundred-fold weaker binding to a Class 3 antibody (REGN10987), and the decreased antibody binding is determined by the L452R mutation. These results indicate that the immune escape from neutralizing antibodies, rather than increased receptor binding, is the main biophysical parameter that determined the fitness landscape of the Delta variant RBD.

Effect of pH on Diclofenac-Lysozyme Interaction: Structural and Functional Aspect

As a nonsteroidal antiinflammatory drug, diclofenac (DCF) is used in the treatment of a variety of human ailments. It has already been reported that the use of this class of drugs for a longer duration is associated with numerous side effects such as cardiovascular implications, reno-medullary complications, etc. In the present study, the effect of DCF on the structure, stability, and function of lysozyme was studied. The study was designed to examine the effect of DCF only at various pH values. Heat-induced denaturation of lysozyme was analyzed in the presence and absence of various molar concentrations of DCF at different pH values. The values of thermodynamic parameters, the midpoint of denaturation (T m), enthalpy change at T m (ΔH m), constant pressure heat capacity change (ΔC p), and Gibbs energy change at 25°C (ΔG D o), thus obtained under a given set of conditions (pH and molar concentration of DCF), demonstrated the following 1) DCF destabilized lysozyme with respect of T m and ΔG D o at all the pH values, 2) the magnitude of protein destabilization is lesser at acidic pH than at physiological pH, 3) structural changes in lysozyme are less projecting at pH 2.0 than at pH 7.0, and 4) quenching is observed at both pH values. Furthermore, the process of protein destabilization in the presence of DCF is entropically driven.

Conversion of the Native N-Terminal Domain of TDP-43 into a Monomeric Alternative Fold with Lower Aggregation Propensity

TAR DNA-binding protein 43 (TDP-43) forms intraneuronal cytoplasmic inclusions associated with amyotrophic lateral sclerosis and ubiquitin-positive frontotemporal lobar degeneration. Its N-terminal domain (NTD) can dimerise/oligomerise with the head-to-tail arrangement, which is essential for function but also favours liquid-liquid phase separation and inclusion formation of full-length TDP-43. Using various biophysical approaches, we identified an alternative conformational state of NTD in the presence of Sulfobetaine 3-10 (SB3-10), with higher content of α-helical structure and tryptophan solvent exposure. NMR shows a highly mobile structure, with partially folded regions and β-sheet content decrease, with a concomitant increase of α-helical structure. It is monomeric and reverts to native oligomeric NTD upon SB3-10 dilution. The equilibrium GdnHCl-induced denaturation shows a cooperative folding and a somewhat lower conformational stability. When the aggregation processes were compared with and without pre-incubation with SB3-10, but at the identical final SB3-10 concentration, a slower aggregation was found in the former case, despite the reversible attainment of the native conformation in both cases. This was attributed to protein monomerization and oligomeric seeds disruption by the conditions promoting the alternative conformation. Overall, the results show a high plasticity of TDP-43 NTD and identify strategies to monomerise TDP-43 NTD for methodological and biomedical applications.