Exploring the Higher Order Structure and Conformational Transitions in Insulin Microcrystalline Biopharmaceuticals by Proton-Detected Solid-State Nuclear Magnetic Resonance at Natural Abundance.

The integrity of a higher order structure (HOS) is an essential requirement to ensure the efficacy, stability, and safety of protein therapeutics. Solution-state nuclear magnetic resonance (NMR) occupies a unique niche as one of the most promising methods to access atomic-level structural information on soluble biopharmaceutical formulations. Another major class of drugs is poorly soluble, such as microcrystalline suspensions, which poses significant challenges for the characterization of the active ingredient in its native state. Here, we have demonstrated a solid-state NMR method for HOS characterization of biopharmaceutical suspensions employing a selective excitation scheme under fast magic angle spinning (MAS). The applicability of the method is shown on commercial insulin suspensions at natural isotopic abundance. Selective excitation aided with proton detection and non-uniform sampling (NUS) provides improved sensitivity and resolution. The enhanced resolution enabled us to demonstrate the first experimental evidence of a phenol-escaping pathway in insulin, leading to conformational transitions to different hexameric states. This approach has the potential to serve as a valuable means for meticulously examining microcrystalline biopharmaceutical suspensions, which was previously not attainable in their native formulation states and can be seamlessly extended to other classes of biopharmaceuticals such as mAbs and other microcrystalline proteins.

Mapping of Methyl Epitopes of a Peptide-Drug with Its Receptor by 2D STDD-Methyl TROSY NMR Spectroscopy.

The current trend in the biopharmaceutical market has boosted the development and production of biological drugs with high efficacy and fidelity for receptor binding. While high-resolution structural insights into binding epitopes of the receptor are indispensable for better therapeutic design, it is tedious and costly. In this work, we develop a protocol by integrating two well-known NMR-based solution-state methods. Saturation transfer double-difference with methyl-TROSY (STDD-Methyl TROSY NMR) was used to probe methyl binding epitopes of the ligand in a label-free environment. This study was carried out with Human insulin as a model peptide drug, with the insulin growth factor receptor (IGFR), which is an off-target receptor for insulin. Methyl epitopes identified from STDD-Methyl TROSY NMR spectroscopy were validated through the HADDOCK platform to generate a drug-receptor model. Since this method can be applied at natural abundance, it has the potential to screen a large set of peptide-drug interactions for optimum receptor binding. Thus, we propose STDD-Methyl TROSY NMR spectroscopy as a technique for rapid screening of biologics for the development of optimized biopharmaceutics.

Physicochemical and functional characterization of MYL-1501D, a proposed biosimilar to insulin glargine.

Insulin glargine is a long-acting analogue of human insulin that has been used to manage hyperglycemia in patients with diabetes mellitus (DM) for nearly 20 years. Insulin glargine has a relatively constant concentration-time profile that mimics basal levels of insulin and allows for once-daily administration. MYL-1501D is a biosimilar insulin glargine designed to offer greater access of insulin glargine to patients, with comparable efficacy and safety to the marketed reference product. We conducted a comprehensive panel of studies based on a formal analysis of critical quality attributes to characterize the structural and functional properties of MYL-1501D and reference insulin glargine products available in the United States and European Union. MYL-1501D was comprehensively shown to have high similarity to the reference products in terms of protein structure, metabolic activity (both in vitro cell-based assays and in vivo rabbit bioassays), and in vitro cell-based assays for mitogenic activity. The structural analyses demonstrated that the primary protein sequence was identical, and secondary and tertiary structures are similar between the proposed biosimilar and the reference products. Insulin receptor binding affinity and phosphorylation studies also established analytical similarity. MYL-1501D demonstrated high similarity in different metabolic assays of glucose uptake, adipogenesis activity, and inhibition of stimulated lipolysis. Rabbit bioassay studies showed MYL-1501D and EU-approved insulin glargine are highly similar to US-licensed insulin glargine. These product quality studies show high similarity between MYL-1501D and licensed or approved insulin glargine products and suggest the potential of MYL-1501D as an alternative cost-effective treatment option for patients and clinicians.

Simultaneous characterization of insulin HMWP and protamine sulphate in complex formulations through SEC-coupled mass spectrometry.

