Thermal unfolding and conformational stability of the recombinant domain II of glutamate dehydrogenase from the hyperthermophile Thermotoga maritima

Protein mimics of fusion core from SARS-CoV-1 can inhibit SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a member of the genus Betacoronavirus (subgenus Sarbecovirus) and shares significant genomic and phylogenetic similarities with severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1). SARS-CoV-2 infection occurs through membrane fusion between the virus and host cell membranes, which is facilitated by the spike glycoprotein subunit 2 (S2). The folding of three heptad-repeat regions 1 (HR1) into a central trimeric core structure, along with the binding of three heptad-repeat regions 2 (HR2) in an antiparallel manner within the groove formed between the HR1 regions, which provides the driving force for membrane fusion. In this study, trimeric and monomeric six-helix bundles (6HB) were created by combining various truncations of the sequences from SARS-CoV-2 HR1 and HR2. In addition, monomeric five-helix bundles (5HB) were constructed using a similar method. Finally, we demonstrated a protein mimic, 5HB_V1 (from SARS-CoV-1), that exhibits activity in inhibiting SARS-CoV-2. These findings suggest a strategy to design monomeric 6HB and 5HB based on the SARS-CoV-2 fusion core: maintain the flanking sequences outside the α-helix region in HR2 and introduce point mutations to enhance hydrogen bonding between the helix bundles. The 5HB could be a target for designing new inhibitors against SARS-CoV-1 and SARS-CoV-2.

Differential scanning calorimetry of proteins and Zimm-Bragg model in water

Differential Scanning Calorimetry (DSC) is a regular and powerful tool to measure the specific heat profile of various materials. Hydrogen bonds play a crucial role in stabilizing the three-dimensional structure of proteins. Naturally, information about the strength of hydrogen bonds is contained in the measured DSC profiles. Despite its obvious importance, there is no approach that would allow the extraction of such information from the heat capacity measurements. In order to connect the measured profile to microscopic properties of a polypeptide chain, a proper model is required to fit. Using recent advances in the Zimm-Bragg (ZB) theory of protein folding in water, we propose a new and efficient algorithm to process the DSC experimental data and to extract the H-bonding energy among other relevant constants. Thus, for the randomly picked set of 33 proteins, we have found a quite narrow distribution of hydrogen bonding energies from 1 to 8 kJ/mol with the average energy of intra-protein hydrogen bonds h¯=4.2±1.5 kJ/mol and the average energy of water-protein bonds as hps¯=3.8±1.5 kJ/mol. This is an important illustration of a tiny disbalance between the water-protein and intraprotein hydrogen bonds. Fitted values of the nucleation parameter σ belong to the range from 0.001 to 0.01, as expected. The reported method can be considered as complementary to the classical two-state approach and together with other parameters provides the protein-water and intraprotein H-bonding energies, not accessible within the two-state paradigm.

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.

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.

Protein Engineering and HDX Identify Structural Regions of G-CSF Critical to Its Stability and Aggregation

The protein engineering and formulation of therapeutic proteins for prolonged shelf-life remain a major challenge in the biopharmaceutical industry. Understanding the influence of mutations and formulations on the protein structure and dynamics could lead to more predictive approaches to their improvement. Previous intrinsic fluorescence analysis of the chemically denatured granulocyte colony-stimulating factor (G-CSF) suggested that loop AB could subtly reorganize to form an aggregation-prone intermediate state. Hydrogen deuterium exchange mass spectrometry (HDX-MS) has also revealed that excipient binding increased the thermal unfolding transition midpoint (T_m) by stabilizing loop AB. Here, we have combined protein engineering with biophysical analyses and HDX-MS to reveal that increased exchange in a core region of the G-CSF comprising loop AB (ABI, a small helix, ABII) and loop CD packed onto helix B and the beginning of loop BC leads to a decrease in T_m and higher aggregation rates. Furthermore, some mutations can increase the population of the aggregation-prone conformation within the native ensemble, as measured by the greater local exchange within this core region.

Thermostable designed ankyrin repeat proteins (DARPins) as building blocks for innovative drugs

Designed ankyrin repeat proteins (DARPins) are antibody mimetics with high and mostly unexplored potential in drug development. By using in silico analysis and a rationally guided Ala scanning, we identified position 17 of the N-terminal capping repeat to play a key role in overall protein thermostability. The melting temperature of a DARPin domain with a single full-consensus internal repeat was increased by 8 °C to 10 °C when Asp17 was replaced by Leu, Val, Ile, Met, Ala, or Thr. We then transferred the Asp17Leu mutation to various backgrounds, including clinically validated DARPin domains, such as the vascular endothelial growth factor-binding domain of the DARPin abicipar pegol. In all cases, these proteins showed improvements in the thermostability on the order of 8 °C to 16 °C, suggesting the replacement of Asp17 could be generically applicable to this drug class. Molecular dynamics simulations showed that the Asp17Leu mutation reduces electrostatic repulsion and improves van-der-Waals packing, rendering the DARPin domain less flexible and more stable. Interestingly, this beneficial Asp17Leu mutation is present in the N-terminal caps of three of the five DARPin domains of ensovibep, a SARS-CoV-2 entry inhibitor currently in clinical development, indicating this mutation could be partly responsible for the very high melting temperature (>90 °C) of this promising anti-COVID-19 drug. Overall, such N-terminal capping repeats with increased thermostability seem to be beneficial for the development of innovative drugs based on DARPins.

A structure-function study of ZraP and ZraS provides new insights into the two-component system Zra

BACKGROUND

