Thermodynamic characterization of monomeric and dimeric forms of CcdB (controller of cell division or death B protein)

Improved Prediction of Stabilizing Mutations in Proteins by Incorporation of Mutational Effects on Ligand Binding

While many computational methods accurately predict destabilizing mutations, identifying stabilizing mutations has remained a challenge, because of their relative rarity. We tested ΔΔG 0 predictions from computational predictors such as Rosetta, ThermoMPNN, RaSP, and DeepDDG, using 82 mutants of the bacterial toxin CcdB as a test case. On this dataset, the best computational predictor is ThermoMPNN, which identifies stabilizing mutations with a precision of 68%. However, the average increase in T m for these predicted mutations was only 1°C for CcdB, and predictions were poorer for a more challenging target, influenza neuraminidase. Using data from multiple previously described yeast surface display libraries and in vitro thermal stability measurements, we trained logistic regression models to identify stabilizing mutations with a precision of 90% and an average increase in T m of 3°C for CcdB. When such libraries contain a population of mutants with significantly enhanced binding relative to the corresponding wild type, there is no benefit in using computational predictors. It is then possible to predict stabilizing mutations without any training, simply by examining the distribution of mutational binding scores. This avoids laborious steps of in vitro expression, purification, and stability characterization. When this is not the case, combining data from computational predictors with high-throughput experimental binding data enhances the prediction of stabilizing mutations. However, this requires training on stability data measured in vitro with known stabilized mutants. It is thus feasible to predict stabilizing mutations rapidly and accurately for any system of interest that can be subjected to a binding selection or screen.

A CcdB toxin-derived peptide acts as a broad-spectrum antibacterial therapeutic in infected mice

The bacterial toxin CcdB (Controller of Cell death or division B) targets DNA Gyrase, an essential bacterial topoisomerase, which is also the molecular target for fluoroquinolones. Here, we present a short cell-penetrating 24-mer peptide, CP1-WT, derived from the Gyrase-binding region of CcdB and examine its effect on growth of Escherichia coli, Salmonella Typhimurium, Staphylococcus aureus and a carbapenem- and tigecycline-resistant strain of Acinetobacter baumannii in both axenic cultures and mouse models of infection. The CP1-WT peptide shows significant improvement over ciprofloxacin in terms of its in vivo therapeutic efficacy in treating established infections of S. Typhimurium, S. aureus and A. baumannii. The molecular mechanism likely involves inhibition of Gyrase or Topoisomerase IV, depending on the strain used. The study validates the CcdB binding site on bacterial DNA Gyrase as a viable and alternative target to the fluoroquinolone binding site.

Ter-Seq: A high-throughput method to stabilize transient ternary complexes and measure associated kinetics

Regulation of biological processes by proteins often involves the formation of transient, multimeric complexes whose characterization is mechanistically important but challenging. The bacterial toxin CcdB binds and poisons DNA Gyrase. The corresponding antitoxin CcdA extracts CcdB from its complex with Gyrase through the formation of a transient ternary complex, thus rejuvenating Gyrase. We describe a high throughput methodology called Ter-Seq to stabilize probable ternary complexes and measure associated kinetics using the CcdA-CcdB-GyrA14 ternary complex as a model system. The method involves screening a yeast surface display (YSD) saturation mutagenesis library of one partner (CcdB) for mutants that show enhanced ternary complex formation. We also isolated CcdB mutants that were either resistant or sensitive to rejuvenation, and used surface plasmon resonance (SPR) with purified proteins to validate the kinetics measured using the surface display. Positions, where CcdB mutations lead to slower rejuvenation rates, are largely involved in CcdA-binding, though there were several notable exceptions suggesting allostery. Mutations at these positions reduce the affinity towards CcdA, thereby slowing down the rejuvenation process. Mutations at GyrA14-interacting positions significantly enhanced rejuvenation rates, either due to reduced affinity or complete loss of CcdB binding to GyrA14. We examined the effect of different parameters (CcdA affinity, GyrA14 affinity, surface accessibilities, evolutionary conservation) on the rate of rejuvenation. Finally, we further validated the Ter-Seq results by monitoring the kinetics of ternary complex formation for individual CcdB mutants in solution by fluorescence resonance energy transfer (FRET) studies.

Mechanistic insights into global suppressors of protein folding defects