High molecular weight protein aggregates present in a recombinant human insulin and analogs are conventionally quantified by SEC-HPLC and identified by SEC-MALS as oligomeric population which lacks precise identification of species. The limitation of these two techniques is vanquished though simultaneous separation and identification by SEC coupled with MS. The identification was established with organic solvent based isocratic elution followed by MS for parallel separation and identification of HMWP species. The developed SEC-MS method is validated to establish the method capability and variability. Further investigations under stress conditions of Insulin analogues revealed the capability of the method to separate higher order oligomeric (Trimeric, and Tetrameric) species alongside covalent dimeric species. Additionally, the method holds good in separating and sequencing protamine peptides used in suspension (Neutral Protamine Hagedorn) and biphasic/mixed (70/30) formulations of Human insulin using ETD-MSMS. The data presented here shows insight towards utilization of state-of-the-art SEC-MS technique in the biopharmaceutical field as a tool to establish the structural comparability of higher order species in biosimilars and to investigate the lot to lot batch variability for protamine sulphate in-terms of sequence confirmation.

Large scale purification and characterization of A21 deamidated variant-most prominent post translational modification (PTM) for insulins which is also widely observed in insulins pharmaceutical manufacturing and storage.

Biopharmaceutical development demands appropriate understanding of product related variants, which are formed due to post-translational modification and during downstream processing. These variants can lead to low yield, reduced biological activity, and suboptimal product quality. In addition, these variants may undergo immune reactions, henceforth need to be appropriately controlled to ensure consistent product quality and patient safety. Deamidation of insulin is the most common post-translational modification occurring in insulin and insulin analogues. AsnA21 desamido variant is also the most prominent product variant formed during human insulin manufacturing process and/or during the storage. Often, this deamidated variant is used as an impurity standard during in-process and final product analysis in the QC system. However, purification of large quantity of purified deamidated material is always being challenging due to highly similar mass, ionic, hydrophobic properties, and high structural similarity of the variant compared to the parent product. Present work demonstrates the simplified and efficient scalable process for generation of AsnA21 deamidated variant in powder form with ~96% purity. The mixed-mode property of anion exchange resin PolyQuat was utilized to purify the deamidated impurity with high recovery. Subsequent reversed-phase high performance liquid chromatography (RP-HPLC) step was introduced for concentration of product in bind elute mode. Elution pool undergone isoelectric precipitation and lyophilisation. The lyophilized product allows users for convenient use of the deamidated impurity for intended purposes. Detailed characterization by Mass spectrometry revealed deamidation is at AsnA21 and further confirmed that, structural and functional characterization as well as the biological activity of isolated variant is equivalent to insulin.

Dynamic studies of H-Ras•GTPγS interactions with nucleotide exchange factor Sos reveal a transient ternary complex formation in solution.

The cycling between GDP- and GTP- bound forms of the Ras protein is partly regulated by the binding of Sos. The structural/dynamic behavior of the complex formed between activated Sos and Ras at the point of the functional cycle where the nucleotide exchange is completed has not been described to date. Here we show that solution NMR spectra of H-Ras∙GTPγS mixed with a functional fragment of Sos (Sos(Cat)) at a 2:1 ratio are consistent with the formation of a rather dynamic assembly. H-Ras∙GTPγS binding was in fast exchange on the NMR timescale and retained a significant degree of molecular tumbling independent of Sos(Cat), while Sos(Cat) also tumbled largely independently of H-Ras. Estimates of apparent molecular weight from both NMR data and SEC-MALS revealed that, at most, only one H-Ras∙GTPγS molecule appears stably bound to Sos. The weak transient interaction between Sos and the second H-Ras∙GTPγS may provide a necessary mechanism for complex dissociation upon the completion of the native GDP → GTP exchange reaction, but also explains measurable GTP → GTP exchange activity of Sos routinely observed in in vitro assays that use fluorescently-labelled analogs of GTP. Overall, the data presents the first dynamic snapshot of Ras functional cycle as controlled by Sos.

Monitoring Ras Interactions with the Nucleotide Exchange Factor Son of Sevenless (Sos) Using Site-specific NMR Reporter Signals and Intrinsic Fluorescence.