Zra belongs to the envelope stress response (ESR) two-component systems (TCS). It is atypical because of its third periplasmic repressor partner (ZraP), in addition to its histidine kinase sensor protein (ZraS) and its response regulator (ZraR) components. Furthermore, although it is activated by Zn²⁺, it is not involved in zinc homeostasis or protection against zinc toxicity. Here, we mainly focus on ZraS but also provide information on ZraP.

METHODS

The purified periplasmic domain of ZraS and ZraP were characterized using biophysical and biochemical technics: multi-angle laser light scattering (MALLS), circular dichroism (CD), differential scanning fluorescence (DSF), inductively coupled plasma atomic emission spectroscopy (ICP-AES), cross-linking and small-angle X-ray scattering (SAXS). In-vivo experiments were carried out to determine the redox state of the cysteine residue in ZraP and the consequences for the cell of an over-activation of the Zra system.

RESULTS

We show that ZraS binds one Zn²⁺ molecule with high affinity resulting in conformational changes of the periplasmic domain, consistent with a triggering function of the metal ion. We also demonstrate that, in the periplasm, the only cysteine residue of ZraP is at least partially reduced. Using SAXS, we conclude that the previously determined X-ray structure is different from the structure in solution.

CONCLUSION

Our results allow us to propose a general mechanism for the Zra system activation and to compare it to the homologous Cpx system.

GENERAL SIGNIFICANCE

We bring new input on the so far poorly described Zra system and notably on ZraS.

HDX and In Silico Docking Reveal that Excipients Stabilize G-CSF via a Combination of Preferential Exclusion and Specific Hotspot Interactions

Assuring the stability of therapeutic proteins is a major challenge in the biopharmaceutical industry, and a better molecular understanding of the mechanisms through which formulations influence their stability is an ongoing priority. While the preferential exclusion effects of excipients are well known, the additional presence and impact of specific protein-excipient interactions have proven to be more elusive to identify and characterize. We have taken a combined approach of in silico molecular docking and hydrogen deuterium exchange-mass spectrometry (HDX-MS) to characterize the interactions between granulocyte colony-stimulating factor (G-CSF), and some common excipients. These interactions were related to their influence on the thermal-melting temperatures (T_m) for the nonreversible unfolding of G-CSF in liquid formulations. The residue-level interaction sites predicted in silico correlated well with those identified experimentally and highlighted the potential impact of specific excipient interactions on the T_m of G-CSF.

nanoDSF: In vitro Label-Free Method to Monitor Picornavirus Uncoating and Test Compounds Affecting Particle Stability

Thermal shift assays measure the stability of macromolecules and macromolecular assemblies as a function of temperature. The Particle Stability Thermal Release Assay (PaSTRy) of picornaviruses is based on probes becoming strongly fluorescent upon binding to hydrophobic patches of the protein capsid (e.g., SYPRO Orange) or to the viral RNA genome (e.g., SYTO-82) that become exposed upon heating virus particles. PaSTRy has been exploited for studying the stability of viral mutants, viral uncoating, and the effect of capsid-stabilizing compounds. While the results were usually robust, the thermal shift assay with SYPRO Orange is sensitive to surfactants and EDTA and failed at least to correctly report the effect of excipients on an inactivated poliovirus 3 vaccine. Furthermore, interactions between the probe and capsid-binding antivirals as well as mutual competition for binding sites cannot be excluded. To overcome these caveats, we assessed differential scanning fluorimetry with a nanoDSF device as a label-free alternative. NanoDSF monitors the changes in the intrinsic tryptophan fluorescence (ITF) resulting from alterations of the 3D-structure of proteins as a function of the temperature. Using rhinovirus A2 as a model, we demonstrate that nanoDFS is well suited for recording the temperature-dependence of conformational changes associated with viral uncoating with minute amounts of sample. We compare it with orthogonal methods and correlate the increase in viral RNA exposure with PaSTRy measurements. Importantly, nanoDSF correctly identified the thermal stabilization of RV-A2 by pleconaril, a prototypic pocket-binding antiviral compound. NanoDFS is thus a label-free, high throughput-customizable, attractive alternative for the discovery of capsid-binding compounds impacting on viral stability.

Thermal unfolding and refolding of a lytic polysaccharide monooxygenase fromThermoascus aurantiacus

TaLPMO9A regains its catalytic power after a thermal unfolding and refolding cycle.

Characterisation of sensor kinase by CD spectroscopy: golden rules and tips