Most amino acid substitutions in a protein either lead to partial loss-of-function or are near neutral. Several studies have shown the existence of second-site mutations that can rescue defects caused by diverse loss-of-function mutations. Such global suppressor mutations are key drivers of protein evolution. However, the mechanisms responsible for such suppression remain poorly understood. To address this, we characterized multiple suppressor mutations both in isolation and in combination with inactive mutants. We examined six global suppressors of the bacterial toxin CcdB, the known M182T global suppressor of TEM-1 β-lactamase, the N239Y global suppressor of p53-DBD and three suppressors of the SARS-CoV-2 spike Receptor Binding Domain. When coupled to inactive mutants, they promote increased in-vivo solubilities as well as regain-of-function phenotypes. In the case of CcdB, where novel suppressors were isolated, we determined the crystal structures of three such suppressors to obtain insight into the specific molecular interactions responsible for the observed effects. While most individual suppressors result in small stability enhancements relative to wildtype, which can be combined to yield significant stability increments, thermodynamic stabilisation is neither necessary nor sufficient for suppressor action. Instead, in diverse systems, we observe that individual global suppressors greatly enhance the foldability of buried site mutants, primarily through increase in refolding rate parameters measured in vitro. In the crowded intracellular environment, mutations that slow down folding likely facilitate off-pathway aggregation. We suggest that suppressor mutations that accelerate refolding can counteract this, enhancing the yield of properly folded, functional protein in vivo.

Low-pH Molten Globule-Like Form of CIA17, a Chitooligosaccharide-Specific Lectin from the Phloem Exudate of Coccinia indica, Retains Carbohydrate-Binding Ability

pH-induced changes in the conformation, structural dynamics, and carbohydrate-binding activity of Coccinia indica agglutinin (CIA17), a PP2-type lectin, were investigated employing biophysical approaches. The secondary structure of CIA17 remains nearly unaltered over a wide pH range (2.0-8.5), while the tertiary structure of the protein exhibits considerable changes. A decrease in the fluorescence intensity and excited-state lifetime at low pH indicated perturbation in the local conformation (near Trp residues) of CIA17, which was further supported by enhancement in the Trp accessibility toward charged quenchers under acidic conditions. Fluorescence correlation spectroscopic studies indicated that at pH 2.0, CIA17 exists as a monomer over the concentration range of 10-200 nM and forms dimers at higher concentrations (K_D ∼ 387 nM) but could not form higher oligomers even at ∼150-fold higher concentrations, unlike under native conditions at pH 7.4. Thermal unfolding of the low pH intermediate involves two distinct steps: dissociation of a dimer to a monomer, followed by the unfolding of the monomer. These results strongly suggest that the acid-induced unfolding pathway of CIA17 involves the formation of a monomeric molten globule-like intermediate, which retains appreciable carbohydrate-binding ability. These observations are of great physiological significance since the PP2 proteins are involved in plant defense responses.

Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries

Although there have been recent transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding. However, reliable screens for mutants that improve protein stability do not yet exist, especially for proteins that are well folded and relatively stable. Here, we demonstrate that incorporation of a single, specific, destabilizing mutation termed parent inactivating mutation into each member of a single-site saturation mutagenesis library, followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. We carried out fluorescence-activated cell sorting of such a yeast surface display, saturation suppressor library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations. We found that multiple stabilizing mutations could be identified after a single round of sorting. In addition, multiple libraries with different parent inactivating mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Finally, we show that individual stabilizing mutations could be combined to result in a multi-mutant that demonstrated an increase in thermal melting temperature of about 20 °C, and that displayed enhanced tolerance to high temperature exposure. We conclude that as this method is robust and employs small library sizes, it can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.

Facile measurement of protein stability and folding kinetics using a nano differential scanning fluorimeter

With advancements in high-throughput generation of phenotypic data on mutant proteins, it has become important to individually characterize different proteins or their variants rapidly and with minimal sample consumption. We have made use of a nano differential scanning fluorimetric device, from NanoTemper technologies, to rapidly carry out isothermal chemical denaturation and measure folding/unfolding kinetics of proteins and compared these to corresponding data obtained from conventional spectrofluorimetry. We show that using sample volumes 10-50-fold lower than with conventional fluorimetric techniques, one can rapidly and accurately measure thermodynamic and kinetic stability, as well as folding/unfolding kinetics. This method also facilitates characterization of proteins that are difficult to express and purify.

