Nonspecific Binding of Adenosine Tripolyphosphate and Tripolyphosphate Modulates the Phase Behavior of Lysozyme.

Adenosine tripolyphosphate (ATP) is a small polyvalent anion that has recently been shown to interact with proteins and have a major impact on assembly processes involved in biomolecular condensate formation and protein aggregation. However, the nature of non-specific protein-ATP interactions and their effects on protein solubility are largely unknown. Here, the binding of ATP to the globular model protein is characterized in detail using X-ray crystallography and nuclear magnetic resonance (NMR). Using NMR, we identified six ATP binding sites on the lysozyme surface, with one known high-affinity nucleic acid binding site and five non-specific previously unknown sites with millimolar affinities that also bind tripolyphosphate (TPP). ATP binding occurs primarily through the polyphosphate moiety, which was confirmed by the X-ray structure of the lysozyme-ATP complex. Importantly, ATP binds preferentially to arginine over lysine in non-specific binding sites. ATP and TPP have similar effects on solution-phase protein-protein interactions. At low salt concentrations, ion binding to lysozyme causes precipitation, while at higher salt concentrations, redissolution occurs. The addition of an equimolar concentration of magnesium to ATP does not alter ATP binding affinities but prevents lysozyme precipitation. These findings have important implications for both protein crystallization and cell biology. Crystallization occurs readily in ATP solutions outside the well-established crystallization window. In the context of cell biology, the findings suggest that ATP binds non-specifically to folded proteins in physiological conditions. Based on the nature of the binding sites identified by NMR, we propose several mechanisms for how ATP binding can prevent the aggregation of natively folded proteins.

The Relationship between Protein-Protein Interactions and Liquid-Liquid Phase Separation for Monoclonal Antibodies.

Being able to predict and control concentrated solution properties for solutions of monoclonal antibodies (mAbs) is critical for developing therapeutic formulations. At higher protein concentrations, undesirable solution properties include high viscosities, opalescence, particle formation, and precipitation. The overall aim of this work is to understand the relationship between commonly measured dilute solution parameters, the reduced osmotic second virial coefficient b22 and the diffusion interaction parameter kD and liquid-liquid phase separation, which occurs at higher protein concentrations. For globular proteins such as lysozyme or γB-crystallin, the location of the liquid-liquid coexistence curve is controlled by the net protein-protein interaction, which is related to b22. Because many mAbs undergo reversible self-association due to forming highly directional interactions, it is not known if b22 can be used as a reliable predictor for LLPS since increasing the anisotropy in the interaction potential causes phase separation to occur at much stonger net protein-protein attractions or lower values of b22. Here, we map the coexistence curves for three mAbs, referred to as COE-01, COE-07, and COE-19, in terms of b22 and kD values. The measurements are carried out at a low salt condition near the pI, where protein-protein interactions are expected to be anisotropic due to the presence of electrostatic attractions, and under salting-out conditions at high ammonium sulfate concentrations, which is expected to reduce the anisotropy by screening electrostatic interactions. We also show that deviations from a linear correlation between b22 and kD can be used as an indicator of reversible self-association. Each of the mAbs under salting-out conditions follows the correlation supporting the hypothesis that protein-protein interactions are nonspecific, while deviations from the correlation occur for COE-01 and COE-19 under low salt conditions indicating the mAbs undergo reversible self-association. For five out of the six conditions, the onset of phase separation, as reflected by the reduced virial coefficient at the critical point b22c occurs in a small window -1.6 > b22c > -2.3, which is similar to what has been observed for lysozyme and for bovine γB-crystallin. Under low salt conditions, b22c ≈ -5.1 for COE-19, which we previously showed to self-associate into small oligomers. Our findings suggest that under conditions where mAb interactions are weakly anisotropic, such as occur at high salt conditions, phase separation will begin to occur in a small window of b22. Deviations from the window can occur when mAbs undergo reversible self-association, although this is not always the case and likely depends upon whether or not highly directional interactions are passivated in the oligomer formation. We expect fitting LLPS measurements to simplified interaction models for mAbs will provide additional insight into the nature of the protein-protein interactions and guide their development for calculating concentrated solution properties.

