Strategy for Stabilization of CutA1 Proteins Due to Ion-Ion Interactions at Temperatures of over 100 °C.

In order to elucidate the contribution of charged residues to protein stabilization at temperatures of over 100 °C, we constructed many mutants of the CutA1 protein ( EcCutA1) from Escherichia coli. The goal was to see if one can achieve the same stability as for a CutA1 from hyperthermophile Pyrococcus horikoshii that has the denaturation temperature near 150 °C. The hydrophobic mutant of EcCutA1 ( Ec0VV) with denaturation temperature ( Td) of 113.2 °C was used as a template for mutations. The highest Td of Ec0VV mutants substituted by a single charged residue was 118.4 °C. Multiple ion mutants were also constructed by combination of single mutants and found to have an increased thermostability. The highest stability of multiple mutants was a mutant substituted by nine charged residues that had a Td of 142.2 °C. To evaluate the energy of ion-ion interactions of mutant proteins, we used the structural ensemble obtained by a molecular dynamics simulation at 300 K. The Td of ionic mutants linearly increases with the increments of the computed energy of ion-ion interactions for ionic mutant proteins even up to the temperatures near 140 °C, suggesting that ion-ion interactions cumulatively contribute to the stabilization of a protein at high temperatures.

Human ER Oxidoreductin-1α (Ero1α) Undergoes Dual Regulation through Complementary Redox Interactions with Protein-Disulfide Isomerase.

In the mammalian endoplasmic reticulum, oxidoreductin-1α (Ero1α) generates protein disulfide bonds and transfers them specifically to canonical protein-disulfide isomerase (PDI) to sustain oxidative protein folding. This oxidative process is coupled to the reduction of O2 to H2O2 on the bound flavin adenine dinucleotide cofactor. Because excessive thiol oxidation and H2O2 generation cause cell death, Ero1α activity must be properly regulated. In addition to the four catalytic cysteines (Cys94, Cys99, Cys104, and Cys131) that are located in the flexible active site region, the Cys208-Cys241 pair located at the base of another flexible loop is necessary for Ero1α regulation, although the mechanistic basis is not fully understood. The present study revealed that the Cys208-Cys241 disulfide was reduced by PDI and other PDI family members during PDI oxidation. Differential scanning calorimetry and small angle X-ray scattering showed that mutation of Cys208 and Cys241 did not grossly affect the thermal stability or overall shape of Ero1α, suggesting that redox regulation of this cysteine pair serves a functional role. Moreover, the flexible loop flanked by Cys208 and Cys241 provides a platform for functional interaction with PDI, which in turn enhances the oxidative activity of Ero1α through reduction of the Cys208-Cys241 disulfide. We propose a mechanism of dual Ero1α regulation by dynamic redox interactions between PDI and the two Ero1α flexible loops that harbor the regulatory cysteines.

Inhibition of the functional interplay between endoplasmic reticulum (ER) oxidoreduclin-1α (Ero1α) and protein-disulfide isomerase (PDI) by the endocrine disruptor bisphenol A.

Bisphenol A (BPA) is an endocrine disruptor that may have adverse effects on human health. We recently isolated protein-disulfide isomerase (PDI) as a BPA-binding protein from rat brain homogenates and found that BPA markedly inhibited PDI activity. To elucidate mechanisms of this inhibition, detailed structural, biophysical, and functional analyses of PDI were performed in the presence of BPA. BPA binding to PDI induced significant rearrangement of the N-terminal thioredoxin domain of PDI, resulting in more compact overall structure. This conformational change led to closure of the substrate-binding pocket in b' domain, preventing PDI from binding to unfolded proteins. The b' domain also plays an essential role in the interplay between PDI and ER oxidoreduclin 1α (Ero1α), a flavoenzyme responsible for reoxidation of PDI. We show that BPA inhibited Ero1α-catalyzed PDI oxidation presumably by inhibiting the interaction between the b' domain of PDI and Ero1α; the phenol groups of BPA probably compete with a highly conserved tryptophan residue, located in the protruding ß-hairpin of Ero1α, for binding to PDI. Consistently, BPA slowed down the reoxidation of PDI and caused the reduction of PDI in HeLa cells, indicating that BPA has a great impact on the redox homeostasis of PDI within cells. However, BPA had no effect on the interaction between PDI and peroxiredoxin-4 (Prx4), another PDI family oxidase, suggesting that the interaction between Prx4 and PDI is different from that of Ero1α and PDI. These results indicate that BPA, a widely distributed and potentially harmful chemical, inhibits Ero1-PDI-mediated disulfide bond formation.