Molecular Determinants of Temperature-Sensitive Phenotypes

Temperature-sensitive (Ts) mutants are important tools for understanding the role of essential gene(s), but their molecular basis is not well understood. We use CcdB ( Controller of Cell Death protein B) as a model system to explore the effects of Ts mutations on protein stability, folding, and ligand binding. Previously isolated Ts CcdB mutants fall broadly into two categories, namely, buried site (<5% accessibility) and active site (involved in DNA gyrase binding). Several mutants from each category were characterized. It was found that buried-site Ts mutants had decreased stability and foldability, higher aggregation propensity, and, in most cases, reduced affinity for gyrase at both permissive and restrictive temperatures. In contrast, exposed, active-site Ts mutants of CcdB exhibited stability either higher than or similar to that of the wild type and weakened inhibition of DNA gyrase function and/or reduced affinity for gyrase at a higher temperature. At all temperatures, Ts mutations at exposed, active-site residues primarily decrease specific activity without affecting protein levels, while Ts mutations at most buried residues decrease both specific activity and protein levels. Ts phenotypes in both cases arise because total activity is decreased below the threshold required for survival at the restrictive temperature but remains above it at the permissive temperatures. For several mutants, Ts phenotypes were ameliorated upon overexpression of the trigger factor chaperone, suggesting that Ts phenotypes may result from mutational effects on in vivo protein folding rather than on protein stability. This study delineates the diverse factors that contribute to Ts phenotypes. These insights can facilitate rational design of Ts mutants.

CcdB at pH 4 Forms a Partially Unfolded State with a Dry Core

pH is an important factor that affects the protein structure, stability, and activity. Here, we probe the nature of the low-pH structural form of the homodimeric CcdB (controller of cell death B) protein. Characterization of CcdB protein at pH 4 and 300 K using circular dichroism spectroscopy, 8-anilino-1-naphthalene-sulphonate binding, and Trp solvation studies suggests that it forms a partially unfolded state with a dry core at equilibrium under these conditions. CcdB remains dimeric at pH 4 as shown by multiple techniques, such as size-exclusion chromatography coupled to multiangle light scattering, analytical ultracentrifugation, and electron paramagnetic resonance. Comparative analysis using two-dimensional ¹⁵N-¹H heteronuclear single-quantum coherence NMR spectra of CcdB at pH 4 and 7 suggests that the pH 4 and native state have similar but nonidentical structures. Hydrogen-exchange-mass-spectrometry studies demonstrate that the pH 4 state has substantial but anisotropic changes in local stability with core regions close to the dimer interface showing lower protection but some other regions showing higher protection relative to pH 7.

Surface plasmon resonance (SPR) based binding studies of refolded single chain antibody fragments

Recent advances in Recombinant antibody technology / Antibody Engineering has given impetus to the genetic manipulation of antibody fragments that has paved the way for better understanding of the structure and functions of immunoglobulins and also has escalated their use in immunotherapy. Bacterial expression system such as Escherichia coli has complemented this technique through the expression of recombinant antibodies. Present communication has attempted to optimize the expression and refolding protocol of single chain fragment variable (ScFv) and single chain antigen binding fragment (ScFab) using E.coli expression system. Efficiency of refolding protocol was validated by structural analysis by CD, native folding by fluorescence and functional analysis by its binding with full length HIV-1 gp120 via SPR. Results show the predominant β–sheet (CD) as secondary structural content and native folding via red shift (tryptophan fluorescence). The single chain fragments have shown good binding with HIV-1 gp120 thus validating the expression and refolding strategy and also reinstating E.coli as model expression system for recombinant antibody engineering. SPR based binding analysis coupled with E.coli based expression and purification will have implication for HIV therapeutics and will set a benchmark for future studies of similar kind.

•

A scFv having V_H and V_L chains joined through a peptide linker and expressed in E.coli.

•

Functional analysis by SPR show good bindng with full length HIV-1 gp120.

•

Protein engineering facilitated improvised scFv with enhanced affinity and specificity.

Mapping Protein Binding Sites and Conformational Epitopes Using Cysteine Labeling and Yeast Surface Display

