Protein engineering to improve the thermostability of glucoamylase from Aspergillus awamori based on molecular dynamics simulations

Twelve mutations were constructed to improve the thermostability of glucoamylase from Aspergillus awamori based on the results of molecular dynamics simulations. The thermal unfolding of the catalytic domain followed a putative hierarchical behavior. In addition, the unfolding of the 13 alpha-helices obeyed the random ordered mechanism, in which the alpha-helices 8, 1 and 11 unfolded more rapidly than the others. The catalytic center was well protected by the (alpha/alpha)(6)-barrel at simulation temperatures up to 600 K, whereas the catalytic base, E400, migrated from its original interior pocket to the surface of the catalytic domain by surmounting the hydrophobic barrier provided by alpha-helices 12 and 13 at 800 K. The disulfide bonds engineered to 'lock' the alpha-helix 11 on the surface of the catalytic domain dramatically increased the thermostability. Substituting G396 and G407 with Ala residues slightly increased the thermostability, whereas their specific activity and catalytic efficiency were reduced. This indicates that the introduced residues with higher hydrophobicity were favorable in the loop between alpha-helices 12 and 13, whereas they partially destroyed the hydrogen bond and salt linkage network in the catalytic center. Alpha-helices 12 and 13 can be stabilized by introducing residues with higher hydrophobicity, except for the H391M mutation.     

Solution structure of human BCL-w: modulation of ligand binding by the C-terminal helix

The structure of human BCL-w, an anti-apoptotic member of the BCL-2 family, was determined by triple-resonance NMR spectroscopy and molecular modeling. Introduction of a single amino acid substitution (P117V) significantly improved the quality of the NMR spectra obtained. The cytosolic domain of BCL-w consists of 8 alpha-helices, which adopt a fold similar to that of BCL-xL, BCL-2, and BAX proteins. Pairwise root meant square deviation values were less than 3 A for backbone atoms of structurally equivalent regions. Interestingly, the C-terminal helix alpha8 of BCL-w folds into the BH3-binding hydrophobic cleft of the protein, in a fashion similar to the C-terminal transmembrane helix of BAX. A peptide corresponding to the BH3 region of the pro-apoptotic protein, BID, could displace helix alpha8 from the BCL-w cleft, resulting in helix unfolding. Deletion of helix alpha8 increased binding affinities of BCL-w for BAK and BID BH3-peptides, indicating that this helix competes for peptide binding to the hydrophobic cleft. These results suggest that although the cytosolic domain of BCL-w exhibits an overall structure similar to that of BCL-xL and BCL-2, the unique organization of its C-terminal helix may modulate BCL-w interactions with pro-apoptotic binding partners.     

Molecular dynamics simulations of the unfolding mechanism of the catalytic domain from Aspergillus awamori var. X100 glucoamylase

In this study, 200 ps molecular dynamics simulations were conducted to investigate the unfolding mechanism of the catalytic domain of glucoamylase from Aspergillus awamori var. X100. The unfolding of this domain was suggested to follow a putative hierarchical manner, in which the heavily O-glycosylated belt region from residues T440 to A471 acted as the initiation site, followed by the alpha-helix secondary structure destruction, and then the collapse of the catalytic center pocket. The O-glycosylated belt region surrounded the surface of the catalytic domain in its native state at low temperature, whereas it was extended and is more suitable to be classified as part of the subsequent linker domain at high temperatures due to its high flexibility. The inner set helices of the (alpha/alpha)(6)-barrel seemed to exhibit higher helical content than the outer set ones at all temperatures examined. The distances between the C(alpha) of the three Cys residue pairs fluctuated rapidly at higher temperatures, indicating that these disulfide bonds have little effect on the structural stabilization. The melting temperature, at which the residual total helicity of the catalytic domain is 50%, is much lower than the critical temperature, at which the catalytic center pocket has lost its structural integrity.     

Replacement and deletion mutations in the catalytic domain and belt region of Aspergillus awamori glucoamylase to enhance thermostability

