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.

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.

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.

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.

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.

Thermodynamics of protein denaturation at temperatures over 100 °C: CutA1 mutant proteins substituted with hydrophobic and charged residues.

Although the thermodynamics of protein denaturation at temperatures over 100 °C is essential for the rational design of highly stable proteins, it is not understood well because of the associated technical difficulties. We designed certain hydrophobic mutant proteins of CutA1 from Escherichia coli, which have denaturation temperatures (Td) ranging from 101 to 113 °C and show a reversible heat denaturation. Using a hydrophobic mutant as a template, we successfully designed a hyperthermostable mutant protein (Td = 137 °C) by substituting six residues with charged ones. Thermodynamic analyses of these mutant proteins indicated that the hydrophobic mutants were stabilized by the accumulation of denaturation enthalpy (ΔH) with no entropic gain from hydrophobic solvation around 100 °C, and that the stabilization due to salt bridges resulted from both the increase in ΔH from ion-ion interactions and the entropic effect of the electrostatic solvation over 113 °C. This is the first experimental evidence that has successfully overcome the typical technical difficulties.

Thermodynamic analysis of unusually thermostable CutA1 protein from human brain and its protease susceptibility.

Unusually stable proteins are a disadvantage for the metabolic turnover of proteins in cells. The CutA1 proteins from Pyrococcus horikoshii and from Oryza sativa (OsCutA1) have unusually high denaturation temperatures (Td) of nearly 150 and 100 °C, respectively, at pH 7.0. It seemed that the CutA1 protein from the human brain (HsCutA1) also has a remarkably high stability. Therefore, the thermodynamic stabilities of HsCutA1 and its protease susceptibility were examined. The Td was remarkably high, being over 95 °C at pH 7.0. The unfolding Gibbs energy (ΔG(0)H2O) was 174 kJ/mol at 37 °C from the denaturant denaturation. The thermodynamic analysis showed that the unfolding enthalpy and entropy values of HsCutA1 were considerably lower than those of OsCutA1 with a similar stability to HsCutA1, which should be related to flexibility of the unstructured properties in both N- and C-terminals of HsCutA1. HsCutA1 was almost completely digested after 1-day incubation at 37 °C by subtilisin, although OsCutA1 was hardly digested at the same conditions. These results indicate that easily available fragmentation of HsCutA1 with remarkably high thermodynamic stability at the body temperature should be important for its protein catabolism in the human cells.

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.

Strategy for cold adaptation of the tryptophan synthase α subunit from the psychrophile Shewanella frigidimarina K14-2: crystal structure and physicochemical properties.

To investigate the molecular basis of cold adaptation of enzymes, we determined the crystal structure of the tryptophan synthase α subunit (SfTSA) from the psychrophile Shewanella frigidimarina K14-2 by X-ray analysis at 2.6-Å resolution and also examined its physicochemical properties. SfTSA was found to have the following characteristics: (i) The stabilities against heat and denaturant of SfTSA were lower than those of an α subunit (EcTSA) from Escherichia coli. This lower equilibrium stability originated from both a faster unfolding rate and a slower refolding rate; (ii) the heat denaturation of SfTSA was completely reversible at pH 7.0 and the solubility of denatured SfTSA was higher than that of denatured EcTSA. The two-state transition of denaturation for SfTSA was highly cooperative, whereas the denaturation process of EcTSA was considerably more complex and (iii) the global structure of SfTSA was quite similar to those of α subunits from other species. Relative to those other proteins, SfTSA exhibited an increase in cavity volume and a decrease in the number of ion pairs. SfTSA also lacks a hydrogen bond near loop B, related to catalytic function. These characteristics of SfTSA might provide the conformational flexibility required for catalytic activity at low temperatures.

Effect of heavy atoms on the thermal stability of α-amylase from Aspergillus oryzae.

Currently, there are no versatile and established methods for improving stability of proteins. In an entirely different approach from conventional techniques such as mutagenesis, we attempted to enhance enzyme stability of α-amylase from Aspergillus oryzae using a heavy-atom derivatization technique. We evaluated changes in stability using differential scanning calorimetry (DSC). Candidate heavy atoms were identified using the Heavy-Atom Database System HATODAS, a Web-based tool designed to assist in heavy-atom derivatization of proteins for X-ray crystallography. The denaturation temperature of α-amylase derivatized with gadolinium (Gd) or samarium (Sm) ions increased by 6.2 or 5.7°C, respectively, compared to that of the native protein (60.6°C). The binding of six Gd ions was confirmed by X-ray crystallography of the enzyme at 1.5 Å resolution. DSC and dynamic light-scattering data revealed a correlation between stability and the aggregation state upon addition of Gd ions. These results show that HATODAS search is an effective tool for selecting heavy atoms for stabilization of this protein.

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.

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.

