Predicted unfolding order of the 13 alpha-helices in the catalytic domain of glucoamylase from Aspergillus awamori var. X100 by molecular dynamics simulations

The unfolding mechanism of the 13 alpha-helices in the catalytic domain of Aspergillus awamori var. X100 glucoamylase was investigated by 200 ps molecular dynamics simulations in explicit water with temperature jump technique. Rather than a simultaneous event, the unfolding of these 13 alpha-helices followed a random ordered mechanism as alpha8-->alpha1-->alpha11-->alpha7-->alpha10-->alpha3-->alpha12-->alpha13-->alpha4-->alpha5-->alpha9-->alpha6-->alpha2. No significant relationships were found between the unfolding order and the length and the hydrophobicity of the helix. alpha-Helix 8 located in the inner region of the catalytic domain was predicted to be the first helix to unfold, indicating that the destruction of the secondary structure motif was initiated from the inner region of the catalytic domain. The dynamic behavior of these alpha-helices induced by increased kinetic energy during the unfolding process is considered to be similar to the expansion and compression of a series of springs under the influence of mechanical stress.     

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.     

Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a robust monomeric protein under kinetic control, which possesses some proline residues at the N-terminal of alpha-helices. Proline residue at the N-terminal of an alpha-helix is thought to stabilize a protein. In this work, the thermostability and folding kinetics of Tk-RNase HII were measured for mutant proteins in which a proline residue is introduced (Xaa to Pro) or removed (Pro to Ala) at the N-terminal of alpha-helices. In the folding experiments, the mutant proteins examined exhibit little influence on the remarkably slow unfolding of Tk-RNase HII. In contrast, E111P and K199P exhibit some thermostabilization, whereas P46A, P70A and P174A have some thermodestabilization. E111P/K199P and P46A/P70A double mutations cause cumulative changes in stability. We conclude that the proline effect on protein thermostability is observed in a hyperthermophilic protein, but each proline residue at the N-terminal of an alpha-helix slightly contributes to the thermostability. The present results also mean that even a natural hyperthermophilic protein can acquire improved thermostability.     

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.     

Effects of helix and fingertip mutations on the thermostability of xyn11A investigated by molecular dynamics simulations and enzyme activity assays

Local conformational changes and global unfolding pathways of wildtype xyn11A recombinant and its mutated structures were studied through a series of atomistic molecular dynamics (MD) simulations, along with enzyme activity assays at three incubation temperatures to investigate the effects of mutations at three different sites to the thermostability. The first mutation was to replace an unstable negatively charged residue at a surface beta turn near the active site (D32G) by a hydrophobic residue. The second mutation was to create a disulphide bond (S100C/N147C) establishing a strong connection between an alpha helix and a distal beta hairpin associated with the thermally sensitive Thumb loop, and the third mutation add an extra hydrogen bond (A155S) to the same alpha helix. From the MD simulations performed, MM/PBSA energy calculations of the unfolding energy were in a good agreement with the enzyme activities measured from the experiment, as all mutated structures demonstrated the improved thermostability, especially the S100C/N147C proved to be the most stable mutant both by the simulations and the experiment. Local conformational analysis at the catalytic sites and the xylan access region also suggested that mutated xyn11A structures could accommodate xylan binding. However, the analysis of global unfolding pathways showed that structural disruptions at the beta sheet regions near the N-terminal were still imminent. These findings could provide the insight on the molecular mechanisms underlying the enhanced thermostability due to mutagenesis and changes in the protein unfolding pathways for further protein engineering of the GH11 family xylanase enzymes.     

Thermal unfolding used as a probe to characterize the intra- and intersubunit stabilizing interactions in phosphorylating D-glyceraldehyde-3-phosphate dehydrogenase from Bacillus stearothermophilus