The activity of Ras is controlled by the interconversion between GTP- and GDP-bound forms partly regulated by the binding of the guanine nucleotide exchange factor Son of Sevenless (Sos). The details of Sos binding, leading to nucleotide exchange and subsequent dissociation of the complex, are not completely understood. Here, we used uniformly (15)N-labeled Ras as well as [(13)C]methyl-Met,Ile-labeled Sos for observing site-specific details of Ras-Sos interactions in solution. Binding of various forms of Ras (loaded with GDP and mimics of GTP or nucleotide-free) at the allosteric and catalytic sites of Sos was comprehensively characterized by monitoring signal perturbations in the NMR spectra. The overall affinity of binding between these protein variants as well as their selected functional mutants was also investigated using intrinsic fluorescence. The data support a positive feedback activation of Sos by Ras·GTP with Ras·GTP binding as a substrate for the catalytic site of activated Sos more weakly than Ras·GDP, suggesting that Sos should actively promote unidirectional GDP → GTP exchange on Ras in preference of passive homonucleotide exchange. Ras·GDP weakly binds to the catalytic but not to the allosteric site of Sos. This confirms that Ras·GDP cannot properly activate Sos at the allosteric site. The novel site-specific assay described may be useful for design of drugs aimed at perturbing Ras-Sos interactions.

Structural and dynamic insights into the energetics of activation loop rearrangement in FGFR1 kinase.

Protein tyrosine kinases differ widely in their propensity to undergo rearrangements of the N-terminal Asp-Phe-Gly (DFG) motif of the activation loop, with some, including FGFR1 kinase, appearing refractory to this so-called 'DFG flip'. Recent inhibitor-bound structures have unexpectedly revealed FGFR1 for the first time in a 'DFG-out' state. Here we use conformationally selective inhibitors as chemical probes for interrogation of the structural and dynamic features that appear to govern the DFG flip in FGFR1. Our detailed structural and biophysical insights identify contributions from altered dynamics in distal elements, including the αH helix, towards the outstanding stability of the DFG-out complex with the inhibitor ponatinib. We conclude that the αC-ß4 loop and 'molecular brake' regions together impose a high energy barrier for this conformational rearrangement, and that this may have significance for maintaining autoinhibition in the non-phosphorylated basal state of FGFR1.

NMR backbone assignments of the tyrosine kinase domain of human fibroblast growth factor receptor 1.

Members of the fibroblast growth factor receptor tyrosine kinase family (FGFR1-4) play an important role in many signalling cascades. Although tightly regulated, aberrant activity of these enzymes may lead to, or become features of, disease pathologies including cancer. FGFR isoforms have been the subject of drug discovery programmes, with a number of kinase-domain inhibitors in pre-clinical and clinical development. Here, we present the first (83% complete) backbone resonance assignments of apo-FGFR1 kinase.

Kinase in motion: insights into the dynamic nature of p38α by high-pressure NMR spectroscopic studies.

Protein kinases are highly dynamic and complex molecules. Here we present high-pressure and relaxation studies of the activated p38α mitogen-activated protein kinase (MAPK). p38α plays a central role in inflammatory diseases such as rheumatoid arthritis and is therefore a highly attractive pharmaceutical target. The combination of high pressure and NMR spectroscopy allowed for a detailed per-residue based assessment of the structural plasticity of p38α and the accessibility of low-lying excited-energy conformations throughout the kinase structure. Such information is uniquely accessible through the combination of liquid-state NMR and high pressure and is of considerable value for the drug discovery process. The interactions of p38α and DFG-in and DFG-out ligands were studied under the application of high pressure, and we demonstrate how we can alter kinase dynamics by pressure in a similar way to what has previously only been observed by ligand binding. Pressure is shown to be a mild and efficient tool for manipulation of intermediate-timescale dynamics.

High-pressure NMR reveals close similarity between cold and alcohol protein denaturation in ubiquitin.

Proteins denature not only at high, but also at low temperature as well as high pressure. These denatured states are not easily accessible for experiment, because usually heat denaturation causes aggregation, whereas cold or pressure denaturation occurs at temperatures well below the freezing point of water or pressures above 5 kbar, respectively. Here we have obtained atomic details of the pressure-assisted, cold-denatured state of ubiquitin at 2,500 bar and 258 K by high-resolution NMR techniques. Under these conditions, a folded, native-like and a disordered state exist in slow exchange. Secondary chemical shifts show that the disordered state has structural propensities for a native-like N-terminal ß-hairpin and α-helix and a nonnative C-terminal α-helix. These propensities are very similar to the previously described alcohol-denatured (A-)state. Similar to the A-state, (15)N relaxation data indicate that the secondary structure elements move as independent segments. The close similarity of pressure-assisted, cold-denatured, and alcohol-denatured states with native and nonnative secondary elements supports a hierarchical mechanism of folding and supports the notion that similar to alcohol, pressure and cold reduce the hydrophobic effect. Indeed, at nondenaturing concentrations of methanol, a complete transition from the native to the A-state can be achieved at ambient temperature by varying the pressure from 1 to 2,500 bar. The methanol-assisted pressure transition is completely reversible and can also be induced in protein G. This method should allow highly detailed studies of protein-folding transitions in a continuous and reversible manner.