Enhanced Thermal Stability and Reduced Aggregation in an Antibody Fab Fragment at Elevated Concentrations.

The aggregation of protein therapeutics such as antibodies remains a major challenge in the biopharmaceutical industry. The present study aimed to characterize the impact of the protein concentration on the mechanisms and potential pathways for aggregation, using the antibody Fab fragment A33 as the model protein. Aggregation kinetics were determined for 0.05 to 100 mg/mL Fab A33, at 65 °C. A surprising trend was observed whereby increasing the concentration decreased the relative aggregation rate, ln(v) (% day-1), from 8.5 at 0.05 mg/mL to 4.4 at 100 mg/mL. The absolute aggregation rate (mol L-1 h-1) increased with the concentration following a rate order of approximately 1 up to a concentration of 25 mg/mL. Above this concentration, there was a transition to an apparently negative rate order of -1.1 up to 100 mg/mL. Several potential mechanisms were examined as possible explanations. A greater apparent conformational stability at 100 mg/mL was observed from an increase in the thermal transition midpoint (Tm) by 7-9 °C, relative to those at 1-4 mg/mL. The associated change in unfolding entropy (△Svh) also increased by 14-18% at 25-100 mg/mL, relative to those at 1-4 mg/mL, indicating reduced conformational flexibility in the native ensemble. Addition of Tween or the crowding agents Ficoll and dextran, showed that neither surface adsorption, diffusion limitations nor simple volume crowding affected the aggregation rate. Fitting of kinetic data to a wide range of mechanistic models implied a reversible two-state conformational switch mechanism from aggregation-prone monomers (N*) into non-aggregating native forms (N) at higher concentrations. kD measurements from DLS data also suggested a weak self-attraction while remaining colloidally stable, consistent with macromolecular self-crowding within weakly associated reversible oligomers. Such a model is also consistent with compaction of the native ensemble observed through changes in Tm and △Svh.

Chemometrics in Protein Formulation: Stability Governed by Repulsion and Protein Unfolding.

Therapeutic proteins can be challenging to develop due to their complexity and the requirement of an acceptable formulation to ensure patient safety and efficacy. To date, there is no universal formulation development strategy that can identify optimal formulation conditions for all types of proteins in a fast and reliable manner. In this work, high-throughput characterization, employing a toolbox of five techniques, was performed on 14 structurally different proteins formulated in 6 different buffer conditions and in the presence of 4 different excipients. Multivariate data analysis and chemometrics were used to analyze the data in an unbiased way. First, observed changes in stability were primarily determined by the individual protein. Second, pH and ionic strength are the two most important factors determining the physical stability of proteins, where there exists a significant statistical interaction between protein and pH/ionic strength. Additionally, we developed prediction methods by partial least-squares regression. Colloidal stability indicators are important for prediction of real-time stability, while conformational stability indicators are important for prediction of stability under accelerated stress conditions at 40 °C. In order to predict real-time storage stability, protein-protein repulsion and the initial monomer fraction are the most important properties to monitor.

Cluster Percolation Causes Shear Thinning Behavior in Concentrated Solutions of Monoclonal Antibodies.

High-concentration (>100 g/L) solutions of monoclonal antibodies (mAbs) are typically characterized by anomalously large solution viscosity and shear thinning behavior for strain rates ≥103 s-1. Here, the link between protein-protein interactions (PPIs) and the rheology of concentrated solutions of COE-03 and COE-19 mAbs is studied by means of static and dynamic light scattering and microfluidic rheometry. By comparing the experimental data with predictions based on the Baxter sticky hard-sphere model, we surprisingly find a connection between the observed shear thinning and the predicted percolation threshold. The longest shear relaxation time of mAbs was much larger than that of model sticky hard spheres within the same region of the phase diagram, which is attributed to the anisotropy of the mAb PPIs. Our results suggest that not only the strength but also the patchiness of short-range attractive PPIs should be explicitly accounted for by theoretical approaches aimed at predicting the shear rate-dependent viscosity of dense mAb solutions.

Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization.

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Arginine to Lysine Mutations Increase the Aggregation Stability of a Single-Chain Variable Fragment through Unfolded-State Interactions.