Three single-residue mutations, Asp71-->Asn, Gln409-->Pro and Gly447-->Ser, two long-to-short loop replacement mutations, Gly23-Ala24-Asp25-Gly26-Ala27-Trp28- Val29-Ser30-->Asn-Pro-Pro (23-30 replacement) and Asp297-Ser298-Glu299-Ala300-Val301-->Ala-G ly-Ala (297-301 replacement) and one deletion mutation removing Glu439, Thr440 and Ser441 (Delta439-441), all based on amino acid sequence alignments, were made to improve Aspergillus awamori glucoamylase thermostability. The first and second single-residue mutations were designed to introduce a potential N:-glycosylation site and to restrict backbone bond rotation, respectively, and therefore to decrease entropy during protein unfolding. The third single-residue mutation was made to decrease flexibility and increase O:-glycosylation in the already highly O:-glycosylated belt region that extends around the globular catalytic domain. The 23-30 replacement mutation was designed to eliminate a very thermolabile extended loop on the catalytic domain surface and to bring the remainder of this region closer to the rest of the catalytic domain, therefore preventing it from unfolding. The 297-301 replacement mutant GA was made to understand the function of the random coil region between alpha-helices 9 and 10. Delta439-441 was constructed to decrease belt flexibility. All six mutations increased glucoamylase thermostability without significantly changing enzyme kinetic properties, with the 23-30 replacement mutation increasing the activation free energy for thermoinactivation by about 4 kJ/mol, which leads to a 4 degrees C increase in operating temperature at constant thermostability.     

Influence of pH on NMR structure and stability of the human prion protein globular domain

The NMR structure of the globular domain of the human prion protein (hPrP) with residues 121-230 at pH 7.0 shows the same global fold as the previously published structure determined at pH 4.5. It contains three alpha-helices, comprising residues 144-156, 174-194, and 200-228, and a short anti-parallel beta-sheet, comprising residues 128-131 and 161-164. There are slight, strictly localized, conformational changes at neutral pH when compared with acidic solution conditions: helix alpha1 is elongated at the C-terminal end with residues 153-156 forming a 310-helix, and the population of helical structure in the C-terminal two turns of helix alpha 2 is increased. The protonation of His155 and His187 presumably contributes to these structural changes. Thermal unfolding monitored by far UV CD indicates that hPrP-(121-230) is significantly more stable at neutral pH. Measurements of amide proton protection factors map local differences in protein stability within residues 154-157 at the C-terminal end of helix alpha 1 and residues 161-164 of beta-strand 2. These two segments appear to form a separate domain that at acidic pH has a larger tendency to unfold than the overall protein structure. This domain could provide a "starting point" for pH-induced unfolding and thus may be implicated in endosomic PrPC to PrPSc conformational transition resulting in transmissible spongiform encephalopathies.     

Molecular Dynamics Simulations of the Unfolding Mechanism of the Catalytic Domain from Aspergillus awamori var. X100 Glucoamylase

Abstract

In this study, 200 ps molecular dynamics simulations were conducted to investigate the unfolding mechanism of the catalytic domain of glucoamylase from Aspergillus awamori var. X100. The unfolding of this domain was suggested to follow a putative hierarchical manner, in which the heavily O-glycosylated belt region from residues T440 to A471 acted as the initiation site, followed by the a-helix secondary structure destruction, and then the collapse of the catalytic center pocket. The O-glycosylated belt region surrounded the surface of the catalytic domain in its native state at low temperature, whereas it was extended and is more suitable to be classified as part of the subsequent linker domain at high temperatures due to its high flexibility. The inner set helices of the (α/α)₆-barrel seemed to exhibit higher helical content than the outer set ones at all temperatures examined. The distances between the C_α of the three Cys residue pairs fluctuated rapidly at higher temperatures, indicating that these disulfide bonds have little effect on the structural stabilization. The melting temperature, at which the residual total helicity of the catalytic domain is 50%, is much lower than the critical temperature, at which the catalytic center pocket has lost its structural integrity.     

Structural characterization of the N-terminal oligomerization domain of the bacterial chromatin-structuring protein, H-NS