We describe a facile method for mapping protein:ligand binding sites and conformational epitopes. The method uses a combination of Cys scanning mutagenesis, chemical labeling, and yeast surface display. While Ala scanning is widely used for similar purposes, often mutation to Ala (or other amino acids) has little effect on binding, except at hotspot residues. Many residues in physical contact with a binding partner are insensitive to substitution with Ala. In contrast, we show that labeling of Cys residues in a binding site consistently abrogates binding. We couple this methodology to yeast surface display and deep sequencing to map conformational epitopes targeted by both monoclonal antibodies and polyclonal sera as well as a protein:ligand binding site. The method does not require purified protein, can distinguish buried and exposed residues, and can be extended to other display formats, including mammalian cells and viruses, emphasizing its wide applicability.

Homodimeric Escherichia coli Toxin CcdB (Controller of Cell Division or Death B Protein) Folds via Parallel Pathways

The existence of parallel pathways in the folding of proteins seems intuitive, yet remains controversial. We explore the folding kinetics of the homodimeric Escherichia coli toxin CcdB (Controller of Cell Division or Death B protein) using multiple optical probes and approaches. Kinetic studies performed as a function of protein and denaturant concentrations demonstrate that the folding of CcdB is a four-state process. The two intermediates populated during folding are present on parallel pathways. Both form by rapid association of the monomers in a diffusion limited manner and appear to be largely unstructured, as they are silent to the optical probes employed in the current study. The existence of parallel pathways is supported by the insensitivity of the amplitudes of the refolding kinetic phases to the different probes used in the study. More importantly, interrupted refolding studies and ligand binding studies clearly demonstrate that the native state forms in a biexponential manner, implying the presence of at least two pathways. Our studies indicate that the CcdA antitoxin binds only to the folded CcdB dimer and not to any earlier folding intermediates. Thus, despite being part of the same operon, the antitoxin does not appear to modulate the folding pathway of the toxin encoded by the downstream cistron. This study highlights the utility of ligand binding in distinguishing between sequential and parallel pathways in protein folding studies, while also providing insights into molecular interactions during folding in Type II toxin-antitoxin systems.

Molecular Determinants of Mutant Phenotypes, Inferred from Saturation Mutagenesis Data

Understanding how mutations affect protein activity and organismal fitness is a major challenge. We used saturation mutagenesis combined with deep sequencing to determine mutational sensitivity scores for 1,664 single-site mutants of the 101 residue Escherichia coli cytotoxin, CcdB at seven different expression levels. Active-site residues could be distinguished from buried ones, based on their differential tolerance to aliphatic and charged amino acid substitutions. At nonactive-site positions, the average mutational tolerance correlated better with depth from the protein surface than with accessibility. Remarkably, similar results were observed for two other small proteins, PDZ domain (PSD95^pdz3) and IgG-binding domain of protein G (GB1). Mutational sensitivity data obtained with CcdB were used to derive a procedure for predicting functional effects of mutations. Results compared favorably with those of two widely used computational predictors. In vitro characterization of 80 single, nonactive-site mutants of CcdB showed that activity in vivo correlates moderately with thermal stability and solubility. The inability to refold reversibly, as well as a decreased folding rate in vitro, is associated with decreased activity in vivo. Upon probing the effect of modulating expression of various proteases and chaperones on mutant phenotypes, most deleterious mutants showed an increased in vivo activity and solubility only upon over-expression of either Trigger factor or SecB ATP-independent chaperones. Collectively, these data suggest that folding kinetics rather than protein stability is the primary determinant of activity in vivo. This study enhances our understanding of how mutations affect phenotype, as well as the ability to predict fitness effects of point mutations.

Use of Protein Folding Reagents

The reagents and methods for purification and use of the most commonly used denaturants, guanidine hydrochloride (guanidine-HCl) and urea, are described. Other protein denaturants and reagents used to fold proteins are briefly mentioned. Sulfhydryl reagents (reducing agents) and "oxido-shuffling" (or oxidative regeneration) systems are also described.

Thermal, chemical and pH induced unfolding of turmeric root lectin: modes of denaturation

