Methionine oxidation of monomeric lambda repressor: the denatured state ensemble under nondenaturing conditions

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔGu°'(H2O), of UBA(1) is 2.4 kcal mol-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, kf, in the absence of a denaturant of 13,000 s-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone 15NH, 13CO, and 13Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

High-throughput oxidation screen of antibody-drug conjugates by analytical protein A chromatography following IdeS digest

OBJECTIVES

Oxidation of protein therapeutics is a major chemical degradation pathway which may impact bioactivity, serum half-life and stability. Therefore, oxidation is a relevant parameter which has to be monitored throughout formulation development. Methods such as HIC, RPLC and LC/MS achieve a separation of oxidized and non-oxidized species by differences in hydrophobicity. Antibody-drug conjugates (ADC) although are highly more complex due to the heterogeneity in linker, drug, drug-to-antibody ratio (DAR) and conjugation site. The analytical protein A chromatography can provide a simple and fast alternative to these common methods.

METHODS

A miniature analytical protein A chromatography method in combination with an IdeS digest was developed to analyse ADCs. The IdeS digest efficiency of an IgG1 was monitored using SEC-HPLC and non-reducing SDS-PAGE. An antibody-fluorescent dye conjugate was conjugated at different dye-to-antibody ratios as model construct to mimic an ADC.

KEY FINDINGS

With IdeS, an almost complete digest of a model IgG1 can be achieved (digested protein amount >98%). This enables subsequent analytical protein A chromatography, which consequently eliminates any interference of payload with the stationary phase.

CONCLUSION

A novel high-throughput method for an interchain cysteine-linked ADC oxidation screens during formulation development was developed.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Can Local Probes Go Global? A Joint Experiment-Simulation Analysis of λ(6-85) Folding

The process of protein folding is known to involve global motions in a cooperative affair; the structure of most of the protein sequences is gained or lost over a narrow range of temperature, denaturant, or pressure perturbations. At the same time, recent simulations and experiments reveal a complex structural landscape with a rich set of local motions and conformational changes. We couple experimental kinetic and thermodynamic measurements with specifically tailored analysis of simulation data to isolate local versus global folding probes. We find that local probes exhibit lower melting temperatures, smaller surface area changes, and faster kinetics compared to global ones. We also see that certain local probes of folding match the global behavior more closely than others. Our work highlights the importance of using multiple probes to fully characterize protein folding dynamics by theory and experiment.

Cooperative folding near the downhill limit determined with amino acid resolution by hydrogen exchange

The relationship between folding cooperativity and downhill, or barrier-free, folding of proteins under highly stabilizing conditions remains an unresolved topic, especially for proteins such as λ-repressor that fold on the microsecond timescale. Under aqueous conditions where downhill folding is most likely to occur, we measure the stability of multiple H bonds, using hydrogen exchange (HX) in a λYA variant that is suggested to be an incipient downhill folder having an extrapolated folding rate constant of 2 × 10(5) s(-1) and a stability of 7.4 kcal·mol(-1) at 298 K. At least one H bond on each of the three largest helices (α1, α3, and α4) breaks during a common unfolding event that reflects global denaturation. The use of HX enables us to both examine folding under highly stabilizing, native-like conditions and probe the pretransition state region for stable species without the need to initiate the folding reaction. The equivalence of the stability determined at zero and high denaturant indicates that any residual denatured state structure minimally affects the stability even under native conditions. Using our ψ analysis method along with mutational ϕ analysis, we find that the three aforementioned helices are all present in the folding transition state. Hence, the free energy surface has a sufficiently high barrier separating the denatured and native states that folding appears cooperative even under extremely stable and fast folding conditions.

Simplified protein models: predicting folding pathways and structure using amino acid sequences

We demonstrate the ability of simultaneously determining a protein's folding pathway and structure using a properly formulated model without prior knowledge of the native structure. Our model employs a natural coordinate system for describing proteins and a search strategy inspired by the observation that real proteins fold in a sequential fashion by incrementally stabilizing nativelike substructures or "foldons." Comparable folding pathways and structures are obtained for the twelve proteins recently studied using atomistic molecular dynamics simulations [K. Lindorff-Larsen, S. Piana, R. O. Dror, D. E. Shaw, Science 334, 517 (2011)], with our calculations running several orders of magnitude faster. We find that nativelike propensities in the unfolded state do not necessarily determine the order of structure formation, a departure from a major conclusion of the molecular dynamics study. Instead, our results support a more expansive view wherein intrinsic local structural propensities may be enhanced or overridden in the folding process by environmental context. The success of our search strategy validates it as an expedient mechanism for folding both in silico and in vivo.

