A model of the molten globule state from molecular dynamics simulations

Improvement the thermostability and specific activity of acidic xylanase PjxA from Penicillium janthinellum via rigid flexible sites

Acidic xylanase PjxA from Penicillium janthinellum MA21601, with good eosinophilic and enzymatic activity, is an excellent candidate for xylan degradation to achieve effective utilization of biomass materials. However, the low thermal stability of PjxA has become a major bottleneck in its application. In this study, the flexible sites of PjxA were identified and rigidified through computational simulations of structure and sequence analysis combined with folding free energy calculations. Finally, a combined mutase PjxA-DS was constructed by rational integration of the two single mutants S82N and D45N. Compared to PjxA, PjxA-DS showed a 115.11-fold longer half-life at 50 °C and a 2.02-fold higher specific enzyme activity. Computer simulation analysis showed that S82N and D45N acted synergistically to improve the thermostability of PjxA. The stabilization of the N-terminus and the active center of PjxA, the increase in surface positive charge and hydrophilicity are the main reasons for the improved thermostability and catalytic activity of PjxA. Rigidification of the flexible site is an effective method for improving the thermostability of enzymes, S82N and D45N can be used as effective targets for the thermostability engineering modification of GH11 acidic xylanase.

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Insights into Fluctuations of Structure of Proteins: Significance of Intermediary States in Regulating Biological Functions

Proteins are indispensable to cellular communication and metabolism. The structure on which cells and tissues are developed is deciphered from proteins. To perform functions, proteins fold into a three-dimensional structural design, which is specific and fundamentally determined by their characteristic sequence of amino acids. Few of them have structural versatility, allowing them to adapt their shape to the task at hand. The intermediate states appear momentarily, while protein folds from denatured (D) ⇔ native (N), which plays significant roles in cellular functions. Prolific effort needs to be taken in characterizing these intermediate species if detected during the folding process. Protein folds into its native structure through definite pathways, which involve a limited number of transitory intermediates. Intermediates may be essential in protein folding pathways and assembly in some cases, as well as misfolding and aggregation folding pathways. These intermediate states help to understand the machinery of proper folding in proteins. In this review article, we highlight the various intermediate states observed and characterized so far under in vitro conditions. Moreover, the role and significance of intermediates in regulating the biological function of cells are discussed clearly.

ATP Can Efficiently Stabilize Protein through a Unique Mechanism

Recent experiments suggested that ATP can effectively stabilize protein structure and inhibit protein aggregation when its concentration is less than 10 mM, which is significantly lower than cosolvent concentrations required in conventional mechanisms. The ultrahigh efficiency of ATP suggests a unique mechanism that is fundamentally different from previous models of cosolvents. In this work, we used molecular dynamics simulation and experiments to study the interactions of ATPs with three proteins: lysozyme, ubiquitin, and malate dehydrogenase. ATP tends to bind to the surface regions with high flexibility and high degree of hydration. These regions are also vulnerable to thermal perturbations. The bound ATPs further assemble into ATP clusters mediated by Mg2+ and Na+ ions. More interestingly, in Mg2+-free ATP solution, Na+ at higher concentration (150 mM under physiological conditions) can similarly mediate the formation of the ATP cluster on protein. The ATP cluster can effectively reduce the fluctuations of the vulnerable region and thus stabilize the protein against thermal perturbations. Both ATP binding and the considerable improvement of thermal stability of ATP-bound protein were verified by experiments.

Size-Selective VAILase Proteolysis Provides Dynamic Insights into Protein Structures

Monitoring the dynamic alterations of protein structures within an aqueous solution remains enormously challenging. In this study, we describe a size-selective VAILase proteolysis (SVP)-mass spectrometry (MS) strategy to probe the protein structure changes without strict control of the proteolysis kinetics. The unique conformation selectivity of SVP depends on the uniform nano-sized entrance pores of the VAILase hexameric cage as well as the six inherent molecular rulers in the VAILase-substrate recognition and cleavage. The dynamic insights into subtle conformation alterations of both myoglobin unfolding transition and Aurora kinase A-inhibitor binding are successfully captured using the SVP strategy, which matches well with the results in the molecular dynamics simulation. Our work provides a new paradigm of size-selective native proteolysis for exploring the aqueous protein structure-function relationships.

Polyethylene Glycol Crowder's Effect on Enzyme Aggregation, Thermal Stability, and Residual Catalytic Activity

