HDX and In Silico Docking Reveal that Excipients Stabilize G-CSF via a Combination of Preferential Exclusion and Specific Hotspot Interactions

Investigation of the Solid-State Interactions in Lyophilized Human G-CSF Using Hydrogen-Deuterium Exchange Mass Spectrometry

Hydrogen/deuterium exchange mass spectrometry (HDX-MS) previously elucidated the interactions between excipients and proteins for liquid granulocyte colony stimulating factor (G-CSF) formulations, confirming predictions made using computational structure docking. More recently, solid-state HDX mass spectrometry (ssHDX-MS) was developed for proteins in the lyophilized state. Deuterium uptake in ssHDX-MS has been shown for various proteins, including monoclonal antibodies, to be highly correlated with storage stability, as measured by protein aggregation and chemical degradation. As G-CSF is known to lose activity through aggregation upon lyophilization, we applied the ssHDX-MS method with peptide mapping to four different lyophilized formulations of G-CSF to compare the impact of three excipients on local structure and exchange dynamics. HDX at 22 °C was confirmed to correlate well with the monomer content remaining after lyophilization and storage at -20 °C, with sucrose providing the greatest protection, and then phenylalanine, mannitol, and no excipient leading to progressively less protection. Storage at 45 °C led to little difference in final monomer content among the formulations, and so there was no discernible relationship with total deuterium uptake on ssHDX. Incubation at 45 °C may have led to a structural conformation and/or aggregation mechanism no longer probed by HDX at 22 °C. Such a conformational change was observed previously at 37 °C for liquid-formulated G-CSF using NMR. Peptide mapping revealed that tolerance to lyophilization and -20 °C storage was linked to increased stability in the small helix, loop AB, helix C, and loop CD. LC-MS HDX and NMR had previously linked loop AB and loop CD to the formation of a native-like state (N*) prior to aggregation in liquid formulations, suggesting a similar structural basis for G-CSF aggregation in the liquid and solid states.

NMR Reveals Functionally Relevant Thermally Induced Structural Changes within the Native Ensemble of G-CSF

Structure–function relationships in proteins refer to a trade-off between stability and bioactivity, molded by evolution of the molecule. Identifying which protein amino acid residues jeopardize global or local stability for the benefit of bioactivity would reveal residues pivotal to this structure–function trade-off. Here, we use ¹⁵N–¹H heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectroscopy to probe the microenvironment and dynamics of residues in granulocyte colony-stimulating factor (G-CSF) through thermal perturbation. From this analysis, we identified four residues (G4, A6, T133, and Q134) that we classed as significant to global stability, given that they all experienced large environmental and dynamic changes and were closely correlated to each other in their NMR characteristics. Additionally, we observe that roughly four structural clusters are subject to localized conformational changes or partial unfolding prior to global unfolding at higher temperature. Combining NMR observables with structure relaxation methods reveals that these structural clusters concentrate around loop AB (binding site III inclusive). This loop has been previously implicated in conformational changes that result in an aggregation prone state of G-CSF. Residues H43, V48, and S63 appear to be pivotal to an opening motion of loop AB, a change that is possibly also important for function. Hence, we present here an approach to profiling residues in order to highlight their potential roles in the two vital characteristics of proteins: stability and bioactivity.

Biological activities of drug inactive ingredients

In a drug formulation (DFM), the major components by mass are not Active Pharmaceutical Ingredient (API) but rather Drug Inactive Ingredients (DIGs). DIGs can reach much higher concentrations than that achieved by API, which raises great concerns about their clinical toxicities. Therefore, the biological activities of DIG on physiologically relevant target are widely demanded by both clinical investigation and pharmaceutical industry. However, such activity data are not available in any existing pharmaceutical knowledge base, and their potentials in predicting the DIG-target interaction have not been evaluated yet. In this study, the comprehensive assessment and analysis on the biological activities of DIGs were therefore conducted. First, the largest number of DIGs and DFMs were systematically curated and confirmed based on all drugs approved by US Food and Drug Administration. Second, comprehensive activities for both DIGs and DFMs were provided for the first time to pharmaceutical community. Third, the biological targets of each DIG and formulation were fully referenced to available databases that described their pharmaceutical/biological characteristics. Finally, a variety of popular artificial intelligence techniques were used to assess the predictive potential of DIGs' activity data, which was the first evaluation on the possibility to predict DIG's activity. As the activities of DIGs are critical for current pharmaceutical studies, this work is expected to have significant implications for the future practice of drug discovery and precision medicine.

Investigation of the pH-dependent aggregation mechanisms of GCSF using low resolution protein characterization techniques and advanced molecular dynamics simulations