Structure and function of a TetR family transcriptional regulator, SbtR, from thermus thermophilus HB8.

SbtR is one of the four TetR family transcriptional regulators present in the extremely thermophilic bacterium, Thermus thermophilus HB8. We identified 10 genes controlled by four promoters with negative regulation by SbtR in vitro. The SbtR-regulated gene products include probable transporters, probable enzymes for sugar or amino acid metabolism, and nucleic acid-related enzymes. SbtR binds pseudopalindromic sequences, with the consensus sequence of 5'-TGACCCNNKGGTCA-3' surrounding the promoters, and has a proposed 1:1 dimer binding stoichiometry. The X-ray crystal structure analysis revealed that SbtR comprises either nine or 10 α-helices and forms a dimer, as in the typical TetR family proteins. Similar to many characterized TetR family regulators, SbtR has a predicted ligand-binding pocket at the center of each monomer. Interestingly, the SbtR dimer contains an intermolecular disulfide bridge, formed between the Cys164 residues at the entrance of the pocket. The Cys164Ser and Cys164Ala mutant SbtR proteins formed homodimers similar to that of the wild type, but their thermal stabilities were lower by about 8°C, indicating that the disulfide bridge contributes to the thermal stability of the protein. However, altered repression activity of the mutants was not observed in vitro. From these results, we propose that ligand-binding is essential for SbtR to disengage from DNA, in a similar manner to the other characterized TetR family regulators. The formation and reduction of the disulfide bond might function in controlling the ligand-binding affinity of this transcriptional regulator.

Site-specific relaxation of peptide bond planarity induced by electrically attracted proton/deuteron observed by neutron crystallography.

In structural biology, peptide bonds, fundamental linkages between hundreds of amino acids, of which a protein molecule is composed, have been commonly treated as a plane structure just as Linus Pauling et al. proposed. In this paper, a site-specific peptide bond relaxation mechanism by deuterons whose localization has been suggested by neutron crystallography is proposed. Such deuteron was observed as an arm of neutron scattering length density protruding from the carbonyl oxygen atoms in the main chain in the omit map drawn by neutron crystallography of human lysozyme. Our comprehensive study using x-ray and neutron diffraction and 15 N chemical shifts of individual amide nitrogen atoms within the same peptide bond strongly suggests the relaxation of the electronic resonance structure because of site-specific modulation by protons/deuterons localized on the electron orbital of the carbonyl oxygen. All experimental data used in this examination were obtained at room temperature, which is preferable for enzymatic activity. Such a close interaction between the electron resonance structure of a peptide bond and the exchangeable protons/deuterons well agreed with that observed in an intermediate state in an amide hydrolytic reaction simulated by the ab-initio calculation including water molecules.

Structural characterization of a trapped folding intermediate of pyrrolidone carboxyl peptidase from a hyperthermophile.

The refolding of cysteine-free pyrrolidone carboxyl peptidase (PCP-0SH) from a hyperthermophile is unusually slow. PCP-0SH is trapped in the denatured (D1) state at 4 °C and pH 2.3, which is different from the highly denatured state in the presence of concentrated denaturant. In order to elucidate the mechanism of the unusually slow folding, we investigated the structure of the D1 state using NMR techniques with amino acid selectively labeled PCP-0SH. The HSQC spectrum of the D1 state showed that most of the resonances arising from the 114-208 residues are broadened, indicating that conformations of the 114-208 residues are in intermediate exchange on the microsecond to millisecond time scale. Paramagnetic relaxation enhancement data indicated the lack of long-range interactions between the 1-113 and the 114-208 segments in the D1 state. Furthermore, proline scanning mutagenesis showed that the 114-208 segment in the D1 state forms a loosely packed hydrophobic core composed of α4- and α6-helices. From these findings, we conclude that the 114-208 segment of PCP-0SH folds into a stable compact structure with non-native helix-helix association in the D1 state. Therefore, in the folding process from the D1 state to the native state, the α4- and α6-helices become separated and the central ß-sheet is folded between these helices. That is, the non-native interaction between the α4- and α6-helices may be responsible for the unusually slow folding of PCP-0SH.