Tetrameric phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from Bacillus stearothermophilus can be described as a dimer of dimers with three nonequivalent interfaces. To investigate the contribution of intra- and intersubunit interactions to GAPDH thermostability, 10 residues located either at the cofactor domain (amino acids 1-148 and 313-333) or at the catalytic domain (amino acids 149-312) were mutated and the thermal unfolding of the mutants was studied by differential scanning calorimetry in the absence and presence of saturating concentrations of NAD. Disruptions of intrasubunit interactions lead to a drastic decrease in thermostability of the N313T, Y283V, and W310F mutants. Moreover, for the N313T mutant, a weakening of cooperative interactions between the catalytic and the cofactor domains and an inefficient binding of NAD are observed. This is likely the consequences of modification or loss of the hydrogen bonding network associating N313 and residues 236-238 and N313 and the nicotinamide carboxyamide of NAD, respectively. For the residues Y283 and W310, which are involved in stacking hydrophobic interactions, mutating both positions does not affect the efficiency of NAD binding. This shows that the factors involved in the thermostability of the tetrameric apo GAPDH are then different from those induced by NAD binding. Disruption of intersubunit hydrogen bonds between the catalytic domain and the NAD-binding domain of a neighboring subunit also leads to a significant destabilization of the apo tetrameric form as observed for the D282G mutant. Moreover, no efficient binding of NAD is observed. Both results are likely the consequence of a loss of hydrogen bonds across the P-axis and the Q-axis between D282 and R197 and between D282 and R52, respectively. Similar results, i.e., a destabilizing effect and inefficient NAD binding, are observed with the T34Q/T39S/L43Q mutant in which steric hindrance is introduced at the S-loop of the R-axis-related subunit via mutations at the adenosine subsite. The dimeric form of the D282G mutant exhibits a single partial heat absorption peak, whereas the Y46G/R52G mutant which exists only as a dimer shows two peaks. Taking into account the recent small-angle X-ray scattering studies which suggested that the dimeric form of the D282G mutant and of the dimeric Y46G/R52G mutant are of the O-R and O-P types, respectively (Vachette, unpublished results), we propose that the presence of one or two peaks in thermal unfolding of dimers is a signature of the dimer type.     

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.     

First crystallographic structure of a xylanase from glycoside hydrolase family 5: implications for catalysis

The room-temperature structure of xylanase (EC 3.2.1.8) from the bacterial plant pathogen Erwinia chrysanthemi expressed in Escherichia coli, a 45 kDa, 413-amino acid protein belonging to glycoside hydrolase family 5, has been determined by multiple isomorphous replacement and refined to a resolution of 1.42 A. This represents the first structure of a xylanase not belonging to either glycoside hydrolase family 10 or family 11. The enzyme is composed of two domains similar to most family 10 xylanases and the alpha-amylases. The catalytic domain (residues 46-315) has a (beta/alpha)(8)-barrel motif with a binding cleft along the C-terminal side of the beta-barrel. The catalytic residues, Glu165 and Glu253, determined by correspondence to other family 5 and family 10 glycoside hydrolases, lie inside this cleft on the C-terminal ends of beta-strands 4 and 7, respectively, with an O(epsilon)2...O(epsilon)1 distance of 4.22 A. The smaller domain (residues 31-43 and 323-413) has a beta(9)-barrel motif with five of the strands interfacing with alpha-helices 7 and 8 of the catalytic domain. The first 13 N-terminal residues form one beta-strand of this domain. Residues 44, 45, and 316-322 form the linkers between this domain and the catalytic domain.     

The crystallographic structure of Na,K-ATPase N-domain at 2.6A resolution

The structure of the N-domain of porcine alpha(2) Na,K-ATPase was determined crystallographically to 3.2A resolution by isomorphous heavy-atom replacement using a single mercury derivative. The structure was finally refined against 2.6A resolution synchrotron data. The domain forms a seven-stranded antiparallel beta-sheet with two additional beta-strands forming a hairpin and five alpha-helices. Approximately 75% of the residues were superimposable with residues from the structure of Ca-ATPase N-domain, and a structure-based sequence alignment is presented. The positions of key residues are discussed in relation to the pattern of hydrophobicity, charge and sequence conservation of the molecular surface. The structure of a hexahistidine tag binding to nickel ions is presented.     

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.     

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800?K, for 25?ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500?K has several conformations reaching into the protein core, whereas at 600?K some of the alpha-helices were lost. At 800?K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.     