Increased protein solubility is known to correlate with an increase in the proportion of lysine over arginine residues. Previous work has shown that the aggregation propensity of a single-chain variable fragment (scFv) does not correlate with its conformational stability or native-state protein-protein interactions. Here, we test the hypothesis that aggregation is driven by the colloidal stability of partially unfolded states, studying the behavior of scFv mutants harboring single or multiple site-specific arginine to lysine mutations in denaturing buffers. In 6 M guanidine hydrochloride (GdmCl) or 8 M urea, repulsive protein-protein interactions were measured for the wild-type and lysine-enriched (4RK) scFvs reflecting weakened short-range attractions and increased excluded volume. In contrast to the arginine-enriched mutant (7KR) scFv exhibited strong reversible association. In 3 M GdmCl, the minimum concentration at which the scFvs were unfolded, the hydrodynamic radius of 4RK remained constant but increased for the wild type and especially for 7KR. Studies of single-point arginine to lysine scFv mutants indicated that the observed aggregation propensity of arginine under denaturing conditions was nonspecific. Interestingly, one such swap generated a scFv with especially low aggregation rates under low/high ionic strengths and denaturing buffers; molecular modeling identified hydrogen bonding between the arginine side chain and main chain peptide groups, stabilizing the structure. The arginine/lysine ratio is not routinely considered in biopharmaceutical scaffold design or current amyloid prediction methods. This work therefore suggests a simple method for increasing the stability of a biopharmaceutical protein against aggregation.

The Macromolecular Femtosecond Crystallography Instrument at the Linac Coherent Light Source.

The Macromolecular Femtosecond Crystallography (MFX) instrument at the Linac Coherent Light Source (LCLS) is the seventh and newest instrument at the world's first hard X-ray free-electron laser. It was designed with a primary focus on structural biology, employing the ultrafast pulses of X-rays from LCLS at atmospheric conditions to overcome radiation damage limitations in biological measurements. It is also capable of performing various time-resolved measurements. The MFX design consists of a versatile base system capable of supporting multiple methods, techniques and experimental endstations. The primary techniques supported are forward scattering and crystallography, with capabilities for various spectroscopic methods and time-resolved measurements. The location of the MFX instrument allows for utilization of multiplexing methods, increasing user access to LCLS by running multiple experiments simultaneously.

Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies.

The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.

Protein-Sol: a web tool for predicting protein solubility from sequence.

MOTIVATION: Protein solubility is an important property in industrial and therapeutic applications. Prediction is a challenge, despite a growing understanding of the relevant physicochemical properties. RESULTS: Protein-Sol is a web server for predicting protein solubility. Using available data for Escherichia coli protein solubility in a cell-free expression system, 35 sequence-based properties are calculated. Feature weights are determined from separation of low and high solubility subsets. The model returns a predicted solubility and an indication of the features which deviate most from average values. Two other properties are profiled in windowed calculation along the sequence: fold propensity, and net segment charge. The utility of these additional features is demonstrated with the example of thioredoxin. AVAILABILITY AND IMPLEMENTATION: The Protein-Sol webserver is available at http://protein-sol.manchester.ac.uk. CONTACT: jim.warwicker@manchester.ac.uk.

Proofreading of substrate structure by the Twin-Arginine Translocase is highly dependent on substrate conformational flexibility but surprisingly tolerant of surface charge and hydrophobicity changes.

