Completely Buried, Non-Ion-Paired Glutamic Acid Contributes Favorably to the Conformational Stability of Pyrrolidone Carboxyl Peptidases from Hyperthermophiles†,‡

Identification of candidate genes and residues for improving nitrogen use efficiency in the N-sensitive medicinal plant Panax notoginseng

BACKGROUND

Nitrogen (N) metabolism-related key genes and conserved amino acid sites in key enzymes play a crucial role in improving N use efficiency (NUE) under N stress. However, it is not clearly known about the molecular mechanism of N deficiency-induced improvement of NUE in the N-sensitive rhizomatous medicinal plant Panax notoginseng (Burk.) F. H. Chen. To explore the potential regulatory mechanism, the transcriptome and proteome were analyzed and the three-dimensional (3D) information and molecular docking models of key genes were compared in the roots of P. notoginseng grown under N regimes.

RESULTS

Total N uptake and the proportion of N distribution to roots were significantly reduced, but the NUE, N use efficiency in biomass production (NUEb), the recovery of N fertilizer (RNF) and the proportion of N distribution to shoot were increased in the N0-treated (without N addition) plants. The expression of N uptake- and transport-related genes NPF1.2, NRT2.4, NPF8.1, NPF4.6, AVP, proteins AMT and NRT2 were obviously up-regulated in the N0-grown plants. Meanwhile, the expression of CIPK23, PLC2, NLP6, TCP20, and BT1 related to the nitrate signal-sensing and transduction were up-regulated under the N0 condition. Glutamine synthetase (GS) activity was decreased in the N-deficient plants, while the activity of glutamate dehydrogenase (GDH) increased. The expression of genes GS1-1 and GDH1, and proteins GDH1 and GDH2 were up-regulated in the N0-grown plants, there was a significantly positive correlation between the expression of protein GDH1 and of gene GDH1. Glu192, Glu199 and Glu400 in PnGS1 and PnGDH1were the key amino acid residues that affect the NUE and lead to the differences in GDH enzyme activity. The 3D structure, docking model, and residues of Solanum tuberosum and P. notoginseng was similar.

CONCLUSIONS

N deficiency might promote the expression of key genes for N uptake (genes NPF8.1, NPF4.6, AMT, AVP and NRT2), transport (NPF1.2 and NRT2.4), assimilation (proteins GS1 and GDH1), signaling and transduction (genes CIPK23, PLC2, NLP6, TCP20, and BT1) to enhance NUE in the rhizomatous species. N deficiency might induce Glu192, Glu199 and Glu400 to improve the biological activity of GS1 and GDH, this has been hypothesized to be the main reason for the enhanced ability of N assimilation in N-deficient rhizomatous species. The key genes and residues involved in improving NUE provide excellent candidates for the breeding of medicinal plants.

Structure of oxidized pyrrolidone carboxypeptidase from Fervidobacterium islandicum AW-1 reveals unique structural features for thermostability and keratinolysis

Pyrrolidone carboxypeptidases (Pcps) (E.C. 3.4.19.3) can cleave the peptide bond adjacent to pyro-glutamic acid (pGlu), an N-terminal modification observed in some proteins that provides protection against common proteases. Pcp derived from extremely thermophilic Fervidobacterium islandicum AW-1 (FiPcp), that belongs to the cysteine protease family, is involved in keratin utilization under stress conditions. Although an irreversible oxidative modification of active cysteine to its sulfonic acid derivative (Cys-SO₃H) renders the enzyme inactive, the molecular details for the sulfonic acid modification in inactive Pcp remain unclear. Here, we determined the crystal structure of FiPcp at 1.85 Å, revealing the oxidized form of cysteine sulfonic acid (C156-SO₃H) in the catalytic triad (His-Cys-Glu), which participates in the hydrolysis of pGlu residue containing peptide bond. The three oxygen atoms of cysteine sulfonic acid were stabilized by hydrogen bonds with H180, carbonyl backbone of Q83, and water molecules, resulting in inactivation of FiPcp. Furthermore, FiPcp demonstrated a unique ¹³⁹KKKK¹⁴² motif involved in inter-subunit electrostatic interactions whose mutation significantly affects the thermostability of tetrameric FiPcp. Thus, our high-resolution structure of the first inactive FiPcp with irreversible oxidative modification of active cysteine provides not only the molecular basis of the redox-dependent catalysis of Pcp, but also the structural features of its thermostability.