A stabilizing alpha/beta-hydrophobic core greatly contributes to hyperthermostability of archaeal [P62A]Ssh10b

The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which has been postulated to play a role in chromosomal organization in Archaea. Ssh10b is capable of significantly constraining negative DNA supercoils at elevated temperatures. In this study, the solution structure of the dimeric P62A mutant Ssh10b ([P62A]Ssh10b) was determined by multidimensional NMR spectroscopy. The backbone 15N dynamics, H/D exchange with and without the denaturant GdmSCN, and chemical and thermal denaturation experiments were performed to investigate the molecular basis of high thermostability of [P62A]Ssh10b. Data analysis has revealed an alpha/beta-hydrophobic core consisting of two alpha-helices and one beta-sheet which are stabilized by cooperative hydrophobic and hydrogen-bonding interactions. This stabilizing alpha/beta-hydrophobic core of [P62A]Ssh10b exhibiting highly restricted internal motions is composed of residues having highly protected amide protons which exchange with solvent mostly by means of a global unfolding process. The K40N mutation greatly destabilizes the mutant [P62A]Ssh10b because this mutation disturbs the packing of alpha-helix against the beta-sheet reducing the stability of the alpha/beta-hydrophobic core in the mutant protein. In comparison with homologous mesophilic and thermophilic proteins, it can be presumed that the stabilizing alpha/beta-hydrophobic core in the [P62A]Ssh10b structure greatly contributes to the high thermostability of the protein.     

Mapping of the RNA recognition site of Escherichia coli ribosomal protein S7

Bacterial ribosomal protein S7 initiates the folding of the 3' major domain of 16S ribosomal RNA by binding to its lower half. The X-ray structure of protein S7 from thermophilic bacteria was recently solved and found to be a modular structure, consisting of an alpha-helical domain with a beta-ribbon extension. To gain further insights into its interaction with rRNA, we cloned the S7 gene from Escherichia coli K12 into a pET expression vector and introduced 4 deletions and 12 amino acid substitutions in the protein sequence. The binding of each mutant to the lower half of the 3' major domain of 16S rRNA was assessed by filtration on nitrocellulose membranes. Deletion of the N-terminal 17 residues or deletion of the B hairpins (residues 72-89) severely decreased S7 affinity for the rRNA. Truncation of the C-terminal portion (residues 138-178), which includes part of the terminal alpha-helix, significantly affected S7 binding, whereas a shorter truncation (residues 148-178) only marginally influenced its binding. Severe effects were also observed with several strategic point mutations located throughout the protein, including Q8A and F17G in the N-terminal region, and K35Q, G54S, K113Q, and M115G in loops connecting the alpha-helices. Our results are consistent with the occurrence of several sites of contact between S7 and the 16S rRNA, in line with its role in the folding of the 3' major domain.     

G alpha12 interaction with alphaSNAP induces VE-cadherin localization at endothelial junctions and regulates barrier function

