The HPr proteins from the thermophile Bacillus stearothermophilus can form domain-swapped dimers

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

Differentiating stable and unstable protein using convolution neural network and molecular dynamics simulations

Protein stability is a critical aspect of molecular biology and biochemistry, hinges on an intricate balance of thermodynamic and structural factors. Determining protein stability is crucial for understanding and manipulating biological machineries, as it directly correlated with the protein function. Thus, this study delves into the intricacies of protein stability, highlighting its dependence on various factors, including thermodynamics, thermal conditions, and structural properties. Moreover, a notable focus is placed on the free energy change of unfolding (ΔGunfolding), change in heat capacity (ΔCp) with protein structural transition, melting temperature (Tm) and number of disulfide bonds, which are critical parameters in understanding protein stability. In this study, a machine learning (ML) predictive model was developed to estimate these four parameters using the primary sequence of the protein. The shortfall of available tools for protein stability prediction based on multiple parameters propelled the completion of this study. Convolutional Neural Network (CNN) with multiple layers was adopted to develop a more reliable ML model. Individual predictive models were prepared for each property, and all the prepared models showed results with high accuracy. The R2 (coefficient of determination) of these models were 0.79, 0.78, 0.92 and 0.92, respectively, for ΔG, ΔCp, Tm and disulfide bonds. A case study on stability analysis of two homologous proteins was presented to validate the results predicted through the developed model. The case study included in silico analysis of protein stability using molecular docking and molecular dynamic simulations. This validation study assured the accuracy of each model in predicting the stability associated properties. The alignment of physics-based principles with ML models has provided an opportunity to develop a fast machine learning solution to replace the computationally demanding physics-based calculations used to determine protein stability. Furthermore, this work provided valuable insights into the impact of mutation on protein stability, which has implications for the field of protein engineering. The source codes are available at https://github.com/Growdeatechnology.

Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study

The histidine-containing phosphocarrier (HPr) is a monomeric protein conserved in Gram-positive bacteria, which may be of mesophilic or thermophilic nature. In particular, the HPr protein from the thermophilic organism B. stearothermophilus is a good model system for thermostability studies, since experimental data, such as crystal structure and thermal stability curves, are available. However, its unfolding mechanism at higher temperatures is yet unclear at a molecular level. Therefore, in this work, we researched the thermal stability of this protein using molecular dynamics simulations, subjecting it to five different temperatures during a time span of 1 μs. The analyses of the structural parameters and molecular interactions were compared with those of the mesophilic homologue HPr protein from B. subtilis. Each simulation was run in triplicate using identical conditions for both proteins. The results showed that the two proteins lose stability as the temperature increases, but the mesophilic structure is more affected. We found that the salt bridge network formed by the triad of Glu3-Lys62-Glu36 residues and the salt bridge made up of Asp79-Lys83 ion pair are key factors to keep stable the thermophilic protein, maintaining the hydrophobic core protected and the structure packed. In addition, these molecular interactions neutralize the negative surface charge, acting as "natural molecular staples".

Two putative zinc metalloproteases contribute to the virulence of Clostridium perfringens strains that cause avian necrotic enteritis

Two putative zinc metalloproteases encoded by Clostridium perfringens have been implicated in the pathogenesis of necrotic enteritis, an economically significant poultry disease that is caused by this anaerobic bacterium. These proteases have ~64% amino acid identity and are encoded by the zmpA and zmpB genes. We screened 83 C. perfringens isolates by PCR for the presence of these genes. The first gene, zmpB, is chromosomally located and was present in all screened strains of C. perfringens, regardless of their origin and virulence. The second gene, zmpA, is plasmid-borne and was only found in isolates derived from chickens with necrotic enteritis. We describe the generation of insertionally inactivated mutants of both zmpA and zmpB in a virulent C. perfringens isolate. For each mutant, a significant (p < 0.001) reduction in virulence was observed in a chicken necrotic enteritis disease model. Examples of each mutant strain were characterized by whole genome sequencing, which showed that there were a few off-site mutations with the potential to affect the virulence of these strains. To confirm the importance of these genes, independently derived zmpA and zmpB mutants were constructed in different virulent C. perfringens isolates and shown to have reduced virulence in the experimental disease induction model. A zmpA-zmpB double mutant also was generated and shown to have significantly reduced virulence, to the same extent as the respective single mutants. Our results provide evidence that both putative zinc metalloproteases play an important role in disease pathogenesis.