Active site metals mediate an oligomeric equilibrium in Plasmodium M17 aminopeptidases

M17 leucyl aminopeptidases are metal-dependent exopeptidases that rely on oligomerization to diversify their functional roles. The M17 aminopeptidases from Plasmodium falciparum (PfA-M17) and Plasmodium vivax (Pv-M17) function as catalytically active hexamers to generate free amino acids from human hemoglobin and are drug targets for the design of novel antimalarial agents. However, the molecular basis for oligomeric assembly is not fully understood. In this study, we found that the active site metal ions essential for catalytic activity have a secondary structural role mediating the formation of active hexamers. We found that PfA-M17 and Pv-M17 exist in a metal-dependent dynamic equilibrium between active hexameric species and smaller inactive species that can be controlled by manipulating the identity and concentration of metals available. Mutation of residues involved in metal ion binding impaired catalytic activity and the formation of active hexamers. Structural resolution of Pv-M17 by cryoelectron microscopy and X-ray crystallography together with solution studies revealed that PfA-M17 and Pv-M17 bind metal ions and substrates in a conserved fashion, although Pv-M17 forms the active hexamer more readily and processes substrates faster than PfA-M17. On the basis of these studies, we propose a dynamic equilibrium between monomer ↔ dimer ↔ tetramer ↔ hexamer, which becomes directional toward the large oligomeric states with the addition of metal ions. This sophisticated metal-dependent dynamic equilibrium may apply to other M17 aminopeptidases and underpin the moonlighting capabilities of this enzyme family.

Structural and functional characterization of buffalo oviduct-specific glycoprotein (OVGP1) expressed during estrous cycle

Oviduct-specific glycoprotein (OVGP1) is a high molecular weight chitinase-like protein belonging to GH18 family. It is secreted by non-ciliated epithelial cells of oviduct during estrous cycle providing an essential milieu for fertilization and embryo development. The present study reports the characterization of buffalo OVGP1 through structural modeling, carbohydrate-binding properties and evolutionary analysis. Structural model displayed the typical fold of GH18 family members till the boundary of chitinase-like domain further consisting of a large (β/α)8 TIM barrel sub-domain and a small (α+β) sub-domain. Two critical catalytic residues were found substituted in the catalytic centre (Asp to Phe118, Glu to Leu120) compared with the active chitinase. The carbohydrate-binding groove in TIM barrel was lined with various conserved aromatic residues. Molecular docking with different sugars revealed the involvement of various residues in hydrogen-bonding and non-bonded contacts. Most of the substrate-binding residues were conserved except for a few replacements (Ser13, Lys48, Asp49, Pro50, Asp167, Glu199, Gln272 and Phe275) in comparison with other GH18 members. The residues Trp10, Trp79, Asn80, Gln272, Phe275 and Trp334 were involved in recognition of all six ligands. The α+β sub-domain participated in sugar-binding through Thr270, Gln272, Tyr242 and Phe275. The binding assays revealed significant sugar-binding with purified native and recombinant OVGP1. Phylogenetic analysis revealed that OVGP1 was closely related to AMCases followed by other CLPs and evolution of OVGP1 occurred through several gene duplications. This is the first study describing the structural characteristics of OVGP1 that will further help to understand its interaction with gametes to perform crucial reproductive functions.

Functional analysis of recombinant buffalo lactoferrin and monoferric lobes and their cytotoxic effect on buffalo mammary epithelial cells