The involvement of heterotrimeric G proteins in the regulation of adherens junction function is unclear. We identified alphaSNAP as an interactive partner of G alpha12 using yeast two-hybrid screening. Glutathione S-transferase pull-down assays showed the selective interaction of alphaSNAP with G alpha12 in COS-7 as well as in human umbilical vein endothelial cells. Using domain swapping experiments, we demonstrated that the N-terminal region of G alpha12 (1-37 amino acids) was necessary and sufficient for its interaction with alphaSNAP. G alpha13 with its N-terminal extension replaced by that of G alpha12 acquired the ability to bind to alphaSNAP, whereas G alpha12 with its N terminus replaced by that of G alpha13 lost this ability. Using four point mutants of alphaSNAP, which alter its ability to bind to the SNARE complex, we determined that the convex rather than the concave surface of alphaSNAP was involved in its interaction with G alpha12. Co-transfection of human umbilical vein endothelial cells with G alpha12 and alphaSNAP stabilized VE-cadherin at the plasma membrane, whereas down-regulation of alphaSNAP with siRNA resulted in the loss of VE-cadherin from the cell surface and, when used in conjunction with G alpha12 overexpression, decreased endothelial barrier function. Our results demonstrate a direct link between the alpha subunit of G12 and alphaSNAP, an essential component of the membrane fusion machinery, and implicate a role for this interaction in regulating the membrane localization of VE-cadherin and endothelial barrier function.     

NMR hydrogen exchange of the OB-fold protein LysN as a function of denaturant: the most conserved elements of structure are the most stable to unfolding.

The structure of LysN contains an OB-fold motif composed of a structurally conserved five-stranded beta-barrel capped by a poorly conserved alpha-helix between strands beta3 and beta4. Two additional alpha-helices, unique to the LysN structure, flank the N terminus of the OB-fold. The stability of LysN to unfolding has been investigated with NMR native state hydrogen exchange measurements as a function of guanidinium hydrochloride concentration, and equilibrium unfolding transitions monitored by ellipticity at 222 nm and fluorescence at 350 nm. The spectrophotometric measurements suggest an apparent two-state unfolding transition with DeltaGu(0) approximately 6 kcal/mol and m approximately 3 kcal/(molM). By contrast, NMR hydrogen exchange measurements manifest a distribution of DeltaGu(0) and m values which indicate that the protein can undergo subglobal unfolding. The largest DeltaGu(0) values from hydrogen exchange are for residues in the beta-sheet of the protein. These values, which reflect complete unfolding of the protein, are between 3 and 4 kcal/mol higher than those obtained from circular dichroism or fluorescence. This discrepancy may be due to the comparison of NMR hydrogen exchange parameters measured at residue-level resolution, with spectrophotometric parameters that reflect an unresolved super position of unfolding transitions of the alpha-helices and beta-strands. The largest DeltaGu(0) values obtained from hydrogen exchange for the subset of residues in the alpha-helices of the protein, agree with the DeltaGu(0) values obtained from circular dichroism or fluorescence. Based on the hydrogen exchange data, however, the three alpha-helices of LysN are on average 3 kcal/mol less stable than the beta-sheet. Consistent with the subglobal unfolding of LysN evinced by hydrogen exchange, a deletion mutant that lacks the first alpha-helix of the protein retains a cooperatively folded structure. Taken together with previous results on the OB-fold proteins SN and CspA, the present results for LysN suggest that the most conserved elements of structure in the OB-fold motif are the most resistant to denaturation. In all three proteins, stability to denaturation correlates with sequence hydrophobicity.     

Structure of a hyperthermophilic archaeal homing endonuclease, I-Tsp061I: contribution of cross-domain polar networks to thermostability

A novel LAGLIDADG-type homing endonuclease (HEase), I-Tsp061I, from the hyperthermophilic archaeon Thermoproteus sp. IC-061 16 S rRNA gene (rDNA) intron was characterized with respect to its structure, catalytic properties and thermostability. It was found that I-Tsp061I is a HEase isoschizomer of the previously described I-PogI and exhibits the highest thermostability among the known LAGLIDADG-type HEases. Determination of the crystal structure of I-Tsp061I at 2.1 A resolution using the multiple isomorphous replacement and anomalous scattering method revealed that the overall fold is similar to that of other known LAGLIDADG-type HEases, despite little sequence similarity between I-Tsp061I and those HEases. However, I-Tsp061I contains important cross-domain polar networks, unlike its mesophilic counterparts. Notably, the polar network Tyr6-Asp104-His180-107O-HOH12-104O-Asn177 exists across the two packed alpha-helices containing both the LAGLIDADG catalytic motif and the GxxxG hydrophobic helix bundle motif. Another important structural feature is the salt-bridge network Asp29-Arg31-Glu182 across N and C-terminal domain interface, which appears to contribute to the stability of the domain/domain packing. On the basis of these structural analyses and extensive mutational studies, we conclude that such cross-domain polar networks play key roles in stabilizing the catalytic center and domain packing, and underlie the hyperthermostability of I-Tsp061I.     