De novo prediction of protein folding pathways and structure using the principle of sequential stabilization

Motivated by the relationship between the folding mechanism and the native structure, we develop a unified approach for predicting folding pathways and tertiary structure using only the primary sequence as input. Simulations begin from a realistic unfolded state devoid of secondary structure and use a chain representation lacking explicit side chains, rendering the simulations many orders of magnitude faster than molecular dynamics simulations. The multiple round nature of the algorithm mimics the authentic folding process and tests the effectiveness of sequential stabilization (SS) as a search strategy wherein 2° structural elements add onto existing structures in a process of progressive learning and stabilization of structure found in prior rounds of folding. Because no a priori knowledge is used, we can identify kinetically significant non-native interactions and intermediates, sometimes generated by only two mutations, while the evolution of contact matrices is often consistent with experiments. Moreover, structure prediction improves substantially by incorporating information from prior rounds. The success of our simple, homology-free approach affirms the validity of our description of the primary determinants of folding pathways and structure, and the effectiveness of SS as a search strategy.

HIC resolution of an IgG1 with an oxidized Trp in a complementarity determining region

Susceptibility of tryptophan (Trp) in a complementarity-determining region (CDR) to oxidation is a significant issue for recombinant monoclonal antibody (mAb) therapeutics due to the clinical efficacy and stability concerns. Here we present a case study using hydrophobic interaction chromatography (HIC) to separate an oxidized Trp containing population of an IgG1. The best separation was achieved using dual Dionex ProPac HIC-10 columns, and the oxidized Trp population was isolated as a separated pre-peak. Peptide map analysis revealed that the oxidized Trp is located in a heavy chain CDR. In addition, the HIC method was capable of monitoring the oxidation status of the CDR Trp, as the oxidation rate of the CDR Trp measured by HIC directly correlated with the results of the peptide maps. The same method conditions were also capable of separating oxidized methionine (Met) and isomerization/deamidation products, which co-elute as another pre-peak at a different retention time from the oxidized Trp species. These observations indicate that the HIC procedure can be utilized to monitor the oxidative status of the CDR Trp in the IgG1.

Non-sequence-specific interactions can account for the compaction of proteins unfolded under "native" conditions

Proteins unfolded by high concentrations of chemical denaturants adopt expanded, largely structure-free ensembles of conformations that are well approximated as random coils. In contrast, globular proteins unfolded under less denaturing conditions (via mutations, or transiently unfolded after a rapid jump to native conditions) and molten globules (arising due to mutations or cosolvents) are often compact. Here we explore the origins of this compaction using a truncated equilibrium-unfolded variant of the 57-residue FynSH3 domain. As monitored by far-UV circular dichroism, NMR spectroscopy, and hydrogen-exchange kinetics, CDelta4 (a 4-residue carboxy-terminal deletion variant of FynSH3) appears to be largely unfolded even in the absence of denaturant. Nevertheless, CDelta4 is quite compact under these conditions, with a hydrodynamic radius only slightly larger than that of the native protein. In order to understand the origins of this molten-globule-like compaction, we have characterized a random sequence polypeptide of identical amino acid composition to CDelta4. Notably, we find that the hydrodynamic radius of this random sequence polypeptide also approaches that of the native protein. Thus, while native-like interactions may contribute to the formation of compact "unfolded" states, it appears that non-sequence-specific monomer-monomer interactions can also account for the dramatic compaction observed for molten globules and the "physiological" unfolded state.

Methionine oxidation in human IgG2 Fc decreases binding affinities to protein A and FcRn

Susceptibility of methionine residues to oxidation is a significant issue of protein therapeutics. Methionine oxidation may limit the product's clinical efficacy or stability. We have studied kinetics of methionine oxidation in the Fc portion of the human IgG2 and its impact on the interaction with FcRn and Protein A. Our results confirm previously published observations for IgG1 that two analogous solvent-exposed methionine residues in IgG2, Met 252 and Met 428, oxidize more readily than the other methionine residue, Met 358, which is buried inside the Fc. Met 397, which is not present in IgG1 but in IgG2, oxidizes at similar rate as Met 358. Oxidation of two labile methionines, Met 252 and Met 428, weakens the binding of the intact antibody with Protein A and FcRn, two natural protein binding partners. Both of these binding partners share the same binding site on the Fc. Additionally, our results shows that Protein A may serve as a convenient and inexpensive surrogate for FcRn binding measurements.

