The native-state ensemble of proteins provides clues for folding, misfolding and function

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.

Human aldose reductase unfolds through an intermediate

Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway.  Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR.

Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. 

Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR.

Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state.

Discrimination power of knowledge-based potential dictated by the dominant energies in native protein structures

Extracting a well-designed energy function is important for protein structure evaluation. Knowledge-based potential functions are one type of the energy functions which can be obtained from known protein structures. The pairwise potential between atom types is approximated using Boltzmann's law which relates the frequency of atom types to its potential. The total energy is approximated as a summation of pairwise potential between the atomic pairs. In the present study, the performance of knowledge-based potential function was assessed based on the strength of interaction between groups of amino acids. The dominant energies involved in the pairwise potentials were revealed by eigenvalue analysis of the matrix, the elements of which represent the energy between amino acids. For this purpose, the matrix including the mean of the energies of residue-residue interaction types was constructed using 500 native protein structures. The matrix has a dominant eigenvalue and amino acids, with LEU, VAL, ILE, PHE, TYR, ALA and TRP having high values along the dominant eigenvector. The results show that the ranking of amino acids is consistent with the power of amino acids in discriminating native structures using K-alphabet reduced model. In the reduced interactions, only amino acids from a subset of all 20 amino acids, along with their interactions are considered to assess the energy. In the K-alphabet reduced model, the reduced structures are constructed based on only the K-amino acid types. The dominant K-alphabet reduced model derived for the k-first amino acids in the list [LEU, VAL, PHE, ILE, TYR, ALA, TRP] of amino acids has the best discrimination of native structure among all possible K-alphabet reduced models. Knowledge-based potentials might be improved with a new strategy.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

The mechanism of water/ion exchange at a protein surface: a weakly bound chloride in Helicobacter pylori apoflavodoxin

Binding/unbinding of small ligands, such as ions, to/from proteins influences biochemical processes such as protein folding, enzyme catalysis or protein/ligand recognition. We have investigated the mechanism of chloride/water exchange at a protein surface (that of the apoflavodoxin from Helicobacter pylori) using classical all-atom molecular dynamics simulations. They reveal a variety of chloride exit routes and residence times; the latter is related to specific coordination modes of the anion. The role of solvent molecules in the mechanism of chloride unbinding has been studied in detail. We see no temporary increase in chloride coordination along the release process. Instead, the coordination of new water molecules takes place in most cases after the chloride/protein atom release event has begun. Moreover, the distribution function of water entrance events into the first chloride solvation shell peaks after chloride protein atom dissociation events. All these observations together seem to indicate that water molecules simply fill the vacancies left by the previously coordinating protein residues. We thus propose a step-by-step dissociation pathway in which protein/chloride interactions gradually break down before new water molecules progressively fill the vacant positions left by protein atoms. As observed for other systems, water molecules associated with bound chloride or with protein atoms have longer residence times than those bound to the free anion. The implications of the exchange mechanism proposed for the binding of the FMN (Flavin Mononucleotide) protein cofactor are discussed.

Experimental and computational evidence on conformational fluctuations as a source of catalytic defects in genetic diseases

Catalytic mutants causing inherited type III galactosemia alter active site structural dynamics and shift the native conformational equilibrium towards inactive conformations.

Folding of a cyclin box: linking multitarget binding to marginal stability, oligomerization, and aggregation of the retinoblastoma tumor suppressor AB pocket domain

The retinoblastoma tumor suppressor (Rb) controls the proliferation, differentiation, and survival of cells in most eukaryotes with a role in the fate of stem cells. Its inactivation by mutation or oncogenic viruses is required for cellular transformation and eventually carcinogenesis. The high conservation of the Rb cyclin fold prompted us to investigate the link between conformational stability and ligand binding properties of the RbAB pocket domain. RbAB unfolding presents a three-state transition involving cooperative secondary and tertiary structure changes and a partially folded intermediate that can oligomerize. The first transition corresponds to unfolding of the metastable B subdomain containing the binding site for the LXCXE motif present in cellular and viral targets, and the second transition corresponds to the stable A subdomain. The low thermodynamic stability of RbAB translates into a propensity to rapidly oligomerize and aggregate at 37 °C (T50 = 28 min) that is suppressed by human papillomavirus E7 and E2F peptide ligands, suggesting that Rb is likely stabilized in vivo through binding to target proteins. We propose that marginal stability and associated oligomerization may be conserved for function as a "hub" protein, allowing the formation of multiprotein complexes, which could constitute a robust mechanism to retain its cell cycle regulatory role throughout evolution. Decreased stability and oligomerization are shared with the p53 tumor suppressor, suggesting a link between folding and function in these two essential cell regulators that are inactivated in most cancers and operate within multitarget signaling pathways.

Electrostatics of hydrogen exchange for analyzing protein flexibility