Residual dipolar couplings measured in unfolded proteins are sensitive to amino-acid-specific geometries as well as local conformational sampling.

Many functional proteins do not have well defined folded structures. In recent years, both experimental and computational approaches have been developed to study the conformational behaviour of this type of protein. It has been shown previously that experimental RDCs (residual dipolar couplings) can be used to study the backbone sampling of disordered proteins in some detail. In these studies, the backbone structure was modelled using a common geometry for all amino acids. In the present paper, we demonstrate that experimental RDCs are also sensitive to the specific geometry of each amino acid as defined by energy-minimized internal co-ordinates. We have modified the FM (flexible-Meccano) algorithm that constructs conformational ensembles on the basis of a statistical coil model, to account for these differences. The modified algorithm inherits the advantages of the FM algorithm to efficiently sample the potential energy landscape for coil conformations. The specific geometries incorporated in the new algorithm result in a better reproduction of experimental RDCs and are generally applicable for further studies to characterize the conformational properties of intrinsically disordered proteins. In addition, the internal-co-ordinate-based algorithm is an order of magnitude more efficient, and facilitates side-chain construction, surface osmolyte simulation, spin-label distribution sampling and proline cis/trans isomer simulation.

Facile backbone structure determination of human membrane proteins by NMR spectroscopy.

Although nearly half of today's major pharmaceutical drugs target human integral membrane proteins (hIMPs), only 30 hIMP structures are currently available in the Protein Data Bank, largely owing to inefficiencies in protein production. Here we describe a strategy for the rapid structure determination of hIMPs, using solution NMR spectroscopy with systematically labeled proteins produced via cell-free expression. We report new backbone structures of six hIMPs, solved in only 18 months from 15 initial targets. Application of our protocols to an additional 135 hIMPs with molecular weight <30 kDa yielded 38 hIMPs suitable for structural characterization by solution NMR spectroscopy without additional optimization.

Side-chain chi(1) conformations in urea-denatured ubiquitin and protein G from (3)J coupling constants and residual dipolar couplings.