The H-NS protein plays a key role in condensing DNA and modulating gene expression in bacterial nucleoids. The mechanism by which this is achieved is dependent, at least in part, on the oligomerization of the protein. H-NS consists of two distinct domains; the N-terminal domain responsible for protein oligomerization, and the C-terminal DNA binding domain, which are separated by a flexible linker region. We present a multidimensional NMR study of the amino-terminal 64 residues of H-NS (denoted H-NS1-64) from Salmonella typhimurium, which constitute the oligomerization domain. This domain exists as a homotrimer, which is predicted to be self-associated through a coiled-coil configuration. NMR spectra show an equivalent magnetic environment for each monomer indicating that the polypeptide chains are arranged in parallel with complete 3-fold symmetry. Despite the limited resonance dispersion, an almost complete backbone assignment for 1H(N), 1H(alpha), 15N, 13CO and 13C(alpha) NMR resonances was obtained using a suite of triple resonance experiments applied to uniformly 15N-, 13C/15N- and 2H/13C/15N-labelled H-NS1-64 samples. The secondary structure of H-NS1-64 has been identified on the basis of the analysis of 1H(alpha), 13C(alpha), 13Cbeta and 13CO chemical shifts, NH/solvent exchange rates, intra-chain H(N)-H(N) and medium-range nuclear Overhauser enhancements (NOEs). Within the context of the homotrimer, each H-NS1-64 protomer consists of three alpha-helices spanning residues 2-8, 12-20 and 22-53, respectively. A topological model is presented for the symmetric H-NS1-64 trimer based upon the combined analysis of the helical elements and the pattern of backbone amide group 15N nuclear relaxation rates within the context of axially asymmetric diffusion tensor. In this model, the longest of the three helices (helix 3, residues 22-53) forms a coiled-coil interface with the other chains in the homotrimer. The two shorter N-terminal helices fold back onto the outer surface of the coiled-coil core and potentially act to stabilise this configuration.     

Inherent asymmetry of the structure of F1-ATPase from bovine heart mitochondria at 6.5 A resolution.

ATP synthase, the assembly which makes ATP in mitochondria, chloroplasts and bacteria, uses transmembrane proton gradients generated by respiration or photosynthesis to drive the phosphorylation of ADP. Its membrane domain is joined by a slender stalk to a peripheral catalytic domain, F1-ATPase. This domain is made of five subunits with stoichiometries of 3 alpha: 3 beta: 1 gamma: 1 delta: 1 epsilon, and in bovine mitochondria has a molecular mass of 371,000. We have determined the 3-dimensional structure of bovine mitochondrial F1-ATPase to 6.5 A resolution by X-ray crystallography. It is an approximately spherical globule 110 A in diameter, on a 40 A stem which contains two alpha-helices in a coiled-coil. This stem is presumed to be part of the stalk that connects F1 with the membrane domain in the intact ATP synthase. A pit next to the stem penetrates approximately 35 A into the F1 particle. The stem and the pit are two examples of the many asymmetric features of the structure. The central element in the asymmetry is the longer of the two alpha-helices in the stem, which extends for 90 A through the centre of the assembly and emerges on top into a dimple 15 A deep. Features with threefold and sixfold symmetry, presumed to be parts of homologous alpha and beta subunits, are arranged around the central rod and pit, but the overall structure is asymmetric. The central helix provides a possible mechanism for transmission of conformational changes induced by the proton gradient from the stalk to the catalytic sites of the enzyme.     

Disruption of a calmodulin central helix-like region of 10-formyltetrahydrofolate dehydrogenase impairs its dehydrogenase activity by uncoupling the functional domains

10-Formyltetrahydrofolate dehydrogenase (FDH) is composed of three domains and possesses three catalytic activities but has only two catalytic centers. The amino-terminal domain (residue 1-310) bears 10-formyltetrahydrofolate hydrolase activity, the carboxyl-terminal domain (residue 420-902) bears an aldehyde dehydrogenase activity, and the full-length FDH produces 10-formyltetrahydrofolate dehydrogenase activity. The intermediate linker (residues 311-419) connecting the two catalytic domains does not contribute directly to the enzyme catalytic centers but is crucial for 10-formyltetrahydrofolate dehydrogenase activity. We have identified a region within the intermediate domain (residues 384-405) that shows sequence similarity to the central helix of calmodulin. Deletion of either the entire putative helix or the central part of the helix or replacement of the six residues within the central part with alanines resulted in total loss of the 10-formyltetrahydrofolate dehydrogenase activity, whereas the full hydrolase and aldehyde dehydrogenase activities were retained. Alanine-scanning mutagenesis revealed that neither of the six residues alone is required for FDH activity. Analysis of the predicted secondary structures and circular dichroic and fluorescence spectroscopy studies of the intermediate domain expressed as a separate protein showed that this region is likely to consist of two alpha-helices connected by a flexible loop. Our results suggest that flexibility within the putative helix is important for FDH function and could be a point for regulation of the enzyme.     

