Native Ensemble Research Articles

Molecular Dynamics Simulations Reveal How Competing Protein-Surface Interactions for Glycine, Citrate, and Water Modulate Stability in Antibody Fragment Formulations.

The design of stable formulations remains a major challenge for protein therapeutics, particularly the need to minimize aggregation. Experimental formulation screens are typically based on thermal transition midpoints (Tm), and forced degradation studies at elevated temperatures. Both approaches give limited predictions of long-term storage stability, particularly at low temperatures. Better understanding of the mechanisms of action for formulation of excipients and buffers could lead to improved strategies for formulation design. Here, we identified a complex impact of glycine concentration on the experimentally determined stability of an antibody Fab fragment and then used molecular dynamics simulations to reveal mechanisms that underpin these complex behaviors. Tm values increased monotonically with glycine concentration, but associated ΔSvh measurements revealed more complex changes in the native ensemble dynamics, which reached a maximum at 30 mg/mL. The aggregation kinetics at 65 °C were similar at 0 and 20 mg/mL glycine, but then significantly slower at 50 mg/mL. These complex behaviors indicated changes in the dominant stabilizing mechanisms as the glycine concentration was increased. MD revealed a complex balance of glycine self-interaction, and differentially preferred interactions of glycine with the Fab as it displaced hydration-shell water, and surface-bound water and citrate buffer molecules. As a result, glycine binding to the Fab surface had different effects at different concentrations, and led from preferential interactions at low concentrations to preferential exclusion at higher concentrations. During preferential interaction, glycine displaced water from the Fab hydration shell, and a small number of water and citrate molecules from the Fab surface, which reduced the protein dynamics as measured by root-mean-square fluctuation (RMSF) on the short time scales of MD. By contrast, the native ensemble dynamics increased according to ΔSvh, suggesting increased conformational changes on longer time scales. The aggregation kinetics did not change at low glycine concentrations, and so the opposing dynamics effects either canceled out or were not directly relevant to aggregation. During preferential exclusion at higher glycine concentrations, glycine could only bind to the Fab surface through the displacement of citrate buffer molecules already favorably bound on the Fab surface. Displacement of citrate increased the flexibility (RMSF) of the Fab, as glycine formed fewer bridging hydrogen bonds to the Fab surface. Overall, the slowing of aggregation kinetics coincided with reduced flexibility in the Fab ensemble at the very highest glycine concentrations, as determined by both RMSF and ΔSvh, and occurred at a point where glycine binding displaced neither water nor citrate. These final interactions with the Fab surface were driven by mass action and were the least favorable, leading to a macromolecular crowding effect under the regime of preferential exclusion that stabilized the dynamics of Fab.

Determinants of Unfolding Cooperativity and Binding Are Decoupled in a DNA Binding Domain.

The relative magnitudes of noncovalent stabilization energies or the coupling free energies in folded proteins are anisotropically distributed, uniquely influencing folding and functional behaviors. In this regard, the fructose repressor (FruR) DBD belonging to the LacR repressor family harbors a three-residue insertion─KQY─between the canonical second and third helices. This sequence insertion promotes a strong Tyr-Tyr stacking interaction that is not observed in related homologues. Combining experiments with simulations, we show that the Tyr-Tyr stacking contributes to a decoupled unfolding due to the localization of a large part of the stabilization energy in this specific structural region. This leads to melting temperatures from different probes spanning nearly 10 K, while concomitantly stabilizing a partially structured intermediate state. Disruption of the aromatic stacking interaction via an alanine mutation promotes a molten-globular state whose native ensemble is replete with non-native interactions while displaying enhanced thermodynamic fluctuations and minimal calorimetric cooperativity. Surprisingly, the molten-globular variant of FruR DBD binds to the operator site on DNA with an affinity similar to that of the wild-type but with altered secondary-structure characteristics in the bound state, underscoring the chaperone-like role of DNA through its large negative electrostatic potential. FruR DBD thus appears to be at the verge of disorder as expected of an entropically destabilizing three-residue insertion but is rescued by the aromatic stacking interaction that distinctly dictates the finer details of stability, cooperativity, and binding.

A finely balanced order–disorder equilibrium sculpts the folding–binding landscape of an antibiotic sequestering protein

TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

Molecular and thermodynamic determinants of self-assembly and hetero-oligomerization in the enterobacterial thermo-osmo-regulatory protein H-NS.