Structure and function of the mitochondrial bc1 complex. A mutational analysis of the yeast Rieske iron-sulfur protein

Respiratory-defective mutants of Saccharomyces cerevisiae assigned to a single complementation group (G12) have been determined to have lesions in the iron-sulfur protein (Rieske protein) of ubiquinol: cytochrome c reductase. Mutants capable of expressing the protein were chosen for further studies. The genes from 13 independent isolates were cloned and their mutations sequenced. Twelve mutations were ascertained to cause single amino acid substitutions in the carboxyl-terminal regions of the protein between residues 127 and 173. This region is proposed to be part of the catalytic domain with the ligands responsible for co-ordinating the two irons of the 2Fe-2S cluster. Based on the catalytic properties of the ubiquinol: cytochrome c reductase complex and the electron paramagnetic resonance (e.p.r.) signals of the iron-sulfur protein, the mutants describe two different phenotypes. A subset of mutants have no detectable iron-sulfur cluster and are completely deficient in ubiquinol: cytochrome c reductase activity. These strains identify mutations in residues considered to be essential for binding of the iron or for maintaining a proper tertiary structure of the catalytic domain. A second group of mutants have reduced levels of enzymatic activity and exhibit e.p.r. spectra characteristic of the Rieske iron-sulfur cluster. The mutations in the latter strains have been ascribed to residues that influence the redox properties of the cluster by distorting the iron-binding pocket. A secondary and tertiary structure model is presented of the carboxyl-terminal 65 residues constituting the catalytic domain of the iron-sulfur protein. It is postulated that the two irons of the cluster are co-ordinated by three cysteine and a single histidine residue located in a loop structure. The catalytic domain also contains two short alpha-helices and three beta-strands that form a partial beta-barrel. Most of the hydrophilic amino acids are present in turns that map to one pole of the domain. When viewed in the context of the model, mutations that abolish the iron-sulfur cluster are mostly in residues defining the boundaries of the alpha-helices and beta-strands. The notable exception is a cysteine residue that has been assigned to the loop with the iron ligands. This cysteine residue is proposed to co-ordinate one iron of the cluster. Mutations that reduce ubiquinol: cytochrome c reductase activity and alter the redox potential of the cluster occur in residues located in the loop that contains the ligands of the cluster.     

Crystal structure of D-Hydantoinase from Burkholderia pickettii at a resolution of 2.7 Angstroms: insights into the molecular basis of enzyme thermostability

D-Hydantoinase (D-HYD) is an industrial enzyme that is widely used in the production of D-amino acids which are precursors for semisynthesis of antibiotics, peptides, and pesticides. This report describes the crystal structure of D-hydantoinase from Burkholderia pickettii (HYD(Bp)) at a 2.7-A resolution. The structure of HYD(Bp) consists of a core (alpha/beta)(8) triose phosphate isomerase barrel fold and a beta-sheet domain, and the catalytic active site consists of two metal ions and six highly conserved amino acid residues. Although HYD(Bp) shares only moderate sequence similarity with D-HYDs from Thermus sp. (HYD(Tsp)) and Bacillus stearothermophilus (HYD(Bst)), whose structures have recently been solved, the overall structure and the structure of the catalytic active site are strikingly similar. Nevertheless, the amino acids that compose the substrate-binding site are less conserved and have different properties, which might dictate the substrate specificity. Structural comparison has revealed insights into the molecular basis of the differential thermostability of D-HYDs. The more thermostable HYD(Tsp) contains more aromatic residues in the interior of the structure than HYD(Bp) and HYD(Bst). Changes of large aromatic residues in HYD(Tsp) to smaller residues in HYD(Bp) or HYD(Bst) decrease the hydrophobicity and create cavities inside the structure. HYD(Tsp) has more salt bridges and hydrogen-bonding interactions and less oxidation susceptible Met and Cys residues on the protein surface than HYD(Bp) and HYD(Bst). Besides, HYD(Tsp) also contains more rigid Pro residues. These factors are likely to make major contributions to the varying thermostability of these enzymes. This information could be exploited in helping to engineer more thermostable mesophilic enzymes.     

