Synthesis of aromatic monothiols and aromatic dithiols to increase the folding rate and yield of disulfide containing proteins

Most pharmaceutically relevant proteins and many extracellular proteins contain disulfide bonds. Formation of the correct disulfide bonds is essential for stability in almost all cases. Disulfide containing proteins can be rapidly and inexpensively overexpressed in bacteria. However, the overexpressed proteins usually form aggregates inside the bacteria, called inclusion bodies, which contains inactive and non-native protein. To obtain native protein, inclusion bodies need to be isolated and resolubilized, and then the resulting protein refolded in vitro. In vitro protein folding is aided by the addition of a redox buffer, which is composed of a small molecule disulfide and/or a small molecule thiol. The most commonly used redox buffer contains reduced and oxidized glutathione. Recently, aliphatic dithiols and aromatic monothiols have been employed as redox buffers. Aliphatic dithiols improved the yield of native protein as compared to the aliphatic thiol, glutathione. Dithiols mimic the in vivo protein folding catalyst, protein disulfide isomerase, which has two thiols per active site. Furthermore, aromatic monothiols increased the folding rate and yield of lysozyme and RNase A relative to glutathione. By combining the beneficial properties of aliphatic dithiols and aromatic monothiols, aromatic dithiols were designed and were expected to increase in vitro protein folding rates and yields. Aromatic monothiols (1-4) and their corresponding disulfides (5-8), two series of ortho- and para-substituted ethylene glycol dithiols (9-15), and a series of aromatic quaternary ammonium salt dithiols (16-17) were synthesized on a multigram scale. Monothiols and disulfides (1-8) were utilized to fold lysozyme and bovine pancreatic trypsin inhibitor. Dithiols (11-17) were tested for their ability to fold lysozyme. At pH 7.0 and pH 8.0, and high protein concentration (1 mg/mL), aromatic dithiols (16, 17) and a monothiol (3) significantly enhanced the in vitro folding rate and yield of lysozyme relative to the aliphatic thiol, glutathione. Additionally, aromatic dithiols (16, 17) significantly enhance the folding yield as compared to the corresponding aromatic monothiol (3). Thus, the folding rate and yield enhancements achieved in in vitro protein folding at high protein concentration will decrease the volume of renaturation solution required for large scale processes and consequently reduce processing time and cost.

Synthesis of aromatic monothiols to increase the folding rate and yields of disulfide-containing proteins

Almost all pharmaceutically relevant proteins and many extracellular proteins contain disulfide bonds, which are essential for protein functions. In many cases, disulfidecontaining proteins are produced via in vitro protein folding that involves the oxidation of reduced protein to native protein, a complex process. The in vitro folding of reduced lysozyme has been extensively studied as a model system because native lysozyme is small, inexpensive, and has only four disulfide bonds. The folding of reduced lysozyme is conducted with the aid of a redox buffer consisting of a small molecule disulfide and a small molecule thiol, such as oxidized and reduced glutathione. Herein, in vitro folding rates and yields of lysozyme obtained in the presence of a series of aromatic thiols and oxidized glutathione are compared to those obtained with reduced and oxidized glutathione. Results showed that aromatic thiols significantly increase the folding rate of lysozyme compared to glutathione.

An investigation into the folding and assembly of engineered antibodies and novel antibody formats

Monoclonal antibodies and novel antibody formats are currently one of the principal therapeutic in the biopharmaceutical industry worldwide and are widely used in the treatment of autoimmune diseases and cancer. It is for this reason that the productivity and quality of antibody production requires improvement; specifically investigations into the engineering of antibodies and any issues that may arise from the production of these therapeutics. The work presented in this thesis describes an investigation into the folding and assembly of seven antibodies plus the novel antibody format FabFv. IgG is comprised of two identical HCs and two identical LCs. The folding process of immunoglobulin is controlled by the CH1 domain within the HC. The CH1 domain remains in a disordered state and is sequestered by BiP in the endoplasmic reticulum. Upon the addition of a folded CL domain, BiP is displaced, the CH1 domain is able to fold and the complete IgG protein can then be secreted from the cell. The results presented in this thesis however, have outlined an additional mechanism for the folding of the CH1 domain. We have shown that the CH1 domain is able to fold in the absence of LC resulting in the secretion of HC dimers in a VH dependent manner. The proposed mechanism for the secretion of HC dimers suggests that some VH domains can interact with each other in order to bring the CH1 domains in close proximity to enable folding to occur. As HC dimer secretion is a hindrance in antibody production, this result has highlighted an engineering target to improve antibody yield. Examination of the folding of IgG4 with the variable region A33 has revealed the inability to secrete LC dimers, cleavage of the HC during expression and secretion of HC dimers in the Fab, FabFv and full length forms. The attributes described have also been shown to be variable region dependent. This has introduced a new concept that the variable domain is important in determining the expression and secretion of antibodies and their individual chains. Pulse chase and 2D gel electrophoresis analysis of the novel antibody format FabFv has revealed that the folding and expression of the LC and HC causes multimeric species of FabFv to be secreted, as opposed to the monomeric form which is the desired therapeutic. Our hypothesis is that this process occurs via a LC dependent mechanism. The proposed hypothesis suggests that further engineering to the LC could diminish the formation and secretion of FabFv multimers. The results from these investigations can be applied to increase the productivity of therapeutics and increase the biological understanding of the domain interactions of IgG during folding, assembly and secretion.