Lactoferrin (Lf) has been involved in diverse type of cellular activities and its biochemical properties are species specific. Lf is a bilobal molecule in which each lobe binds with one Fe²⁺/Fe³⁺ ion. A lot of physiological effects of Lf are regulated by its iron binding and release properties; however these properties are species-specific. To understand the iron-binding, thermal stability and cytotoxic effect of buffalo Lf (buLf) and contribution of individual N- and C-terminal lobes therein, buLf and the truncated monoferric lobes were expressed in Kluyveromyces lactis or Pichia pastoris yeast expression systems. The iron-uptake/release behavior and thermal stability of recombinant buLf was observed similar to the Lf purified from buffalo milk. Supplementation of recombinant buLf to the buffalo mammary epithelial cells (BuMEC) culture decreased their proliferation and the cell viability in a dose dependent manner. The cell growth decreased by 37% at 1.0 mg/ml Lf. C-lobe decreased the viability of BuMEC by 15% at 1 mg/ml. The C-lobe showed greater cytotoxic effect against BuMEC in comparison to N-lobe. buLf caused a reduced expression of the casein in BuMEC. At 1.0 mg/ml of buLf, CSN2 transcript level was reduced by 74% and 78% in the normal and hormone free media, respectively. The expression of IL-1β gene in BuMEC increased by 4-5 fold in the presence of 1.0 mg/ml of Lf. The effect was similar to that observed in the involutory mammary gland, suggesting the role of elevated level of Lf in remodeling of buffalo mammary tissue during involution.

Structural, Functional and Phylogenetic Analysis of Sperm Lysozyme-Like Proteins

Sperm lysozyme-like proteins belonging to c-type lysozyme family evolved in multiple forms. Lysozyme-like proteins, viz., LYZL2, LYZL3 or SLLP1, LYZL4, LYZL5 and LYZL6 are expressed in the testis of mammals. Not all members of LYZL family have been uniformly and unambiguously identified in the genome and proteome of mammals. Some studies suggested a role of SLLP1 and LYZL4 in fertilization; however, the function of other LYZL proteins is unknown. We identified all known forms of LYZL proteins in buffalo sperm by LC-MS/MS. Cloning and sequence analysis of the Lyzl cDNA showed 38-50% identity at amino acid level among the buffalo LYZL paralogs, complete conservation of eight cysteines and other signature sequences of c-type lysozyme family. Catalytic residues in SLLP1, LYZL4 and LYZL5 have undergone replacement. The substrate binding residues showed significant variation in LYZL proteins. Residues at sites 62, 101, 114 in LYZL4; 101 in SLLP1; 37, 62, and 101 in LYZL6 were more variable among diverse species. Sites 63 and 108 occupied by tryptophan were least tolerant to variation. Site 37 also showed lower tolerance to substitution in SLLP1, LYZL4 and LYZL5, but more variable in non-testicular lysozymes. Models of LYZL proteins were created by homology modeling and the substrate binding pockets were analyzed in term of binding energies and contacting residues of LYZL proteins with tri-N-acetylglucosamine (NAG)3 in the A-B-C and B-C-D binding mode. Except LYZL6, LYZL proteins did not show significant difference in binding energies in comparison to hen egg white lysozyme in the A-B-C mode. (NAG)3 binding energy in the B-C-D mode was higher by 1.3-2.2 kcal/mol than in A-B-C mode. Structural analysis indicated that (NAG)3 was involved in making more extensive interactions including hydrogen bonding with LYZL proteins in B-C-D mode than in A-B-C mode. Despite large sequence divergence among themselves and with respect to c-type lysozymes, substrate binding residues as well as hydrogen bonding network between (NAG)3 and proteins were mostly conserved. LYZL5 in buffalo and other mammalian species contained additional 10-12 amino acid sequence at c-terminal that matched with ankyrin repeat domain-containing protein 27. Phylogenetic analysis indicated LYZL2 to be most ancient among all the LYZL proteins and that the evolution of LYZL proteins occurred through several gene duplications preceding the speciation of mammals from other vertebrates as distant as reptiles and amphibians.

Design, challenge, and promise of stimuli-responsive nanoantibiotics

Over the past few years, there have been calls for novel antimicrobials to combat the rise of drug-resistant bacteria. While some promising new discoveries have met this call, it is not nearly enough. The major problem is that although these new promising antimicrobials serve as a short-term solution, they lack the potential to provide a long-term solution. The conventional method of creating new antibiotics relies heavily on the discovery of an antimicrobial compound from another microbe. This paradigm of development is flawed due to the fact that microbes can easily transfer a resistant mechanism if faced with an environmental pressure. Furthermore, there has been some evidence to indicate that the environment of the microbe can provide a hint as to their virulence. Because of this, the use of materials with antimicrobial properties has been garnering interest. Nanoantibiotics, (nAbts), provide a new way to circumvent the current paradigm of antimicrobial discovery and presents a novel mechanism of attack not found in microbes yet; which may lead to a longer-term solution against drug-resistance formation. This allows for environment-specific activation and efficacy of the nAbts but may also open up and create new design methods for various applications. These nAbts provide promise, but there is still ample work to be done in their development. This review looks at possible ways of improving and optimizing nAbts by making them stimuli-responsive, then consider the challenges ahead, and industrial applications.Graphical abstractA graphic detailing how the current paradigm of antibiotic discovery can be circumvented by the use of nanoantibiotics.