Curcuma longa rhizome lectin, of non-seed origin having antifungal, antibacterial and α-glucosidase inhibitory activities, forms a homodimer with high thermal stability as well as acid tolerance. Size exclusion chromatography and dynamic light scattering show it to be a dimer at pH 7, but it converts to a monomer near pH 2. Circular dichroism spectra and fluorescence emission maxima are virtually indistinguishable from pH 7 to 2, indicating secondary and tertiary structures remain the same in dimer and monomer within experimental error. The tryptophan environment as probed by acrylamide quenching data yielded very similar data at pH 2 and pH 7, implying very similar folding for monomer and dimer. Differential scanning calorimetry shows a transition at 350.3 K for dimer and at 327.0 K for monomer. Thermal unfolding and chemical unfolding induced by guanidinium chloride for dimer are both reversible and can be described by two-state models. The temperatures and the denaturant concentrations at which one-half of the protein molecules are unfolded, are protein concentration-dependent for dimer but protein concentration-independent for monomer. The free energy of unfolding at 298 K was found to be 5.23 Kcal mol⁻¹ and 14.90 Kcal mol⁻¹ for the monomer and dimer respectively. The value of change in excess heat capacity upon protein denaturation (ΔC_p) is 3.42 Kcal mol⁻¹ K⁻¹ for dimer. The small ΔC_p for unfolding of CLA reflects a buried hydrophobic core in the folded dimeric protein. These unfolding experiments, temperature dependent circular dichroism and dynamic light scattering for the dimer at pH 7 indicate its higher stability than for the monomer at pH 2. This difference in stability of dimeric and monomeric forms highlights the contribution of inter-subunit interactions in the former.

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Additional role for the ccd operon of F-plasmid as a transmissible persistence factor

Toxin-antitoxin (TA) systems are found on both bacterial plasmids and chromosomes, but in most cases their functional role is unclear. Gene knockouts often yield limited insights into functions of individual TA systems because of their redundancy. The well-characterized F-plasmid-based CcdAB TA system is important for F-plasmid maintenance. We have isolated several point mutants of the toxin CcdB that fail to bind to its cellular target, DNA gyrase, but retain binding to the antitoxin, CcdA. Expression of such mutants is shown to result in release of the WT toxin from a functional preexisting TA complex as well as derepression of the TA operon. One such inactive, active-site mutant of CcdB was used to demonstrate the contribution of CcdB to antibiotic persistence. Transient activation of WT CcdB either by coexpression of the mutant or by antibiotic/heat stress was shown to enhance the generation of drug-tolerant persisters in a process dependent on Lon protease and RecA. An F-plasmid containing a ccd locus can, therefore, function as a transmissible persistence factor.

Structural and thermodynamic characterization of Vibrio fischeri CcdB

CcdB(Vfi) from Vibrio fischeri is a member of the CcdB family of toxins that poison covalent gyrase-DNA complexes. In solution CcdB(Vfi) is a dimer that unfolds to the corresponding monomeric components in a two-state fashion. In the unfolded state, the monomer retains a partial secondary structure. This observation correlates well with the crystal and NMR structures of the protein, which show a dimer with a hydrophobic core crossing the dimer interface. In contrast to its F plasmid homologue, CcdB(Vfi) possesses a rigid dimer interface, and the apparent relative rotations of the two subunits are due to structural plasticity of the monomer. CcdB(Vfi) shows a number of non-conservative substitutions compared with the F plasmid protein in both the CcdA and the gyrase binding sites. Although variation in the CcdA interaction site likely determines toxin-antitoxin specificity, substitutions in the gyrase-interacting region may have more profound functional implications.

A decision tree model for the prediction of homodimer folding mechanism

The formation of protein homodimer complexes for molecular catalysis and regulation is fascinating. The homodimer formation through 2S (2 state), 3SMI (3 state with monomer intermediate) and 3SDI (3 state with dimer intermediate) folding mechanism is known for 47 homodimer structures. Our dataset of forty-seven homodimers consists of twenty-eight 2S, twelve 3SMI and seven 3SDI. The dataset is characterized using monomer length, interface area and interface/total (I/T) residue ratio. It is found that 2S are often small in size with large I/T ratio and 3SDI are frequently large in size with small I/T ratio. Nonetheless, 3SMI have a mixture of these features. Hence, we used these parameters to develop a decision tree model. The decision tree model produced positive predictive values (PPV) of 72% for 2S, 58% for 3SMI and 57% for 3SDI in cross validation. Thus, the method finds application in assigning homodimers with folding mechanism.

Driving forces of gyrase recognition by the addiction toxin CcdB