Study of the dynamics of protein folding through minimalistic models

Starting with the Levinthal paradox, a brief introduction to the protein folding problem is presented. The existing theories of protein folding, including the folding funnel scenario, are discussed. After briefly discussing different simulation studies of model proteins, we discuss our recent work on the dynamics of folding of the model HP-36 (the chicken villin headpiece) protein by using a simplified hydropathy scale. Special attention has been paid to the statics and dynamics of contact formation among the hydrophobic residues. The results obtained from this simple model appear to be surprisingly similar to several features observed in the folding of real proteins. The account concludes with a discussion of future problems.

Stretching single polymer chains of donor-acceptor foldamers: toward the quantitative study on the extent of folding

Single-molecule force spectroscopy has proven to be an efficient tool for the quantitative characterization of flexible foldamers on the single-molecule level in this study. The extent of folding has been estimated quantitatively for the first time to the best of our knowledge, which is crucial for a better understanding of the ``folding-process'' on single-molecule level. Therefore, this study may provide a guidance to regulate folding for realizing rational control over the functions of bulk materials.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.

The complex kinetics of protein folding in wide temperature ranges

The complex protein folding kinetics in wide temperature ranges is studied through diffusive dynamics on the underlying energy landscape. The well-known kinetic chevron rollover behavior is recovered from the mean first passage time, with the U-shape dependence on temperature. The fastest folding temperature T-0 is found to be smaller than the folding transition temperature T-f. We found that the fluctuations of the kinetics through the distribution of first passage time show rather universal behavior, from high-temperature exponential Poissonian kinetics to the relatively low-temperature highly nonexponential kinetics. The transition temperature is at T-k and T-0, T-k, T-f. In certain low-temperature regimes, a power law behavior at long time emerges. At very low temperatures ( lower than trapping transition temperature T< T-0/(4&SIM;6)), the kinetics is an exponential Poissonian process again.

Probing the kinetics of single molecule protein folding

We propose an approach to integrate the theory, simulations, and experiments in protein-folding kinetics. This is realized by measuring the mean and high-order moments of the first-passage time and its associated distribution. The full kinetics is revealed in the current theoretical framework through these measurements. In the experiments, information about the statistical properties of first-passage times can be obtained from the kinetic folding trajectories of single molecule experiments ( for example, fluorescence). Theoretical/simulation and experimental approaches can be directly related. We study in particular the temperature-varying kinetics to probe the underlying structure of the folding energy landscape. At high temperatures, exponential kinetics is observed; there are multiple parallel kinetic paths leading to the native state. At intermediate temperatures, nonexponential kinetics appears, revealing the nature of the distribution of local traps on the landscape and, as a result, discrete kinetic paths emerge. At very low temperatures, exponential kinetics is again observed; the dynamics on the underlying landscape is dominated by a single barrier.

Diffusion dynamics, moments, and distribution of first-passage time on the protein-folding energy landscape, with applications to single molecules

We study the dynamics of protein folding via statistical energy-landscape theory. In particular, we concentrate on the local-connectivity case with the folding progress described by the fraction of native conformations. We found that the first passage-time (FPT) distribution undergoes a dynamic transition at a temperature below which the FPT distribution develops a power-law tail, a signature of the intermittent nonexponential kinetic phenomena for the folding dynamics. Possible applications to single-molecule dynamics experiments are discussed.

Osmolyte-induced folding of an intrinsically disordered protein: folding mechanism in the absence of ligand.

Understanding the interconversion between thermodynamically distinguishable states present in a protein folding pathway provides not only the kinetics and energetics of protein folding but also insights into the functional roles of these states in biological systems. The protein component of the bacterial RNase P holoenzyme from Bacillus subtilis (P protein) was previously shown to be unfolded in the absence of its cognate RNA or other anionic ligands. P protein was used in this study as a model system to explore general features of intrinsically disordered protein (IDP) folding mechanisms. The use of trimethylamine N-oxide (TMAO), an osmolyte that stabilizes the unliganded folded form of the protein, enabled us to study the folding process of P protein in the absence of ligand. Transient stopped-flow kinetic traces at various final TMAO concentrations exhibited multiphasic kinetics. Equilibrium "cotitration" experiments were performed using both TMAO and urea during the titration to produce a urea-TMAO titration surface of P protein. Both kinetic and equilibrium studies show evidence of a previously undetected intermediate state in the P protein folding process. The intermediate state is significantly populated, and the folding rate constants are relatively slow compared to those of intrinsically folded proteins similar in size and topology. The experiments and analysis described serve as a useful example for mechanistic folding studies of other IDPs.

The Origin of Nonmonotonic Complex Behavior and the Effects of Nonnative Interactions on the Diffusive Properties of Protein Folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The water factor in the protein-folding problem

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion. Including the configurational dependence will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory also will have to be modified quantitatively.