Protein stability and performance in various natural and artificial systems incorporating many other macromolecules for therapeutic, diagnostic, sensor, and biotechnological applications attract increasing interest with the expansion of these technologies. Here we address the catalytic activity of lysozyme protein (LYZ) in the presence of a polyethylene glycol (PEG) crowder in a broad range of concentrations and temperatures in aqueous solutions of two different molecular mass PEG samples (M_w = 3350 and 10000 g/mol). The phase behavior of PEG-protein solutions is examined by using dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS), while the enzyme denaturing is monitored by using an activity assay (AS) and circular dichroism (CD) spectroscopy. Molecular dynamic (MD) simulations are used to illustrate the effect of PEG concentration on protein stability at high temperatures. The results demonstrate that LYZ residual activity after 1 h incubation at 80 °C is improved from 15% up to 55% with the addition of PEG. The improvement is attributed to two underlying mechanisms. (i) Primarily, the stabilizing effect is due to the suppression of the enzyme aggregation because of the stronger PEG-protein interactions caused by the increased hydrophobicity of PEG and lysozyme at elevated temperatures. (ii) The MD simulations showed that the addition of PEG to some degree stabilizes the secondary structures of the enzyme by delaying unfolding at elevated temperatures. The more pronounced effect is observed with an increase in PEG concentration. This trend is consistent with CD and AS experimental results, where the thermal stability is strengthened with increasing of PEG concentration and molecular mass. The results show that the highest stabilizing effect is approached at the critical overlap concentration of PEG.

Experimentally-driven protein structure modeling

Revolutions in natural and exact sciences started at the dawn of last century have led to the explosion of theoretical, experimental, and computational approaches to determine structures of molecules, complexes, as well as their rich conformational dynamics. Since different experimental methods produce information that is attributed to specific time and length scales, corresponding computational methods have to be tailored to these scales and experiments. These methods can be then combined and integrated in scales, hence producing a fuller picture of molecular structure and motion from the "puzzle pieces" offered by various experiments. Here, we describe a number of computational approaches to utilize experimental data to glance into structure of proteins and understand their dynamics. We will also discuss the limitations and the resolution of the constraints-based modeling approaches. SIGNIFICANCE: Experimentally-driven computational structure modeling and determination is a rapidly evolving alternative to traditional approaches for molecular structure determination. These new hybrid experimental-computational approaches are proving to be a powerful microscope to glance into the structural features of intrinsically or partially disordered proteins, dynamics of molecules and complexes. In this review, we describe various approaches in the field of experimentally-driven computational structure modeling.

Protein aggregation and mitigation strategy in low pH viral inactivation for monoclonal antibody purification

Significant amounts of soluble product aggregates were observed during low-pH viral inactivation (VI) scale-up for an IgG4 monoclonal antibody (mAb IgG4-N1), while small-scale experiments in the same condition showed negligible aggregation. Poor mixing and product exposure to low pH were identified as the root cause. To gain a mechanistic understanding of the problem, protein aggregation properties were studied by varying critical parameters including pH, hold time and protein concentration. Comprehensive biophysical characterization of product monomers and aggregates was performed using fluorescence-size-exclusion chromatography, differential scanning fluorimetry, fluorescence spectroscopy, and dynamic light scattering. Results showed IgG4-N1 partially unfolds at about pH 3.3 where the product molecules still exist largely as monomers owing to strong inter-molecular repulsions and favorable colloidal stability. In the subsequent neutralization step, however, the conformationally changed monomers are prone to aggregation due to weaker inter-molecular repulsions following the pH transition from 3.3 to 5.5. Surface charge calculations using homology modeling suggested that intra-molecular repulsions, especially between CH2 domains, may contribute to the IgG4-N1 unfolding at ≤ pH 3.3. Computational fluid dynamics (CFD) modeling was employed to simulate the conditions of pH titration to reduce the risk of aggregate formation. The low-pH zones during acid addition were characterized using CFD modeling and correlated to the condition causing severe product aggregation. The CFD tool integrated with the mAb solution properties was used to optimize the VI operating parameters for successful scale-up demonstration. Our research revealed the governing aggregation mechanism for IgG4-N1 under acidic conditions by linking its molecular properties and various process-related parameters to macroscopic aggregation phenomena. This study also provides useful insights into the cause and mitigation of low-pH-induced IgG4 aggregation in downstream VI operation.

Determination of the structural ensemble of the molten globule state of a protein by computer simulations

We have used computer simulations to investigate the structural nature of the molten globule (MG) state of canine milk lysozyme. To sample the conformational space efficiently, we performed replica-exchange umbrella sampling simulations with the radius of gyration as a reaction coordinate. We applied the Weighted Histogram Analysis Method to the trajectory of the simulations to obtain the potential of mean force, from which we identified representative structures corresponding to local minima in the free energy surface. The representative structures obtained in this way are in accord with the characteristics of the MG state reported previously by experimental studies. We conjecture that the MG state comprises a series of partially structured states undergoing relatively fast conformational interchange.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Thermal unfolding simulations of NBD1 domain variants reveal structural motifs associated with the impaired folding of F508del-CFTR

