Protein thermal denaturation is modulated by central residues in the protein structure network

Trikafta rescues F508del-CFTR by tightening specific phosphorylation-dependent interdomain interactions

Trikafta is well-known for correcting thermal and gating defects caused by the most common cystic fibrosis mutation F508del in the human cystic fibrosis transmembrane conductance regulator even at a physiological temperature. However, the exact correction pathway is still unclear. Here, noncovalent interactions among two transmembrane domains (TMD1 and TMD2), the regulatory (R) domain and two nucleotide binding domains (NBD1 and NBD2) were analyzed. The thermal stability of NBD1 was also evaluated through its tertiary constrained noncovalent interaction networks or thermoring structures. The results demonstrated that Trikafta binding to flexible TMD1 and TMD2 rearranged their interactions with the R domain upon phosphorylation, coupling tightened cytoplasmic TMD1-TMD2 interactions to tightened Mg/ATP-dependent NBD1-NBD2 dimerization, which stabilized NBD1 above human body temperature. Overall, although the deletion of F508 induces the primary thermal defect in NBD1 and then the gating defect at the TMD1-TMD2 interface, Trikafta rescued them in a reverse manner allosterically. These mechanistic insights into the precise correction pathway of this misfolded channel facilitate optimizing cystic fibrosis treatment.

Tris inhibits a GH1 β-glucosidase by a linear mixed inhibition mechanism

Here we demonstrate that Tris (2-amino-2-(hydroxymethyl)-1,3-propanediol), largely used as a buffering agent, is a linear mixed inhibitor (Ki =  12 ±  2 mM and α =  3 ±  1) of the GH1 β-glucosidase from the insect Spodoptera frugiperda (Sfβgly). Such an inhibition mechanism implies the formation of a non-productive ESI complex involving Sfβgly, substrate, and Tris. In addition, Tris binding reduces by 3 fold the enzyme affinity for the substrate. Hence, at concentrations higher than the Ki, Tris can completely abolish Sfβgly activity, whereas even at lower concentrations the presence of Tris causes underestimation of β-glucosidase kinetic parameters (Km and kcat). In agreement with the inhibition mechanism, computational docking showed that Tris could bind to a pocket placed at the lateral of the active site opening in the Sfβgly-substrate complex, hence leading to the formation of an ESI complex. In agreement with the crystallographic data available, computational docking also showed that Tris may find binding spots in the interior of the active site of the Sfβgly and several GH1 β-glucosidases. Moreover, the variety of their active site shapes results in a multiplicity of binding profiles, foreseeing different inhibition mechanisms. Thus, Tris inhibition may affect other GH1 β-glucosidases. This remark should be taken into account in their study, highlighting the importance of the appropriate buffer for accurate enzyme characterization.

Tailoring industrial enzymes for thermostability and activity evolution by the machine learning-based iCASE strategy

The pursuit of obtaining enzymes with high activity and stability remains a grail in enzyme evolution due to the stability-activity trade-off. Here, we develop an isothermal compressibility-assisted dynamic squeezing index perturbation engineering (iCASE) strategy to construct hierarchical modular networks for enzymes of varying complexity. Molecular mechanism analysis elucidates that the peak of adaptive evolution is reached through a structural response mechanism among variants. Furthermore, this dynamic response predictive model using structure-based supervised machine learning is established to predict enzyme function and fitness, demonstrating robust performance across different datasets and reliable prediction for epistasis. The universality of the iCASE strategy is validated by four sorts of enzymes with different structures and catalytic types. This machine learning-based iCASE strategy provides guidance for future research on the fitness evolution of enzymes.

Characterization and synergistic activity of heterologously expressed microbial-derived endoglucanase and bifunctional cellulase on wheat straw

Cellulases are divided into endoglucanase, exoglucanase, and β-glucosidase based on their catalytic activity. Eight cellulases were recombinantly expressed in Escherichia coli BL21 to investigate their effects on the enzymatic hydrolysis of wheat straw. Among them, cellulase 2006 exhibited the highest endoglucanase activity (432.25 U/mg), while Bf1 displayed superior exoglucanase and β-glucosidase activities (577.46 and 1991.63 U/mg respectively). Bioinformatic and enzymatic analyses revealed that both cellulases displayed notable thermal and pH stability. The enzyme kinetics parameters revealed that Km values for cellulases 2006 and Bf1 were 15.98 and 14.19 mg/mL, respectively, with Vmax values of 20.78 and 16.38 μmol/min/mg. In a prokaryotic co-expression system, the mixed cellulase Bf2006 exhibited endoglucanase, exoglucanase, and β-glucosidase activities (130.78, 1406.36, and 1119.25 U/mg, respectively). The enzymatic hydrolysis assay revealed that these three cellulases acted on the cellulose macromolecules in wheat straw, increasing reducing sugar content and decreasing crystallinity. Endoglucanase 2006 acted on various organic compounds rich in phenols and aromatic heterocycles, while Bf1 primarily acted on compounds containing glucose units. Bf2006 significantly affected the content of lignin, neutral detergent fiber, acid detergent fiber, and the microstructure of wheat straw, with degradation products primarily consisting of disaccharides, oligosaccharides, polysaccharides, glycosides, and other carbohydrates. This study provides theoretical guidance for the production and application of mixed cellulase Bf2006, serving as a reference for the industrial use of different cellulase types.