Structural and functional insights into the catalytic inactivity of the major fraction of buffalo milk xanthine oxidoreductase

Background

Xanthine oxidoreductase (XOR) existing in two interconvertible forms, xanthine dehydrogenase (XDH) and xanthine oxidase (XO), catabolises xanthine to uric acid that is further broken down to antioxidative agent allantoin. XOR also produces free radicals serving as second messenger and microbicidal agent. Large variation in the XO activity has been observed among various species. Both hypo and hyper activity of XOR leads to pathophysiological conditions. Given the important nutritional role of buffalo milk in human health especially in south Asia, it is crucial to understand the functional properties of buffalo XOR and the underlying structural basis of variations in comparison to other species.

Methods and Findings

Buffalo XO activity of 0.75 U/mg was almost half of cattle XO activity. Enzymatic efficiency (k_cat/K_m) of 0.11 sec⁻¹ µM⁻¹ of buffalo XO was 8–10 times smaller than that of cattle XO. Buffalo XOR also showed lower antibacterial activity than cattle XOR. A CD value (Δε_{430 nm}) of 46,000 M⁻¹ cm⁻¹ suggested occupancy of 77.4% at Fe/S I centre. Buffalo XOR contained 0.31 molybdenum atom/subunit of which 48% existed in active sulfo form. The active form of XO in buffalo was only 16% in comparison to ∼30% in cattle. Sequencing revealed 97.4% similarity between buffalo and cattle XOR. FAD domain was least conserved, while metal binding domains (Fe/S and Molybdenum) were highly conserved. Homology modelling of buffalo XOR showed several variations occurring in clusters, especially close to FAD binding pocket which could affect NAD⁺ entry in the FAD centre. The difference in XO activity seems to be originating from cofactor deficiency, especially molybdenum.

Conclusion

A major fraction of buffalo milk XOR exists in a catalytically inactive form due to high content of demolybdo and desulfo forms. Lower Fe/S content and structural factors might be contributing to lower enzymatic efficiency of buffalo XOR in a minor way.

Slow unfolding pathway of hyperthermophilic Tk-RNase H2 examined by pulse proteolysis using the stable protease Tk-subtilisin

The unfolding speed of some hyperthermophilic proteins is significantly slower than those of their mesostable homologues. Ribonuclease H2 from the hyperthermophilic archaeon Thermococcus kodakarensis (Tk-RNase H2) is stabilized by its remarkably slow unfolding rate. In this work, we examined the slow unfolding pathway of Tk-RNase H2 by pulse proteolysis using a superstable subtilisin-like serine protease from T. kodakarensis (Tk-subtilisin). Tk-subtilisin has enzymatic activity in highly concentrated guanidine hydrochloride (GdnHCl), in which Tk-RNase H2 unfolds slowly. The native state of Tk-RNase H2 was completely resistant to Tk-subtilisin, whereas the unfolded state (induced by 4 M GdnHCl) was degraded by Tk-subtilisin. Degradation products of Tk-RNase H2 created from pulse proteolysis during its unfolding were detected by tricine-sodium dodecyl sulfate-polyacrylamide gel electrophoresis. We identified the cleavage sites in Tk-RNase H2 by N-terminal sequencing and mass spectrometry and constructed mimics of the unfolding intermediate of Tk-RNase H2 by protein engineering. The mimics were biophysically characterized. We found that the native state of Tk-RNase H2 (N-state) changed to the I(A)-state that was digested by Tk-subtilisin in the early stage of unfolding. In the slow unfolding pathway, the I(A)-state shifted to two intermediate forms, I(B)-state and I(C)-state. The I(B)-state was digested by Tk-subtilisin in the C-terminal region, but the I(C)-state was a Tk-subtilisin resistant form. These states gradually unfolded through the I(D)-state, in which the N-terminal region was digested. The results indicate that pulse proteolysis, by a superstable protease, was a suitable strategy and an effective tool for analyzing intermediate structures of proteins with slow unfolding properties. We also showed that the N-terminal region contributes to the slow unfolding of Tk-RNase H2, and the C-terminal region is important for folding and stability.

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

