Strategies for predictive modeling of overloaded oligonucleotide elution profiles in ion-pair chromatography

Due to their potential for gene regulation, oligonucleotides have moved into focus as one of the preferred modalities modulating currently undruggable disease-associated targets. In the course of synthesis and storage of oligonucleotides a significant number of compound-related impurities can be generated. Purification protocols and analytical methods have become crucial for the therapeutic application of any oligonucleotides, be they antisense oligonucleotides (ASOs), small interfering ribonucleic acids (siRNAs) or conjugates. Ion-pair chromatography is currently the standard method for separating and analyzing therapeutic oligonucleotides. Although mathematical modeling can improve the accuracy and efficiency of ion-pair chromatography, its application remains challenging. Simple models may not be suitable to treat advanced single molecules, while complex models are still inefficient for industrial oligonucleotide optimization processes. Therefore, fundamental research to improve the accuracy and simplicity of mathematical models in ion-pair chromatography is still a necessity. In this study, we predict overloaded concentration profiles of oligonucleotides in ion-pair chromatography and compare relatively simple and more advanced predictive models. The experimental system consists of a traditional C18 column using (dibutyl)amine as the ion-pair reagent and acetonitrile as organic modifier. The models were built and tested based on three crude 16-mer oligonucleotides with varying degrees of phosphorothioation, as well as their respective n - 1 and (P = O)1 impurities. In short, the proposed models were suitable to predict the overloaded concentration profiles for different slopes of the organic modifier gradient and column load.

Analyzing proteolytic stability and metabolic hotspots of therapeutic peptides in two rodent pulmonary fluids

Peptides and peptide drug conjugates are emerging modalities to treat pulmonary diseases. Peptides are susceptible to proteolytic cleavage. Expression levels of specific proteases in the lung can be significantly increased in disease state and may lead to exaggerated peptide proteolysis. To support optimization of peptides for inhaled administration, we have recently reported a streamlined high-throughput LC-HRMS protocol to determine enzymatic protease stability of peptides. This method has now been complemented with profiling of peptide metabolic stability in two respiratory fluids, a lung supernatant (lung S9) and a bronchioalveolar lavage fluid (BALF) taken from rats. We have tested a set of 28 peptides with high structural diversity, analyzed the whole data set for formed metabolites, and identified the differences of cleavage pattern in the two test fluids. Comparison of our experimental results and literature-derived cleavage site estimates based on e.g. MEROPS show significant differences for a number of peptides. This indicates the need for an experimental workflow using both protease panels and testing of metabolic stability in lung fluid (BALF) to guide peptide optimization and selection of peptides for inhaled in vivo PK/PD studies in our drug discovery projects.

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP, TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) 1H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

Estimation of Nonlinear Adsorption Isotherms in Gradient Elution RP-LC of Peptides in the Presence of an Adsorbing Additive

ABSTRACT

In electrostatic repulsive interaction chromatography, using a charged surface hybrid sorbent carrying positive charges can improve the peak shape of peptides in reversed-phase liquid chromatography (RP-LC), especially in overloaded conditions, compared with standard C18 sorbents. However, the positive surface charges can interact with anionic additives commonly used in peptide separations, e.g., trifluoroacetic acid (TFA), complicating adsorption isotherm estimation. We investigated how the competition for available adsorption sites between TFA and two peptides influenced the adsorption isotherm in gradient elution. A model accounting for the competition with TFA was compared with a model neglecting TFA adsorption. We found that the two models predicted elution profiles with the same accuracy. We also found that the adsorption isotherms were extremely similar in shape, leading to the conclusion that neglecting the competition with TFA is a valid approximation enabling faster and more robust adsorption isotherm estimation for the studied type of sorbent.

GRAPHICAL ABSTRACT

The importance of ion-pairing in peptide purification by reversed-phase liquid chromatography

The adsorption mechanism for three peptides was studied under overloaded conditions through adsorption isotherm measurements in the presence of an ion-pairing reagent, trifluoroacetic acid (TFA), on an end-capped C₁₈-bonded stationary phase. The overall aim of the study was to obtain a better understanding of how the acetonitrile and the TFA fractions in the eluent affected the overloaded elution profiles and the selectivity between peptides using mechanistic modelling and multivariate design of experiments. When studying the effect of TFA, direct evidence for ion pair formation between a peptide and TFA in acetonitrile-water solutions was provided by fluorine-proton nuclear Overhauser NMR enhancement experiments and the adsorption of TFA on the stationary phase was measured by frontal analysis. The adsorption isotherms for each peptide were then determined by the inverse method at eight TFA concentrations ranging from 2.6mM to 37.3mM (0.02-0.29vol-%) in isocratic elution. The equilibrium between the peptide ion and the peptide-TFA complex was modelled by coupling the mass-balance to reaction kinetics and determining separate adsorption isotherms for the two species. We found that a Langmuir isotherm described the elution profile of peptide-TFA complex well while the peptide ion was described by a bi-Langmuir adsorption isotherm since it exhibited strong secondary interactions. The elution profiles had an unfavorable shape at low TFA concentrations consisting of a spike in their front and a long tailing rear due to the secondary interactions for the peptide ion having very low saturation capacity. The acetonitrile dependence on the adsorption isotherms was studied by determination of adsorption isotherms directly from elution profiles obtained in gradient elution which enabled a broad acetonitrile interval to be studied. Here, it was found that the column saturation capacity was quickly reached at very low acetonitrile fractions and that there were significant variations in adsorption with the molecular weight. Finally, practical implications for method development are discussed based on an experimental design where gradient slope and TFA concentrations are used as factors.

Improved method for the simultaneous determination of proguanil and its metabolites by high-performance liquid chromatography and solid-phase extraction of 100-microl capillary blood samples dried on sampling paper

An improved method is presented for the determination of proguanil, cycloguanil and 4-chlorophenylbiguanide in 100-microl capillary blood samples applied to sampling paper. This method also utilises a solid-phase extraction technique and high-performance liquid chromatography. Different kinds of sampling paper, such as ion-exchange and cellulose sampling paper were tested. The best elution recovery (70-80%) was obtained after treatment of cellulose sampling paper with a quaternary ammonium compound. The limit of determination was 50 nmol/l for cycloguanil and 4-chlorophenylbiguanide and 125 nmol/l for proguanil using 100 microl capillary blood. The stability of the analytes and elution performance from sampling paper was validated at different temperature and storage time. Venous blood and capillary blood concentrations of proguanil and metabolites were found to be similar.

Improved validated assay for the determination of proguanil and its metabolites in plasma, whole blood, and urine using solid-phase extraction and high-performance liquid chromatography

An improved and validated method is presented for the determination of proguanil, cycloguanil, and 4-chlorophenylbiguanide in plasma, whole blood, and urine using solid-phase extraction (SPE) technique and reversed-phase high-performance liquid chromatography (HPLC). The HPLC method uses isocratic elution with acetonitrile:phosphate buffer 0.1 mol/l, pH 2.6 (21.5:78.5 vol/vol) at a flow rate of 1.0 ml/min for the separation. The recovery of proguanil and metabolites ranged from 82% to 104%. The limit of determination was 20 nmol/l for proguanil and its metabolites in plasma and approximately 50 nmol/l for proguanil and metabolites in whole blood. Different stationary phases for HPLC and SPE were tested and the best chromatographic separation from endogenous constituents and other antimalarial drugs was achieved with cyanopropyl stationary phases.