Mechanism of imidazole inhibition of a GH1 β-glucosidase

Imidazole is largely employed in recombinant protein purification, including GH1 β-glucosidases, but its effect on the enzyme activity is rarely taken into consideration. Computational docking suggested that imidazole interacts with residues forming the active site of the GH1 β-glucosidase from Spodoptera frugiperda (Sfβgly). We confirmed this interaction by showing that imidazole reduces the activity of Sfβgly, which does not result from enzyme covalent modification or promotion of transglycosylation reactions. Instead, this inhibition occurs through a partial competitive mechanism. Imidazole binds to the Sfβgly active site, reducing the substrate affinity by about threefold, whereas the rate constant of product formation remains unchanged. The binding of imidazole within the active site was further confirmed by enzyme kinetic experiments in which imidazole and cellobiose competed to inhibit the hydrolysis of p-nitrophenyl β-glucoside. Finally, imidazole interaction in the active site was also demonstrated by showing that it hinders access of carbodiimide to the Sfβgly catalytic residues, protecting them from chemical inactivation. In conclusion, imidazole binds in the Sfβgly active site, generating a partial competitive inhibition. Considering that GH1 β-glucosidases share conserved active sites, this inhibition phenomenon is probably widespread among these enzymes, and this should be taken into account when considering the characterization of their recombinant forms.

Role of a high centrality residue in protein dynamics and thermal stability

Centralities determined from Residue Interaction Networks (RIN) in proteins have been used to predict aspects of their structure and dynamics. Here, we correlate the Eigenvector Centrality (Ec) with the rate constant for thermal denaturation (kden) of the HisF protein from Thermotoga maritima based on 12 single alanine substitution mutants. The molecular basis for this correlation was further explored by studying a mutant containing a replacement of a high Ec residue, Y182A, which displayed increased kden at 80 °C. The crystallographic structure of this mutant showed few changes, mostly in two flexible loops. The 1H-15N -HSQC showed only subtle changes of cross peak positions for residues located near the mutation site and scattered throughout the structure. However, the comparison of the RIN showed that Y182 is the vertex of a set of high centrality residues that spreads throughout the HisF structure, which is lacking in the mutant. Cross-correlation displacements of Cα calculated from a molecular dynamics simulation at different temperatures showed that the Y182A mutation reduced the correlated movements in the HisF structure above 70 °C. 1H-15N NMR chemical shift covariance using temperature as perturbation were consistent with these results. In conclusion the increase in temperature drives the structure of the mutant HisF-Y182A into a less connected state, richer in non-concerted motions, located predominantly in the C-terminal half of the protein where Y182 is placed. Conversely, wild-type HisF responds to increased temperature as a single unit. Hence the replacement of a high Ec residue alters the distribution of thermal energy through HisF structure.

Homodimerization of a glycoside hydrolase family GH1 β-glucosidase suggests distinct activity of enzyme different states

In this work, we investigated how activity and oligomeric state are related in a purified GH1 β-glucosidase from Spodoptera frugiperda (Sfβgly). Gel filtration chromatography coupled to a multiple angle light scattering detector allowed separation of the homodimer and monomer states and determination of the dimer dissociation constant (K_D ), which was in the micromolar range. Enzyme kinetic parameters showed that the dimer is on average 2.5-fold more active. Later, we evaluated the kinetics of homodimerization, scanning the changes in the Sfβgly intrinsic fluorescence over time when the dimer dissociates into the monomer after a large dilution. We described how the rate constant of monomerization (k_off ) is affected by temperature, revealing the enthalpic and entropic contributions to the process. We also evaluated how the rate constant (k_obs ) by which equilibrium is reached after dimer dilution behaves when varying the initial Sfβgly concentration. These data indicated that Sfβgly dimerizes through the conformational selection mechanism, in which the monomer undergoes a conformational exchange and then binds to a similar monomer, forming a more active homodimer. Finally, we noted that conformational selection reports and experiments usually rely on a ligand whose concentration is in excess, but for homodimerization, this approach does not hold. Hence, since our approach overcomes this limitation, this study not only is a new contribution to the comprehension of GH1 β-glucosidases, but it can also help to elucidate protein interaction pathways.

