Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding.

The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.

The influence of pH on the structure and stability of the Grb2 dimer reveals changes in the inter-domain and molecular interaction: Could it be a modulation mechanism?

Cancer cells present an increased replicative potential as a hallmark. The increased replication leads to a higher intracellular pH. Grb2, an adapter protein, is mainly involved in several types of cancers due to its role in signaling pathways responsible for cell growth and proliferation. At pH 7, we observed a more compact structure, as seen by DLS and 1H NMR relaxation experiments, with high cooperativity within domains. On the other hand, we observed an increase in disordered structures at pH 8, with relative independence between domains characterized by higher melting temperatures and enthalpy of unfolding. CD and DLS corroborate with these observations at pH 8, conferring more flexibility among the domains, followed by lower unfolding cooperativity and increased hydrodynamic diameter at higher pH. In addition, 15N-HSQC chemical shift perturbations experiments showed significant differences in the positions of several amino acids spread on the Grb2 structure when pH was changed, which agrees with the previous results. Finally, the molecular dynamic analysis demonstrates that Grb2 presents a movement pattern where both SH3 domains move toward the center of the protein at pH 7. On the contrary, the pattern changes its direction at pH 8, where domains move outside the center of the protein, conferring a more elongated structure at higher pH. So, Grb2 presents significant structural and dynamic changes modulated by pH. If considering the role of Grb2 in cell signaling upstream, these conformational changes could be a critical mechanistic behavior of this protein, preventing/disrupting the stability of the cell signaling pathways related to cancer.

Information and redundancy in the burial folding code of globular proteins within a wide range of shapes and sizes.

Recent ab initio folding simulations for a limited number of small proteins have corroborated a previous suggestion that atomic burial information obtainable from sequence could be sufficient for tertiary structure determination when combined to sequence-independent geometrical constraints. Here, we use simulations parameterized by native burials to investigate the required amount of information in a diverse set of globular proteins comprising different structural classes and a wide size range. Burial information is provided by a potential term pushing each atom towards one among a small number L of equiprobable concentric layers. An upper bound for the required information is provided by the minimal number of layers L(min) still compatible with correct folding behavior. We obtain L(min) between 3 and 5 for seven small to medium proteins with 50 ≤ Nr ≤ 110 residues while for a larger protein with Nr = 141 we find that L ≥ 6 is required to maintain native stability. We additionally estimate the usable redundancy for a given L ≥ L(min) from the burial entropy associated to the largest folding-compatible fraction of "superfluous" atoms, for which the burial term can be turned off or target layers can be chosen randomly. The estimated redundancy for small proteins with L = 4 is close to 0.8. Our results are consistent with the above-average quality of burial predictions used in previous simulations and indicate that the fraction of approachable proteins could increase significantly with even a mild, plausible, improvement on sequence-dependent burial prediction or on sequence-independent constraints that augment the detectable redundancy during simulations.

Conformational changes in a hyperthermostable glycoside hydrolase: enzymatic activity is a consequence of the loop dynamics and protonation balance.

Endo-ß-1, 4-mannanase from Thermotoga petrophila (TpMan) is a modular hyperthermostable enzyme involved in the degradation of mannan-containing polysaccharides. The degradation of these polysaccharides represents a key step for several industrial applications. Here, as part of a continuing investigation of TpMan, the region corresponding to the GH5 domain (TpManGH5) was characterized as a function of pH and temperature. The results indicated that the enzymatic activity of the TpManGH5 is pH-dependent, with its optimum activity occurring at pH 6. At pH 8, the studies demonstrated that TpManGH5 is a molecule with a nearly spherical tightly packed core displaying negligible flexibility in solution, and with size and shape very similar to crystal structure. However, TpManGH5 experiences an increase in radius of gyration in acidic conditions suggesting expansion of the molecule. Furthermore, at acidic pH values, TpManGH5 showed a less globular shape, probably due to a loop region slightly more expanded and flexible in solution (residues Y88 to A105). In addition, molecular dynamics simulations indicated that conformational changes caused by pH variation did not change the core of the TpManGH5, which means that only the above mentioned loop region presents high degree of fluctuations. The results also suggested that conformational changes of the loop region may facilitate polysaccharide and enzyme interaction. Finally, at pH 6 the results indicated that TpManGH5 is slightly more flexible at 65°C when compared to the same enzyme at 20°C. The biophysical characterization presented here is well correlated with the enzymatic activity and provide new insight into the structural basis for the temperature and pH-dependent activity of the TpManGH5. Also, the data suggest a loop region that provides a starting point for a rational design of biotechnological desired features.

Oligomeric state and structural stability of two hyperthermophilic ß-glucosidases from Thermotoga petrophila.

The ß-glucosidases are enzymes essential for several industrial applications, especially in the field of plant structural polysaccharides conversion into bioenergy and bioproducts. In a recent study, we have provided a biochemical characterization of two hyperthermostable ß-glucosidases from Thermotoga petrophila belonging to the families GH1 (TpBGL1) and GH3 (TpBGL3). Here, as part of a continuing investigation, the oligomeric state, the net charge, and the structural stability, at acidic pH, of the TpBGL1 and TpBGL3 were characterized and compared. Enzymatic activity is directly related to the balance between protonation and conformational changes. Interestingly, our results indicated that there were no significant changes in the secondary, tertiary and quaternary structures of the ß-glucosidases at temperatures below 80 °C. Furthermore, the results indicated that both the enzymes are stable homodimers in solution. Therefore, the observed changes in the enzymatic activities are due to variations in pH that modify protonation of the enzymes residues and the net charge, directly affecting the interactions with ligands. Finally, the results showed that the two ß-glucosidases displayed different pH dependence of thermostability at temperatures above 80 °C. TpBGL1 showed higher stability at pH 6 than at pH 4, while TpBGL3 showed similar stability at both pH values. This study provides a useful comparison of the structural stability, at acidic pH, of two different hyperthermostable ß-glucosidases and how it correlates with the activity of the enzymes. The information described here can be useful for biotechnological applications in the biofuel and food industries.