Solution structure and backbone dynamics of long-[Arg(3)]insulin-like growth factor-I

Long-[Arg(3)]insulin-like growth factor-I (IGF-I) is a potent analog of insulin-like growth factor-I that has been modified by a Glu(3) --> Arg mutation and a 13-amino acid extension appended to the N terminus. We have determined the solution structure of (15)N-labeled Long-[Arg(3)]-IGF-I using high resolution NMR and restrained molecular dynamics techniques to a precision of 0.82 +/- 0.28 A root mean square deviation for the backbone heavy atoms in the three alpha-helices and 3.5 +/- 0.9 A root mean square deviation for all backbone heavy atoms excluding the 8 N-terminal residues and the 8 C-terminal eight residues. Overall, the structure of the IGF-I domain is consistent with earlier studies of IGF-I with some minor changes remote from the N terminus. The major variations in the structure, compared with IGF-I, occur at the N terminus with a substantial reorientation of the N-terminal three residues of the IGF-I domain. These results are interpreted in terms of the lower binding affinity for insulin-like growth factor-binding proteins. The backbone dynamics of Long-[Arg(3)]IGF-I were investigated using (15)N nuclear spin relaxation and the heteronuclear nuclear Overhauser enhancement (NOE). There is a considerable degree of flexibility in Long-[Arg(3)]IGF-I, even in the alpha-helices, as indicated by an average ((1)H)(15)N NOE of 0.55 for the regions. The largest heteronuclear NOEs are observed in the helical regions, lower heteronuclear NOEs are observed in the C-domain loop separating helix 1 from helix 2, and negative heteronuclear NOEs are observed in the N-terminal extension and at the C terminus. Despite these data indicating conformational flexibility for the N-terminal extension, slow amide proton exchange was observed for some residues in this region, suggesting some transitory structure does exist, possibly a molten helix. A certain degree of flexibility may be necessary in all insulin-like growth factors to enable association with various receptors and binding proteins.     

Recombinant collagen studies link the severe conformational changes induced by osteogenesis imperfecta mutations to the disruption of a set of interchain salt bridges

The clinical severity of Osteogenesis Imperfecta (OI), also known as the brittle bone disease, relates to the extent of conformational changes in the collagen triple helix induced by Gly substitution mutations. The lingering question is why Gly substitutions at different locations of collagen cause different disruptions of the triple helix. Here, we describe markedly different conformational changes of the triple helix induced by two Gly substitution mutations placed only 12 residues apart. The effects of the Gly substitutions were characterized using a recombinant collagen fragment modeling the 63-residue segment of the alpha1 chain of type I collagen containing no Hyp (residues 877-939) obtained from Escherichia coli. Two Gly --> Ser substitutions at Gly-901 and Gly-913 associated with, respectively, mild and severe OI variants were introduced by site-directed mutagenesis. Biophysical characterization and limited protease digestion experiments revealed that while the substitution at Gly-901 causes relatively minor destabilization of the triple helix, the substitution at Gly-913 induces large scale unfolding of an unstable region C-terminal to the mutation site. This extensive unfolding is caused by the intrinsic low stability of the C-terminal region of the helix and the mutation induced disruption of a set of salt bridges, which functions to lock this unstable region into the triple helical conformation. The extensive conformational changes associated with the loss of the salt bridges highlight the long range impact of the local interactions of triple helix and suggest a new mechanism by which OI mutations cause severe conformational damages in collagen.     

Stabilities and activities of the N- and C-domains of FKBP22 from a psychrotrophic bacterium overproduced in Escherichia coli