The catalytic efficiency (kcat/Km) of the class A beta-lactamase Toho-1 correlates with the thermal stability of its catalytic intermediate analog.

The extended-spectrum beta-lactamases are associated with antibiotic resistance. Toho-1 R274N/R276N, a Class A beta-lactamase of CTX-M-type, efficiently hydrolyzes first generationcephalosporins (for example, cephalothin), in addition to cefotaxime, a third generation cephalosporin. However, this enzyme only marginally hydrolyzes the third generation cephalosporin ceftazidime, and the monobactam aztreonam. The deacylation defectiveness of the mutant Toho-1 E166A/R274N/R276N, which lacks the deacylation activity, results in the accumulation of the complex of an acylated-enzyme intermediate analog. For drug design, it would be useful if a quantitative prediction of a catalytic property were available without the need of enzymatic measurements. Therefore, we examined whether there is a correlation between the thermal stability of a catalytic intermediate (analog) and its kinetic parameters. First we measured the hydrolytic kinetics of the 14 species of beta-lactam antibiotics by Toho-1 R274N/R276N, and also measured the thermal stability of the accumulated acyl-intermediates of Toho-1 E166A/R274N/R276 by differential scanning calorimetry. Here we report the correlation of these parameters. The logarithm of the catalytic efficiency for Toho-1 R274N/R276N, log(k(cat)/K(m)) exhibited the best linear correlation with T(m,) which is the heat-denaturation temperature midpoint of the corresponding acylated complex of Toho-1 E166A/R274N/R276N. The correlation coefficient was 0.947, indicating that a relationship exists between the kinetic parameters and the stability of the intermediates. The results demonstrate a new method for investigating the catalytic properties of enzymes against any substrates, and a new approach to designing enzymes.

Structure of indole-3-glycerol phosphate synthase from Thermus thermophilus HB8: implications for thermal stability.

The three-dimensional structure of indole-3-glycerol phosphate synthase (IGPS) from the thermophilic bacterium Thermus thermophilus HB8 (TtIGPS) has been determined at 1.8 Å resolution. The structure adopts a typical (ß/α)(8)-barrel fold with an additional N-terminal extension of 46 residues. A detailed comparison of the crystal structure of TtIGPS with available structures of IGPS from the archaeon Sulfolobus solfataricus (SsIGPS) and the bacteria Thermotoga maritima (TmIGPS) and Escherichia coli (EcIGPS) has been performed. Although the overall folds of the proteins are the same, there are differences in amino-acid composition, structural rigidity, ionic features and stability clusters which may account for the high thermostability of the hyperthermophilic (SsIGPS and TmIGPS) and thermophilic (TtIGPS) proteins when compared with the mesophilic EcIGPS. The thermostability of IGPS seems to be established mainly by favourable interactions of charged residues, salt bridges and the spatial distribution of relatively rigid clusters of extensively interacting residues.

Crystal structure of TTHA0061, an uncharacterized protein from Thermus thermophilus HB8, reveals a novel fold.

The crystal structure of an uncharacterized protein TTHA0061 from Thermus thermophilus HB8, was determined and refined to 1.8 A by a single wavelength anomalous dispersion (SAD) method. The structural analysis and comparison of TTHA0061 with other existing structures in the Protein Data Bank (PDB) revealed a novel fold, suggesting that this protein may belong to a translation initiation factor or ribosomal protein family. Differential scanning calorimetry analysis suggested that the thermostability of TTHA0061 increased at pH ranges of 5.8-6.2, perhaps due to the abundance of glutamic acid residues.

Structural basis of the initial binding of tRNA(Ile) lysidine synthetase TilS with ATP and L-lysine.