This is a review that describes the golden rules and tips on how to characterise the molecular interactions of membrane sensor kinase proteins with ligands using mainly circular dichroism (CD) spectroscopy. CD spectroscopy is essential for this task as any conformational change observed in the far-UV (secondary structures (α-helix, β-strands, poly-proline of type II, β-turns, irregular and folding) and near-UV regions [local environment of the aromatic side-chains of amino acid residues (Phe, Tyr and Trp) and ligands (drugs) and prosthetic groups (porphyrins, cofactors and coenzymes (FMN, FAD, NAD))] upon ligand addition to the protein can be used to determine qualitatively and quantitatively ligand-binding interactions. Advantages of using CD versus other techniques will be discussed. The difference CD spectra of the protein-ligand mixtures calculated subtracting the spectra of the ligand at various molar ratios can be used to determine the type of conformational changes induced by the ligand in terms of the estimated content of the various elements of protein secondary structure. The highly collimated microbeam and high photon flux of Diamond Light Source B23 beamline for synchrotron radiation circular dichroism (SRCD) enable the use of minimal amount of membrane proteins (7.5 µg for a 0.5 mg/ml solution) for high-throughput screening. Several examples of CD titrations of membrane proteins with a variety of ligands are described herein including the protocol tips that would guide the choice of the appropriate parameters to conduct these titrations by CD/SRCD in the best possible way.

Coupled molecular dynamics mediate long- and short-range epistasis between mutations that affect stability and aggregation kinetics

Multiple mutations are typically required to significantly improve protein stability or aggregation kinetics. However, when several substitutions are made in a single protein, the mutations can potentially interact in a nonadditive manner, resulting in epistatic effects, which can hamper protein-engineering strategies to improve thermostability or aggregation kinetics. Here, we have examined the role of protein dynamics in mediating epistasis between pairs of mutations. With Escherichia coli transketolase (TK) as a model, we explored the epistatic interactions between two single variants H192P and A282P, and also between the double-mutant H192P/A282P and two single variants, I365L or G506A. Epistasis was determined for several measures of protein stability, including the following: the free-energy barrier to kinetic inactivation, ∆∆G ^‡; thermal transition midpoint temperatures, T _m; and aggregation onset temperatures, T _agg Nonadditive epistasis was observed between neighboring mutations as expected, but also for distant mutations located in the surface and core regions of different domains. Surprisingly, the epistatic behaviors for each measure of stability were often different for any given pairwise recombination, highlighting that kinetic and thermodynamic stabilities do not always depend on the same structural features. Molecular-dynamics simulations and a pairwise cross-correlation analysis revealed that mutations influence the dynamics of their local environment, but also in some cases the dynamics of regions distant in the structure. This effect was found to mediate epistatic interactions between distant mutations and could therefore be exploited in future protein-engineering strategies.

Computational Design To Reduce Conformational Flexibility and Aggregation Rates of an Antibody Fab Fragment

Computationally guided semirational design has significant potential for improving the aggregation kinetics of protein biopharmaceuticals. While improvement in the global conformational stability can stabilize proteins to aggregation under some conditions, previous studies suggest that such an approach is limited, because thermal transition temperatures ( T_m) and the fraction of protein unfolded ( f_T) tend to only correlate with aggregation kinetics where the protein is incubated at temperatures approaching the T_m. This is because under these conditions, aggregation from globally unfolded protein becomes dominant. However, under native conditions, the aggregation kinetics are presumed to be dependent on local structural fluctuations or partial unfolding of the native state, which reveal regions of high propensity to form protein-protein interactions that lead to aggregation. In this work, we have targeted the design of stabilizing mutations to regions of the A33 Fab surface structure, which were predicted to be more flexible. This Fab already has high global stability, and global unfolding is not the main cause of aggregation under most conditions. Therefore, the aim was to reduce the conformational flexibility and entropy of the native protein at various locations and thus identify which of those regions has the greatest influence on the aggregation kinetics. Highly dynamic regions of structure were identified through both molecular dynamics simulation and B-factor analysis of related X-ray crystal structures. The most flexible residues were mutated into more stable variants, as predicted by Rosetta, which evaluates the ΔΔ G_ND for each potential point mutation. Additional destabilizing variants were prepared as controls to evaluate the prediction accuracy and also to assess the general influence of conformational stability on aggregation kinetics. The thermal conformational stability, and aggregation rates of 18 variants at 65 °C, were each examined at pH 4, 200 mM ionic strength, under which conditions the initial wild-type protein was <5% unfolded. Variants with decreased T_m values led to more rapid aggregation due to an increase in the fraction of protein unfolded under the conditions studied. As expected, no significant improvements were observed in the global conformational stability as measured by T_m. However, 6 of the 12 stable variants led to an increase in the cooperativity of unfolding, consistent with lower conformational flexibility and entropy in the native ensemble. Three of these had 5-11% lower aggregation rates, and their structural clustering indicated that the local dynamics of the C-terminus of the heavy chain had a role in influencing the aggregation rate.

Influence of solubilization and AD-mutations on stability and structure of human presenilins

Presenilin (PS1 or PS2) functions as the catalytic subunit of γ-secretase, which produces the toxic amyloid beta peptides in Alzheimer’s disease (AD). The dependence of folding and structural stability of PSs on the lipophilic environment and mutation were investigated by far UV CD spectroscopy. The secondary structure content and stability of PS2 depended on the lipophilic environment. PS2 undergoes a temperature-dependent structural transition from α-helical to β-structure at 331 K. The restructured protein formed structures which tested positive in spectroscopic amyloid fibrils assays. The AD mutant PS1L266F, PS1L424V and PS1ΔE9 displayed reduced stability which supports a proposed ‘loss of function’ mechanism of AD based on protein instability. The exon 9 coded sequence in the inhibitory loop of the zymogen was found to be required for the modulation of the thermal stability of PS1 by the lipophilic environment.

Cold adaptation of tRNA nucleotidyltransferases: A tradeoff in activity, stability and fidelity

Cold adaptation is an evolutionary process that has dramatic impact on enzymatic activity. Increased flexibility of the protein structure represents the main evolutionary strategy for efficient catalysis and reaction rates in the cold, but is achieved at the expense of structural stability. This results in a significant activity-stability tradeoff, as it was observed for several metabolic enzymes. In polymerases, however, not only reaction rates, but also fidelity plays an important role, as these enzymes have to synthesize copies of DNA and RNA as exact as possible. Here, we investigate the effects of cold adaptation on the highly accurate CCA-adding enzyme, an RNA polymerase that uses an internal amino acid motif within the flexible catalytic core as a template to synthesize the CCA triplet at tRNA 3'-ends. As the relative orientation of these residues determines nucleotide selection, we characterized how cold adaptation impacts template reading and fidelity. In a comparative analysis of closely related psychro-, meso-, and thermophilic enzymes, the cold-adapted polymerase shows a remarkable error rate during CCA synthesis in vitro as well as in vivo. Accordingly, CCA-adding activity at low temperatures is not only achieved at the expense of structural stability, but also results in a reduced polymerization fidelity.

DNA and IκBα Both Induce Long-Range Conformational Changes in NFκB

We recently discovered that IκBα enhances the rate of release of nuclear factor kappa B (NFκB) from DNA target sites in a process we have termed molecular stripping. Coarse-grained molecular dynamics simulations of the stripping pathway revealed two mechanisms for the enhanced release rate: the negatively charged PEST region of IκBα electrostatically repels the DNA, and the binding of IκBα appears to twist the NFκB heterodimer so that the DNA can no longer bind. Here, we report amide hydrogen/deuterium exchange data that reveal long-range allosteric changes in the NFκB (RelA-p50) heterodimer induced by DNA or IκBα binding. The data suggest that the two Ig-like subdomains of each Rel-homology region, which are connected by a flexible linker in the heterodimer, communicate in such a way that when DNA binds to the N-terminal DNA-binding domains, the nuclear localization signal becomes more highly exchanging. Conversely, when IκBα binds to the dimerization domains, amide exchange throughout the DNA-binding domains is decreased as if the entire domain is becoming globally stabilized. The results help understand how the subtle mechanism of molecular stripping actually occurs.

Linkage between Fitness of Yeast Cells and Adenylate Kinase Catalysis

Enzymes have evolved with highly specific values of their catalytic parameters k_cat and K_M. This poses fundamental biological questions about the selection pressures responsible for evolutionary tuning of these parameters. Here we are address these questions for the enzyme adenylate kinase (Adk) in eukaryotic yeast cells. A plasmid shuffling system was developed to allow quantification of relative fitness (calculated from growth rates) of yeast in response to perturbations of Adk activity introduced through mutations. Biophysical characterization verified that all variants studied were properly folded and that the mutations did not cause any substantial differences to thermal stability. We found that cytosolic Adk is essential for yeast viability in our strain background and that viability could not be restored with a catalytically dead, although properly folded Adk variant. There exist a massive overcapacity of Adk catalytic activity and only 12% of the wild type k_cat is required for optimal growth at the stress condition 20°C. In summary, the approach developed here has provided new insights into the evolutionary tuning of k_cat for Adk in a eukaryotic organism. The developed methodology may also become useful for uncovering new aspects of active site dynamics and also in enzyme design since a large library of enzyme variants can be screened rapidly by identifying viable colonies.

Macromolecular crowding effects on two homologs of ribosomal protein s16: protein-dependent structural changes and local interactions

Proteins function in cellular environments that are crowded with biomolecules, and in this reduced available space, their biophysical properties may differ from those observed in dilute solutions in vitro. Here, we investigated the effects of a synthetic macromolecular crowding agent, dextran 20, on the folded states of hyperthermophilic (S16Thermo) and mesophilic (S16Meso) homologs of the ribosomal protein S16. As expected for an excluded-volume effect, the resistance of the mesophilic protein to heat-induced unfolding increased in the presence of dextran 20, and chemical denaturation experiments at different fixed temperatures showed the macromolecular crowding effect to be temperature-independent. Förster resonance energy transfer experiments show that intramolecular distances between an intrinsic Trp residue and BODIPY-labeled S16Meso depend on the level of the crowding agent. The BODIPY group was attached at three specific positions in S16Meso, allowing measurements of three intraprotein distances. All S16Meso variants exhibited a decrease in the average Trp-BODIPY distance at up to 100 mg/mL dextran 20, whereas the changes in distance became anisotropic (one distance increased, two distances decreased) at higher dextran concentrations. In contrast, the two S16Thermo mutants did not show any changes in Trp-BODIPY distances upon increase of dextran 20 concentrations. It should be noted that the fluorescence quantum yields and lifetimes of BODIPY attached to the two S16 homologs decreased gradually in the presence of dextran 20. To investigate the origin of this decrease, we studied the BODIPY quantum yield in three protein variants in the presence of a tyrosine-labeled dextran. The experiments revealed distinct tyrosine quenching behaviors of BODIPY in the three variants, suggesting a dynamic local interaction between dextran and one particular S16 variant.

Structural basis of Lewis(b) antigen binding by the Helicobacter pylori adhesin BabA

Helicobacter pylori is a leading cause of peptic ulceration and gastric cancer worldwide. To achieve colonization of the stomach, this Gram-negative bacterium adheres to Lewis(b) (Le(b)) antigens in the gastric mucosa using its outer membrane protein BabA. Structural information for BabA has been elusive, and thus, its molecular mechanism for recognizing Le(b) antigens remains unknown. We present the crystal structure of the extracellular domain of BabA, from H. pylori strain J99, in the absence and presence of Le(b) at 2.0- and 2.1-Å resolutions, respectively. BabA is a predominantly α-helical molecule with a markedly kinked tertiary structure containing a single, shallow Le(b) binding site at its tip within a β-strand motif. No conformational change occurs in BabA upon binding of Le(b), which is characterized by low affinity under acidic [K D (dissociation constant) of ~227 μM] and neutral (K D of ~252 μM) conditions. Binding is mediated by a network of hydrogen bonds between Le(b) Fuc1, GlcNAc3, Fuc4, and Gal5 residues and a total of eight BabA amino acids (C189, G191, N194, N206, D233, S234, S244, and T246) through both carbonyl backbone and side-chain interactions. The structural model was validated through the generation of two BabA variants containing N206A and combined D233A/S244A substitutions, which result in a reduction and complete loss of binding affinity to Le(b), respectively. Knowledge of the molecular basis of Le(b) recognition by BabA provides a platform for the development of therapeutics targeted at inhibiting H. pylori adherence to the gastric mucosa.

Completing the structural family portrait of the human EphB tyrosine kinase domains

The EphB receptors have key roles in cell morphology, adhesion, migration and invasion, and their aberrant action has been linked with the development and progression of many different tumor types. Their conflicting expression patterns in cancer tissues, combined with their high sequence and structural identity, present interesting challenges to those seeking to develop selective therapeutic molecules targeting this large receptor family. Here, we present the first structure of the EphB1 tyrosine kinase domain determined by X-ray crystallography to 2.5Å. Our comparative crystalisation analysis of the human EphB family kinases has also yielded new crystal forms of the human EphB2 and EphB4 catalytic domains. Unable to crystallize the wild-type EphB3 kinase domain, we used rational engineering (based on our new structures of EphB1, EphB2, and EphB4) to identify a single point mutation which facilitated its crystallization and structure determination to 2.2 Å. This mutation also improved the soluble recombinant yield of this kinase within Escherichia coli, and increased both its intrinsic stability and catalytic turnover, without affecting its ligand-binding profile. The partial ordering of the activation loop in the EphB3 structure alludes to a potential cis-phosphorylation mechanism for the EphB kinases. With the kinase domain structures of all four catalytically competent human EphB receptors now determined, a picture begins to emerge of possible opportunities to produce EphB isozyme-selective kinase inhibitors for mechanistic studies and therapeutic applications.

The proline-rich motif of the proDer p 3 allergen propeptide is crucial for protease-protease interaction

The majority of proteases are synthesized in an inactive form, termed zymogen, which consists of a propeptide and a protease domain. The propeptide is commonly involved in the correct folding and specific inhibition of the enzyme. The propeptide of the house dust mite allergen Der p 3, NPILPASPNAT, contains a proline-rich motif (PRM), which is unusual for a trypsin-like protease. By truncating the propeptide or replacing one or all of the prolines in the non-glycosylated zymogen with alanine(s), we demonstrated that the full-length propeptide is not required for correct folding and thermal stability and that the PRM is important for the resistance of proDer p 3 to undesired proteolysis when the protein is expressed in Pichia pastoris. Additionally, we followed the maturation time course of proDer p 3 by coupling a quenched-flow assay to mass spectrometry analysis. This approach allowed to monitor the evolution of the different species and to determine the steady-state kinetic parameters for activation of the zymogen by the major allergen Der p 1. This experiment demonstrated that prolines 5 and 8 are crucial for proDer p 3-Der p 1 interaction and for activation of the zymogen.

Refolding of a thermostable glyceraldehyde dehydrogenase for application in synthetic cascade biomanufacturing

The production of chemicals from renewable resources is gaining importance in the light of limited fossil resources. One promising alternative to widespread fermentation based methods used here is Synthetic Cascade Biomanufacturing, the application of minimized biocatalytic reaction cascades in cell free processes. One recent example is the development of the phosphorylation independent conversion of glucose to ethanol and isobutanol using only 6 and 8 enzymes, respectively. A key enzyme for this pathway is aldehyde dehydrogenase from Thermoplasma acidophilum, which catalyzes the highly substrate specific oxidation of d-glyceraldehyde to d-glycerate. In this work the enzyme was recombinantly expressed in Escherichia coli. Using matrix-assisted refolding of inclusion bodies the yield of enzyme production was enhanced 43-fold and thus for the first time the enzyme was provided in substantial amounts. Characterization of structural stability verified correct refolding of the protein. The stability of the enzyme was determined by guanidinium chloride as well as isobutanol induced denaturation to be ca. -8 kJ/mol both at 25°C and 40°C. The aldehyde dehydrogenase is active at high temperatures and in the presence of small amounts of organic solvents. In contrast to previous publications, the enzyme was found to accept NAD(+) as cofactor making it suitable for application in the artificial glycolysis.

Biochemical and biophysical characterization of four EphB kinase domains reveals contrasting thermodynamic, kinetic and inhibition profiles

The Eph (erythropoietin-producing hepatocellular carcinoma) B receptors are important in a variety of cellular processes through their roles in cell-to-cell contact and signalling; their up-regulation and down-regulation has been shown to have implications in a variety of cancers. A greater understanding of the similarities and differences within this small, highly conserved family of tyrosine kinases will be essential to the identification of effective therapeutic opportunities for disease intervention. In this study, we have developed a route to production of multi-milligram quantities of highly purified, homogeneous, recombinant protein for the kinase domain of these human receptors in Escherichia coli. Analyses of these isolated catalytic fragments have revealed stark contrasts in their amenability to recombinant expression and their physical properties: e.g., a >16°C variance in thermal stability, a 3-fold difference in catalytic activity and disparities in their inhibitor binding profiles. We find EphB3 to be an outlier in terms of both its intrinsic stability, and more importantly its ligand-binding properties. Our findings have led us to speculate about both their biological significance and potential routes for generating EphB isozyme-selective small-molecule inhibitors. Our comprehensive methodologies provide a template for similar in-depth studies of other kinase superfamily members.

Autoproteolysis and intramolecular dissociation of Yersinia YscU precedes secretion of its C-terminal polypeptide YscU(CC)

Type III secretion system mediated secretion and translocation of Yop-effector proteins across the eukaryotic target cell membrane by pathogenic Yersinia is highly organized and is dependent on a switching event from secretion of early structural substrates to late effector substrates (Yops). Substrate switching can be mimicked in vitro by modulating the calcium levels in the growth medium. YscU that is essential for regulation of this switch undergoes autoproteolysis at a conserved N↑PTH motif, resulting in a 10 kDa C-terminal polypeptide fragment denoted YscU(CC). Here we show that depletion of calcium induces intramolecular dissociation of YscU(CC) from YscU followed by secretion of the YscU(CC) polypeptide. Thus, YscU(CC) behaved in vivo as a Yop protein with respect to secretion properties. Further, destabilized yscU mutants displayed increased rates of dissociation of YscU(CC)in vitro resulting in enhanced Yop secretion in vivo at 30°C relative to the wild-type strain.These findings provide strong support to the relevance of YscU(CC) dissociation for Yop secretion. We propose that YscU(CC) orchestrates a block in the secretion channel that is eliminated by calcium depletion. Further, the striking homology between different members of the YscU/FlhB family suggests that this protein family possess regulatory functions also in other bacteria using comparable mechanisms.

Modulation of a pre-existing conformational equilibrium tunes adenylate kinase activity

Structural plasticity is often required for distinct microscopic steps during enzymatic reaction cycles. Adenylate kinase from Escherichia coli (AK(eco)) populates two major conformations in solution; the open (inactive) and closed (active) state, and the overall turnover rate is inversely proportional to the lifetime of the active conformation. Therefore, structural plasticity is intimately coupled to enzymatic turnover in AK(eco). Here, we probe the open to closed conformational equilibrium in the absence of bound substrate with NMR spectroscopy and molecular dynamics simulations. The conformational equilibrium in absence of substrate and, in turn, the turnover number can be modulated with mutational- and osmolyte-driven perturbations. Removal of one hydrogen bond between the ATP and AMP binding subdomains results in a population shift toward the open conformation and a resulting increase of k(cat). Addition of the osmolyte TMAO to AK(eco) results in population shift toward the closed conformation and a significant reduction of k(cat). The Michaelis constants (K(M)) scale with the change in k(cat), which follows from the influence of the population of the closed conformation for substrate binding affinity. Hence, k(cat) and K(M) are mutually dependent, and in the case of AK(eco), any perturbation that modulates k(cat) is mirrored with a proportional response in K(M). Thus, our results demonstrate that the equilibrium constant of a pre-existing conformational equilibrium directly affects enzymatic catalysis. From an evolutionary perspective, our findings suggest that, for AK(eco), there exists ample flexibility to obtain a specificity constant (k(cat)/K(M)) that commensurate with the exerted cellular selective pressure.

Extraordinary μs-ms backbone dynamics in Arabidopsis thaliana peroxiredoxin Q

Peroxiredoxin Q (PrxQ) isolated from Arabidopsis thaliana belongs to a family of redox enzymes called peroxiredoxins, which are thioredoxin- or glutaredoxin-dependent peroxidases acting to reduce peroxides and in particular hydrogen peroxide. PrxQ cycles between an active reduced state and an inactive oxidized state during its catalytic cycle. The catalytic mechanism involves a nucleophilic attack of the catalytic cysteine on hydrogen peroxide to generate a sulfonic acid intermediate with a concerted release of a water molecule. This intermediate is subsequently relaxed by the reaction of a second cysteine, denoted the resolving cysteine, generating an intramolecular disulfide bond and release of a second water molecule. PrxQ is recycled to the active state by a thioredoxin-dependent reduction. Previous structural studies of PrxQ homologues have provided the structural basis for the switch between reduced and oxidized conformations. Here, we have performed a detailed study of the activity, structure and dynamics of PrxQ in both the oxidized and reduced states. Reliable and experimentally validated structural models of PrxQ in both oxidation states were generated using homology based modeling. Analysis of NMR spin relaxation rates shows that PrxQ is monomeric in both oxidized and reduced states. As evident from R(2) relaxation rates the reduced form of PrxQ undergoes unprecedented dynamics on the slow μs-ms timescale. The ground state of this conformational dynamics is likely the stably folded reduced state as implied by circular dichroism spectroscopy. We speculate that the extensive dynamics is intimately related to the catalytic function of PrxQ.

A second look at mini-protein stability: analysis of FSD-1 using circular dichroism, differential scanning calorimetry, and simulations

Mini-proteins that contain <50 amino acids often serve as model systems for studying protein folding because their small size makes long timescale simulations possible. However, not all mini-proteins are created equal. The stability and structure of FSD-1, a 28-residue mini-protein that adopted the betabetaalpha zinc-finger motif independent of zinc binding, was investigated using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The broad melting transition of FSD-1, similar to that of a helix-to-coil transition, was observed by using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The N-terminal beta-hairpin was found to be flexible. The FSD-1 apparent melting temperature of 41 degrees C may be a reflection of the melting of its alpha-helical segment instead of the entire protein. Thus, despite its attractiveness due to small size and purposefully designed helix, sheet, and turn structures, the status of FSD-1 as a model system for studying protein folding should be reconsidered.

The zinc center influences the redox and thermodynamic properties of Escherichia coli thioredoxin 2

Thioredoxins are small, ubiquitous redox enzymes that reduce protein disulfide bonds by using a pair of cysteine residues present in a strictly conserved WCGPC catalytic motif. The Escherichia coli cytoplasm contains two thioredoxins, Trx1 and Trx2. Trx2 is special because it is induced under oxidative stress conditions and it has an additional N-terminal zinc-binding domain. We have determined the redox potential of Trx2, the pK(a) of the active site nucleophilic cysteine, as well as the stability of the oxidized and reduced form of the protein. Trx2 is more oxidizing than Trx1 (-221 mV versus -284 mV, respectively), which is in good agreement with the decreased value of the pK(a) of the nucleophilic cysteine (5.1 versus 7.1, respectively). The difference in stability between the oxidized and reduced forms of an oxidoreductase is the driving force to reduce substrate proteins. This difference is smaller for Trx2 (DeltaDeltaG degrees(H2O)=9 kJ/mol and DeltaT(m)=7. 4 degrees C) than for Trx1 (DeltaDeltaG degrees(H2O)=15 kJ/mol and DeltaT(m)=13 degrees C). Altogether, our data indicate that Trx2 is a significantly less reducing enzyme than Trx1, which suggests that Trx2 has a distinctive function. We disrupted the zinc center by mutating the four Zn(2+)-binding cysteines to serine. This mutant has a more reducing redox potential (-254 mV) and the pK(a) of its nucleophilic cysteine shifts from 5.1 to 7.1. The removal of Zn(2+) also decreases the overall stability of the reduced and oxidized forms by 3.2 kJ/mol and 5.8 kJ/mol, respectively. In conclusion, our data show that the Zn(2+)-center of Trx2 fine-tunes the properties of this unique thioredoxin.

Extreme temperature tolerance of a hyperthermophilic protein coupled to residual structure in the unfolded state

Understanding the mechanisms that dictate protein stability is of large relevance, for instance, to enable design of temperature-tolerant enzymes with high enzymatic activity over a broad temperature interval. In an effort to identify such mechanisms, we have performed a detailed comparative study of the folding thermodynamics and kinetics of the ribosomal protein S16 isolated from a mesophilic (S16(meso)) and hyperthermophilic (S16(thermo)) bacterium by using a variety of biophysical methods. As basis for the study, the 2.0 A X-ray structure of S16(thermo) was solved using single wavelength anomalous dispersion phasing. Thermal unfolding experiments yielded midpoints of 59 and 111 degrees C with associated changes in heat capacity upon unfolding (DeltaC(p)(0)) of 6.4 and 3.3 kJ mol(-1) K(-1), respectively. A strong linear correlation between DeltaC(p)(0) and melting temperature (T(m)) was observed for the wild-type proteins and mutated variants, suggesting that these variables are intimately connected. Stopped-flow fluorescence spectroscopy shows that S16(meso) folds through an apparent two-state model, whereas S16(thermo) folds through a more complex mechanism with a marked curvature in the refolding limb indicating the presence of a folding intermediate. Time-resolved energy transfer between Trp and N-(4,4-difluoro-5,7-dimethyl-4-bora-3a,4a-diaza-s-indacene-3-yl)methyl iodoacetamide of proteins mutated at selected positions shows that the denatured state ensemble of S16(thermo) is more compact relative to S16(meso). Taken together, our results suggest the presence of residual structure in the denatured state ensemble of S16(thermo) that appears to account for the large difference in quantified DeltaC(p)(0) values and, in turn, parts of the observed extreme thermal stability of S16(thermo). These observations may be of general importance in the design of robust enzymes that are highly active over a wide temperature span.

Lessons in stability from thermophilic proteins

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.

Conformational stability and thermodynamic characterization of the lipoic acid bearing domain of human mitochondrial branched chain alpha-ketoacid dehydrogenase

The lipoic acid bearing domain (hbLBD) of human mitochondrial branched chain alpha-ketoacid dehydrogenase (BCKD) plays important role of substrate channeling in oxidative decarboxylation of the branched chain alpha-ketoacids. Recently hbLBD has been found to follow two-step folding mechanism without detectable presence of stable or kinetic intermediates. The present study describes the conformational stability underlying the folding of this small beta-barrel domain. Thermal denaturation in presence of urea and isothermal urea denaturation titrations are used to evaluate various thermodynamic parameters defining the equilibrium unfolding. The linear extrapolation model successfully describes the two-step; native state <-->denatured state unfolding transition of hbLBD. The average temperature of maximum stability of hbLBD is estimated as 295.6 +/- 0.9 K. Cold denaturation of hbLBD is also predicted and discussed.

Thermal and conformational stability of Ssh10b protein from archaeon Sulfolobus shibattae

The secondary structure of the DNA binding protein Ssh10b is largely unaffected by change in temperature between 25 degrees C and 85 degrees C, indicating that the protein is highly thermostable. Here, we report the temperature-dependent equilibrium denaturation of Ssh10b in the presence of guanidine hydrochloride (GdnHCl). It was found that the transition midpoint values of the temperature (T(m)), and changes of enthalpy (DeltaH(m)) and entropy (DeltaS(m)) of Ssh10b unfolding were linearly decreasing with increasing GdnHCl concentration. The true values of the thermodynamic parameters, T(m)=402 K, DeltaH(m)=590+/-40 kJ x mol(-1) and DeltaS(m)=1.4+/-0.15 kJ x T(-1) x mol(-1), were obtained by linear extrapolation to 0 M GdnHCl. The value of the heat capacity change of Ssh10b unfolding, DeltaC(p)=3.8+/-0.2 kJ x T(-1) x mol(-1) (approx. 19 J T(-1) x mol residue(-1)), was obtained from the measured thermodynamic parameters. This is significantly smaller than that of the average value for mesophilic proteins (50 J.K(-1) x mol residue(-1)) or the value calculated from the Ssh10b structural data (64 J T(-1) x mol residue(-1)). A consequence of the small DeltaC(p) is that the DeltaG of Ssh10b is larger than that of mesophilic proteins, while the values of DeltaH and T*DeltaS are smaller. The small DeltaC(p) of Ssh10b appears to result mainly from the presence of compactness in the denatured state.

Experiment-guided thermodynamic simulations on reversible two-state proteins: implications for protein thermostability

Here, we perform protein thermodynamic simulations within a set of boundary conditions, effectively blanketing the experimental data. The thermodynamic parameters, melting temperature (TG), enthalpy change at the melting temperature (DeltaHG) and heat capacity change (DeltaCp) were systematically varied over the experimentally observed ranges for small single domain reversible two-state proteins. Parameter sets that satisfy the Gibbs-Helmholtz equation and yield a temperature of maximal stability (TS) around room temperature were selected. The results were divided into three categories by arbitrarily chosen TG ranges. The TG ranges in these categories correspond to typical values of the melting temperatures observed for the majority of the proteins from mesophilic, thermophilic and hyperthermophilic organisms. As expected, DeltaCp values tend to be high in mesophiles and low in hyperthermophiles. An increase in TG is accompanied by an up-shift and broadening of the protein stability curves, however, with a large scatter. Furthermore, the simulations reveal that the average DeltaHG increases with TG up to approximately 360 K and becomes constant thereafter. DeltaCp decreases with TG with different rates before and after approximately 360 K. This provides further justification for the separate grouping of proteins into thermophiles and hyperthermophiles to assess their thermodynamic differences. This analysis of the Gibbs-Helmholtz equation has allowed us to study the interdependence of the thermodynamic parameters TG, DeltaHG and DeltaCp and their derivatives in a more rigorous way than possible by the limited experimental protein thermodynamics data available in the literature. The results provide new insights into protein thermostability and suggest potential strategies for its manipulation.

A flexible approach for understanding protein stability

A distance constraint model (DCM) is presented that identifies flexible regions within protein structure consistent with specified thermodynamic condition. The DCM is based on a rigorous free energy decomposition scheme representing structure as fluctuating constraint topologies. Entropy non-additivity is problematic for naive decompositions, limiting the success of heat capacity predictions. The DCM resolves non-additivity by summing over independent entropic components determined by an efficient network-rigidity algorithm. A minimal 3-parameter DCM is demonstrated to accurately reproduce experimental heat capacity curves. Free energy landscapes and quantitative stability-flexibility relationships are obtained in terms of global flexibility. Several connections to experiment are made.

Why is trehalose an exceptional protein stabilizer? An analysis of the thermal stability of proteins in the presence of the compatible osmolyte trehalose

Trehalose, a naturally occurring osmolyte, is known to be an exceptional stabilizer of proteins and helps retain the activity of enzymes in solution as well as in the freeze-dried state. To understand the mechanism of action of trehalose in detail, we have conducted a thorough investigation of its effect on the thermal stability in aqueous solutions of five well characterized proteins differing in their various physico-chemical properties. Among them, RNase A has been used as a model enzyme to investigate the effect of trehalose on the retention of enzymatic activity upon incubation at high temperatures. 2 m trehalose was observed to raise the transition temperature, Tm of RNase A by as much as 18 degrees C and Gibbs free energy by 4.8 kcal mol-1 at pH 2.5. There is a decrease in the heat capacity of protein denaturation (DeltaCp) in trehalose solutions for all the studied proteins. An increase in the DeltaG and a decrease in the DeltaCp values for all the proteins points toward a general mechanism of stabilization due to the elevation and broadening of the stability curve (DeltaG versus T). A direct correlation of the surface tension of trehalose solutions and the thermal stability of various proteins has been observed. Wyman linkage analysis indicates that at 1.5 m concentration 4-7 molecules of trehalose are excluded from the vicinity of protein molecules upon denaturation. We further show that an increase in the stability of proteins in the presence of trehalose depends upon the length of the polypeptide chain. The pH dependence data suggest that even though the charge status of a protein contributes significantly, trehalose can be expected to work as a universal stabilizer of protein conformation due to its exceptional effect on the structure and properties of solvent water compared with other sugars and polyols.

The unusually slow relaxation kinetics of the folding-unfolding of pyrrolidone carboxyl peptidase from a hyperthermophile, Pyrococcus furiosus

In order to understand the thermodynamic and kinetic basis of the intrinsic stability of proteins from hyperthermophiles, the folding-unfolding reactions of cysteine-free pyrrolidone carboxyl peptidase (Cys142/188Ser) (PCP-0SH) from Pyrococcus furiosus were examined using circular dichroism (CD) and differential scanning calorimetry (DSC) at pH 2.3, where PCP-0SH exists in monomeric form. DSC showed a strong dependence of the shape and position of the unfolding profiles on the scan rate, suggesting the stability of PCP-0SH under kinetic control. On DSC timescales, even at a scan rate of 1 deg. C/hour, heat denaturation of PCP-0SH was non-equilibrium. However, over a long period of incubation of the heat-denatured PCP-0SH at pre-transition temperatures, it refolded completely, indicating reversibility with very slow relaxation kinetics. The rates of refolding of the heat-denatured PCP-0SH determined from the time-resolved DSC and CD spectroscopic progress curves were found to be similar within experimental error, confirming the mechanism of refolding to be a two-state process. The equilibrium established with a relaxation time of 5080 seconds (at t(m)=46.5 degrees C), which is unusually higher than the relaxation times observed for mesophilic and hyperthermophilic proteins. The long relaxation time may lead to the apparent irreversibility of an unfolding process occurring on the DSC experiment timescale. The refolding rate (9.8 x 10(-5) s(-1)) peaked near the t(m) (=46.5 degrees C), whereas the stability profile reached maxima (11.8 kJ mol(-1)) at 17 degrees C. The results clearly indicate the unusual mode of protein destabilization via a drastic decrease in the rate of folding at low pH and still maintaining a high activation energy barrier (284 kJ mol(-1)) for unfolding, which provides an effective kinetic advantage to unusually stable proteins from hyperthermophiles.