Electrostatic interactions at the protein-aqueous interface modulate the reactivity of solvent-exposed backbone amides by a factor of at least a billion fold. The brief (∼10 ps) lifetime of the peptide anion formed during the hydroxide-catalyzed exchange reaction helps enable the experimental rates to be robustly predictable by continuum dielectric methods. Since this ability to predict the structural dependence of exchange reactivity also applies to the protein amide hydrogens that are only rarely exposed to the bulk solvent phase, electrostatic analysis of the experimental exchange rates provides an effective assessment of whether a given model ensemble is consistent with the properly weighted Boltzmann conformational distribution of the protein native state.

Design and structure of an equilibrium protein folding intermediate: a hint into dynamical regions of proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.

Energetics of nucleotide-induced DnaK conformational states

Hsp70 chaperones are molecular switches that use the free energy of ATP binding and hydrolysis to modulate their affinity for protein substrates and, most likely, to remodel non-native interactions allowing proper substrate folding. By means of isothermal titration calorimetry, we have measured the thermodynamics of ATP and ADP binding to (i) wild-type DnaK, the main bacterial Hsp70; (ii) two single-point mutants, DnaK(T199A), which lacks ATPase activity but maintains conformational changes similar to those observed in the wild-type protein, and DnaK(R151A), defective in interdomain communication; and iii) two deletion mutants, the isolated nucleotide binding domain (K-NBD) and a DeltaLid construct [DnaK(1-507)]. At 25 degrees C, ATP binding to DnaK results in a fast endothermic and a slow exothermic process due to ATP hydrolysis. We demonstrate that the endothermic event is due to the allosteric coupling between ATP binding to the nucleotide binding domain and the conformational rearrangement of the substrate binding domain. The interpretation of our data is compatible with domain docking upon ATP binding and shows that this conformational change carries an energy penalty of ca. 1 kcal/mol. The conformational energy stored in the ATP-bound DnaK state, together with the free energy of ATP hydrolysis, can be used in remodeling bound substrates.

Rescue of misfolded proteins and stabilization by small molecules

Increasing stability of functional proteins by binding small compounds and ions has long been used to extend shelf-life of protein formulations in the pharmacological and biotechnological industry. Likewise, the therapeutic application of small molecules for in vivo recovery and maintenance of structure and function of proteins is steadily increasing. Compounds that can rescue misfolded proteins by stimulating their correct folding and/or the stabilization of native-like conformations in vivo are referred to as pharmacological chaperones. Here we present thermal-shift and isothermal methods for the high-throughput screening of stabilizing pharmacological chaperones for soluble and membrane proteins. The effect of selected hit compounds on the kinetics of protein synthesis is further evaluated by an in vitro transcription-translation rapid translation system. These procedures can be integrated in an interdisciplinary and translational approach for the search of personalized pharmacological chaperones in genetic misfolding diseases.

Native state dynamics drive the unfolding of the SH3 domain of PI3 kinase at high denaturant concentration

Little is known about the role of protein dynamics in directing protein unfolding along a specific pathway and about the role played by chemical denaturants in modulating the dynamics and the initiation of unfolding. In this study, deuterium-hydrogen exchange (HX) detected by electrospray ionization mass spectrometry (ESI-MS) was used to study the unfolding of the SH3 domain of the PI3 kinase. Unfolding on the principal unfolding pathway occurs in 2 steps, both in the absence and in the presence of 1.8 M guanidine hydrochloride (GdnHCl). In both cases, the first step leads to the formation of an intermediate, I(N), with 5 fewer protected amide hydrogen sites than in N. In the second step, I(N) loses the structure protecting the remaining 14 amide hydrogen sites from HX as it unfolds completely. ESI-MS analysis of fragments of the protein created by proteolytic digestion, after completion of the HX reaction, shows that I(N) has lost protection against HX in the same segments of native structure during unfolding in the absence and presence of 1.8 M GdnHCl. Hence, GdnHCl does not appear to play a direct active role in the initiation of unfolding. However, at higher GdnHCl concentrations, a second unfolding pathway is shown to compete effectively with the N <--> I(N) <--> U pathway. In this way, the denaturant modulates the energy landscape of unfolding.

Peptide conformer acidity analysis of protein flexibility monitored by hydrogen exchange