The Tat system transports folded proteins across the bacterial plasma membrane, and in Escherichia coli preferentially transports correctly-folded proteins. Little is known of the mechanism by which Tat proofreads a substrate's conformational state, and in this study we have addressed this question using a heterologous single-chain variable fragment (scFv) with a defined structure. We introduced mutations to surface residues while leaving the folded structure intact, and also tested the importance of conformational flexibility. We show that while the scFv is stably folded and active in the reduced form, formation of the 2 intra-domain disulphide bonds enhances Tat-dependent export 10-fold, indicating Tat senses the conformational flexibility and preferentially exports the more rigid structure. We further show that a 26-residue unstructured tail at the C-terminus blocks export, suggesting that even this short sequence can be sensed by the proofreading system. In contrast, the Tat system can tolerate significant changes in charge or hydrophobicity on the scFv surface; substitution of uncharged residues by up to 3 Lys-Glu pairs has little effect, as has the introduction of up to 5 Lys or Glu residues in a confined domain, or the introduction of a patch of 4 to 6 Leu residues in a hydrophilic region. We propose that the proofreading system has evolved to sense conformational flexibility and detect even very transiently-exposed internal regions, or the presence of unfolded peptide sections. In contrast, it tolerates major changes in surface charge or hydrophobicity.

The X-ray Correlation Spectroscopy instrument at the Linac Coherent Light Source.

The X-ray Correlation Spectroscopy instrument is dedicated to the study of dynamics in condensed matter systems using the unique coherence properties of free-electron lasers. It covers a photon energy range of 4-25 keV. The intrinsic temporal characteristics of the Linac Coherent Light Source, in particular the 120 Hz repetition rate, allow for the investigation of slow dynamics (milliseconds) by means of X-ray photon correlation spectroscopy. Double-pulse schemes could probe dynamics on the picosecond timescale. A description of the instrument capabilities and recent achievements is presented.

The Coherent X-ray Imaging instrument at the Linac Coherent Light Source.

The Coherent X-ray Imaging (CXI) instrument specializes in hard X-ray, in-vacuum, high power density experiments in all areas of science. Two main sample chambers, one containing a 100 nm focus and one a 1 µm focus, are available, each with multiple diagnostics, sample injection, pump-probe and detector capabilities. The flexibility of CXI has enabled it to host a diverse range of experiments, from biological to extreme matter.

Lysine and arginine content of proteins: computational analysis suggests a new tool for solubility design.

Prediction and engineering of protein solubility is an important but imprecise area. While some features are routinely used, such as the avoidance of extensive non-polar surface area, scope remains for benchmarking of sequence and structural features with experimental data. We study properties in the context of experimental solubilities, protein gene expression levels, and families of abundant proteins (serum albumin and myoglobin) and their less abundant paralogues. A common feature that emerges for proteins with elevated solubility and at higher expression and abundance levels is an increased ratio of lysine content to arginine content. We suggest that the same properties of arginine that give rise to its recorded propensity for specific interaction surfaces also lead to favorable interactions at nonspecific contacts, and thus lysine is favored for proteins at relatively high concentration. A survey of protein therapeutics shows that a significant subset possesses a relatively low lysine to arginine ratio, and therefore may not be favored for high protein concentration. We conclude that modulation of lysine and arginine content could prove a useful and relatively simple addition to the toolkit available for engineering protein solubility in biotechnological applications.

Quantifying the impact of an invasive Hornet on Bombus terrestris Colonies.

The invasive hornet Vespa velutina nigrithorax is considered a proliferating threat to pollinators in Europe and Asia. While the impact of this species on managed honey bees is well-documented, effects upon other pollinator populations remain poorly understood. Nonetheless, dietary analyses indicate that the hornets consume a diversity of prey, fuelling concerns for at-risk taxa. Here, we quantify the impact of V. velutina upon standardised commercially-reared colonies of the European bumblebee, Bombus terrestris terrestris. Using a landscape-scale experimental design, we deploy colonies across a gradient of local V. velutina densities, utilising automated tracking to non-invasively observe bee and hornet behaviour, and quantify subsequent effects upon colony outcomes. Our results demonstrate that hornets frequently hunt at B. terrestris colonies, being preferentially attracted to those with high foraging traffic, and engaging in repeated-yet entirely unsuccessful-predation attempts at nest entrances. Notably however, we show that B. terrestris colony weights are negatively associated with local V. velutina densities, indicating potential indirect effects upon colony growth. Taken together, these findings provide the first empirical insight into impacts on bumblebees at the colony level, and inform future mitigation efforts for wild and managed pollinators.

