A program to calculate the isoelectric point of 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.

Structure-based charge calculations for predicting isoelectric point, viscosity, clearance, and profiling antibody therapeutics

The effect of hydrophobicity on antibody aggregation is well understood, and it has been shown that charge calculations can be useful for high-concentration viscosity and pharmacokinetic (PK) clearance predictions. In this work, structure-based charge descriptors are evaluated for their predictive performance on recently published antibody pI, viscosity, and clearance data. From this, we devised four rules for therapeutic antibody profiling which address developability issues arising from hydrophobicity and charged-based solution behavior, PK, and the ability to enrich for those that are approved by the U.S. Food and Drug Administration. Differences in strategy for optimizing the solution behavior of human IgG1 antibodies versus the IgG2 and IgG4 isotypes and the impact of pH alterations in formulation are discussed.

Highland games: A benchmarking exercise in predicting biophysical and drug properties of monoclonal antibodies from amino acid sequences

Biopharmaceutical product and process development do not yet take advantage of predictive computational modeling to nearly the degree seen in industries based on smaller molecules. To assess and advance progress in this area, spirited coopetition (mutually beneficial collaboration between competitors) was successfully used to motivate industrial scientists to develop, share, and compare data and methods which would normally have remained confidential. The first "Highland Games" competition was held in conjunction with the October 2018 Recovery of Biological Products Conference in Ashville, NC, with the goal of benchmarking and assessment of the ability to predict development-related properties of six antibodies from their amino acid sequences alone. Predictions included purification-influencing properties such as isoelectric point and protein A elution pH, and biophysical properties such as stability and viscosity at very high concentrations. Essential contributions were made by a large variety of individuals, including companies which consented to provide antibody amino acid sequences and test materials, volunteers who undertook the preparation and experimental characterization of these materials, and prediction teams who attempted to predict antibody properties from sequence alone. Best practices were identified and shared, and areas in which the community excels at making predictions were identified, as well as areas presenting opportunities for considerable improvement. Predictions of isoelectric point and protein A elution pH were especially good with all-prediction average errors of 0.2 and 1.6 pH unit, respectively, while predictions of some other properties were notably less good. This manuscript presents the events, methods, and results of the competition, and can serve as a tutorial and as a reference for in-house benchmarking by others. Organizations vary in their policies concerning disclosure of methods, but most managements were very cooperative with the Highland Games exercise, and considerable insight into common and best practices is available from the contributed methods. The accumulated data set will serve as a benchmarking tool for further development of in silico prediction tools.

Design of a software for calculating isoelectric point of a polypeptide according to their net charge using the graphical programming language LabVIEW

A software to calculate the net charge and to predict the isoelectric point (pI) of a polypeptide is developed in this work using the graphical programming language LabVIEW. Through this instrument the net charges of the ionizable residues of the polypeptide chains of the proteins are calculated at different pH values, tabulated, pI is predicted and an Excel (-xls) type file is generated. In this work, the experimental values of the pIs (pI) of different proteins are compared with the values of the pIs (pI) calculated graphically, achieving a correlation coefficient (R) of 0.934746 which represents a good reliability for a p < 0.01. In this way the generated program can constitute an instrument applicable in the laboratory, facilitating the calculation to graduate students and junior researchers. © 2017 by The International Union of Biochemistry and Molecular Biology, 46(1):39-46, 2018.

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.

Refolding and characterization of methionine adenosyltransferase from Euglena gracilis

Methionine adenosyltransferase from Euglena gracilis (MATX) is a recently discovered member of the MAT family of proteins that synthesize S-adenosylmethionine. Heterologous overexpression of MATX in Escherichia coli rendered the protein mostly in inclusion bodies under all conditions tested. Therefore, a refolding and purification procedure from these aggregates was developed to characterize the enzyme. Maximal recovery was obtained using inclusion bodies devoid of extraneous proteins by washing under mild urea (2M) and detergent (5%) concentrations. Refolding was achieved in two steps following solubilization in the presence of Mg(2+); chaotrope dilution to <1M and dialysis under reducing conditions. Purified MATX is a homodimer that exhibits Michaelis kinetics with a V(max) of 1.46 μmol/min/mg and K(m) values of approximately 85 and 260 μM for methionine and ATP, respectively. The activity is dependent on Mg(2+) and K(+) ions, but is not stimulated by dimethylsulfoxide. MATX exhibits tripolyphosphatase activity that is stimulated in the presence of S-adenosylmethionine. Far-UV circular dichroism revealed β-sheet and random coil as the main secondary structure elements of the protein. The high level of sequence conservation allowed construction of a structural model that preserved the main features of the MAT family, the major changes involving the N-terminal domain.

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.

Refolding and characterization of rat liver methionine adenosyltransferase from Escherichia coli inclusion bodies

Methionine adenosyltransferase (MAT) catalyzes the synthesis of S-adenosylmethionine, the major methyl donor for transmethylation reactions. Attempts to perform structural studies using rat liver MAT have met with problems because the protein purified from cellular extracts is heterogeneous. Overexpression of the enzyme in Escherichia coli rendered most of the protein as inclusion bodies. These aggregates were purified by specific washes using urea and Triton X-100 and used for refolding. Maximal activity was obtained when chaotropic solubilization included the structural cation Mg(2+), the protein concentration was kept below 0.1 mg/ml, and denaturant removal was carried out in a two-step process, namely, a fast dilution followed by dialysis in the presence of 10 mM DTT or GSH/GSSG redox buffers. Refolding by this procedure generated the oligomeric forms, MAT I and III, which were basically indistinguishable from the purified rat liver forms in secondary structure and catalytic properties.

H2-forming methylenetetrahydromethanopterin dehydrogenase, a novel type of hydrogenase without iron-sulfur clusters in methanogenic archaea

A novel hydrogenase has recently been found in methanogenic archaea. It catalyzes the reversible dehydrogenation of methylenetetrahydromethanopterin (CH2 = H4MPT) to methenyltetrahydromethanopterin (CH identical to H4MPT+) and H2 and was therefore named H2-forming methylenetetrahydromethanopterin dehydrogenase. The hydrogenase, which is composed of only one polypeptide with an apparent molecular mass of 43 kDa, does not mediate the reduction of viologen dyes with either H2 or CH2 = H4MPT. We report here that the purified enzyme from Methanobacterium thermoautotrophicum exhibits the following other unique properties: (a) the colorless protein with a specific activity of 2000 U/mg (Vmax) did not contain iron-sulfur clusters, nickel, or flavins; (b) the activity was not inhibited by carbon monoxide, acetylene, nitrite, cyanide, or azide; (c) the enzyme did not catalyze an isotopic exchange between 3H2 and 1H+; (d) the enzyme catalyzed the reduction of CH identical to H4MPT+ with 3H2 generating [methylene-3H]CH2 = H4MPT; and (e) the primary structure contained at most four conserved cysteines as revealed by a comparison of the DNA-deduced amino acid sequence of the proteins from M. thermoautotrophicum and Methanopyrus kandleri. None of the four cysteines were closely spaced as would be indicative for a (NiFe) hydrogenase or a ferredoxin-type iron-sulfur protein. Properties of the H2-forming methylenetetrahydromethanopterin dehydrogenase from Methanobacterium wolfei are also described indicating that the enzyme from this methanogenic archaeon is very similar to the enzyme from M. thermoautotrophicum with respect both to molecular and catalytic properties.