Small angle X-ray scattering and molecular dynamic simulations provide molecular insight for stability of recombinant human transferrin

•

Combination of SAXS and MD simulations can explain the molecular background for the stability studies.

•

Destabilizing excipients in general induce opening of the N-lobe by binding to specific residues in the protein.

•

The fully open conformation of transferrin is not seen without aggregation.

Transferrin is an attractive candidate for drug delivery due to its ability to cross the blood brain barrier. However, in order to be able to use it for therapeutic purposes, it is important to investigate how its stability depends on different formulation conditions. Combining high-throughput thermal and chemical denaturation studies with small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations, it was possible to connect the stability of transferrin with its conformational changes. Lowering pH induces opening of the transferrin N-lobe, which results in a negative effect on the stability. Presence of NaCl or arginine at low pH enhances the opening and has a negative impact on the overall protein stability.

Statement of Significance

Protein-based therapeutics have become an essential part of medical treatment. They are highly specific, have high affinity and fewer off-target effects. However, stabilization of proteins is critical, time-consuming, and expensive, and it is not yet possible to predict the behavior of proteins under different conditions. The current work is focused on a molecular understanding of the stability of human serum transferrin; a protein which is abundant in blood serum, may pass the blood brain barrier and therefore with high potential in drug delivery. Combination of high throughput unfolding techniques and structural studies, using small angle X-ray scattering and molecular dynamic simulations, allows us to understand the behavior of transferrin on a molecular level.

Concentrated protein solutions investigated using acoustic levitation and small-angle X-ray scattering

An acoustically levitated droplet has been used to collect synchrotron SAXS data on human serum albumin protein solutions up to a protein concentration of 400 mg ml^-1. A careful selection of experiments allows for fast data collection of a large amount of data, spanning a protein concentration/solvent concentration space with limited sample consumption (down to 3 µL per experiment) and few measurements. The data analysis shows data of high quality that are reproducible and comparable with data from standard flow-through capillary-based experiments. Furthermore, using this methodology, it is possible to achieve concentrations that would not be accessible by conventional cells. The protein concentration and ionic strength parameter space diagram may be covered easily and the amount of protein sample is significantly reduced (by a factor of 100 in this work). Used in routine measurements, the benefits in terms of protein cost and time spent are very significant.

Solution structures of long-acting insulin analogues and their complexes with albumin

The lipidation of peptide drugs is one strategy to obtain extended half-lives, enabling once-daily or even less frequent injections for patients. The half-life extension results from a combination of self-association and association with human serum albumin (albumin). The self-association and association with albumin of two insulin analogues, insulin detemir and insulin degludec, were investigated by small-angle X-ray scattering (SAXS) and dynamic light scattering (DLS) in phenolic buffers. Detemir shows concentration-dependent self-association, with an equilibrium between hexamer, dihexamer, trihexamer and larger species, while degludec appears as a dihexamer independent of concentration. The solution structure of the detemir trihexamer has a bent shape. The stoichiometry of the association with albumin was studied using DLS. For albumin-detemir the molar stoichiometry was determined to be 1:6 (albumin:detemir ratio) and for albumin-degludec it was between 1:6 and 1:12 (albumin:degludec ratio). Batch SAXS measurements of a 1:6 albumin:detemir concentration series revealed a concentration dependence of complex formation. The data allowed the modelling of a complex between albumin and a detemir hexamer and a complex consisting of two albumins binding to opposite ends of a detemir dihexamer. Measurements of size-exclusion chromatography coupled to SAXS revealed a complex between a degludec dihexamer and albumin. Based on the results, equilibria for the albumin-detemir and albumin-degludec mixtures are proposed.

Self-Interaction of Human Serum Albumin: A Formulation Perspective

In the present study, small-angle X-ray scattering (SAXS) and static light scattering (SLS) have been used to study the solution properties and self-interaction of recombinant human serum albumin (rHSA) molecules in three pharmaceutically relevant buffer systems. Measurements are carried out up to high protein concentrations and as a function of ionic strength by adding sodium chloride to probe the role of electrostatic interactions. The effective structure factors (S _eff) as a function of the scattering vector magnitude q have been extracted from the scattering profiles and fit to the solution of the Ornstein-Zernike equation using a screened Yukawa potential to describe the double-layer force. Although only a limited q range is used, accurate fits required including an electrostatic repulsion element in the model at low ionic strength, while only a hard sphere model with a tunable diameter is necessary for fitting to high-ionic-strength data. The fit values of net charge agree with available data from potentiometric titrations. Osmotic compressibility data obtained by extrapolating the SAXS profiles or directly from SLS measurements has been fit to a 10-term virial expansion for hard spheres and an equation of state for hard biaxial ellipsoids. We show that modeling rHSA as an ellipsoid, rather than a sphere, provides a much more accurate fit for the thermodynamic data over the entire concentration range. Osmotic virial coefficient data, derived at low protein concentration, can be used to parameterize the model for predicting the behavior up to concentrations as high as 450 g/L. The findings are especially important for the biopharmaceutical sector, which require approaches for predicting concentrated protein solution behavior using minimal sample consumption.

Glucagon-like Peptide 1 Conjugated to Recombinant Human Serum Albumin Variants with Modified Neonatal Fc Receptor Binding Properties. Impact on Molecular Structure and Half-Life