The deletion of phenylalanine 508 reshapes the conformational space of the NBD1 domain that populates unique intermediate states that provide insights into the molecular events that underlie the impaired folding of F508del-NBD1.

Structural characterization of MG and pre-MG states of proteins by MD simulations, NMR, and other techniques

Almost all proteins fold via a number of partially structured intermediates such as molten globule (MG) and pre-molten globule states. Understanding the structure of these intermediates at atomic level is often a challenge, as these states are observed under extreme conditions of pH, temperature, and chemical denaturants. Furthermore, several other processes such as chemical modification, site-directed mutagenesis (or point mutation), and cleavage of covalent bond of natural proteins often lead to MG like partially unfolded conformation. However, the dynamic nature of proteins in these states makes them unsuitable for most structure determination at atomic level. Intermediate states studied so far have been characterized mostly by circular dichroism, fluorescence, viscosity, dynamic light scattering measurements, dye binding, infrared techniques, molecular dynamics simulations, etc. There is a limited amount of structural data available on these intermediate states by nuclear magnetic resonance (NMR) and hence there is a need to characterize these states at the molecular level. In this review, we present characterization of equilibrium intermediates by biophysical techniques with special reference to NMR.

Protein thermal stability enhancement by designing salt bridges: a combined computational and experimental study

Protein thermal stability is an important factor considered in medical and industrial applications. Many structural characteristics related to protein thermal stability have been elucidated, and increasing salt bridges is considered as one of the most efficient strategies to increase protein thermal stability. However, the accurate simulation of salt bridges remains difficult. In this study, a novel method for salt-bridge design was proposed based on the statistical analysis of 10,556 surface salt bridges on 6,493 X-ray protein structures. These salt bridges were first categorized based on pairing residues, secondary structure locations, and Cα–Cα distances. Pairing preferences generalized from statistical analysis were used to construct a salt-bridge pair index and utilized in a weighted electrostatic attraction model to find the effective pairings for designing salt bridges. The model was also coupled with B-factor, weighted contact number, relative solvent accessibility, and conservation prescreening to determine the residues appropriate for the thermal adaptive design of salt bridges. According to our method, eight putative salt-bridges were designed on a mesophilic β-glucosidase and 24 variants were constructed to verify the predictions. Six putative salt-bridges leaded to the increase of the enzyme thermal stability. A significant increase in melting temperature of 8.8, 4.8, 3.7, 1.3, 1.2, and 0.7°C of the putative salt-bridges N437K–D49, E96R–D28, E96K–D28, S440K–E70, T231K–D388, and Q277E–D282 was detected, respectively. Reversing the polarity of T231K–D388 to T231D–D388K resulted in a further increase in melting temperatures by 3.6°C, which may be caused by the transformation of an intra-subunit electrostatic interaction into an inter-subunit one depending on the local environment. The combination of the thermostable variants (N437K, E96R, T231D and D388K) generated a melting temperature increase of 15.7°C. Thus, this study demonstrated a novel method for the thermal adaptive design of salt bridges through inference of suitable positions and substitutions.

Native electrospray ionization and electron-capture dissociation for comparison of protein structure in solution and the gas phase

The importance of protein and protein-complex structure motivates improvements in speed and sensitivity of structure determination in the gas phase and comparison with that in solution or solid state. An opportunity for the gas phase measurement is mass spectrometry (MS) combined with native electrospray ionization (ESI), which delivers large proteins and protein complexes in their near-native states to the gas phase. In this communication, we describe the combination of native ESI, electron-capture dissociation (ECD), and top-down MS for exploring the structures of ubiquitin and cytochrome c in the gas phase and their relation to those in the solid-state and solution. We probe structure by comparing the protein's flexible regions, as predicted by the B-factor in X-ray crystallography, with the ECD fragments. The underlying hypothesis is that maintenance of structure gives fragments that can be predicted from B-factors. This strategy may be applicable in general when X-ray structures are available and extendable to the study of intrinsically disordered proteins.

The mechanism by which a propeptide-encoded pH sensor regulates spatiotemporal activation of furin

The proprotein convertase furin requires the pH gradient of the secretory pathway to regulate its multistep, compartment-specific autocatalytic activation. Although His-69 within the furin prodomain serves as the pH sensor that detects transport of the propeptide-enzyme complex to the trans-Golgi network, where it promotes cleavage and release of the inhibitory propeptide, a mechanistic understanding of how His-69 protonation mediates furin activation remains unclear. Here we employ biophysical, biochemical, and computational approaches to elucidate the mechanism underlying the pH-dependent activation of furin. Structural analyses and binding experiments comparing the wild-type furin propeptide with a nonprotonatable His-69 → Leu mutant that blocks furin activation in vivo revealed protonation of His-69 reduces both the thermodynamic stability of the propeptide as well as its affinity for furin at pH 6.0. Structural modeling combined with mathematical modeling and molecular dynamic simulations suggested that His-69 does not directly contribute to the propeptide-enzyme interface but, rather, triggers movement of a loop region in the propeptide that modulates access to the cleavage site and, thus, allows for the tight pH regulation of furin activation. Our work establishes a mechanism by which His-69 functions as a pH sensor that regulates compartment-specific furin activation and provides insights into how other convertases and proteases may regulate their precise spatiotemporal activation.