FKBP22 from a psychrotrophic bacterium Shewanella sp. SIB1, is a dimeric protein with peptidyl prolyl cis-trans isomerase (PPIase) activity. According to homology modeling, it consists of an N-terminal domain, which is involved in dimerization of the protein, and a C-terminal catalytic domain. A long alpha3 helix spans these domains. An N-domain with the entire alpha3 helix (N-domain+) and a C-domain with the entire alpha3 helix (C-domain+) were overproduced in Escherichia coli in a His-tagged form, purified, and their biochemical properties were compared with those of the intact protein. C-domain+ was shown to be a monomer and enzymatically active. Its optimum temperature for activity (10 degrees C) was identical to that of the intact protein. Determination of the PPIase activity using peptide and protein substrates suggests that dimerization is required to make the protein fully active for the protein substrate or that the N-domain is involved in substrate-binding. The differential scanning calorimetry studies revealed two distinct heat absorption peaks at 32.5 degrees C and 46.6 degrees C for the intact protein, and single heat absorption peaks at 44.7 degrees C for N-domain+ and 35.6 degrees C for C-domain+. These results indicate that the thermal unfolding transitions of the intact protein at lower and higher temperatures represent those of C- and N-domains, respectively. Because the unfolding temperature of C-domain+ is much higher than its optimum temperature for activity, SIB1 FKBP22 may adapt to low temperatures by increasing a local flexibility around the active site. This study revealed the relationship between the stability and the activity of a psychrotrophic FKBP22.     

Src kinase conformational activation: thermodynamics, pathways, and mechanisms

Tyrosine kinases of the Src-family are large allosteric enzymes that play a key role in cellular signaling. Conversion of the kinase from an inactive to an active state is accompanied by substantial structural changes. Here, we construct a coarse-grained model of the catalytic domain incorporating experimental structures for the two stable states, and simulate the dynamics of conformational transitions in kinase activation. We explore the transition energy landscapes by constructing a structural network among clusters of conformations from the simulations. From the structural network, two major ensembles of pathways for the activation are identified. In the first transition pathway, we find a coordinated switching mechanism of interactions among the αC helix, the activation-loop, and the β strands in the N-lobe of the catalytic domain. In a second pathway, the conformational change is coupled to a partial unfolding of the N-lobe region of the catalytic domain. We also characterize the switching mechanism for the αC helix and the activation-loop in detail. Finally, we test the performance of a Markov model and its ability to account for the structural kinetics in the context of Src conformational changes. Taken together, these results provide a broad framework for understanding the main features of the conformational transition taking place upon Src activation.     

Thermal‐ and urea‐induced unfolding processes of glutathione S‐transferase by molecular dynamics simulation

The Schistosoma juponicum 26 kDa glutathione S‐transferase (sj26GST) consists of the N‐terminal domain (N‐domain), containing three alpha‐helices (named H1‐H3) and four anti‐parallel beta‐strands (S1‐S4), and the C‐terminal domain (C‐domain), comprising five alpha‐helices (named H4‐H8). In present work, molecular dynamics simulations and fluorescence spectroscopic were used to gain insights into the unfolding process of sj26GST. The molecular dynamics simulations on sj26GST subunit both in water and in 8 M urea were carried out at 300 K, 400 K and 500 K, respectively. Spectroscopic measurements were employed to monitor structural changes. Molecular dynamics simulations of sj26GST subunit induced by urea and temperature showed that the initial unfolding step of sj26GST both in water and urea occurred on N‐domain, involving the disruption of helices H2, H3 and strands S3 and S4, whereas H6 was the last region exposed to solution and was the last helix to unfold. Moreover, simulations analyses combining with fluorescence and circular dichroism spectra indicated that N‐domain could not fold independent, suggesting that correct folding of N‐domain depended on its interactions with C‐domain. We further proposed that the folding of GSTs could begin with the hydrophobic collapse of C‐domain whose H4, H5, H6 and H7 could move close to each other and form a hydrophobic core, especially H6 wrapped in the hydrophobic center and beginning spontaneous formation of the helix. S3, S4, H3, and H2 could form in the wake of the interaction between C‐domain and N‐domain. The paper can offer insights into the molecular mechanism of GSTs unfolding. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 247–259, 2015.     