Glucagon-like peptide 1 (GLP-1) is a small incretin hormone stimulated by food intake, resulting in an amplification of the insulin response. Though GLP-1 is interesting as a drug candidate for the treatment of type 2 diabetes mellitus, its short plasma half-life of <3 min limits its clinical use. A strategy for extending the half-life of GLP-1 utilizes the long half-life of human serum albumin (HSA) by combining the two via chemical conjugation or genetic fusion. HSA has a plasma half-life of around 21 days because of its interaction with the neonatal Fc receptor (FcRn) expressed in endothelial cells of blood vessels, which rescues circulating HSA from lysosomal degradation. We have conjugated GLP-1 to C34 of native sequence recombinant HSA (rHSA) and two rHSA variants, one with increased and one with decreased binding affinity for human FcRn. We have investigated the impact of conjugation on FcRn binding affinities, GLP-1 potency, and pharmacokinetics, combined with the solution structure of the rHSA variants and GLP-1-albumin conjugates. The solution structures, determined by small-angle X-ray scattering, show the GLP-1 pointing away from the surface of rHSA. Combining the solution structures with the available structural information about the FcRn and GLP-1 receptor obtained from X-ray crystallography, we can explain the observed in vitro and in vivo behavior. We conclude that the conjugation of GLP-1 to rHSA does not affect the interaction between rHSA and FcRn, while the observed decrease in the potency of GLP-1 can be explained by a steric hindrance of binding of GLP-1 to its receptor.

Oligomerization of a Glucagon-like Peptide 1 Analog: Bridging Experiment and Simulations

The glucagon-like peptide 1 (GLP-1) analog, liraglutide, is a GLP-1 agonist and is used in the treatment of type-2 diabetes mellitus and obesity. From a pharmaceutical perspective, it is important to know the oligomerization state of liraglutide with respect to stability. Compared to GLP-1, liraglutide has an added fatty acid (FA) moiety that causes oligomerization of liraglutide as suggested by small-angle x-ray scattering (SAXS) and multiangle static light scattering (MALS) results. SAXS data suggested a global shape of a hollow elliptical cylinder of size hexa-, hepta-, or octamer, whereas MALS data indicate a hexamer. To elaborate further on the stability of these oligomers and the role of the FA chains, a series of molecular-dynamics simulations were carried out on 11 different hexa-, hepta-, and octameric systems. Our results indicate that interactions of the fatty acid chains contribute noticeably to the stabilization. The simulation results indicate that the heptamer with paired FA chains is the most stable oligomer when compared to the 10 other investigated structures. Theoretical SAXS curves extracted from the simulations qualitatively agree with the experimentally determined SAXS curves supporting the view that liraglutide forms heptamers in solution. In agreement with the SAXS data, the heptamer forms a water-filled oligomer of elliptical cylindrical shape.

Combining SAXS with Rosetta for Identification of Native-like Protein Complexes

It is well-known that recombinant human serum albumin (rHSA) has the ability to stabilize proteins in solution preventing protein adsorption, aggregation and oxidation. For this reason, rHSA is used as an excipient in the formulation of protein pharmaceuticals. To shed light on the molecular interactions, we have studied a variety of protein drugs that are known to bind to rHSA and thereby being stabilized. We observe that the interactions depend on protein concentrations and differ significantly between protein drugs. One approach to study these systems on a molecular level is the combination of small angle X-ray scattering (SAXS) and in-silico modelling. SAXS can be used to identify the overall shape of proteins and protein complexes and ab initio models can be derived from the scattering profile using programs such as Dammif [1]. These programs allow us to assess the overall conformation of the macromolecular structures, but cannot provide detailed information on the molecular level regarding protein-protein interfaces of the complexes. Here, the Rosetta modelling suite, a multipurpose software suite, can be utilized to perform protein-protein docking and to study the complexes. The challenge in using the Rosetta docking tool [2] is the difficulties in efficiently identifying the native-like structure. For better identification we apply SAXS constrains during the docking procedure. Although the method has been applied previously [3], no benchmarking has been published regarding the relative success of using SAXS constrains. We therefore have conducted an elaborate benchmarking, where we have used SAXS constrains for determination of complexes of non-identical components using the Rosetta docking protocol. A pool of complex structures has been chosen to evaluate the difference between conventional docking and docking performed using SAXS constrains. This allows us to optimize the parameters in the protocol and pave the way to study unknown complex structures.

The intrinsically disordered RNR inhibitor Sml1 is a dynamic dimer

Sml1 is a small ribonucleotide reductase (RNR) regulatory protein in Saccharomyces cerevisiae that binds to and inhibits RNR activation. NMR studies of 15N-labeled Sml1 (104 residues), as well as of a truncated variant (residues 50-104), have allowed characterization of their molecular properties. Sml1 belongs to the class of intrinsically disordered proteins with a high degree of dynamics and very little stable structure. Earlier suggestions for a dimeric structure of Sml1 were confirmed, and from translation diffusion NMR measurements, a dimerization dissociation constant of 0.1 mM at 4 degreesC could be determined. The hydrodynamic radius for the monomeric form of Sml1 was determined to be 23.4 A, corresponding to a protein size between those of a globular protein and a coil. Formation of a dimer results in a hydrodynamic radius of 34.4 A. The observed chemical shifts showed in agreement with previous studies two segments with transient helical structure, residues 4-20 and 60-86, and relaxation studies clearly showed restricted motion in these segments. A spin-label attached to C14 showed long-range interactions with residues 60-70 and 85-95, suggesting that the N-terminal domain folds onto the C-terminal domain. Importantly, protease degradation studies combined with mass spectrometry indicated that the N-terminal domain is degraded before the C-terminal region and thus may serve as a protection against proteolysis of the functionally important C-terminal region. Dimer formation was not associated with significant induction of structure but was found to provide further protection against proteolysis. We propose that this molecular shielding and protection of vital functional structures from degradation by functionally unimportant sites may be a general attribute of other natively disordered proteins.