Analysis of the free-energy surface of proteins from reversible folding simulations

Computer generated trajectories can, in principle, reveal the folding pathways of a protein at atomic resolution and possibly suggest general and simple rules for predicting the folded structure of a given sequence. While such reversible folding trajectories can only be determined ab initio using all-atom transferable force-fields for a few small proteins, they can be determined for a large number of proteins using coarse-grained and structure-based force-fields, in which a known folded structure is by construction the absolute energy and free-energy minimum. Here we use a model of the fast folding helical λ-repressor protein to generate trajectories in which native and non-native states are in equilibrium and transitions are accurately sampled. Yet, representation of the free-energy surface, which underlies the thermodynamic and dynamic properties of the protein model, from such a trajectory remains a challenge. Projections over one or a small number of arbitrarily chosen progress variables often hide the most important features of such surfaces. The results unequivocally show that an unprojected representation of the free-energy surface provides important and unbiased information and allows a simple and meaningful description of many-dimensional, heterogeneous trajectories, providing new insight into the possible mechanisms of fast-folding proteins.

The process of protein folding is a complex transition from a disordered to an ordered state. Here, we simulate a specific fast-folding protein at the point at which the native and denatured states are at equilibrium and show that obtaining an accurate description of the mechanisms of folding and unfolding is far from trivial. Using simple quantities which quantify the degree of native order is, in the case of this protein, clearly misleading. We show that an unbiased representation of the free-energy surface can be obtained; using such a representation we are able to redesign the landscape and thus modify, upon site-specific “mutations”, the folding and unfolding rates. This leads us to formulate a hypothesis to explain the very fast folding of many proteins.

Structural determinants of protein folding

The last several decades have seen an explosion of knowledge in the field of structural biology. With critical advances in spectroscopic techniques in examining structures of biomacromolecules, in maturation of molecular biology techniques, as well as vast improvements in computation prowess, protein structures are now being elucidated at an unprecedented rate. In spite of all the recent advances, the protein folding puzzle remains as one of the fundamental biochemical challenges. A facet to this empiric problem is the structural determinants of protein folding. What are the driving forces that pivot a polypeptide chain to a specific conformation amongst the vast conformation space? In this review, we shall discuss some of the structural determinants to protein folding that have been identified in the recent decades.

Small-angle X-ray scattering of reduced ribonuclease A: effects of solution conditions and comparisons with a computational model of unfolded proteins

The disulfide-reduced form of bovine ribonuclease A, with the Cys thiols irreversibly blocked, was characterized by small-angle x-ray scattering. To help resolve the conflicting results and interpretations from previous studies of this model unfolded protein, we measured scattering profiles using a range of solution conditions and compared them with the profiles predicted by a computational model for a random-coil polypeptide. Analysis of the simulated and experimental profiles reveals that scattering intensities at intermediate angles, corresponding to interatomic distances in the range of 5-20 A, are particularly sensitive to changes in solvation and can be used to assess the internal scaling behavior of the polypeptide chain, expressed as a mass fractal dimension, D(m). This region of the scattering curve is also much less sensitive to experimental artifacts than is the very small angle regime (the Guinier region) that has been more typically used to characterize unfolded proteins. The experimental small-angle x-ray scattering profiles closely matched those predicted by the computational model assuming relatively small solvation energies. The scaling behavior of the polypeptide approaches that of a well-solvated polymer under conditions where it has a large net charge and at high urea concentrations. At lower urea concentrations and neutral pH, the behavior of the chain approaches that expected for theta-conditions, where the effects of slightly unfavorable interactions with solvent balance those of excluded volume, leading to scaling behavior comparable to that of an idealized random walk chain. Though detectable, the shift toward more compact conformations at lower urea concentrations does not correspond to a transition to a globule state and is associated with little or no reduction in conformational entropy. This type of collapse, therefore, is unlikely to greatly reduce the conformational search for the native state.

Structural effect of deglycosylation and methionine oxidation on a recombinant monoclonal antibody

