An algorithm for the computer calculation of the coefficients of a polynomial that allows determination of isoelectric points of proteins and other macromolecules

Mechanistic modeling of cation exchange chromatography scale-up considering packing inhomogeneities

In process development and characterization, the scale-up of chromatographic steps is a crucial part and brings a number of challenges. Usually, scale-down models are used to represent the process step, and constant column properties are assumed. The scaling is then typically based on the concept of linear scale-up. In this work, a mechanistic model describing an anti-Langmuirian to Langmuirian elution behavior of a polypeptide, calibrated with a pre-packed 1 ml column, is applied to demonstrate the scalability to larger column volumes up to 28.2 ml. Using individual column parameters for each column size, scaling to similar eluting salt concentrations, peak heights, and shapes is experimentally demonstrated by considering the model's relationship between the normalized gradient slope and the eluting salt concentration. Further scale-up simulations show improved model predictions when radial inhomogeneities in packing quality are considered.

Anti-Langmuir elution behavior of a bispecific monoclonal antibody in cation exchange chromatography: Mechanistic modeling using a pH-dependent Self-Association Steric Mass Action isotherm

The objective of this scientific work was to model and simulate the complex anti-Langmuir elution behavior of a bispecific monoclonal antibody (bsAb) under high loading conditions on the strong cation exchange resin POROS™ XS. The bsAb exhibited anti-Langmuirian elution behavior as a consequence of self-association expressed both in uncommon retentions and peak shapes highly atypical for antibodies. The widely applied Steric Mass Action (SMA) model was unsuitable here because it can only describe Langmuirian elution behavior and is not able to describe protein-protein interactions in the form of self-association. For this reason, a Self-Association SMA (SAS-SMA) model was applied, which was extended by two activity coefficients for the salt and protein in solution. This model is able to describe protein-protein interactions in the form of self-dimerization and thus can describe anti-Langmuir elution behavior. Linear gradient elution (LGE) experiments were carried out to obtain a broad dataset ranging from pH 4.5 to 7.3 and from 50 to 375 mmol/L Na⁺ for model parameter determination. High loading LGE experiments were conducted with an increasing load from 0.5 up to 75.0 mg_bsAb/mL_resin. Thereby, pH-dependent empirical correlations for the activity coefficient of the solute protein, for the equilibrium constant of the self-dimerization process and for the shielding factor could be set up and ultimately incorporated into the SAS-SMA model. This pH-dependent SAS-SMA model was thus able to simulate anti-Langmuir behavior over extended ranges of pH, counterion concentration, and column loading. The model was confirmed by experimental verification of simulated linear pH gradient elutions up to a load of 75.0 mg_bsAb/mL_resin.

Mechanistic modeling and simulation of a complex low and high loading elution behavior of a polypeptide in cation exchange chromatography

The mechanistic modeling of preparative liquid chromatography is still a challenging task. Non-ideal thermodynamic conditions may require activity coefficients for the mechanistic description of preparative chromatography. In this work, a chromatographic cation exchange step with a polypeptide having a complex elution behavior in low and high loading situations is modeled. Model calibration in the linear range of the isotherm is done by applying counterion-induced linear gradient elution experiments between pH 3.3 and pH 4.3. Inverse fitting with column loads up to 25 mg/mL_CV is performed for parameter estimation in the non-linear range. The polypeptide elution peak shows an anti-Langmuirian behavior with fronting under low loading conditions and a switch to a Langmuirian behavior with increasing load. This unusual elution behavior could be described with an extended version of the sigmoidal Self-Association isotherm, including two activity coefficients for the polypeptide and counterion in solution. The activity coefficient of the solute polypeptide shows a strong influence on the model parameters and is crucial in the linear and non-linear range of the isotherm. The modeling procedure results in a unique and robust model parameter set that is sufficient to describe the complex elution behavior and allows modeling over the full isotherm range. This article is protected by copyright. All rights reserved.

Solvent modulated linear pH gradient elution for the purification of conventional and bispecific antibodies: Modeling and application

Classical ion-exchange chromatography using a linear salt gradient to elute the adsorbed protein at fixed pH is the most common method to separate product-related impurities during downstream processing of biopharmaceuticals. Linear pH gradient elution provides a useful alternative by separating proteins in a linear pH gradient at fixed salt concentration. Although linear pH gradient elution provides excellent selectivity, it is rarely encountered in industrial purification processes. Here, a stoichiometric displacement model is used to characterize pH gradient elution based on simple linear gradient elution experiments. Protein retention behavior is described with respect to the pH dependencies of the characteristic binding charge and the equilibrium constant of the ion exchange reaction. Furthermore, the influence of solvent composition using PEG as a mobile phase modifier is investigated. Validity and applicability of the model are demonstrated for the purification of a conventional monoclonal antibody from soluble aggregates and for a novel bispecific antibody format containing a unique product-related impurity profile. pH step elution protocols are derived from model calculations without further optimization experiments necessary.

Isoelectric point, electric charge, and nomenclature of the acid-base residues of proteins

The main object of this work is to present the pedagogical usefulness of the theoretical methods, developed in this laboratory, for the determination of the isoelectric point (pI) and the net electric charge of proteins together with some comments on the naming of the acid-base residues of proteins.

Improved application of the oscillating method for the isoelectric point determination of protein: Potential connection with protein data banks