The histidine-phosphocarrier protein of the phosphoenolpyruvate: sugar phosphotransferase system of Bacillus sphaericus self-associates

The phosphotransferase system (PTS) is involved in the use of carbon sources in bacteria. Bacillus sphaericus, a bacterium with the ability to produce insecticidal proteins, is unable to use hexoses and pentoses as the sole carbon source, but it has ptsHI genes encoding the two general proteins of the PTS: enzyme I (EI) and the histidine phosphocarrier (HPr). In this work, we describe the biophysical and structural properties of HPr from B. sphaericus, HPr^bs, and its affinity towards EI of other species to find out whether there is inter-species binding. Conversely to what happens to other members of the HPr family, HPr^bs forms several self-associated species. The conformational stability of the protein is low, and it unfolds irreversibly during heating. The protein binds to the N-terminal domain of EI from Streptomyces coelicolor, EIN^sc, with a higher affinity than that of the natural partner of EIN^sc, HPr^sc. Modelling of the complex between EIN^sc and HPr^bs suggests that binding occurs similarly to that observed in other HPr species. We discuss the functional implications of the oligomeric states of HPr^bs for the glycolytic activity of B. sphaericus, as well as a strategy to inhibit binding between HPr^sc and EIN^sc.

Is it possible to stabilize a thermophilic protein further using sequences and structures of mesophilic proteins: a theoretical case study concerning DgAS

Incorporating structural elements of thermostable homologs can greatly improve the thermostability of a mesophilic protein. Despite the effectiveness of this method, applying it is often hampered. First, it requires alignment of the target mesophilic protein sequence with those of thermophilic homologs, but not every mesophilic protein has a thermophilic homolog. Second, not all favorable features of a thermophilic protein can be incorporated into the structure of a mesophilic protein. Furthermore, even the most stable native protein is not sufficiently stable for industrial applications. Therefore, creating an industrially applicable protein on the basis of the thermophilic protein could prove advantageous. Amylosucrase (AS) can catalyze the synthesis of an amylose-like polysaccharide composed of only α-1,4-linkages using sucrose as the lone energy source. However, industrial development of AS has been hampered owing to its low thermostability. To facilitate potential industrial applications, the aim of the current study was to improve the thermostability of Deinococcus geothermalis amylosucrase (DgAS) further; this is the most stable AS discovered to date. By integrating ideas from mesophilic AS with well-established protein design protocols, three useful design protocols are proposed, and several promising substitutions were identified using these protocols. The successful application of this hybrid design method indicates that it is possible to stabilize a thermostable protein further by incorporating structural elements of less-stable homologs.

3.8 Protein and Nucleic Acid Folding: Domain Swapping in Proteins

Among thousands of homo-oligomeric protein structures, there is a small but growing subset of ‘domain-swapped’ proteins. The term ‘domain swapping,’ originally coined by D. Eisenberg, describes a scenario in which two or more polypeptide chains exchange identical units for oligomerization. This type of assembly could play a role in disease-related aggregation and amyloid formation or as a specific mechanism for regulating function. This chapter introduces terms and features concerning domain swapping, summarizes ideas about its putative mechanisms, reports on domain-swapped structures collected from the literature, and describes a few notable examples in detail.

Dynamics of Tryptophan in the Histidine-Containing Phosphocarrier Protein of Streptomyces coelicolor: Evidence of Multistate Equilibrium Unfolding

The nanosecond dynamics of the single tryptophan, Trp10, of HPr from Streptomyces coelicolor, HPrsc, has been monitored at different pHs. Time-resolved fluorescence methods and DOSY measurements have been used to map the compactness of the protein. At low pHs, where a molten globule-like species has been described, the correlation times from fluorescence showed an abrupt change as the pH was increased. When the protein was folded (above pH 4), two correlation times were observed, which remained practically constant up to pH 9.5. The long correlation time, around 7.5 ns, corresponds to the global rotational motion of the protein, since this value is in agreement with that determined theoretically from hydrodynamic measurements. The short correlation time, around 1.4 ns, must report on fast movements of the protein segment containing the tryptophan residue. On the other hand, fluorescence lifetimes showed the same abrupt change as the correlation times at low pH, but, in addition, a sigmoidal change with a pKa approximately 4.3 was also observed. On the basis of the modeled structure of HPrsc, this last transition could be due to the proximity of Glu12 to Trp10. The changes monitored by the fluorescence lifetimes agree with those observed previously by steady-state fluorescence, CD, and ANS binding experiments. Taken together, these data suggest a multistate equilibrium during folding of HPrsc starting from low pHs.