Environmentally regulated gene expression is critical for bacterial survival under stress conditions, including extremes in temperature, osmolarity and nutrient availability. Here, we dissect the thermo- and osmo-responsory behavior of the transcriptional repressor H-NS, an archetypal nucleoid-condensing sensory protein, ubiquitous in enterobacteria that infect the mammalian gut. Through experiments and thermodynamic modeling, we show that H-NS exhibits osmolarity, temperature and concentration dependent self-association, with a highly polydisperse native ensemble dominated by monomers, dimers, tetramers and octamers. The relative population of these oligomeric states is determined by an interplay between dimerization and higher-order oligomerization, which in turn drives a competition between weak homo- versus hetero-oligomerization of protein-protein and protein-DNA complexes. A phosphomimetic mutation, Y61E, fully eliminates higher-order self-assembly and preserves only dimerization while weakening DNA binding, highlighting that oligomerization is a prerequisite for strong DNA binding. We further demonstrate the presence of long-distance thermodynamic connectivity between dimerization and oligomerization sites on H-NS which influences the binding of the co-repressor Cnu, and switches the DNA binding mode of the hetero-oligomeric H-NS:Cnu complex. Our work thus uncovers important organizational principles in H-NS including a multi-layered thermodynamic control, and provides a molecular framework broadly applicable to other thermo-osmo sensory proteins that employ similar mechanisms to regulate gene expression.

Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins.

Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.

Functional regulation of an intrinsically disordered protein via a conformationally excited state.

A longstanding goal in the field of intrinsically disordered proteins (IDPs) is to characterize their structural heterogeneity and pinpoint the role of this heterogeneity in IDP function. Here, we use multinuclear chemical exchange saturation (CEST) nuclear magnetic resonance to determine the structure of a thermally accessible globally folded excited state in equilibrium with the intrinsically disordered native ensemble of a bacterial transcriptional regulator CytR. We further provide evidence from double resonance CEST experiments that the excited state, which structurally resembles the DNA-bound form of cytidine repressor (CytR), recognizes DNA by means of a "folding-before-binding" conformational selection pathway. The disorder-to-order regulatory switch in DNA recognition by natively disordered CytR therefore operates through a dynamical variant of the lock-and-key mechanism where the structurally complementary conformation is transiently accessed via thermal fluctuations.

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.

Computational identification of drug-like marine natural products as potential RNA polymerase inhibitors against Nipah virus

Nipah virus (NiV) has been an alarming threat to human populations in southern Asia for more than a decade. It is one of the most deadly viruses in the Mononegavirales order. Despite its high mortality rate and virulence, no chemotherapeutic agent or vaccine is publicly available. Hence, this work was conducted to computationally screen marine natural products database for drug-like potential inhibitors for the viral RNA-dependent RNA polymerase (RdRp). The structural model was subjected to molecular dynamics (MD) simulation to obtain the native ensemble of the protein. The CMNPDB dataset of marine natural products was filtered to retain only compounds following Lipinski's five rules. The molecules were energy minimized and docked into different conformers of the RdRp using AutoDock Vina. The best 35 molecules were rescored by GNINA, a deep learning-based docking software. The resulting nine compounds were evaluated for their pharmacokinetic profiles and medicinal chemistry properties. The best five compounds were subjected to MD simulation for 100 ns, followed by binding free energy estimation via Molecular Mechanics/ Generalized Born Surface Area (MM/GBSA) calculations. The results showed remarkable behavior of five hits as inferred by stable binding pose and orientation to block the exit channel of RNA synthesis products in the RdRp cavity. These hits are promising starting materials for in vitro validation and structural modifications to enhance the pharmacokinetic and medicinal chemistry properties for developing antiviral lead compounds.