Combining Free Energy Simulations and NMR Chemical-Shift Perturbation To Identify Transient Cation-π Contacts in Proteins

Flexible protein regions containing cationic and aromatic side-chains exposed to solvent may form transient cation-π interactions with structural and functional roles. To evaluate their stability and identify important intramolecular cation-π contacts, a combination of free energy profiles estimated from umbrella sampling with molecular dynamics simulations and chemical shift perturbations (CSP) obtained from nuclear magnetic resonance (NMR) experiments is applied here to the complete catalytic domain of human phosphatase Cdc25B. This protein is a good model system for transient cation-π interactions as it contains only one Trp residue (W550) in the disordered C-terminal segment and a total of 17 Arg residues, many exposed to solvent. Eight putative Arg-Trp pairs were simulated here. Only R482 and R544 show bound profiles corresponding to important transient cation-π interactions, while the others have dissociative or almost flat profiles. These results are corroborated by CSP analysis of three Cdc25B point mutants (W550A, R482A, and R544A) disrupting cation-π contacts. The proposed validation of statistically representative molecular simulations by NMR spectroscopy could be applied to identify transient contacts of proteins in general but carefully, as NMR chemical shifts are sensitive to changes in both molecular contacts and conformational distributions.

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Optimum temperature may be a misleading parameter in enzyme characterization and application

The optimum temperature is commonly determined in enzyme characterization. A search in the PubMed database for “optimum temperature” and “enzymes” yielded more than 1,700 manuscripts reporting this parameter over the last five years. Here, we show that the optimum temperature is not a constant. The catalytic activity of the mesophylic β-glucosidase Sfβgly was determined at different temperatures using different assay times and enzyme concentrations. We observed that the optimum temperature for Sfβgly changed by 5°C simply by modifying the assay length, and it was inversely correlated with enzyme concentration. These observations rely on the fact that close to the melting temperature, thermal denaturation continuously decreases the active enzyme concentration as the assay progresses. Thus, as the denaturation rate increases with increasing temperature, the bell-shaped curves observed in “activity versus temperature plots” occur only if the enzyme is denatured at and above the optimum temperature, which was confirmed using the thermostable β-glucosidase bglTm. Thus, the optimum temperature hardly reflects an intrinsic enzyme property and is actually a mere consequence of the assay condition. Thus, adoption of the “optimum temperature” determined under bench conditions for large-scale uses, which differ in assay length and enzyme concentration, may result in lower yields and financial losses.

Allosteric binding sites in Rab11 for potential drug candidates

Rab11 is an important protein subfamily in the RabGTPase family. These proteins physiologically function as key regulators of intracellular membrane trafficking processes. Pathologically, Rab11 proteins are implicated in many diseases including cancers, neurodegenerative diseases and type 2 diabetes. Although they are medically important, no previous study has found Rab11 allosteric binding sites where potential drug candidates can bind to. In this study, by employing multiple clustering approaches integrating principal component analysis, independent component analysis and locally linear embedding, we performed structural analyses of Rab11 and identified eight representative structures. Using these representatives to perform binding site mapping and virtual screening, we identified two novel binding sites in Rab11 and small molecules that can preferentially bind to different conformations of these sites with high affinities. After identifying the binding sites and the residue interaction networks in the representatives, we computationally showed that these binding sites may allosterically regulate Rab11, as these sites communicate with switch 2 region that binds to GTP/GDP. These two allosteric binding sites in Rab11 are also similar to two allosteric pockets in Ras that we discovered previously.

Mutations close to a hub residue affect the distant active site of a GH1 β-glucosidase

The tertiary structure of proteins has been represented as a network, in which residues are nodes and their contacts are edges. Protein structure networks contain residues, called hubs or central, which are essential to form short connection pathways between any pair of nodes. Hence hub residues may effectively spread structural perturbations through the protein. To test whether modifications nearby to hub residues could affect the enzyme active site, mutations were introduced in the β-glycosidase Sfβgly (PDB-ID: 5CG0) directed to residues that form an α-helix (260–265) and a β-strand (335–337) close to one of its main hub residues, F251, which is approximately 14 Å from the Sfβgly active site. Replacement of residues A263 and A264, which side-chains project from the α-helix towards F251, decreased the rate of substrate hydrolysis. Mutation A263F was shown to weaken noncovalent interactions involved in transition state stabilization within the Sfβgly active site. Mutations placed on the opposite side of the same α-helix did not show these effects. Consistently, replacement of V336, which side-chain protrudes from a β-strand face towards F251, inactivated Sfβgly. Next to V336, mutation S337F also caused a decrease in noncovalent interactions involved in transition state stabilization. Therefore, we suggest that mutations A263F, A264F, V336F and S337F may directly perturb the position of the hub F251, which could propagate these perturbations into the Sfβgly active site through short connection pathways along the protein network.