The oscillating method (OM) for the theoretical determination of the pI values, one by one, of proteins and other macromolecules has been previously published [Sillero and Maldonado, Comput. Biol. Med 36 (2006) 157-166]. An improved application of the method, here named as improved oscillating method (IOM), allows the pI determination of group of proteins. This characteristic may be useful to explore the pI value and electric charge of family of enzymes. As an example the pI values of 1630 enzymes collected in a Swiss-Prot data bank (www.expasy.org), as belonging to the enzymes ligases (EC 6. 2. 1. *) is presented. The method also permits the determination of the pI value of any group of proteins stored in data banks provided that they can be supplied to the program in a FASTA format. Free access to IOM can be reached at http://www.bq.uam.es/otros/pical3.zip.

Isoelectric point determination of proteins and other macromolecules: oscillating method

A program written in Visual Basic has been developed to calculate the isoelectric point of proteins and other macromolecules bearing acid-basic residues. The pI value can be theoretically calculated with the precision required. The computer automatically supplies a representation of the charge of the protein versus pH values. The corresponding values can also be obtained, on command, in the form of table.

Identification, cloning, and characterization of a major cat flea salivary allergen (Cte f 1)

An 18 kDa protein isolated from saliva of the cat flea, Ctenocephalides felis, elicits a positive intradermal skin test (IDST) in 100 and 80% of experimental and clinical flea allergic dogs, respectively. Using solid-phase enzyme-linked immuno assay (ELISA), this protein detected IgE in 100 and 80% of experimental and clinical flea allergic dogs, respectively. A cDNA (pFSI) encoding a full-length Cte f 1 protein was isolated from a C. felis salivary gland cDNA library, using a combination of PCR and hybridization screening. This cDNA is 658 bp in length, and contains an open reading frame of 528 bp. The open reading frame encodes a protein of 176 amino acids, consisting of an 18 amino acid signal sequence and a 158 amino acid mature protein. The calculated molecular weight and pI of the mature protein are 18106 Da and 9.3, respectively. The protein, named Cte f 1, is the first novel major allergen described for canine flea allergy. Recombinant Cte f 1 (rCte f 1) was expressed in Escherichia coli, Pichia pastoris and baculovirus infected Trichoplusia ni cells. Approximately, 90% of the rCte f 1 expressed in E. coli accumulated in insoluble inclusion bodies, which could be refolded to a soluble mixture of disulfide isomers with partial IgE binding activity. Small quantities of an apparently correctly refolded form of rCte f 1, which had IgE binding activity equal to the native antigen, was isolated from the soluble fraction of E. coli cells. However, P. pastoris and baculovirus infected insect cells expressed and secreted a fully processed, correctly refolded and fully active form of rCte f 1. Mass spectrometry analysis of the active forms of rCte f 1confirmed that eight intact disulfide bonds were present, matching the number observed in the native allergen. The relative ability of rCte f 1 to bind IgE in the serum of flea allergic animals, produced in these three expression systems, matched that of the native allergen. Competition ELISA demonstrated that approximately 90% of the specific IgE binding to native Cte f 1 could be blocked by the different forms of rCte f 1.

Matrix metalloproteinase 13 (collagenase 3) in human rheumatoid synovium

OBJECTIVE

To show the eventual presence and extent of production of matrix metalloproteinase 13 (MMP-13, or collagenase 3) in rheumatoid synovial tissue samples and extracts, and to assess the inhibition characteristics of recombinant MMP-13.

METHODS

Immunohistochemical avidin-biotin-peroxidase complex staining/morphometry was used to analyze MMP-13-positive cells in situ. Neutral salt extraction of synovial tissue, electrophoresis of the extract in different buffer systems, and Western blotting were also used. The inhibitory properties of doxycycline, clodronate, pamidronate, and D-penicillamine for recombinant enzyme were determined with a soluble type II collagen assay.

RESULTS

MMP-13 was detected in fibroblast- and macrophage-like mononuclear cells in the synovial lining and stroma and in vascular endothelial cells. The overall expression of MMP-13 in these cells in the synovial stroma was high in rheumatoid arthritis (86 +/- 12%) compared with osteoarthritis (17 +/- 5%) patient samples (P = 0.0027). In a high-pH native electrophoresis gel, immunoreactivity to anti-MMP-1 and anti-MMP-13 were clearly separated, with anti-MMP-13-immunoreactive material migrating faster than anti-MMP-1-immunoreactive material. Finally, in contrast to MMP-1 and MMP-8, MMP-13 was found to be relatively resistant to the inhibitory effects of doxycycline and clodronate in vitro.

CONCLUSION

Due to its localization in synovial tissue, its substrate profile, increased expression, and relative resistance to known MMP inhibitors, MMP-13 is suggested to play a major role in the pathogenesis of tissue destruction in rheumatoid arthritis.

A program to calculate the isoelectric point of macromolecules

A program containing 260 sentences, written in BASIC and adapted to be run in personal computers, has been developed to calculate isoelectric point values of macromolecules (proteins, nucleic acids, polysaccharides, etc.). This implies the calculation of the coefficients and roots of a particular kind of polynomial. With common personal computers and using a compiled version of the program, the time required to calculate pI values of macromolecules containing for example, 9, 12, and 18 acid-base residues with different pK values was around 10 s, 2 min, and 2 h respectively.

Theoretical isoelectric points and electric charges of mutated human hemoglobin subunits

Isoelectric point values have been determined from a theoretical standpoint for alpha, beta, gamma and delta human abnormal subunits bearing elimination or addition of 1-2 different aminoacids with acid-base residues. The charge distribution, in relation to pH values, of the four types of subunits are also reported. These values may be of interest both to predict the electrophoretic behaviour and to characterize human abnormal hemoglobins.