Characterization of folding intermediates of a domain-swapped protein by solid-state NMR spectroscopy

We have employed two-dimensional solid-state NMR to study structure and dynamics of insoluble folding states of the domain-swapped protein Crh. Starting from the protein precipitated at its pI, conformational changes due to a modest temperature increase were investigated at the level of individual residues and in real-time. As compared to the crystalline state, Crh pI-precipitates exhibited a higher degree of molecular mobility for several regions of the protein. A rigidly intact center was observed including a subset of residues of the hydrophobic core. Raising the temperature by 13 K to 282 K created a partially unfolded intermediate state that was converted into beta-sheet-rich aggregates that are mostly of spherical character according to electron microscopy. Residue-by-residue analysis indicated that two out of three alpha-helices in aggregated Crh underwent major structural rearrangements while the third helix was preserved. Residues in the hinge region exhibited major chemical-shift changes, indicating that the domain swap was not conserved in the aggregated form. Our study provides direct evidence that protein aggregates of a domain-swapped protein retain a significant fraction of native secondary structure and demonstrates that solid-state NMR can be used to directly monitor slow molecular folding events.

How phosphotransferase system-related protein phosphorylation regulates carbohydrate metabolism in bacteria

The phosphoenolpyruvate(PEP):carbohydrate phosphotransferase system (PTS) is found only in bacteria, where it catalyzes the transport and phosphorylation of numerous monosaccharides, disaccharides, amino sugars, polyols, and other sugar derivatives. To carry out its catalytic function in sugar transport and phosphorylation, the PTS uses PEP as an energy source and phosphoryl donor. The phosphoryl group of PEP is usually transferred via four distinct proteins (domains) to the transported sugar bound to the respective membrane component(s) (EIIC and EIID) of the PTS. The organization of the PTS as a four-step phosphoryl transfer system, in which all P derivatives exhibit similar energy (phosphorylation occurs at histidyl or cysteyl residues), is surprising, as a single protein (or domain) coupling energy transfer and sugar phosphorylation would be sufficient for PTS function. A possible explanation for the complexity of the PTS was provided by the discovery that the PTS also carries out numerous regulatory functions. Depending on their phosphorylation state, the four proteins (domains) forming the PTS phosphorylation cascade (EI, HPr, EIIA, and EIIB) can phosphorylate or interact with numerous non-PTS proteins and thereby regulate their activity. In addition, in certain bacteria, one of the PTS components (HPr) is phosphorylated by ATP at a seryl residue, which increases the complexity of PTS-mediated regulation. In this review, we try to summarize the known protein phosphorylation-related regulatory functions of the PTS. As we shall see, the PTS regulation network not only controls carbohydrate uptake and metabolism but also interferes with the utilization of nitrogen and phosphorus and the virulence of certain pathogens.

Lessons in stability from thermophilic proteins

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.

A thermodynamic comparison of HPr proteins from extremophilic organisms

A thermodynamic stability study of five histidine-containing phosphocarrier protein (HPr) homologues derived from organisms inhabiting diverse environments is described. These HPr homologues are from Bacillus subtilis (Bs), Streptococcus thermophilus (St), Bacillus staerothermophilus (Bst), Bacillus halodurans (Bh), and Oceanobacillus iheyensis (Oi). Analyses of solvent and thermal denaturation experiments provide the cardinal thermodynamic parameters, like deltaG, deltaH, deltaS, T(m), and deltaC(p), that characterize the conformational stability for each homologue. The homologue from Bacillus staerothermophilus (BstHPr) was established as the most thermostable homologue and also the homologue with highest deltaG at all temperatures. A good correlation between habitat temperature of the organism and thermal stability of the protein is also seen. Stability curves (deltaG vs T) for every homologue are also reported; these reveal very similar deltaC(p) and temperature of maximum stability (T(S)) values for all HPr homologues. Stability curves show that the higher thermal stability of some homologues is not a result of change in curvature of the curve or a shift to higher temperature, but rather a displacement of the stability curves to higher deltaG values. Stability curves also allowed estimation of deltaG at habitat temperature of the organisms, and we find good agreement between homologues. Electrostatic contributions to stability of each homologue were investigated by measuring stability as a function of varying pH and NaCl concentration, and our results suggest that most HPr homologues share similar electrostatic contributions to stability.