BACKGROUND

The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).

RESULTS

To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.

CONCLUSIONS

These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.

The pK(a) values of acidic and basic residues buried at the same internal location in a protein are governed by different factors

The pK(a) values of internal ionizable groups are usually very different from the normal pK(a) values of ionizable groups in water. To examine the molecular determinants of pK(a) values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK(a), whereas Asp38 and Glu38 titrate with elevated pK(a) values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK(a) of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK(a) values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK(a) values of internal groups in proteins. This study suggests that improvements to computational methods for pK(a) calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.

Slow unfolding of monomeric proteins from hyperthermophiles with reversible unfolding

Based on the differences in their optimal growth temperatures microorganisms can be classified into psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Proteins from hyperthermophiles generally exhibit greater stability than those from other organisms. In this review, we collect data about the stability and folding of monomeric proteins from hyperthermophilies with reversible unfolding, from the equilibrium and kinetic aspects. The results indicate that slow unfolding is a general strategy by which proteins from hyperthermophiles adapt to higher temperatures. Hydrophobic interaction is one of the factors in the molecular mechanism of the slow unfolding of proteins from hyperthermophiles.

A buried lysine that titrates with a normal pKa: role of conformational flexibility at the protein-water interface as a determinant of pKa values

Previously we reported that Lys, Asp, and Glu residues at positions 66 and 92 in staphylococcal nuclease (SNase) titrate with pK(a) values shifted by up to 5 pK(a) units in the direction that promotes the neutral state. In contrast, the internal Lys-38 in SNase titrates with a normal pK(a). The crystal structure of the L38K variant shows that the side chain of Lys-38 is buried. The ionizable moiety is approximately 7 A from solvent and ion paired with Glu-122. This suggests that the pK(a) value of Lys-38 is normal because the energetic penalty for dehydration is offset by a favorable Coulomb interaction. However, the pK(a) of Lys-38 was also normal when Glu-122 was replaced with Gln or with Ala. Continuum electrostatics calculations were unable to reproduce the pK(a) of Lys-38 unless the protein was treated with an artificially high dielectric constant, consistent with structural reorganization being responsible for the normal pK(a) value of Lys-38. This reorganization must be local because circular dichroism and NMR spectroscopy indicate that the L38K protein is native-like under all conditions studied. In molecular dynamics simulations, the ion pair between Lys-38 and Glu-122 is unstable. The simulations show that a minor rearrangement of a loop is sufficient to allow penetration of water to the amino moiety of Lys-38. This illustrates both the important roles of local flexibility and water penetration as determinants of pK(a) values of ionizable groups buried near the protein-water interface, and the challenges faced by structure-based pK(a) calculations in reproducing these effects.

Genetic selection reveals the role of a buried, conserved polar residue

The burial of nonpolar surface area is known to enhance markedly the conformational stability of proteins. The contribution from the burial of polar surface area is less clear. Here, we report on the tolerance to substitution of Ser75 of bovine pancreatic ribonuclease (RNase A), a residue that has the unusual attributes of being buried, conserved, and polar. To identify variants that retain biological function, we used a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity. Cell growth at 30 degrees C, 37 degrees C, and 44 degrees C correlated with residue size, indicating that the primary attribute of Ser75 is its small size. The side-chain hydroxyl group of Ser75 forms a hydrogen bond with a main-chain nitrogen. The conformational stability of the S75A variant, which lacks this hydrogen bond, was diminished by DeltaDeltaG = 2.5 kcal/mol. Threonine, which can reinstate this hydrogen bond, provided a catalytically active RNase A variant at higher temperatures than did some smaller residues (including aspartate), indicating that a secondary attribute of Ser75 is the ability of its uncharged side chain to accept a hydrogen bond. These results provide insight on the imperatives for the conservation of a buried polar residue.