Methionine (Met) is one of the most susceptible amino acids to oxidation. Met256 (CH2-Met15.1) and Met432 (CH3-Met107) of a recombinant humanized monoclonal IgG1 antibody are located in the CH2 and CH3 domains, respectively. In three-dimensional structure, these two Met residues are close to the CH2-CH3 interface. In close proximity, oligosaccharides on the conserved asparagine (Asn) residues are enclosed in the CH2 domains. The relationship of Met oxidation with oligosaccharides and their effect on the structure of the antibody was investigated. Removal of oligosaccharides did not alter the oxidation rates of Met256 and Met432, however it caused significant structural changes as evidenced by the susceptibility of the deglycosylated antibody to trypsin and chymotrypsin. Oxidation of Met256 and Met432 did not cause significant conformational changes of the antibody with oligosaccharides, however oxidation of these Met residues accelerated degradation of the deglycosylated antibody. Analysis by mass spectrometry indicated that most of the protease cleavage sites were in the CH2 domains, which suggested that conformational changes induced by the removal of oligosaccharides and further by Met oxidation were local to the CH2 domains.

Backbone dynamics of the monomeric lambda repressor denatured state ensemble under nondenaturing conditions

Oxidizing two native methionine residues predominantly populates the denatured state of monomeric lambda repressor (MetO-lambdaLS) under nondenaturing conditions. NMR was used to characterize the secondary structure and dynamics of MetO-lambdaLS in standard phosphate buffer. 13Calpha and 1Halpha chemical shift indices reveal a region of significant helicity between residues 9 and 29. This helical content is further supported by the observation of medium-range amide NOEs. The remaining residues do not exhibit significant helicity as determined by NMR. We determined 15N relaxation parameters for 64 of 85 residues at 600 and 800 MHz. There are two distinct regions of reduced flexibility, residues 8-32 in the N-terminal third and residues 50-83 in the C-terminal third. The middle third, residues 33-50, has greater flexibility. We have analyzed the amplitude of the backbone motions in terms of the physical properties of the amino acids and conclude that conformational restriction of the backbone MetO-lambdaLS is due to nascent helix formation in the region corresponding to native helix 1. The bulkiness of amino acid residues in the C-terminal third leads to the potential for hydrophobic interactions, which is suggested by chemical exchange detected by the difference in spectral density J(0) at the two static magnetic fields. The more flexible middle region is the result of a predominance of small side chains in this region.

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.

Comparison of methionine oxidation in thermal stability and chemically stressed samples of a fully human monoclonal antibody

Methionine (Met) oxidation is a major degradation pathway of protein therapeutics. Met oxidation of a fully human recombinant monoclonal antibody was investigated under both chemically stressed conditions using tert-butylhydroperoxide (tBHP) and thermal stability conditions where the sample was incubated in formulation buffer at 25 degrees C for 12 months. This antibody has one Met residue on each of the light chains and four Met residues on each of the heavy chains. In the thermal stability sample, only Met residues 256 and 432 in the Fc region were oxidized to form methionine sulfoxide, while Met residues in the Fab region were relatively stable. The susceptibility of Met residues 256 and 432 was further confirmed by incubating samples with tBHP, which has been shown to induce Met oxidation. Further analysis revealed that the susceptible Met residues of each heavy chain were randomly oxidized in samples incubated with tBHP, while in the thermal stability sample, the susceptible Met residues of one heavy chain were preferentially oxidized.

An Isolated Helix Persists in a Sparsely Populated Form of KIX under Native Conditions

NMR relaxation dispersion techniques were used to investigate conformational exchange of the three-helix bundle protein KIX under native conditions. These experiments provide site-resolved kinetic information about microsecond-to-millisecond time scale motions along with structural (chemical shift) information without requiring a perturbation of the equilibrium. All kinetic data are consistent with an apparent two-state transition between natively folded KIX and a partially unfolded high-energy state that is populated to 3.0 +/- 0.2% at 27 degrees C. By combining (13)C- and (15)N-based experiments that probe specific structural aspects, we show that the sparsely populated high-energy state displays a strong conformational preference. An isolated secondary structural element, C-terminal helix alpha3, is highly populated, while the hydrophobic core of the domain and the remainder of the protein backbone, including helices alpha1 and alpha2, are disordered and devoid of specific interactions. This high-energy state presumably represents the equilibrium analogue of a folding intermediate that is transiently populated in stopped-flow kinetic experiments [Horng, J. C., Tracz, S. M., Lumb, K. J., and Raleigh, D. P. (2002) Biochemistry 44, 627-634].