Capturing native/native like structures with a physico-chemical metric (pcSM) in protein folding

Specification of the three dimensional structure of a protein from its amino acid sequence, also called a "Grand Challenge" problem, has eluded a solution for over six decades. A modestly successful strategy has evolved over the last couple of decades based on development of scoring functions (e.g. mimicking free energy) that can capture native or native-like structures from an ensemble of decoys generated as plausible candidates for the native structure. A scoring function must be fast enough in discriminating the native from unfolded/misfolded structures, and requires validation on a large data set(s) to generate sufficient confidence in the score. Here we develop a scoring function called pcSM that detects true native structure in the top 5 with 93% accuracy from an ensemble of candidate structures. If we eliminate the native from ensemble of decoys then pcSM is able to capture near native structure (RMSD<=5Ǻ) in top 10 with 86% accuracy. The parameters considered in pcSM are a C-alpha Euclidean metric, secondary structural propensity, surface areas and an intramolecular energy function. pcSM has been tested on 415 systems consisting 142,698 decoys (public and CASP-largest reported hitherto in literature). The average rank for the native is 2.38, a significant improvement over that existing in literature. In-silico protein structure prediction requires robust scoring technique(s). Therefore, pcSM is easily amenable to integration into a successful protein structure prediction strategy. The tool is freely available at http://www.scfbio-iitd.res.in/software/pcsm.jsp.

Using simulations to provide the framework for experimental protein folding studies

Molecular dynamics simulations are a powerful theoretical tool to model the protein folding process in atomistic details under realistic conditions. Combined with a number of experimental techniques, simulations provide a detailed picture of how a protein folds or unfolds in the presence of explicit solvent and other molecular species, such as cosolvents, osmolytes, cofactors, active binding partners or inert crowding agents. The denaturing effects of temperature, pressure and external mechanical forces can also be probed. Qualitative and quantitative agreement with experiment contributes to a comprehensive molecular picture of protein states along the folding/unfolding pathway. The variety of systems examined reveals key features of the protein folding process.

About the structural role of disulfide bridges in serum albumins: evidence from protein simulated unfolding

The role of the 17 disulfide (S-S) bridges in preserving the native conformation of human serum albumin (HSA) is investigated by performing classical molecular dynamics (MD) simulations on protein structures with intact and, respectively, reduced S-S bridges. The thermal unfolding simulations predict a clear destabilization of the protein secondary structure upon reduction of the S-S bridges as well as a significant distortion of the tertiary structure that is revealed by the changes in the protein native contacts fraction. The effect of the S-S bridges reduction on the protein compactness was tested by calculating Gibbs free energy profiles with respect to the protein gyration radius. The theoretical results obtained using the OPLS-AA and the AMBER ff03 force fields are in agreement with the available experimental data. Beyond the validation of the simulation method, the results here reported provide new insights into the mechanism of the protein reductive/oxidative unfolding/folding processes. It is predicted that in the native conformation of the protein, the thiol (-SH) groups belonging to the same reduced S-S bridge are located in potential wells that maintain them in contact. The -SH pairs can be dispatched by specific conformational transitions of the peptide chain located in the neighborhood of the cysteine residues.

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded proteins.

Propeptides are sufficient to regulate organelle-specific pH-dependent activation of furin and proprotein convertase 1/3

The proprotein convertases (PCs) furin and proprotein convertase 1/3 (PC1) cleave substrates at dibasic residues along the eukaryotic secretory/endocytic pathway. PCs are evolutionarily related to bacterial subtilisin and are synthesized as zymogens. They contain N-terminal propeptides (PRO) that function as dedicated catalysts that facilitate folding and regulate activation of cognate proteases through multiple-ordered cleavages. Previous studies identified a histidine residue (His69) that functions as a pH sensor in the propeptide of furin (PRO(FUR)), which regulates furin activation at pH~6.5 within the trans-Golgi network. Although this residue is conserved in the PC1 propeptide (PRO(PC1)), PC1 nonetheless activates at pH~5.5 within the dense core secretory granules. Here, we analyze the mechanism by which PRO(FUR) regulates furin activation and examine why PRO(FUR) and PRO(PC1) differ in their pH-dependent activation. Sequence analyses establish that while both PRO(FUR) and PRO(PC1) are enriched in histidines when compared with cognate catalytic domains and prokaryotic orthologs, histidine content in PRO(FUR) is ~2-fold greater than that in PRO(PC1), which may augment its pH sensitivity. Spectroscopy and molecular dynamics establish that histidine protonation significantly unfolds PRO(FUR) when compared to PRO(PC1) to enhance autoproteolysis. We further demonstrate that PRO(FUR) and PRO(PC1) are sufficient to confer organelle sensing on folding and activation of their cognate proteases. Swapping propeptides between furin and PC1 transfers pH-dependent protease activation in a propeptide-dictated manner in vitro and in cells. Since prokaryotes lack organelles and eukaryotic PCs evolved from propeptide-dependent, not propeptide-independent prokaryotic subtilases, our results suggest that histidine enrichment may have enabled propeptides to evolve to exploit pH gradients to activate within specific organelles.