Applications of hydrogen deuterium exchange ensemble reweighting to model native protein ensembles, macromolecular interactions and protein drug interactions.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a powerful approach to probe the structural dynamics of a protein in solution. It simultaneously provides time-resolved information on the structural environment of the amide hydrogen of every residue of a protein except Prolines. As such, HDX-MS has been successfully applied to provide insight into the characterization of native state ensembles of proteins, macromolecular interactions, protein-small molecule interactions or protein folding and unfolding studies. However, it remains largely limited to peptide-level resolution information. On the other hand, molecular dynamics (MD) simulations provide a high-resolution view of a protein's structural dynamics by simulating the motion of a protein's atoms over time starting from a high-resolution structure. However, the time scales accessible to MD simulations are limited and can lead to inadequate exploration of a protein's conformational landscape. To bridge this gap, efforts have aimed at developing quantitative and integrative approaches that leverage HDX-MS data to assess and refine computationally generated structural ensembles. HDX ensemble reweighting (HDXer) is one such approach which uses a post hoc maximum entropy approach to adjusts the weights of individual frames of a structural ensemble to conform to HDX-MS data. We had previously applied HDXer to model the native ensemble of the cytoplasmic heme binding protein, PhuS. Here, we have validated the proposed model ensemble by using structure-based drug design approach to successfully target novel, cryptic drug binding sites that would have otherwise not been targetable using HDX-MS or MD simulations alone. Finally, we demonstrate the broader applicability of such HDX-MS guided integrative approach to model macromolecular interactions, missing electron density and protein-small molecule drug interactions

HDX-MS guided ensemble reweighting approach reveals cryptic drug binding sites in the cytoplasmic heme binding protein, PhuS.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a valuable approach to probe protein structure and dynamics. However, interpretation of HDX-MS data is often qualitative and typically limited to peptide-level resolution. Ongoing efforts aim to integrate experimental HDX-MS data with molecular dynamics (MD) approaches, to provide a high-resolution interpretation of the HDX-MS data. One such approach is an HDX-MS based maximum-entropy reweighting approach (HDXer) to reweight computationally generated ensembles to HDX-MS data. Here, HDXer is used in conjunction with enhanced MD (eMD) to characterize the structural dynamics of PhuS. PhuS is a cytoplasmic heme binding protein from P. aeruginosa which shuttles exogenous heme to Heme Oxygenase (HemO) for degradation and is essential for virulence. Through the implementation of a workflow that integrates eMD, HDXer, and dimensional reduction, we had previously modeled the native state ensemble of apo-PhuS revealing a locally diverse apo-PhuS ensemble characterized by large rearrangements of helices surrounding the heme binding. Here, the representative structures from the modeled native ensemble were used as starting points for computer aided drug design (CADD). Implementing the Site Identification by Ligand Competitive Saturation (SILCS) CADD method approximately 1.3 million compounds were screened for each representative structure as well as the published apo PhuS crystal structure, yielding a library of 50 chemically representative candidate compounds. The candidate compounds were then screened using a fluorescence quenching assay. From these assays, several hits were confirmed to bind apo PhuS at nanomolar binding affinities. Most significantly, none of the compounds were derived from CADD using the published apo PhuS crystal structure, demonstrating the utility of HDX-MS guiding ensemble reweighting in revealing cryptic drug binding site in CADD approaches.

Thermodynamic architecture and conformational plasticity of GPCRs

G-protein-coupled receptors (GPCRs) are ubiquitous integral membrane proteins involved in diverse cellular signaling processes. Here, we carry out a large-scale ensemble thermodynamic study of 45 ligand-free GPCRs employing a structure-based statistical mechanical framework. We find that multiple partially structured states co-exist in the GPCR native ensemble, with the TM helices 1, 6 and 7 displaying varied folding status, and shaping the conformational landscape. Strongly coupled residues are anisotropically distributed, accounting for only 13% of the residues, illustrating that a large number of residues are inherently dynamic. Active-state GPCRs are characterized by reduced conformational heterogeneity with altered coupling-patterns distributed throughout the structural scaffold. In silico alanine-scanning mutagenesis reveals that extra- and intra-cellular faces of GPCRs are coupled thermodynamically, highlighting an exquisite structural specialization and the fluid nature of the intramolecular interaction network. The ensemble-based perturbation methodology presented here lays the foundation for understanding allosteric mechanisms and the effects of disease-causing mutations in GCPRs.

Glycan Conformation in the Heavily Glycosylated Protein, CEACAM1.

Glycans attached to glycoproteins can contribute to stability, mediate interactions with other proteins, and initiate signal transduction. Glycan conformation, which is critical to these processes, is highly variable and often depicted as sampling a multitude of conformers. These conformers can be generated by molecular dynamics simulations, and more inclusively by accelerated molecular dynamics, as well as other extended sampling methods. However, experimental assessments of the contribution that various conformers make to a native ensemble are rare. Here, we use long-range pseudo-contact shifts (PCSs) of NMR resonances from an isotopically labeled glycoprotein to identify preferred conformations of its glycans. The N-terminal domain from human Carcinoembryonic Antigen Cell Adhesion Molecule 1, hCEACAM1-Ig1, was used as the model glycoprotein in this study. It has been engineered to include a lanthanide-ion-binding loop that generates PCSs, as well as a homogeneous set of three 13C-labeled N-glycans. Analysis of the PCSs indicates that preferred glycan conformers have extensive contacts with the protein surface. Factors leading to this preference appear to include interactions between N-acetyl methyls of GlcNAc residues and hydrophobic surface pockets on the protein surface.