Role of charged residues in stabilization of Pyrococcus horikoshii CutA1, which has a denaturation temperature of nearly 150 °C.

The CutA1 protein from Pyrococcus horikoshii (PhCutA1), a hyperthermophile, has an unusually high content of charged residues and an unusually high denaturation temperature. To elucidate the role of ion-ion interactions in protein stability, mutant proteins of PhCutA1 in which charged residues were substituted by noncharged residues were comprehensively examined. The denaturation temperatures (T(d)) for 13 of 53 examined mutant proteins were higher than that of the wild-type (148.5 °C at pH 7.0), among which E99Q had the highest T(d) at 154.9 °C. R25A had the largest decrease in T(d) among single mutants at ΔT(d) = -12.4 °C. The average decrease in T(d) of Lys or Arg mutants was greater than that of Glu or Asp mutants, and the average change in T(d) (ΔT(d)) of 21 Glu mutants was negligible, at 0.03 ± 2.05 °C. However, the electrostatic energy (-159.3 kJ·mol(-1)) of PhCutA1 was quite high, compared with that of CutA1 from Escherichia coli (-9.7 kJ·mol(-1)), a mesophile. These results indicate that: (a) many Glu and Asp residues of PhCutA1 should be essential for highly efficient interactions with positively charged residues and for generating high electrostatic energy, although they were forced to be partially repulsive to each other; (b) the changes in stability of mutant proteins with a T(d) value of ~140-150 °C were able to be explained by considering factors important for protein stability and the structural features of mutant sites; and (c) these findings are useful for the design of proteins that are stable at temperatures > 100 °C.

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.

Remarkable improvement in the heat stability of CutA1 from Escherichia coli by rational protein design.

To enhance the heat stability of the CutA1 protein from Escherichia coli (EcCutA1) so that it has comparable stability to CutA1 from Pyrococcus horikoshii with a denaturation temperature (T(d)) of 150°C, we used the Stability Profile of Mutant Protein (SPMP) to examine the structure-sequence (3D-1D) compatibility between the conformation of EcCutA1 and its native sequence [J. Mol. Biol., 248, 733-738, (1995)]. We identified seven residues in EcCutA1 that were incompatible in terms of dihedral angles and hydrophobicity. These residues were replaced with appropriate amino acids, and the mutant proteins were evaluated for changes in stability by DSC and denaturant denaturation. The mutations that were introduced at five out of the seven positions improved the stability of EcCutA1. The T(d) values of single (S11A) and triple (S11V/E61V/Q73V) mutants improved by 16.5 and 26.6°C, respectively, compared to that of the wild-type protein (89.9°C). These analyses showed that (1) the stability of EcCutA1 is remarkably improved by slight substitutions, even though the stability of the wild-type protein is considerably high, (2) remarkable improvements in the stability can be quantitatively explained based on the newly solved native structure, and (3) SPMP is a powerful tool to examine substitutions that improve protein stability.

Large conformational changes in the Escherichia coli tryptophan synthase beta(2) subunit upon pyridoxal 5'-phosphate binding.

To understand the basis for the lower activity of the tryptophan synthase beta(2) subunit in comparison to the alpha(2)beta(2) complex, we determined the crystal structures of apo-beta(2) and holo-beta(2) from Escherichia coli at 3.0 and 2.9 A resolutions, respectively. To our knowledge, this is the first report of both beta(2) subunit structures with and without pyridoxal-5'-phosphate. The apo-type molecule retained a dimeric form in solution, as in the case of the holo-beta(2) subunit. The subunit structures of both the apo-beta(2) and the holo-beta(2) forms consisted of two domains, namely the N domain and the C domain. Although there were significant structural differences between the apo- and holo-structures, they could be easily superimposed with a 22 degrees rigid body rotation of the C domain. The pyridoxal-5'-phosphate-bound holo-form had multiple interactions between the two domains and a long loop (residues 260-310), which were missing in the apo-form. Comparison of the structures of holo-Ecbeta(2) and Stbeta(2) in the alpha(2)beta(2) complex from Salmonella typhimurium (Stalpha(2)beta(2)) identified the cause of the lower enzymatic activity of holo-Ecbeta(2) in comparison with Stalpha(2)beta(2). The substrate (indole) gate residues, Tyr279 and Phe280, block entry of the substrate into the beta(2) subunit, although the indole can directly access the active site as a result of a wider cleft between the N and C domains in the holo-Ecbeta(2) subunit. In addition, the structure around betaAsp305 of the holo-Ecbeta(2) subunit was similar to the open state of Stalpha(2)beta(2) with low activity, resulting in lower activity of holo-Ecbeta(2).