Crystal structure of glucoamylase from Aspergillus awamori var. X100 to 2.2-A resolution.

The crystal structure of a catalytically active fragment of glucoamylase-I from Aspergillus awamori var. X100 has been determined to a resolution of 2.2 A. Twelve of its 13 alpha-helices are arranged into an "alpha/alpha-barrel." An inner core of six mutually parallel alpha-helices are connected to each other through a peripheral set of six alpha-helices. The peripheral helices are parallel to each other, but approximately antiparallel to the inner core of alpha-helices. The putative active site lies in the packing void of the inner set of helices. The last 30 residues of the enzyme comprise a separate domain containing 10 sites of O-glycosylation. Each instance of O-glycosylation involves a serine or threonine side chain linked to the alpha-anomer of a single mannosyl residue. The O-glycosylated domain is in an extended conformation, wrapping around the "waist" of the alpha/alpha-barrel. Two additional sites of N-glycosylation contribute well ordered glycosyl chains that lie in proximity to the belt of O-glycosylation. The model developed for glucoamylase is a rare and valuable structural example of a glycoprotein and an exo-acting amylolytic enzyme.     

Role of the C-terminal helix 9 in the stability and ligandin function of class alpha glutathione transferase A1-1

Helix 9 at the C-terminus of class alpha glutathione transferase (GST) polypeptides is a unique structural feature in the GST superfamily. It plays an important structural role in the catalytic cycle. Its contribution toward protein stability/folding as well as the binding of nonsubstrate ligands was investigated by protein engineering, conformational stability, enzyme activity, and ligand-binding methods. The helix9 sequence displays an unfavorable propensity toward helix formation, but tertiary interactions between the amphipathic helix and the GST seem to contribute sufficient stability to populate the helix on the surface of the protein. The helix's stability is enhanced further by the binding of ligands at the active site. The order of ligand-induced stabilization increases from H-site occupation, to G-site occupation, to the simultaneous occupation of H- and G-sites. Ligand-induced stabilization of helix9 reduces solvent accessible hydrophobic surface by facilitating firmer packing at the hydrophobic interface between helix and GST. This stabilized form exhibits enhanced affinity for the binding of nonsubstrate ligands to ligandin sites (i.e., noncatalytic binding sites). Although helix9 contributes very little toward the global stability of hGSTA1-1, its conformational dynamics have significant implications for the protein's equilibrium unfolding/refolding pathway and unfolding kinetics. Considering the high concentration of reduced glutathione in human cells (about 10 mM), the physiological form of hGSTA1-1 is most likely the thiol-complexed protein with a stabilized helix9. The C-terminus region (including helix9) of the class alpha polypeptide appears not to have been optimized for stability but rather for catalytic and ligandin function.     

Molecular dynamics simulations as a tool for improving protein stability

Haloalkane dehalogenase (DhlA) was used as a model protein to explore the possibility to use molecular dynamics (MD) simulations as a tool to identify flexible regions in proteins that can serve as a target for stability enhancement by introduction of a disulfide bond. DhlA consists of two domains: an alpha/beta-hydrolase fold main domain and a cap domain composed of five alpha-helices. MD simulations of DhlA showed high mobility in a helix-loop-helix region in the cap domain, involving residues 184-211. A disulfide cross-link was engineered between residue 201 of this flexible region and residue 16 of the main domain. The mutant enzyme showed substantial changes in both thermal and urea denaturation. The oxidized form of the mutant enzyme showed an increase of the apparent transition temperature from 47.5 to 52.5 degrees C, whereas the T(m,app) of the reduced mutant decreased by more than 8 degrees C compared to the wild-type enzyme. Urea denaturation results showed a similar trend. Measurement of the kinetic stability showed that the introduction of the disulfide bond caused a decrease in activation free energy of unfolding of 0.43 kcal mol(-1) compared to the wild-type enzyme and also indicated that the helix-loop-helix region was involved early in the unfolding process. The results show that MD simulations are capable of identifying mobile protein domains that can successfully be used as a target for stability enhancement by the introduction of a disulfide cross-link.