Gyrase, an essential bacterial topoisomerase, is the target of several antibiotics (e.g. quinolones) as well as of bacterial toxin CcdB. This toxin, encoded by Escherichia coli toxin-antitoxin module ccd, poisons gyrase by causing inhibition of both transcription and replication. Because the molecular driving forces of gyrase unfolding and CcdB-gyrase binding were unknown, the nature of the CcdB-gyrase recognition remained elusive. Therefore, we performed a detailed thermodynamic analysis of CcdB binding to several fragments of gyrase A subunit (GyrA) that contain the CcdB-binding site. Binding of CcdB to the shorter fragments was studied directly by isothermal titration calorimetry. Its binding to the longer GyrA59 fragment in solution is kinetically limited and was therefore investigated via urea induced unfolding of the GyrA59-CcdB complex and unbound GyrA59 and CcdB, monitored by circular dichroism spectroscopy. Model analysis of experimental data, in combination with the relevant structural information, indicates that CcdB binding to gyrase is an enthalpic process driven mainly by specific interactions between CcdB and the highly stable dimerization domain of the GyrA. The dissection of binding energetics indicates that CcdB-gyrase recognition is accompanied by opening of the tower and catalytic domain of GyrA. Such extensive structural rearrangements appear to be crucial driving forces for the functioning of the ccd toxin-antitoxin module.

Refolding and simultaneous purification by three-phase partitioning of recombinant proteins from inclusion bodies

Many recombinant eukaryotic proteins tend to form insoluble aggregates called inclusion bodies, especially when expressed in Escherichia coli. We report the first application of the technique of three-phase partitioning (TPP) to obtain correctly refolded active proteins from solubilized inclusion bodies. TPP was used for refolding 12 different proteins overexpressed in E. coli. In each case, the protein refolded by TPP gave either higher refolding yield than the earlier reported method or succeeded where earlier efforts have failed. TPP-refolded proteins were characterized and compared to conventionally purified proteins in terms of their spectral characteristics and/or biological activity. The methodology is scaleable and parallelizable and does not require subsequent concentration steps. This approach may serve as a useful complement to existing refolding strategies of diverse proteins from inclusion bodies.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Stereochemical criteria for prediction of the effects of proline mutations on protein stability

When incorporated into a polypeptide chain, proline (Pro) differs from all other naturally occurring amino acid residues in two important respects. The φ dihedral angle of Pro is constrained to values close to −65° and Pro lacks an amide hydrogen. Consequently, mutations which result in introduction of Pro can significantly affect protein stability. In the present work, we describe a procedure to accurately predict the effect of Pro introduction on protein thermodynamic stability. Seventy-seven of the 97 non-Pro amino acid residues in the model protein, CcdB, were individually mutated to Pro, and the in vivo activity of each mutant was characterized. A decision tree to classify the mutation as perturbing or nonperturbing was created by correlating stereochemical properties of mutants to activity data. The stereochemical properties including main chain dihedral angle φ and main chain amide H-bonds (hydrogen bonds) were determined from 3D models of the mutant proteins built using MODELLER. We assessed the performance of the decision tree on a large dataset of 163 single-site Pro mutations of T4 lysozyme, 74 nsSNPs, and 52 other Pro substitutions from the literature. The overall accuracy of this algorithm was found to be 81% in the case of CcdB, 77% in the case of lysozyme, 76% in the case of nsSNPs, and 71% in the case of other Pro substitution data. The accuracy of Pro scanning mutagenesis for secondary structure assignment was also assessed and found to be at best 69%. Our prediction procedure will be useful in annotating uncharacterized nsSNPs of disease-associated proteins and for protein engineering and design.

Unlike other amino acids that constitute proteins, Proline is missing a vital hydrogen atom and also bestows local structural rigidity to the three-dimensional (3D) structure of proteins. In some locations, proline can be introduced with little or no detrimental effect to protein function, while at others it is destabilizing and can result in significant degradation or aggregation of the protein. To determine the features of protein 3D structure that tolerate the introduction of prolines, each of the 101 amino acid residues of the protein CcdB were replaced with Proline, and the functional consequence of the mutations were observed. On correlating these data to features of protein 3D structure, a decision tree was generated to predict the functional consequences of proline mutations in proteins of known (or accurately modeled) 3D structure. The performance of the tree was assessed on three different datasets that contained a total of 289 proline mutants in 37 different proteins. The average accuracy of prediction was 75%. The decision tree will be useful in predicting if known but uncharacterized proline mutations in disease-related proteins are likely to have adverse effects. It will also be useful in engineering and designing new proteins and peptides.