Towards an improved prediction of concentrated antibody solution viscosity using the Huggins coefficient.

The viscosity of a monoclonal antibody solution must be monitored and controlled as it can adversely affect product processing, packaging and administration. Engineering low viscosity mAb formulations is challenging as prohibitive amounts of material are required for concentrated solution analysis, and it is difficult to predict viscosity from parameters obtained through low-volume, high-throughput measurements such as the interaction parameter, kD, and the second osmotic virial coefficient, B22. As a measure encompassing the effect of intermolecular interactions on dilute solution viscosity, the Huggins coefficient, kh, is a promising candidate as a parameter measureable at low concentrations, but indicative of concentrated solution viscosity. In this study, a differential viscometry technique is developed to measure the intrinsic viscosity, [η], and the Huggins coefficient, kh, of protein solutions. To understand the effect of colloidal protein-protein interactions on the viscosity of concentrated protein formulations, the viscometric parameters are compared to kD and B22 of two mAbs, tuning the contributions of repulsive and attractive forces to the net protein-protein interaction by adjusting solution pH and ionic strength. We find a strong correlation between the concentrated protein solution viscosity and the kh but this was not observed for the kD or the b22, which have been previously used as indicators of high concentration viscosity. Trends observed in [η] and kh values as a function of pH and ionic strength are rationalised in terms of protein-protein interactions.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics.

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

Flavocytochrome P450 BM3 mutant W1046A is a NADH-dependent fatty acid hydroxylase: implications for the mechanism of electron transfer in the P450 BM3 dimer.

Bacillus megaterium P450 BM3 (BM3) is a P450/P450 reductase fusion enzyme, where the dimer is considered the active form in NADPH-dependent fatty acid hydroxylation. The BM3 W1046A mutant was generated, removing an aromatic "shield" from its FAD isoalloxazine ring. W1046A BM3 is a catalytically active NADH-dependent lauric acid hydroxylase, with product formation slightly superior to the NADPH-driven enzyme. The W1046A BM3 K(m) for NADH is 20-fold lower than wild-type BM3, and catalytic efficiency of W1046A BM3 with NADH and NADPH are similar in lauric acid oxidation. Wild-type BM3 also catalyzes NADH-dependent lauric acid hydroxylation, but less efficiently than W1046A BM3. A hypothesis that W1046A BM3 is inactive [15] helped underpin a model of electron transfer from FAD in one BM3 monomer to FMN in the other in order to drive fatty acid hydroxylation in native BM3. Our data showing W1046A BM3 is a functional fatty acid hydroxylase are consistent instead with a BM3 catalytic model involving electron transfer within a reductase monomer, and from FMN of one monomer to heme of the other [12]. W1046A BM3 is an efficient NADH-utilizing fatty acid hydroxylase with potential biotechnological applications.

Production of recombinant peptides as fusions with SUMO.

Recombinant production of non-native peptides requires using protein fusion technology to prevent peptide degradation by host-cell proteases. In this work, we have used SUMO protein as a fusion partner for the production of difficult-to-express, antimicrobial, self-assembling and amyloidogenic peptides using Escherichia coli. SUMO-peptide fusions were expressed as intracellular products by utilizing pET based expression vectors constructed by Life Sensors Inc., USA. Histidine tagged SUMO-peptide fusions were purified using Ni-NTA affinity chromatography. Complete (100%) cleavage of the SUMO-peptide fusion was achieved using SUMO protease-1. Our findings demonstrate that SUMO fusion technology is a promising alternative for production of peptides in E. coli. The key advantage of this technology is that the enzymatic activity of SUMO protease-1 is specific and efficient leading to inexpensive costs for cleaving the peptide fusion when compared with other fusion systems.