Quantitative prediction of ensemble dynamics, shapes and contact propensities of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) are highly dynamic systems that play an important role in cell signaling processes and their misfunction often causes human disease. Proper understanding of IDP function not only requires the realistic characterization of their three-dimensional conformational ensembles at atomic-level resolution but also of the time scales of interconversion between their conformational substates. Large sets of experimental data are often used in combination with molecular modeling to restrain or bias models to improve agreement with experiment. It is shown here for the N-terminal transactivation domain of p53 (p53TAD) and Pup, which are two IDPs that fold upon binding to their targets, how the latest advancements in molecular dynamics (MD) simulations methodology produces native conformational ensembles by combining replica exchange with series of microsecond MD simulations. They closely reproduce experimental data at the global conformational ensemble level, in terms of the distribution properties of the radius of gyration tensor, and at the local level, in terms of NMR properties including 15N spin relaxation, without the need for reweighting. Further inspection revealed that 10–20% of the individual MD trajectories display the formation of secondary structures not observed in the experimental NMR data. The IDP ensembles were analyzed by graph theory to identify dominant inter-residue contact clusters and characteristic amino-acid contact propensities. These findings indicate that modern MD force fields with residue-specific backbone potentials can produce highly realistic IDP ensembles sampling a hierarchy of nano- and picosecond time scales providing new insights into their biological function.

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 15N–1H 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.

Quantification of Entropic Excluded Volume Effects Driving Crowding-Induced Collapse and Folding of a Disordered Protein.

We investigate the conformational properties of the intrinsically disordered DNA-binding domain of CytR in the presence of the polymeric crowder polyethylene glycol (PEG). Integrating circular dichroism, nuclear magnetic resonance, and single-molecule Förster resonance energy transfer measurements, we demonstrate that disordered CytR populates a well-folded minor conformation in its native ensemble, while the unfolded ensemble collapses and folds with an increase in crowder density independent of the crowder size. Employing a statistical-mechanical model, the effective reduction in the accessible conformational space of a residue in the unfolded state is estimated to be 10% at 300 mg/mL PEG8000, relative to dilute conditions. The experimentally consistent PEG-temperature phase diagram thus constructed reveals that entropic effects can stabilize disordered CytR by 10 kJ mol-1, driving the equilibrium toward folded conformations under physiological conditions. Our work highlights the malleable conformational landscape of CytR, the presence of a folded conformation in the disordered ensemble, and proposes a scaling relation for quantifying excluded volume effects on protein stability.

Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells

SARS-CoV-2 is a betacoronavirus with a single-stranded, positive-sense, 30-kilobase RNA genome responsible for the ongoing COVID-19 pandemic. Although population average structure models of the genome were recently reported, there is little experimental data on native structural ensembles, and most structures lack functional characterization. Here we report secondary structure heterogeneity of the entire SARS-CoV-2 genome in two lines of infected cells at single nucleotide resolution. Our results reveal alternative RNA conformations across the genome and at the critical frameshifting stimulation element (FSE) that are drastically different from prevailing population average models. Importantly, we find that this structural ensemble promotes frameshifting rates much higher than the canonical minimal FSE and similar to ribosome profiling studies. Our results highlight the value of studying RNA in its full length and cellular context. The genomic structures detailed here lay groundwork for coronavirus RNA biology and will guide the design of SARS-CoV-2 RNA-based therapeutics.

Bivariate bias correction of the regional climate model ensemble over the Adriatic region