Transient molten globules and metastable aggregates induced by brief exposure of a monoclonal IgG to low pH

The presence of aggregates in therapeutic protein formulations is of great concern due to quality, safety, and efficacy issues. Nonetheless, the mechanisms and kinetics of protein aggregation are only partly understood. In this study, metastable immunoglobulin G (IgG) aggregates induced by a brief exposure to pH 1 were kept at 4°C and analyzed over time by size-exclusion chromatography (SEC), nanoparticle tracking analysis, light obscuration, dynamic light scattering, fluorescence spectroscopy, and circular dichroism. The results show the formation of polydisperse aggregates (from dimers to 10-μm particles) shortly after the pH-shift stress. These aggregates increased in size and number over time until a pseudo-equilibrium was reached after 5-7 days. The presence of transient, partially unfolded monomers (molten globules) was detected by SEC with online fluorescent dye detection. The molten globules seemed to either refold into the native state or become involved in aggregation pathways. Seeding pH-shift-induced aggregates into unstressed IgG did not accelerate aggregation during incubation for 3 weeks at 55°C. These results reinforce the role of partially unfolded species in the aggregation of therapeutic proteins. We conclude that the formation of pH-shift-induced IgG aggregates is likely driven by downhill polymerization, as a consequence of successive additions of molten globular monomers.

Native electrospray and electron-capture dissociation FTICR mass spectrometry for top-down studies of protein assemblies

The high sensitivity, extended mass range, and fast data acquisition/processing of mass spectrometry and its coupling with native electrospray ionization (ESI) make the combination complementary to other biophysical methods of protein analysis. Protein assemblies with molecular masses up to MDa are now accessible by this approach. Most current approaches have used quadrupole/time-of-flight tandem mass spectrometry, sometimes coupled with ion mobility, to reveal stoichiometry, shape, and dissociation of protein assemblies. The amino-acid sequence of the subunits, however, still relies heavily on independent bottom-up proteomics. We describe here an approach to study protein assemblies that integrates electron-capture dissociation (ECD), native ESI, and FTICR mass spectrometry (12 T). Flexible regions of assembly subunits of yeast alcohol dehydrogenase (147 kDa), concanavalin A (103 kDa), and photosynthetic Fenna-Matthews-Olson antenna protein complex (140 kDa) can be sequenced by ECD or "activated-ion" ECD. Furthermore, noncovalent metal-binding sites can also be determined for the concanavalin A assembly. Most importantly, the regions that undergo fragmentation, either from one of the termini by ECD or from the middle of a protein, as initiated by CID, correlate well with the B-factor from X-ray crystallography of that protein. This factor is a measure of the extent an atom can move from its coordinated position as a function of temperature or crystal imperfections. The approach provides not only top-down proteomics information of the complex subunits but also structural insights complementary to those obtained by ion mobility.

Coarse-grained modeling study of nonpeptide RGD ligand density and PEG molecular weight on the conformation of poly(γ-glutamyl-glutamate) paclitaxel conjugates

Molecular shape, flexibility, and surface hydrophilicity are thought to influence the ability of nanoparticles to cross biological barriers during drug delivery. In this study, coarse-grained (CG) molecular dynamics (MD) simulations were used to study these properties of a polymer-drug construct in potential clinical development: poly(γ-glutamyl-glutamate)-paclitaxel-poly(ethylene glycol) nonpeptide RGD (PGG-PTX-PEG-npRGD), a linear glutamyl-glutamate polymer with paclitaxel and poly(ethylene glycol)-nonpeptide RGD side groups. It was hypothesized that the PEG molecular weight (MW) (500 Da; 1,000 Da; and 2,000 Da) and nonpeptide RGD ligand density (4, 8, 12, and 16 per molecule), respectively, may have advantageous effects on the shape, flexibility, and surface hydrophilicity of PGG-PTX-PEG-npRGD. Circular dichroism spectroscopy was used to suggest initial structures for the all-atom (AA) models of PGG-PTX-PEG-npRGD, which were further converted to CG models using a commercially available mapping algorithm. Due to its semi-flexibility, PGG-PTX-PEG-npRGD is not limited to one specific conformation. Thus, CG MD simulations were run until statistical equilibrium, at which PGG-PTX-PEG-npRGD is represented as an ensemble of statistically similar conformations. The size of a PGG-PTX-PEG-npRGD molecule is not affected by the PEG MW or the nonpeptide RGD density, but higher PEG MW results in increased surface density of a PGG-PTX-PEG-npRGD molecule. Most PGG-PTX-PEG-npRGD shapes are globular, although filamentous shapes were also observed in the PEG500 and PEG1000 molecules. PEG500 and PEG1000 molecules are more flexible than PEG2000 systems. A higher presence of npRGD ligands results in decrease surface hydrophilicity of PGG-PTX-PEG-npRGD. These results indicate that the PGG-PTX-PEG1000-npRGD₄ and PGG-PTX-PEG1000-npRGD₈ molecules are the most efficacious candidates and are further recommended for experimental preclinical studies.