The isolated C-terminal domain of Ring1B is a dimer made of stable, well-structured monomers

The Ring1B is a core subunit protein of the PRC1 (polycomb repressive complex 1), which plays key roles in the regulation of the Homeobox gene expression, X-chromosome inactivation, stem cell self-renewal, and tumorigenesis. The C-terminal region of Ring1B interacts with RYBP, a transcriptional repressor in transiently transfected cells, and also with M33, another transcriptional repressor involved in mesoderm patterning. In this work, we show that the C-terminal domain of Ring1B, C-Ring1B, is a dimer in solution, with a dissociation constant of 200 microM, as shown by NMR, ITC, and analytical gel filtration. Each monomer is stable at physiological conditions in a wide pH range ( approximately 5 kcal mol-1 at 298 K), with a well-formed core and a spherical shape. The dimer has a high content of alpha-helix and beta-sheet, as indicated by FTIR spectra, and it is formed by the mutual docking of the preformed folded monomers. Since the C-terminal region is important for interaction with other proteins of the PRC1, the dimerization and the presence of those well-structured monomers might be a form of regulation.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of DeltaCp of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Role of stimuli-sensitive polymers in protein refolding: alpha-amylase and CcdB (controller of cell division or death B) as model proteins

Alginate, a calcium-sensitive polymer, could carry out simultaneous purification and refolding of 8 M urea/100 mM dithiothreitol (DTT) denatured and thermally denatured alpha-amylase present in a commercial preparation. Activity recoveries of 80 and 70% in the former and the latter cases, respectively, were obtained. The fluorescence spectra showed refolding, and PAGE showed the absence of any aggregates in the refolded preparation. As another example, Eudragit S-100, a pH-sensitive poly(methyl methacrylate), was used to refold CcdB (controller of cell division or death B) protein. Initial experiments with wild-type (WT) CcdB showed that Eudragit bound and precipitated (upon lowering the pH to 4.0) CcdB quantitatively from the latter's aqueous solution. The bioconjugate showed DNA gyrase inhibition activity of CcdB and could be recycled. The inclusion bodies of CcdB mutant CcdB-17P were solubilized in 8 M urea/100 mM dithiothreitol. This preparation could be refolded by precipitation with Eudragit. The fluorescence and CD spectra showed that protein refolding has occurred.

Design of immunogens that present the crown of the HIV-1 V3 loop in a conformation competent to generate 447-52D-like antibodies

gp120 is a subunit of the envelope glycoprotein of HIV-1. The third variable loop region of gp120 (V3 loop) contains multiple immunodominant epitopes and is also functionally important for deciding cell-tropism of the virus. 447-52D is a monoclonal antibody that recognizes the conserved tip of the V3 loop in a beta-turn conformation. This antibody has previously been shown to neutralize diverse strains of the virus. In an attempt to generate an immunogen competent to generate 447-52D-like antibodies, the known epitope of 447-52D was inserted at three different surface loop locations in the small, stable protein Escherichia coli Trx (thioredoxin). At one of the three locations (between residues 74 and 75), the insertion was tolerated, the resulting protein was stable and soluble, and bound 447-52D with an affinity similar to that of intact gp120. Upon immunization, the V3 peptide-inserted Trx scaffold was able to generate anti-V3 antibodies that could compete out 447-52D binding to gp120. Epitope mapping studies demonstrated that these anti-V3 antibodies recognized the same epitope as 447-52D. Although the 447-52D-type antibodies were estimated to be present at concentrations of 50-400 microg/ml of serum, these were not able to effect neutralization of strains like JRFL and BAL but could neutralize the sensitive MN strain. The data suggest that because of the low accessibility of the V3 loop on primary isolates such as JRFL, it will be difficult to elicit a V3-specific, 447-52D-like antibody response to effectively neutralize such isolates.

Mutagenesis-based definitions and probes of residue burial in proteins