In the bacterial genetic-code system, the codon AUA is decoded as isoleucine by tRNA(Ile)(2) with the lysidine residue at the wobble position. Lysidine is derived from cytidine, with ATP and L-lysine, by tRNA(Ile) lysidine synthetase (TilS), which is an N-type ATP pyrophosphatase. In this study, we determined the crystal structure of Aquifex aeolicus TilS, complexed with ATP, Mg2+, and L-lysine, at 2.5 A resolution. The presence of the TilS-specific subdomain causes the active site to have two separate gateways, a large hole and a narrow tunnel on the opposite side. ATP is bound inside the hole, and L-lysine is bound at the entrance of the tunnel. The conserved Asp36 in the PP-motif coordinates Mg2+. In these initial binding modes, the ATP, Mg2+, and L-lysine are held far apart from each other, but they seem to be brought together for the reaction upon cytidine binding, with putative structural changes of the complex.

Confirmation of the formation of salt bridges in the denatured state of CutA1 protein using molecular dynamics simulations.

It remains unclear how the abundant charged residues in proteins from hyperthermophiles contribute to the stabilization of proteins. Previously, based on molecular dynamics (MD) simulations, we proposed that these charged residues decrease the entropic effect by forming salt bridges in the denatured state under physiological conditions (Yutani et al., Sci. Rep. 8, 7613 (2018)). Because the quality of MD results is strongly dependent on the force fields used, in this study we performed the MD simulations using a different force field (AMBER99SB) along with the one we used before (Gromos43a1), at the same temperatures examined previously as well as at higher temperatures. In these experiments, we used the same ionic mutant (Ec0VV6) of CutA1 from Escherichia coli as in the previous study. In MD simulations at 300 K, Lys87 and Arg88 in the loop region of Ec0VV6 formed salt bridges with different favorable pairs in different force fields. Furthermore, the helical content and radius of gyration differed slightly between two force fields. However, at a higher temperature (600 K), the average numbers of salt bridges for the six substituted residues of Ec0VV6 were 0.87 per residue for Gromos43a1 and 0.88 for AMBER99SB in 400-ns MD simulation, indicating that the values were similar despite the use of different force fields. These observations suggest that the charged residues in Ec0VV6 can form a considerable number of salt bridges, even in the denatured state with drastic fluctuation at 600 K. These results corroborate our previous proposal.

Evaluating the strengths of salt bridges in the CutA1 protein using molecular dynamic simulations: a comparison of different force fields.

Ion-ion interactions (salt bridges) between favorable pairs of charged residues are important for the conformational stability of proteins. Molecular dynamic (MD) simulations are useful for elucidating the interactions among charged residues fluctuating in solution. However, the quality of MD results depends strongly on the force fields used. In this study, we compared the strengths of salt bridges among force fields by performing MD simulations using the CutA1 protein (trimer) from the hyperthermophile Pyrococcus horikoshii (PhCutA1), which has an unusually large proportion of charged residues. The force fields Chemistry at HARvard Macromolecular Mechanics (Charmm)27, Assisted Model Building and Energy Refinement (Amber)99sb, Amber14sb, GROningen Molecular Simulation (Gromos)43a1, and Gromos53a6 were used in combination with two different water models, tip3p (for Charmm27, Amber99sb, and Amber14sb) and simple point charge/extended (for Amber99sb, Gromos43a1, and Gromos53a6), yielding a total of six combinations. The RMSDs of all Cα atoms of PhCutA1 were similar among force fields, except for Charmm27, during 400-ns MD simulations at 300 K; however, the radius of gyration (Rg ) was greater for Amber99sb and shorter for Gromos43a1. The average strengths of salt bridges for each positively charged residue did not differ greatly among force fields, but the strengths at specific sites within the structure depended sensitively on the force field used. In the case of the Gromos group, positively charged residues could engage in favorable interactions with many more charged residues than in the other force fields, especially in loop regions; consequently, the apparent strength at each site was lower.

The confirmation of the denatured structure of pyrrolidone carboxyl peptidase under nondenaturing conditions: difference in helix propensity of two synthetic peptides with single amino acid substitution.