Current NMR information on side-chain conformations of unfolded protein states is sparse due to the poor dispersion particularly of side-chain proton resonances. We present here optimized schemes for the detection of (3)J(HalphaHbeta), (3)J(NHbeta), and (3)J(C'Hbeta) scalar and (1)D(CbetaHbeta) residual dipolar couplings (RDCs) in unfolded proteins. For urea-denatured ubiquitin and protein G, up to six (3)J-couplings to (1)H(beta) are detected, which define the chi(1) angle at very high precision. Interpretation of the (3)J couplings by a model of mixed staggered chi(1) rotamers yields excellent agreement and also provides stereoassignments for (1)H(beta) methylene protons. For all observed amino acids with the exception of leucine, the chemical shift of (1)H(beta3) protons was found downfield from (1)H(beta2). For most residues, the precision of individual chi(1) rotamer populations is better than 2%. The experimental chi(1) rotamer populations are in the vicinity of averages obtained from coil regions in folded protein structures. However, individual variations from these averages of up to 40% are highly significant and indicate sequence- and residue-specific interactions. Particularly strong deviations from the coil average are found for serine and threonine residues, an effect that may be explained by a weakening of side-chain to backbone hydrogen bonds in the urea-denatured state. The measured (1)D(CbetaHbeta) RDCs correlate well with predicted RDCs that were calculated from a sterically aligned coil model ensemble and the (3)J-derived chi(1) rotamer populations. This agreement supports the coil model as a good first approximation of the unfolded state. Deviations between measured and predicted values at certain sequence locations indicate that the description of the local backbone conformations can be improved by incorporation of the RDC information. The ease of detection of a large number of highly precise side-chain RDCs opens the possibility for a more rigorous characterization of both side-chain and backbone conformations in unfolded proteins.

Targeting Bcr-Abl by combining allosteric with ATP-binding-site inhibitors.

In an effort to find new pharmacological modalities to overcome resistance to ATP-binding-site inhibitors of Bcr-Abl, we recently reported the discovery of GNF-2, a selective allosteric Bcr-Abl inhibitor. Here, using solution NMR, X-ray crystallography, mutagenesis and hydrogen exchange mass spectrometry, we show that GNF-2 binds to the myristate-binding site of Abl, leading to changes in the structural dynamics of the ATP-binding site. GNF-5, an analogue of GNF-2 with improved pharmacokinetic properties, when used in combination with the ATP-competitive inhibitors imatinib or nilotinib, suppressed the emergence of resistance mutations in vitro, displayed additive inhibitory activity in biochemical and cellular assays against T315I mutant human Bcr-Abl and displayed in vivo efficacy against this recalcitrant mutant in a murine bone-marrow transplantation model. These results show that therapeutically relevant inhibition of Bcr-Abl activity can be achieved with inhibitors that bind to the myristate-binding site and that combining allosteric and ATP-competitive inhibitors can overcome resistance to either agent alone.

Solution conformations and dynamics of ABL kinase-inhibitor complexes determined by NMR substantiate the different binding modes of imatinib/nilotinib and dasatinib.

Current structural understanding of kinases is largely based on x-ray crystallographic studies, whereas very little data exist on the conformations and dynamics that kinases adopt in the solution state. ABL kinase is an important drug target in the treatment of chronic myelogenous leukemia. Here, we present the first characterization of ABL kinase in complex with three clinical inhibitors (imatinib, nilotinib, and dasatinib) by modern solution NMR techniques. Structural and dynamical results were derived from complete backbone resonance assignments, experimental residual dipolar couplings, and (15)N relaxation data. Residual dipolar coupling data on the imatinib and nilotinib complexes show that the activation loop adopts the inactive conformation, whereas the dasatinib complex preserves the active conformation, which does not support contrary predictions based upon molecular modeling. Nanosecond as well as microsecond dynamics can be detected for certain residues in the activation loop in the inactive and active conformation complexes.

Backbone NMR resonance assignment of the Abelson kinase domain in complex with imatinib.

Imatinib (Glivec or Gleevec) potently inhibits the tyrosine kinase activity of BCR-ABL, a constitutively activated kinase, which causes chronic myelogenous leukemia (CML). Here we report the first almost complete backbone assignment of c-ABL kinase domain in complex with imatinib.

Residual dipolar couplings in short peptides reveal systematic conformational preferences of individual amino acids.

Residual dipolar couplings (RDCs) observed by NMR in solution under weak alignment conditions can monitor average net orientations and order parameters of individual bonds. By their simple geometrical dependence, RDCs bear particular promise for the quantitative characterization of conformations in partially folded or unfolded proteins. We have systematically investigated the influence of amino acid substitutions X on the conformation of unfolded model peptides EGAAXAASS as monitored by their (1)H(Nu)-(15)N and (1)H(alpha)-(13)C(alpha) RDCs detected at natural abundance of (15)N and (13)C in strained polyacrylamide gels. In total, 14 single amino acid substitutions were investigated. The RDCs show a specific dependence on the substitution X that correlates to steric or hydrophobic interactions with adjacent amino acids. In particular, the RDCs for the glycine and proline substitutions indicate less or more order, respectively, than the other amino acids. The RDCs for aromatic substitutions tryptophane and tyrosine give evidence of a kink in the peptide backbone. This effect is also observable for orientation by Pf1 phages and corroborated by variations in (13)C(alpha) secondary shifts and (3)J(HNH)(alpha) scalar couplings in isotropic samples. RDCs for a substitution with the beta-turn sequence KNGE differ from single amino acid substitutions. Terminal effects and next neighbor effects could be demonstrated by further specific substitutions. The results were compared to statistical models of unfolded peptide conformations derived from PDB coil subsets, which reproduce overall trends for (1)H(Nu)-(15)N RDCs for most substitutions, but deviate more strongly for (1)H(alpha)-(13)C(alpha) RDCs. The outlined approach opens the possibility to obtain a systematic experimental characterization of the influence of individual amino acid/amino acid interactions on orientational preferences in polypeptides.