AbstractThe main goal of this analysis is to provide an assessment of the fine grid regional climate model (RCM) ensemble over the broader Adriatic region, which is the domain with orographically complex terrain and a developed coastline. The results of this study are a beneficial complement to earlier work on using modified bivariate bias correction methods. This analysis was based on an ensemble of 12 combinations of three RCMs and four global climate models (GCMs) from 1971 to 2004, divided into calibration and validation parts. A comparison of one univariate and three versions of bivariate bias correction methods of summer and winter monthly temperature and precipitation was performed. To examine the effect of various methods, we used different marginal distributions and distributions of dependence (copula). Bivariate bias correction was conducted with parametric and empirical marginal distributions. In the quantile mapping (univariate) and the parametric bivariate method, gamma (precipitation) and normal (temperature) probability distributions were used. We documented the bias and the influence of bias correction methods on precipitation and temperature, including their climate change signals in historical data. The considered methods maintained the spatial distribution of trends from the native ensemble. If the relationship between variables was not distinct and the correlation coefficient was low, bivariate bias correction methods tended to weaken the correlation coefficient from the model. Two experiments on how the bias correction methods influence statistical measures of the considered variables and their relationships showed that the bivariate method with empirical distributions (eeG) performed slightly better than other methods. Although the eeG method performed better than the other methods, none of the considered methods was able to fully preserve the correlation coefficient between variables detected from observations. Considering two copulas (Gaussian and Student's t), we conclude that the effect of different copula types on method performance is small.

Structural-Energetic Basis for Coupling between Equilibrium Fluctuations and Phosphorylation in a Protein Native Ensemble.

The functioning of proteins is intimately tied to their fluctuations in the native ensemble. The structural–energetic features that determine fluctuation amplitudes and hence the shape of the underlying landscape, which in turn determine the magnitude of the functional output, are often confounded by multiple variables. Here, we employ the FF1 domain from human p190A RhoGAP protein as a model system to uncover the molecular basis for phosphorylation of a buried tyrosine, which is crucial to the transcriptional activity associated with transcription factor TFII-I. Combining spectroscopy, calorimetry, statistical–mechanical modeling, molecular simulations, and in vitro phosphorylation assays, we show that the FF1 domain samples a diverse array of conformations in its native ensemble, some of which are phosphorylation-competent. Upon eliminating unfavorable charge–charge interactions through a single charge-reversal (K53E) or charge-neutralizing (K53Q) mutation, we observe proportionately lower phosphorylation extents due to the altered structural coupling, damped equilibrium fluctuations, and a more compact native ensemble. We thus establish a conformational selection mechanism for phosphorylation in the FF1 domain with K53 acting as a “gatekeeper”, modulating the solvent exposure of the buried tyrosine. Our work demonstrates the role of unfavorable charge–charge interactions in governing functional events through the modulation of native ensemble characteristics, a feature that could be prevalent in ordered protein domains.

Multiply Charged Cation Attachment to Facilitate Mass Measurement in Negative-Mode Native Mass Spectrometry.

Native mass spectrometry (MS) is usually conducted in the positive-ion mode; however, in some cases, it is advantageous to use the negative-ion polarity. Challenges associated with native MS using ensemble measurements (i.e., the measurement of many ions at a time as opposed to the measurement of the charge and the mass-to-charge ratio of individual ions) include narrow charge state distributions with the potential for an overlap in neighboring charge states. These issues can either compromise or preclude confident charge state (and hence mass) determination. Charge state determination in challenging instances can be enabled via the attachment of multiply charged ions of opposite polarity. Multiply charged ion attachment facilitates the resolution of charge states and generates mass-to-charge (m/z) information across a broad m/z range. In this work, we demonstrated the attachment of multiply charged cations to anionic complexes generated under native MS conditions. To illustrate the flexibility available in selecting the mass and charge of the reagents, the 15+ and 20+ charge states of horse skeletal muscle apomyoglobin and the 20+ and 30+ charge states of bovine carbonic anhydrase were demonstrated to attach to model complex anions derived from either β-galactosidase or GroEL. The exclusive attachment of reagent ions is observed with no evidence for proton transfer, which is the key for the unambiguous interpretation of the post-ion/ion reaction product ion spectrum. To illustrate the application to mixtures of complex ions, the 10+ charge state of bovine ubiquitin was attached to mixtures of anions generated from the 30S and 50S particles of the Escherichia coli ribosome. Six and five major components were revealed, respectively. In the case of the 50S anion population, it was shown that the attachment of two 30+ cations of carbonic anhydrase revealed the same information as the attachment of six 10+ cations of ubiquitin. In neither case was the intact 50S particle observed. Rather, particles with different combinations of missing components were observed. This work demonstrated the utility of multiply charged cation attachment to facilitate charge state assignments in native MS ensemble measurements of heterogeneous mixtures.

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 (Tm) 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 Tm 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.