The amide hydrogens that are exposed to solvent in the high-resolution X-ray structures of ubiquitin, FK506-binding protein, chymotrypsin inhibitor 2, and rubredoxin span a billion-fold range in hydroxide-catalyzed exchange rates which are predictable by continuum dielectric methods. To facilitate analysis of transiently accessible amides, the hydroxide-catalyzed rate constants for every backbone amide of ubiquitin were determined under near physiological conditions. With the previously reported NMR-restrained molecular dynamics ensembles of ubiquitin (PDB codes 2NR2 and 2K39) used as representations of the Boltzmann-weighted conformational distribution, nearly all of the exchange rates for the highly exposed amides were more accurately predicted than by use of the high-resolution X-ray structure. More strikingly, predictions for the amide hydrogens of the NMR relaxation-restrained ensemble that become exposed to solvent in more than one but less than half of the 144 protein conformations in this ensemble were almost as accurate. In marked contrast, the exchange rates for many of the analogous amides in the residual dipolar coupling-restrained ubiquitin ensemble are substantially overestimated, as was particularly evident for the Ile 44 to Lys 48 segment which constitutes the primary interaction site for the proteasome targeting enzymes involved in polyubiquitylation. For both ensembles, “excited state” conformers in this active site region having markedly elevated peptide acidities are represented at a population level that is 10² to 10³ above what can exist in the Boltzmann distribution of protein conformations. These results indicate how a chemically consistent interpretation of amide hydrogen exchange can provide insight into both the population and the detailed structure of transient protein conformations.

Molten globule and native state ensemble of Helicobacter pylori flavodoxin: can crowding, osmolytes or cofactors stabilize the native conformation relative to the molten globule?

Partly unfolded protein conformations close in energy to the native state may be involved in protein functioning and also be related to folding diseases, but yet their structure and energetics are poorly understood. One such conformation, the monomeric and well-behaved molten globule of Helicobacter pylori apoflavodoxin, is here investigated to provide, in a wide pH interval, a complete thermodynamic description of its unfolding equilibrium and the equilibrium linking molten globule and native state. All thermodynamic and molecular properties of the molten globule here analyzed are characteristic of a partly unfolded conformation, and their differences with those of the native state are typically quantitative rather than qualitative. The stability data depict a native state ensemble where the relative populations of the different intermediates are strongly modulated by pH. Whereas the molten globule is dominant at pH 2.0, at neutral pH it is just the least stable of three partly unfolded intermediates populated by this protein. It is of interest that the energy rank of these intermediates at pH 7.0 is consistent with their likelihood to overcome the native state and become the more stable conformation when the native state protein is subjected to heat or mutation stress. Given the small volume difference between molten globule and native state, neither crowding agents nor osmolytes can drive the molten globule back to the native state. This observation, which is in qualitative accord with predictions of simple excluded volume theory, indicates that molecular crowding in vivo is not an effective mechanism to minimize partial unfolding events leading to equilibrium intermediates.

The flavodoxin from Helicobacter pylori: structural determinants of thermostability and FMN cofactor binding

Flavodoxin has been recently recognized as an essential protein for a number of pathogenic bacteria including Helicobacter pylori, where it has been proposed to constitute a target for antibacterial drug development. One way we are exploring to screen for novel inhibitory compounds is to perform thermal upshift assays, for which a detailed knowledge of protein thermostability and cofactor binding properties is of great help. However, very little is known on the stability and ligand binding properties of H. pylori flavodoxin, and its peculiar FMN binding site together with the variety of behaviors observed within the flavodoxin family preclude extrapolations. We have thus performed a detailed experimental and computational analysis of the thermostability and cofactor binding energetics of H. pylori flavodoxin, and we have found that the thermal unfolding equilibrium is more complex that any other previously described for flavodoxins as it involves the accumulation of two distinct equilibrium intermediates. Fortunately the entire stability and binding data can be satisfactorily fitted to a model, summarized in a simple phase diagram, where the cofactor only binds to the native state. On the other hand, we show how variability of thermal unfolding behavior within the flavodoxin family can be predicted using structure-energetics relationships implemented in the COREX algorithm. The different distribution and ranges of local stabilities of the Anabaena and H. pylori apoflavodoxins explain the essential experimental differences observed: much lower Tm1, greater resistance to global unfolding, and more pronounced cold denaturation in H. pylori. Finally, a new strategy is proposed to identify using COREX structural characteristics of equilibrium intermediate states populated during protein unfolding.

Molecular dynamics simulations of the native and partially folded states of ubiquitin: influence of methanol cosolvent, pH, and temperature on the protein structure and dynamics

A series of explicit-solvent molecular dynamics simulations of the protein ubiquitin are reported, which investigate the effect of environmental factors (presence of methanol cosolvent in the aqueous solution, neutral or low pH value, room or elevated temperature) on the structure, stability, and dynamics of the protein. The simulations are initiated either from the native structure of the protein or from a model of a partially folded state (A-state) that is known to exist at low pH in methanol-water mixtures. The main results of the simulations are: (1) The ubiquitin native structure is remarkably stable at neutral pH in water; (2) the addition of the methanol cosolvent enhances the stability of the secondary structure but weakens tertiary interactions within the protein; (3) this influence of methanol on the protein structure is enhanced at low pH, while the effect of lowering the pH in pure water is limited; and (4) the A-state of ubiquitin can be described as a set of relatively rigid secondary structure elements (a native-like beta-sheet and native-like alpha-helix plus two nonnative alpha-helices) connected by flexible linkers.