Protein stability and structure in HIC: hydrogen exchange experiments and COREX calculations

Hydrogen exchange mass spectrometry (HXMS) coupled to proteolytic digestion has been used to probe the conformation of bovine β-lactoglobulin (BLG), bovine α-lactalbumin (BLA), and human serum albumin (HSA) in solution and while adsorbed to the hydrophobic interaction chromatography media Phenyl Sepharose 6FF. All three proteins show evidence of EX1 exchange kinetics, indicating a loss of stability on the surface. HX protection patterns for all three proteins also indicate that the unfolded form is only partially solvent exposed. The hydrogen-deuterium exchange patterns of BLG and BLA on the surface suggest a structure that resembles each protein's respective solution phase molten globule state. The low stability of Domain II of HSA observed on Phenyl Sepharose 6FF also suggests a link to solution stability because Domain II is frequently cited as the least stable domain in solution unfolding pathways. COREX, an algorithm used to compute protein folding stabilities, correctly predicts solution hydrogen-deuterium exchange patterns for BLG and offers insight into its adsorbed phase stabilities but is unreliable for BLA predictions. The results of this work demonstrate a link between solution-phase local stability patterns and the nature of partially unfolded states that proteins can adopt on HIC surfaces.

Dynameomics: large-scale assessment of native protein flexibility

Structure is only the first step in understanding the interactions and functions of proteins. In this paper, we explore the flexibility of proteins across a broad database of over 250 solvated protein molecular dynamics simulations in water for an aggregate simulation time of approximately 6 micros. These simulations are from our Dynameomics project, and these proteins represent approximately 75% of all known protein structures. We employ principal component analysis of the atomic coordinates over time to determine the primary axis and magnitude of the flexibility of each atom in a simulation. This technique gives us both a database of flexibility for many protein fold families and a compact visual representation of a particular protein's native-state conformational space, neither of which are available using experimental methods alone. These tools allow us to better understand the nature of protein motion and to describe its relationship to other structural and dynamical characteristics. In addition to reporting general properties of protein flexibility and detailing many dynamic motifs, we characterize the relationship between protein native-state flexibility and early events in thermal unfolding and show that flexibility predicts how a protein will begin to unfold. We provide evidence that fold families have conserved flexibility patterns, and family members who deviate from the conserved patterns have very low sequence identity. Finally, we examine novel aspects of highly inflexible loops that are as important to structural integrity as conventional secondary structure. These loops, which are difficult if not impossible to locate without dynamic data, may constitute new structural motifs.

Probing protein fold space with a simplified model

We probe the stability and near-native energy landscape of protein fold space using powerful conformational sampling methods together with simple reduced models and statistical potentials. Fold space is represented by a set of 280 protein domains spanning all topological classes and having a wide range of lengths (33-300 residues) amino acid composition and number of secondary structural elements. The degrees of freedom are taken as the loop torsion angles. This choice preserves the native secondary structure but allows the tertiary structure to change. The proteins are represented by three-point per residue, three-dimensional models with statistical potentials derived from a knowledge-based study of known protein structures. When this space is sampled by a combination of parallel tempering and equi-energy Monte Carlo, we find that the three-point model captures the known stability of protein native structures with stable energy basins that are near-native (all alpha: 4.77 A, all beta: 2.93 A, alpha/beta: 3.09 A, alpha+beta: 4.89 A on average and within 6 A for 71.41%, 92.85%, 94.29% and 64.28% for all-alpha, all-beta, alpha/beta and alpha+beta, classes, respectively). Denatured structures also occur and these have interesting structural properties that shed light on the different landscape characteristics of alpha and beta folds. We find that alpha/beta proteins with alternating alpha and beta segments (such as the beta-barrel) are more stable than proteins in other fold classes.

Characterization of protein-folding pathways by reduced-space modeling