In the denatured state (D(1) state) of cystein-free pyrrolidone carboxyl peptidase (PCP-0SH) from Pyrococcus furiosus, a hyperthermophile under nondenaturing conditions, a fairly stable alpha-helix (alpha6-helix) has been determined from H/D exchange-NMR experiments. On the other hand, the alpha6-helix region of the proline-mutant at position 199 (A199P) was unstructured in the D(1) state unlike that of the wild-type PCP-0SH, although the folded conformations of both proteins were almost identical to each other. This finding has been deduced from the information regarding the remaining amide hydrogens in the HSQC spectra after H/D exchanges in the D(1) state. To confirm this inference, we examined the helical propensities of two synthetic peptides from their NMR structural analysis in the presence of trifluoroethanol (TFE). One is an 18-residue peptide called the wild-type H6-peptide corresponding to the alpha6-helix (from Ser188 to Glu205) of the wild-type PCP-0SH, and the other is the mutant H6-peptide corresponding to the alpha6-helix region of A199P. The NOE-contact information obtained from the 2D-(1)H-NOESY spectra measured for both peptides in the presence of 30% TFE clearly demonstrated that the wild-type H6-peptide had a high helical propensity, but the mutant H6-peptide was almost totally unstructured. The TFE-induced helical propensities for these peptide fragments confirmed the conclusions deduced from the H/D exchange data measured in the D(1) states of two proteins.

Structure of putative CutA1 from Homo sapiens determined at 2.05 A resolution.

The structure of human brain CutA1 (HsCutA1) has been determined using diffraction data to 2.05 A resolution. HsCutA1 has been implicated in the anchoring of acetylcholinesterase in neuronal cell membranes, while its bacterial homologue Escherichia coli CutA1 is involved in copper tolerance. Additionally, the structure of HsCutA1 bears similarity to that of the signal transduction protein PII, which is involved in regulation of nitrogen metabolism. Although several crystal structures of CutA1 from various sources with different rotation angles and degrees of interaction between trimer interfaces have been reported, the specific functional role of CutA1 is still unclear. In this study, the X-ray structure of HsCutA1 was determined in space group P2(1)2(1)2(1), with unit-cell parameters a = 68.69, b = 88.84, c = 125.33 A and six molecules per asymmetric unit. HsCutA1 is a trimeric molecule with intertwined antiparallel beta-strands; each subunit has a molecular weight of 14.6 kDa and contains 135 amino-acid residues. In order to obtain clues to the possible function of HsCutA1, its crystal structure was compared with those of other CutA1 and PII proteins.

Ion-ion interactions in the denatured state contribute to the stabilization of CutA1 proteins.

In order to elucidate features of the denatured state ensembles that exist in equilibrium with the native state under physiological conditions, we performed 1.4-µs molecular dynamics (MD) simulations at 400 K and 450 K using the monomer subunits of three CutA1 mutants from Escherichia coli: an SH-free mutant (Ec0SH) with denaturation temperature (Td) = 85.6 °C, a hydrophobic mutant (Ec0VV) with Td = 113.3 °C, and an ionic mutant (Ec0VV_6) with Td = 136.8 °C. The occupancy of salt bridges by the six substituted charged residues in Ec0VV_6 was 140.1% at 300 K and 89.5% at 450 K, indicating that even in the denatured state, salt bridge occupancy was high, approximately 60% of that at 300 K. From these results, we can infer that proteins from hyperthermophiles with a high ratio of charged residues are stabilized by a decrease in conformational entropy due to ion-ion interactions in the denatured state. The mechanism must be comparable to the stabilization conferred by disulfide bonds within a protein. This suggests that introduction of charged residues, to promote formation of salt bridges in the denatured state, would be a simple way to rationally design stability-enhanced mutants.

Hyper-thermostability of CutA1 protein, with a denaturation temperature of nearly 150 degrees C.