Search for independent (β/α)4 subdomains in a (β/α)8 barrel β-glucosidase

Proteins that fold as (β/α)₈ barrels are thought to have evolved from half-barrels that underwent duplication and fusion events. The evidence is particularly clear for small barrels, which have almost identical halves. Additionally, computational calculations of the thermodynamic stability of these structures in the presence of denaturants have revealed that (β/α)₈ barrels contain two subunits or domains corresponding to half-barrels. Hence, within (β/α)₈ barrels, half-barrels are self-contained units. Here, we tested this hypothesis using β-glucosidase from the bacterium Thermotoga maritima (bglTm), which has a (β/α)₈ barrel structure. Mutations were introduced to disrupt the noncovalent contacts between its halves and reveal the presence of two domains within bglTm, thus resulting in the creation of mutants T1 (containing W12A and I217A mutations) and T2 (containing W12A, H195A, I217A and F404A mutations). Mutants T1 and T2 were properly folded, as indicated by their fluorescence spectra and enzyme kinetic parameters. T1 and wild-type bglTm were equally stable, as shown by the results of thermal inactivation, differential scanning fluorimetry and guanidine hydrochloride denaturation experiments. However, T2 showed a first-order inactivation at 80°C, a single melting temperature of 82°C and only one transition concentration (c₅₀) in 2.4 M guanidine hydrochloride. Additionally, T1 and T2 exhibited a cooperative denaturation process that followed a two-state model (m-values equal to 1.4 and 1.6 kcal/mol/M, respectively), similar to that of wild-type bglTm (1.2 kcal/mol/M). Hence, T1 and T2 each denatured as a single unit, although they contained different degrees of disruption between their halves. In conclusion, bglTm halves are equivalent in terms of their thermal and chemical stability; thus, their separate contributions to (β/α)₈ barrel unfolding cannot be disentangled.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Using the Amino Acid Network to Modulate the Hydrolytic Activity of β-Glycosidases

The active site residues in GH1 β-glycosidases are compartmentalized into 3 functional regions, involved in catalysis or binding of glycone and aglycone motifs from substrate. However, it still remains unclear how residues outside the active site modulate the enzymatic activity. To tackle this question, we solved the crystal structure of the GH1 β-glycosidase from Spodoptera frugiperda (Sfβgly) to systematically map its residue contact network and correlate effects of mutations within and outside the active site. External mutations neighbouring the functional residues involved in catalysis and glycone-binding are deleterious, whereas mutations neighbouring the aglycone-binding site are less detrimental or even beneficial. The large dataset of new and previously characterized Sfβgly mutants supports that external perturbations are coherently transmitted to active site residues possibly through contacts and specifically disturb functional regions they interact to, reproducing the effects observed for direct mutations of functional residues. This allowed us to suggest that positions related to the aglycone-binding site are preferential targets for introduction of mutations aiming to further improve the hydrolytic activity of β-glycosidases.

Molecular mechanism of carbon nanotube to activate Subtilisin Carlsberg in polar and non-polar organic media

In the work, we mainly used molecular dynamics (MD) simulation and protein structure network (PSN) to study subtilisin Carlsberg (SC) immobilized onto carbon nanotube (CNT) in water, acetonitrile and heptane solvents, in order to explore activation mechanism of enzymes in non-aqueous media. The result indicates that the affinity of SC with CNT follows the decreasing order of water > acetonitrile > heptane. The overall structure of SC and the catalytic triad display strong robustness to the change of environments, responsible for the activity retaining. However, the distances between two β-strands of substrate-binding pocket are significantly expanded by the immobilization in the increasing order of water < acetonitrile < heptane, contributing to the highest substrate-binding energy in heptane media. PSN analysis further reveals that the immobilization enhances structural communication paths to the substrate-binding pocket, leading to its larger change than the free-enzymes. Interestingly, the increase in the number of the pathways upon immobilization is not dependent on the absorbed extent but the desorbed one, indicating significant role of shifting process of experimental operations in influencing the functional region. In addition, some conserved and important hot-residues in the paths are identified, providing molecular information for functional modification.

Use of network model to explore dynamic and allosteric properties of three GPCR homodimers

We used an elastic network model and protein structure network to study three class A GPCR homodimers.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.