Chimeric G alpha i2/G alpha 13 proteins reveal the structural requirements for the binding and activation of the RGS-like (RGL)-containing Rho guanine nucleotide exchange factors (GEFs) by G alpha 13

The alpha-subunit of G proteins of the G(12/13) family stimulate Rho by their direct binding to the RGS-like (RGL) domain of a family of Rho guanine nucleotide exchange factors (RGL-RhoGEFs) that includes PDZ-RhoGEF (PRG), p115RhoGEF, and LARG, thereby regulating cellular functions as diverse as shape and movement, gene expression, and normal and aberrant cell growth. The structural features determining the ability of G alpha(12/13) to bind RGL domains and the mechanism by which this association results in the activation of RGL-RhoGEFs are still poorly understood. Here, we explored the structural requirements for the functional interaction between G alpha(13) and RGL-RhoGEFs based on the structure of RGL domains and their similarity with the area by which RGS4 binds the switch region of G alpha(i) proteins. Using G alpha(i2), which does not bind RGL domains, as the backbone in which G alpha(13) sequences were swapped or mutated, we observed that the switch region of G alpha(13) is strictly necessary to bind PRG, and specific residues were identified that are critical for this association, likely by contributing to the binding surface. Surprisingly, the switch region of G alpha(13) was not sufficient to bind RGL domains, but instead most of its GTPase domain is required. Furthermore, membrane localization of G alpha(13) and chimeric G alpha(i2) proteins was also necessary for Rho activation. These findings revealed the structural features by which G alpha(13) interacts with RGL domains and suggest that molecular interactions occurring at the level of the plasma membrane are required for the functional activation of the RGL-containing family of RhoGEFs.     

Locating the stabilizing residues in (alpha/beta)8 barrel proteins based on hydrophobicity, long-range interactions, and sequence conservation

In nature, 1 out of every 10 proteins has an (alpha/beta)(8) (TIM)-barrel fold, and in most cases, pairwise comparisons show no sequence similarity between them. Hence, delineating the key residues that induce very different sequences to share a common fold is important for understanding the folding and stability of TIM-barrel domains. In this work, we propose a new consensus approach for locating these stabilizing residues based on long-range interactions, hydrophobicity, and conservation of amino acid residues. We have identified 957 stabilizing residues in 63 proteins from a nonredundant set of 71 TIM-barrel domains. Most of these residues are located in the 8-stranded beta-sheet, with nearly one half of them oriented toward the interior of the barrel and the other half oriented toward the surrounding alpha-helices. Several stabilizing residues are found in the N- and C-terminal loops, whereas very few appear in the alpha-helices that surround the internal beta-sheet. Further, these 957 residues are placed in 434 stabilizing segments of various sizes, and each domain contains 1-10 of these segments. We found that 8 segments per domain is the most abundant one, and two thirds of the proteins have 7-9 stabilizing segments. Finally, we verified the identified residues with experimental temperature factors and found that these residues are among the ones with less mobility in the considered proteins. We suggest that our new protocol serves as a powerful tool to identify the stabilizing residues in TIM-barrel domains, which can be used as potential candidates for studying protein folding and stability by means of protein engineering experiments.