Every residue of the 101-aa Escherichia coli toxin CcdB was substituted with Ala, Asp, Glu, Lys, and Arg by using site-directed mutagenesis. The activity of each mutant in vivo was characterized as a function of Controller of Cell Division or Death B protein (CcdB) transcriptional level. The mutation data suggest that an accessibility value of 5% is an appropriate cutoff for definition of buried residues. At all buried positions, introduction of Asp results in an inactive phenotype at all CcdB transcriptional levels. The average amount of destabilization upon substitution at buried positions decreases in the order Asp>Glu>Lys>Arg>Ala. Asp substitutions at buried sites in two other proteins, maltose-binding protein and thioredoxin, also were shown to be severely destabilizing. Ala and Asp scanning mutagenesis, in combination with dose-dependent expression phenotypes, was shown to yield important information on protein structure and activity. These results also suggest that such scanning mutagenesis data can be used to rank order sequence alignments and their corresponding homology models, as well as to distinguish between correct and incorrect structural alignments. With continuous reductions in oligonucleotide costs and increasingly efficient site-directed mutagenesis procedures, comprehensive scanning mutagenesis experiments for small proteins/domains are quite feasible.

Structural features differentiate the mechanisms between 2S (2 state) and 3S (3 state) folding homodimers

The formation of homodimer complexes for interface stability, catalysis and regulation is intriguing. The mechanisms of homodimer complexations are even more interesting. Some homodimers form without intermediates (two-state (2S)) and others through the formation of stable intermediates (three-state (3S)). Here, we analyze 41 homodimer (25 2S and 16 3S) structures determined by X-ray crystallography to estimate structural differences between them. The analysis suggests that a combination of structural properties such as monomer length, subunit interface area, ratio of interface to interior hydrophobicity can predominately distinguish 2S and 3S homodimers. These findings are useful in the prediction of homodimer folding and binding mechanisms using structural data.

Oligomerization endows enormous stability to soybean agglutinin: a comparison of the stability of monomer and tetramer of soybean agglutinin

Soybean agglutinin is a tetrameric legume lectin, each of whose subunits are glycosylated. This protein shows a very high degree of stability when compared to the other proteins of the same family. In a previous work, it was shown that the unusual stability of the protein is due to a high degree of subunit interactions. In this study we present the thermodynamic parameters for the stability of soybean agglutinin monomer. The monomeric species is found at pH 2 and below which it is most populated at pH 1.9, as evident from size-exclusion chromatographic and dynamic light scattering studies. The analyses of circular dichroism and fluorescence spectroscopy suggest that the monomer is well folded, and that it has certain characteristic features when compared to its tetrameric counterpart. The conformational stabilities of the tetramer and the monomer at the temperature of their maximum stabilities (310 K) are 59.2 kcal/mol and 9.8 kcal/mol, respectively, indicating that oligomerization contributes significantly to the stability of the native molecule. Also, the T(g) difference for the two forms of the protein is approximately 40 K, whereas the difference in DeltaC(p) is only 1.6 kcal/mol/K. This suggests that the major hydrophobic core is present in the monomer itself, and that oligomerization involves mainly ionic interactions.

Equilibrium unfolding of dimeric and engineered monomeric forms of lambda Cro (F58W) repressor and the effect of added salts: evidence for the formation of folded monomer induced by sodium perchlorate

The equilibrium unfolding transitions of Cro repressor variants, dimeric variant Cro F58W and monomer Cro K56[DGEVK]F58W, have been studied by urea and guanidine hydrochloride to probe the folding mechanism. The unfolding transitions of a dimeric variant are well described by a two state process involving native dimer and unfolded monomer with a free energy of unfolding, DeltaG(0,un)(0), of approximately 10-11 kcal/mol. The midpoint of transition curves is dependent on total protein concentration and DeltaG(0,un)(0) is independent of protein concentration, as expected for this model. Unfolding of Cro monomer is well described by the standard two state model. The stability of both forms of protein increases in the presence of salt but decreases with the decrease in pH. Because of the suggested importance of a N2<-->2F dimerization process in DNA binding, we have also studied the effect of sodium perchlorate, containing the chaotropic perchlorate anion, on the conformational transition of Cro dimer by CD, fluorescence and NMR (in addition to urea and guanidine hydrochloride) in an attempt both to characterize the thermodynamics of the process and to identify conditions that lead to an increase in the population of the folded monomers. Data suggest that sodium perchlorate stabilizes the protein at low concentration (<1.5 M) and destabilizes the protein at higher perchlorate concentration with the formation of a "significantly folded" monomer. The tryptophan residue in the "significantly folded" monomer induced by perchlorate is more exposed to the solvent than in native dimer.