Ab initio simulations of the folding pathways are currently limited to very small proteins. For larger proteins, some approximations or simplifications in protein models need to be introduced. Protein folding and unfolding are among the basic processes in the cell and are very difficult to characterize in detail by experiment or simulation. Chymotrypsin inhibitor 2 (CI2) and barnase are probably the best characterized experimentally in this respect. For these model systems, initial folding stages were simulated by using CA-CB-side chain (CABS), a reduced-space protein-modeling tool. CABS employs knowledge-based potentials that proved to be very successful in protein structure prediction. With the use of isothermal Monte Carlo (MC) dynamics, initiation sites with a residual structure and weak tertiary interactions were identified. Such structures are essential for the initiation of the folding process through a sequential reduction of the protein conformational space, overcoming the Levinthal paradox in this manner. Furthermore, nucleation sites that initiate a tertiary interactions network were located. The MC simulations correspond perfectly to the results of experimental and theoretical research and bring insights into CI2 folding mechanism: unambiguous sequence of folding events was reported as well as cooperative substructures compatible with those obtained in recent molecular dynamics unfolding studies. The correspondence between the simulation and experiment shows that knowledge-based potentials are not only useful in protein structure predictions but are also capable of reproducing the folding pathways. Thus, the results of this work significantly extend the applicability range of reduced models in the theoretical study of proteins.

Dynamics of lysozyme structure network: probing the process of unfolding

Recently we showed that the three-dimensional structure of proteins can be investigated from a network perspective, where the amino acid residues represent the nodes in the network and the noncovalent interactions between them are considered for the edge formation. In this study, the dynamical behavior of such networks is examined by considering the example of T4 lysozyme. The equilibrium dynamics and the process of unfolding are followed by simulating the protein at 300 K and at higher temperatures (400 K and 500 K), respectively. The snapshots of the protein structure from the simulations are represented as protein structure networks in which the strength of the noncovalent interactions is considered an important criterion in the construction of edges. The profiles of the network parameters, such as the degree distribution and the size of the largest cluster (giant component), were examined as a function of interaction strength at different temperatures. Similar profiles are seen at all the temperatures. However, the critical strength of interaction (Icritical) and the size of the largest cluster at all interaction strengths shift to lower values at 500 K. Further, the folding/unfolding transition is correlated with contacts evaluated at Icritical and with the composition of the top large clusters obtained at interaction strengths greater than Icritical. Finally, the results are compared with experiments, and predictions are made about the residues, which are important for stability and folding. To summarize, the network analysis presented in this work provides insights into the details of the changes occurring in the protein tertiary structure at the level of amino acid side-chain interactions, in both the equilibrium and the unfolding simulations. The method can also be employed as a valuable tool in the analysis of molecular dynamics simulation data, since it captures the details at a global level, which may elude conventional pairwise interaction analysis.

Hydration dependence of the mass fractal dimension and anomalous diffusion of vibrational energy in proteins

Vibrational dynamics of proteins and energy flow depend on protein geometry as well as interactions of a protein molecule with the surrounding solvent. We compute the mass fractal dimension D of proteins ranging from 100 to over 10,000 amino acids comparing values for the bare protein with those computed when buried and hydration waters are included in the calculation. Including water in the calculation increases D by about 0.3 to 2.87 on average above D computed for the dehydrated protein. The mass fractal dimension of proteins that are partially unfolded by molecular dynamics (MD) simulation is also computed and found to vary little when the radius of gyration changes within about 10% of that for the Protein Data Bank structure. MD simulations of vibrational energy diffusion in proteins reveal that the exponent characterizing anomalous diffusion of vibrational energy does not change much with hydration, which is seen to be due to an increase in the spectral dimension with hydration by a factor similar to the increase in D.

Multiple folding mechanisms of protein ubiquitin

Based on the C(alpha) Go-type model, the folding kinetics and mechanisms of protein ubiquitin with mixed alpha/beta topology are studied by molecular dynamics simulations. The relaxation kinetics shows that there are three phases, namely the major phase, the intermediate phase and the slowest minor phase. The existence of these three phases are relevant to the phenomenon found in experiments. According to our simulations, the folding at high temperatures around the folding transition temperature T(f) is of a two-state process, and the folding nucleus is consisted of contacts between the front end of alpha-helix and the turn(4). The folding at low temperature (approximately T = 0.8) is also studied, where an A-state like structure is found lying on the major folding pathway. The appearance of this structure is related to the stability of the first part (residue 1-51) of protein ubiquitin. As the temperature decreases, the formation of secondary structures, tertiary structures and collapse of the protein are found to be decoupled gradually and the folding mechanism changes from the nucleation-condensation to the diffusion-collision. This feature indicates a unifying common folding mechanism for proteins. The intermediate phase is also studied and is found to represent a folding process via a long-lived intermediate state which is stabilized by strong interactions between the beta(1) and the beta(5) strand. These strong interactions are important for the function of protein ubiquitin as a molecular chaperone. Thus the intermediate phase is assumed as a byproduct of the requirement of protein function. In addition, the validity of the current Go-model is also investigated, and a lower limited temperature for protein ubiquitin T(limit) = 0.8 is proposed. At temperatures higher than this value, the kinetic traps due to glass dynamics cannot be significantly populated and the intermediate states can be reliably identified although there is slight chevron rollover in the folding rates. At temperature lower than T(limit), however, the traps due to glass dynamics become dominant and may be mistaken for real intermediate states. This limitation of valid temperature range prevents us to reveal the burst phase intermediate in the major folding phase since it might only be stabilized at temperatures lower than T(limit), according to experiments. Our works show that caution must be taken when studying low-temperature intermediate states by using the C(alpha) Go-models.