Granulocyte-colony stimulating factor (GCSF) is a widely used therapeutic protein to treat neutropenia. GCSF has an increased propensity to aggregate if the pH is increased above 5.0. Although GCSF is very well experimentally characterized, the exact pH-dependent aggregation mechanism of GCSF is still under debate. This study aimed to model the complex pH-dependent aggregation behavior of GCSF using state-of-the-art simulation techniques. The conformational stability of GCSF was investigated by performing metadynamics simulations, while the protein-protein interactions were investigated using coarse-grained (CG) simulations of multiple GCSF monomers. The CG simulations were directly compared with small-angle X-ray (SAXS) data. The metadynamics simulations demonstrated that the orientations of Trp residues in GCSF are dependent on pH. The conformational change of Trp residues is due to the loss of Trp-His interactions at the physiological pH, which in turn may increase protein flexibility. The helical structure of GCSF was not affected by the pH conditions of the simulations. Our CG simulations indicate that at pH 4.0, the colloidal stability may be more important than the conformational stability of GCSF. The electrostatic potential surface and CG simulations suggested that the basic residues are mainly responsible for colloidal stability as deprotonation of these residues causes a reduction of the highly positively charged electrostatic barrier close to the aggregation-prone long loop regions.

Protein Engineering and HDX Identify Structural Regions of G-CSF Critical to Its Stability and Aggregation

The protein engineering and formulation of therapeutic proteins for prolonged shelf-life remain a major challenge in the biopharmaceutical industry. Understanding the influence of mutations and formulations on the protein structure and dynamics could lead to more predictive approaches to their improvement. Previous intrinsic fluorescence analysis of the chemically denatured granulocyte colony-stimulating factor (G-CSF) suggested that loop AB could subtly reorganize to form an aggregation-prone intermediate state. Hydrogen deuterium exchange mass spectrometry (HDX-MS) has also revealed that excipient binding increased the thermal unfolding transition midpoint (T_m) by stabilizing loop AB. Here, we have combined protein engineering with biophysical analyses and HDX-MS to reveal that increased exchange in a core region of the G-CSF comprising loop AB (ABI, a small helix, ABII) and loop CD packed onto helix B and the beginning of loop BC leads to a decrease in T_m and higher aggregation rates. Furthermore, some mutations can increase the population of the aggregation-prone conformation within the native ensemble, as measured by the greater local exchange within this core region.

Recent advancements in mass spectrometry for higher order structure characterization of protein therapeutics

Molecular characterization of higher order structure (HOS) in protein therapeutics is crucial to the selection of candidate molecules, understanding of structure-function relationships, formulation development, stability assessment, and comparability studies. Recent advances in mass spectrometry (MS), including native MS, hydrogen/deuterium exchange (HDX)-MS, and fast photochemical oxidation of proteins (FPOP) coupled with MS, have provided orthogonal ways to characterize HOS of protein therapeutics. In this review, we present the utility of native MS, HDX-MS and FPOP-MS in protein therapeutics discovery and development, with a focus on epitope mapping, aggregation assessment, and comparability studies. We also discuss future trends in the application of these MS methods to HOS characterization.

The Role of Cyclodextrins against Interface-Induced Denaturation in Pharmaceutical Formulations: A Molecular Dynamics Approach

Protein-based pharmaceutical products are subject to a variety of environmental stressors, during both production and shelf-life. In order to preserve their structure, and, therefore, functionality, it is necessary to use excipients as stabilizing agents. Among the eligible stabilizers, cyclodextrins (CDs) have recently gained interest in the scientific community thanks to their properties. Here, a computational approach is proposed to clarify the role of β-cyclodextrin (βCD) and 2-hydroxypropyl-β-cyclodextrin (HPβCD) against granulocyte colony-stimulating (GCSF) factor denaturation at the air–water and ice–water interfaces, and also in bulk water at 300 or 260 K. Both traditional molecular dynamics (MD) simulations and enhanced sampling techniques (metadynamics, MetaD) are used to shed light on the underlying molecular mechanisms. Bulk simulations revealed that CDs were preferentially included within the surface hydration layer of GCSF, and even included some peptide residues in their hydrophobic cavity. HPβCD was able to stabilize the protein against surface-induced denaturation in proximity of the air–water interface, while βCD had a destabilizing effect. No remarkable conformational changes of GCSF, or noticeable effect of the CDs, were instead observed at the ice surface. GCSF seemed less stable at low temperature (260 K), which may be attributed to cold-denaturation effects. In this case, CDs did not significantly improve conformational stability. In general, the conformationally altered regions of GCSF seemed not to depend on the presence of excipients that only modulated the extent of destabilization with either a positive or a negative effect.