We found that the CutA1 protein, from Pyrococcus horikoshii (PhCutA1), has an extremely high denaturation temperature (T(d)) of nearly 150 degrees C, which exceeds the highest record determined by DSC by about 30 degrees C. To elucidate the mechanism of the ultra-high stability of PhCutA1, we analyzed the crystal structures of CutA1 proteins from three different sources, P. horikoshii, Thermus thermophilus, and Escherichia coli, with different growth temperatures (98, 75, and 37 degrees C). This analysis revealed that the remarkably increased number of ion pairs in the monomeric structure contributes to the stabilization of the trimeric structure and plays an important role in enhancing the T(d), up to 150 degrees C, for PhCutA1.

The unusually slow relaxation kinetics of the folding-unfolding of pyrrolidone carboxyl peptidase from a hyperthermophile, Pyrococcus furiosus.

In order to understand the thermodynamic and kinetic basis of the intrinsic stability of proteins from hyperthermophiles, the folding-unfolding reactions of cysteine-free pyrrolidone carboxyl peptidase (Cys142/188Ser) (PCP-0SH) from Pyrococcus furiosus were examined using circular dichroism (CD) and differential scanning calorimetry (DSC) at pH 2.3, where PCP-0SH exists in monomeric form. DSC showed a strong dependence of the shape and position of the unfolding profiles on the scan rate, suggesting the stability of PCP-0SH under kinetic control. On DSC timescales, even at a scan rate of 1 deg. C/hour, heat denaturation of PCP-0SH was non-equilibrium. However, over a long period of incubation of the heat-denatured PCP-0SH at pre-transition temperatures, it refolded completely, indicating reversibility with very slow relaxation kinetics. The rates of refolding of the heat-denatured PCP-0SH determined from the time-resolved DSC and CD spectroscopic progress curves were found to be similar within experimental error, confirming the mechanism of refolding to be a two-state process. The equilibrium established with a relaxation time of 5080 seconds (at t(m)=46.5 degrees C), which is unusually higher than the relaxation times observed for mesophilic and hyperthermophilic proteins. The long relaxation time may lead to the apparent irreversibility of an unfolding process occurring on the DSC experiment timescale. The refolding rate (9.8 x 10(-5) s(-1)) peaked near the t(m) (=46.5 degrees C), whereas the stability profile reached maxima (11.8 kJ mol(-1)) at 17 degrees C. The results clearly indicate the unusual mode of protein destabilization via a drastic decrease in the rate of folding at low pH and still maintaining a high activation energy barrier (284 kJ mol(-1)) for unfolding, which provides an effective kinetic advantage to unusually stable proteins from hyperthermophiles.

An Insight into the pathway of the amyloid fibril formation of hen egg white lysozyme obtained from a small-angle X-ray and neutron scattering study.

It is known that hen egg white lysozyme (HEWL) forms amyloid fibrils. Since HEWL is one of the proteins that have been studied most extensively and is closely related to human lysozyme, the variants of which form the amyloid fibrils that are related to hereditary systemic amyloidosis, this protein is an ideal model to study the mechanism of amyloid fibril formation. In order to gain an insight into the mechanism of amyloid fibril formation, systematic and detailed studies to detect and characterize various structural states of HEWL were conducted. Since HEWL forms amyloid fibrils in highly concentrated ethanol solutions, solutions of various concentrations of HEWL in various concentrations of ethanol were prepared, and the structures of HEWL in these solutions were investigated by small-angle X-ray and neutron scattering. It was shown that the structural states of HEWL were distinguished as the monomer state, the state of the dimer formation, the state of the protofilament formation, the protofilament state, and the state towards the formation of amyloid fibrils. A phase diagram of these structural states was obtained as a function of protein, water and ethanol concentrations. It was found that under the monomer state the structural changes of HEWL were not gross changes in shape but local conformational changes, and the dimers, formed by the association at the end of the long axis of HEWL, had an elongated shape. Circular dichroism measurements showed that the large changes in the secondary structures of HEWL occurred during dimer formation. The protofilaments were formed by stacking of the dimers with their long axis (nearly) perpendicular to and rotated around the protofilament axis to form a helical structure. These protofilaments were characterized by their radius of gyration of the cross-section of 2.4nm and the mass per unit length of 16,000(+/-2300)Da/nm. It was shown that the changes of the structural states towards the amyloid fibril formation occurred via lateral association of the protofilaments. A pathway of the amyloid fibril formation of HEWL was proposed from these results.