Terminal ion pairs stabilize the second β-hairpin of the B1 domain of protein G

The effects of terminal ion pairs on the stability of a beta-hairpin peptide corresponding to the C-terminal residues of the B1 domain of protein G were determined using thermal unfolding as monitored by nuclear magnetic resonance and circular dichroism spectroscopy. Molecular dynamics (MD) simulations were also performed to examine the effect of ion pairs on the structures. Eight peptides were studied including the wild type (G41) and the N-terminal modified sequences that had the first residue deleted (E42), replaced with a Lys (K41), or extended by an additional Gly (G40). Acetylated variants were made to examine the effect of removing the positive N-terminal charge on beta-hairpin stability. The rank in stability determined experimentally is K41 > E42 approximately G41 approximately G40 > Ac-K41 > Ac-E42 approximately Ac-G41 > Ac-G40. The Tm of the K41 peptide is 12 degrees C higher than G41, while the Tm values for the acetylated peptides are less than their unacetylated forms by more than 15 degrees C. NOE cross-peaks between side-chain methylene groups at the N- and C-termini and larger CalphaH shifts compared to random values are seen for K41. The addition of 20% methanol increases the stability in K41 and G41. The MD studies complement these results by showing that the charged N-terminus is important to stability. The type of ion pair observed varies with peptide, and when formed the simulations show that the ion pair can prevent fraying of the beta-strands through electrostatic and hydrophobic contacts. Therefore, introducing favorable electrostatic interactions at the N- and C-termini can substantially enhance beta-hairpin stability and help define the structure.

Symmetric connectivity of secondary structure elements enhances the diversity of folding pathways

The influence of native connectivity of secondary structure elements (SSE) on folding is studied using coarse-grained models of proteins with mixed alpha and beta structure and the analysis of the structural database of wild-type proteins. We found that the distribution of SSE along a sequence determines the diversity of folding pathways. If alpha and beta SSE are localized in different parts of a sequence, the diversity of folding pathways is restricted. An even (symmetric) distribution of alpha and beta SSE with respect to sequence midpoint favors multiple folding routes. Simulations are supplemented by the database analysis of the distribution of SSE in wild-type protein sequences. On an average, two-thirds of wild-type proteins with mixed alpha and beta structure have symmetric distribution of alpha and beta SSE. The propensity for symmetric distribution of SSE is especially evident for large proteins with the number of SSE > or = 10. We suggest that symmetric SSE distribution in protein sequences may arise due to nearly random allocation of alpha and beta structure along wild-type sequences. The tendency of long sequences to misfold is perhaps compensated by the enhanced pathway diversity. In addition, folding pathways are shown to progress via hierarchic assembly of SSE in accordance with their proximity along a sequence. We demonstrate that under mild denaturation conditions folding and unfolding pathways are similar. However, the reversibility of folding/unfolding pathways is shown to depend on the distribution of SSE. If alpha and beta SSE are localized in different parts of a sequence, folding and unfolding pathways are likely to coincide.

Low-temperature and high-pressure induced swelling of a hydrophobic polymer-chain in aqueous solution

We report molecular dynamics simulations of a hydrophobic polymer-chain in aqueous solution between 260 K and 420 K at pressures of 1 bar, 3000 bar, and 4500 bar. The simulations reveal a hydrophobically collapsed structure at low pressures and high temperatures. At 3000 bar and about 260 K and at 4500 bar and about 260 K, however, an abrupt transition to a swelled state is observed. The transition is driven by a smaller volume and a remarkably strong lower enthalpy of the swelled state, indicating a steep positive slope of the corresponding transition line. The swelling is strongly stabilized by the energetically favorable state of water in the polymer's hydrophobic first hydration shell at low temperatures. This finding is consistent with the observation of a positive heat capacity of hydrophobic solvation. Moreover, the slope and location of the estimated swelling transition line for the collapsed hydrophobic chain coincides remarkably well with the cold denaturation transition of proteins.