Hydrogen-bonding classes in proteins and their contribution to the unfolding reaction

Amyloids, amorphous aggregates and assemblies of peptides - Assessing aggregation

Amyloid and amorphous aggregates represent the two major categories of aggregates associated with diseases, and although exhibiting distinct features, researchers often treat them as equivalent, which demonstrates the need for more thorough characterization. Here, we compare amyloid and amorphous aggregates based on their biochemical properties, kinetics, and morphological features. To further decipher this issue, we propose the use of peptide self-assemblies as minimalistic models for understanding the aggregation process. Peptide building blocks are significantly smaller than proteins that participate in aggregation, however, they make a plausible means to bridge the gap in discerning the aggregation process at the more complex, protein level. Additionally, we explore the potential use of peptide-inspired models to research the liquid-liquid phase separation as a feasible mechanism preceding amyloid formation. Connecting these concepts can help clarify our understanding of aggregation-related disorders and potentially provide novel drug targets to impede and reverse these serious illnesses.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Fourier Transform Infrared Spectroscopy for Assessing Structural and Enzymatic Reactivity Changes Induced during Feather Hydrolysis

Chicken feathers are major byproducts of the livestock processing industry with high potential in the feed sector. In this study, we present a new approach using Fourier transform infrared (FTIR) spectroscopy to detect the structural changes of feather keratin and its availability for enzymatic hydrolysis (AEH) induced by the thermal pressure hydrolysis (TPH) process. Compared to time-consuming in vitro measurement techniques, the proposed method provides rapid information about the structural changes during TPH which enables quick adaptation of TPH conditions as the quality of the incoming feather changes. By analyzing the FTIR spectra of raw and processed feathers, it was found that AEH negatively relates to the β-sheet content (represented by two IR peaks centered at 1635 and 1689 cm^-1), while it positively relates to a new series of peaks centered around 1700 cm^-1 appearing after the TPH process. The proposed FTIR technique provides a reliable and rapid approach to determine the digestibility indicated by AEH of the processed feather and may be used in process control and optimization.

Fever as an evolutionary agent to select immune complexes interfaces

We herein analyzed all available protein–protein interfaces of the immune complexes from the Protein Data Bank whose antigens belong to pathogens or cancers that are modulated by fever in mammalian hosts. We also included, for comparison, protein interfaces from immune complexes that are not significantly modulated by the fever response. We highlight the distribution of amino acids at these viral, bacterial, protozoan and cancer epitopes, and at their corresponding paratopes that belong strictly to monoclonal antibodies. We identify the “hotspots”, i.e. residues that are highly connected at such interfaces, and assess the structural, kinetic and thermodynamic parameters responsible for complex formation. We argue for an evolutionary pressure for the types of residues at these protein interfaces that may explain the role of fever as a selective force for optimizing antibody binding to antigens.

Effect of the single mutation N9Y on the catalytical properties of xylanase Xyn11A from Cellulomonas uda: a biochemical and molecular dynamic simulation analysis

Cellulomonas uda produces Xyn11A, moderately thermostable xylanase, with optimal activity at 50 °C and pH 6.5. An improvement in the biochemical properties of Xyn11A was achieved by site-directed mutagenesis approach. Wild-type xylanase, Xyn11A-WT, and its mutant Xyn11A-N9Y were expressed in Escherichia coli, and then both enzymes were purified and characterized. Xyn11A-N9Y displayed optimal activity at 60 °C and pH 7.5, an upward shift of 10 ºC in the optimum temperature, and an upward shift of one unit in optimum pH; also, it manifested an 11-fold increase in thermal stability at 60 ºC, compared to that displayed by Xyn11A-WT. Molecular dynamics (MD) simulations of Xyn11A-WT and Xyn11A-N9Y suggest the substitution N9Y leads to an array of secondary structure changes at the N-terminal end and an increase in the number of hydrogen bonds in Xyn11A-N9Y. Based on the significant improvements, Xyn11A-N9Y may be considered as a candidate for several biotechnological applications.

Time-dependent conformational analysis of ALK5-lumican complex in presence of graphene and graphene oxide employing molecular dynamics and MMPBSA calculation

Lumican, an extracellular matrix protein avails wound healing by binding to ALK5 membrane receptor (TGF-beta receptor I). Their interaction enables epithelialization and substantiates rejuvenation of injured tissue. To enrich permanence of ALK5-lumican interaction, we employed graphene and graphene oxide co-factors. Herein, this study explicates concomitancy of graphene and graphene oxide with ALK5-lumican. We performed an in silico approach involving molecular modelling, molecular docking, molecular dynamics for 200 ns, DSSP analysis and MMPBSA calculations. Results of molecular dynamics indicate cofactors influential in altering bioactive site of lumican than ALK5. Similarly, MMPBSA calculations unveiled binding energy of apoenzyme as -108.09 kcal/mol, holoenzyme (G) as -79.20 kcal/mol and holoenzyme (GO) as -114.33 kcal/mol. This concludes graphene oxide lucrative in enhancing binding energy of ALK5-lumican in holoenzyme (GO) via coil formation of Lum C₁₃ domain. In contrast, graphene reduced binding energy of ALK5-lumican in holoenzyme (G) modifying Lum C₁₃ into beta sheets. MMPBSA residual contribution analysis of Lum C₁₃ residues revealed binding energy of -13.9 kcal/mol for apoenzyme, -6.8 kcal/mol for holoenzyme (G) and -19.5 kcal/mol for holoenzyme (GO). This supports coil formation propitious for better ALK5-Lum interaction. Highest SASA energy of -21.05 kcal/mol of holoenzyme (G) assures graphene reasonable for improved ALK5-lumican hydrophobicity. As per the motive of the study, graphene oxide enriches permanence of ALK5-lumican. This provides counsel for plausible exploitation of lumican and graphene oxide as targeted/nano drug delivery system to reinstate acute wounds, chronic wounds, corneal wounds, hypertrophic scars and keloids in near future. Communicated by Ramaswamy H. Sarma.

Biochemical characterization of a thermophilic hyaluronate lyase TcHly8C from Thermasporomyces composti DSM22891

Hyaluronic acid (HA) is an anionic linear polysaccharide abundantly distributed in the extracellular matrix of mammalian connective, growing, and tumor tissues. Hyaluronidase is used as an important drug diffusion promoter and a tool enzyme to produce HA oligosaccharides. However, there is no thermostable hyaluronidase suitable for application to date. In this study, a thermophilic hyaluronate lyase, TcHly8C, from Thermasporomyces composti DSM22891 was expressed in Escherichia coli. The recombinant TcHly8C was most active at 70 °C, and it retained about 30% of initial activity after incubation at 60 °C for 28 days. The half-lives of TcHly8C at 60 °C and 70 °C were 16.1 d and 2.3 h, respectively. The optimum pH of TcHly8C is 5.93, and it was stable at pH 6.15-10.90. The presence of Mg²⁺ could enhance its enzymatic activity significantly. K_m, k_cat, and k_cat/K_m of TcHly8C towards HA were 3.69 mg∙ml^-1, 17.82 s^-1, and 4.82 ml∙mg^-1∙s^-1, respectively. TcHly8C degraded HA in an exolytic mode, and the end product was unsaturated HA disaccharide (ΔUA-GlcNAc). Overall, our results show that TcHly8C is the first reported PL8 exo-type hyaluronate lyase with high thermostability, which provides a potential enzyme used in medicine and production of HA oligosaccharides.

Computational study of the intermolecular interactions and their effect on the UV-visible spectra of the ternary liquid mixture of benzene, ethanol and propylene glycol

Quantum chemical calculations are well-equipped to provide answers to the questions regarding the different aspects of intermolecular interactions. We investigate the benzene, ethanol and 1,2 propanediol ternary mixture with theoretical as well as experimental UV-Vis spectroscopy. An extensive theoretical study on the molecular structure and UV-Vis spectral analysis was undertaken using density functional theory (DFT) method. Structural parameter analysis and the HOMO-LUMO (highest occupied molecular orbital-lowest unoccupied molecular orbital) energy gap help to describe the possible interaction between molecules in dimer and in combination. Interaction energy has been calculated from topological study. Time-dependent density functional theory (TDDFT) calculations on dimer/cluster in gas phase help to understand the effect of the molecular interaction on the overall spectral shift and related intensity variation. Our results show that in the ternary mixture, the interaction energies of the interactions are π-π interaction: 0.52-2.57 kcal/mol, Hp-π interaction: 1.15 kcal/mol and H-bonding: 2.49 to 4.46 kcal/mol. The π-π interaction and H-bonding cause red shift in absorption spectra while Hp-π interaction causes blue shift. In the ternary mixture, the strength of different kinds of interaction depends on the concentration, and as each interaction has its own effect on spectral shift, the overall experimental spectra get broader and distorted from the Gaussian shape.

Changes of Thermostability, Organic Solvent, and pH Stability in Geobacillus zalihae HT1 and Its Mutant by Calcium Ion

Thermostable T1 lipase from Geobacillus zalihae has been crystallized using counter-diffusion method under space and Earth conditions. The comparison of the three-dimensional structures from both crystallized proteins show differences in the formation of hydrogen bond and ion interactions. Hydrogen bond and ion interaction are important in the stabilization of protein structure towards extreme temperature and organic solvents. In this study, the differences of hydrogen bond interactions at position Asp43, Thr118, Glu250, and Asn304 and ion interaction at position Glu226 was chosen to imitate space-grown crystal structure, and the impact of these combined interactions in T1 lipase-mutated structure was studied. Using space-grown T1 lipase structure as a reference, subsequent simultaneous mutation D43E, T118N, E226D, E250L, and N304E was performed on recombinant wild-type T1 lipase (wt-HT1) to generate a quintuple mutant term as 5M mutant lipase. This mutant lipase shared similar characteristics to its wild-type in terms of optimal pH and temperature. The stability of mutant 5M lipase improved significantly in acidic and alkaline pH as compared to wt-HT1. 5M lipase was highly stable in organic solvents such as dimethyl sulfoxide (DMSO), methanol, and n-hexane compared to wt-HT1. Both wild-type and mutant lipases were found highly activated in calcium as compared to other metal ions due to the presence of calcium-binding site for thermostability. The presence of calcium prolonged the half-life of mutant 5M and wt-HT1, and at the same time increased their melting temperature (T_m). The melting temperature of 5M and wt-HT1 lipases increased at 8.4 and 12.1 °C, respectively, in the presence of calcium as compared to those without. Calcium enhanced the stability of mutant 5M in 25% (v/v) DMSO, n-hexane, and n-heptane. The lipase activity of wt-HT1 also increased in 25% (v/v) ethanol, methanol, acetonitrile, n-hexane, and n-heptane in the presence of calcium. The current study showed that the accumulation of amino acid substitutions D43E, T118N, E226D, E250L, and N304E produced highly stable T1 mutant when hydrolyzing oil in selected organic solvents such as DMSO, n-hexane, and n-heptane. It is also believed that calcium ion plays important role in regulating lipase thermostability.

Exploring the Structural Mechanism of Covalently Bound E3 Ubiquitin Ligase: Catalytic or Allosteric Inhibition?

Covalent inhibition has recently gained a resurgence of interest in several drug discovery areas. The expansion of this approach is based on evidence elucidating the selectivity and potency of covalent inhibitors when bound to particular amino acids of a biological target. The Nedd4-1, an E3 ubiquitin ligase, is characterized by two covalent binding sites, of which catalytic Cys^cat and allosteric Cys^allo are enclosed. This enzyme has demonstrated inhibition at both the above-mentioned binding sites; however, a detailed molecular understanding of the structural mechanism of inhibition upon Cys^cat and Cys^allo binding remains vague. This prompted us to provide the first account of investigating the preferential covalent binding mode and the underlying structural and molecular dynamic implications. Based on the molecular dynamic analyses, it was evident that although both catalytic and allosteric covalent binding led to greater stability of the enzyme, a preferential covalent mechanism of inhibition was seen in the allosteric-targeted system. This was supported by a more favorable binding energy in the allosteric site compared to the catalytic site, in addition to the larger number of residue interactions and stabilizing hydrogen bonds occurring in the allosteric covalent bound complex. The fundamental dynamic analysis presented in this report compliments, as well as adds to previous experimental findings, thus leading to a crucial understanding of the structural mechanism by which Nedd4-1 is inhibited. The findings from this study may assist in the design of more target-specific Nedd4-1 covalent inhibitors exploring the surface-exposed cysteine residues.

Identification and nanomechanical characterization of the fundamental single-strand protofilaments of amyloid α-synuclein fibrils

The formation and spreading of amyloid aggregates from the presynaptic protein α-synuclein in the brain play central roles in the pathogenesis of Parkinson's disease. Here, we use high-resolution atomic force microscopy to investigate the early oligomerization events of α-synuclein with single monomer angstrom resolution. We identify, visualize, and characterize directly the smallest elementary unit in the hierarchical assembly of amyloid fibrils, termed here single-strand protofilaments. We show that protofilaments form from the direct molecular assembly of unfolded monomeric α-synuclein polypeptide chains. To unravel protofilaments' internal structure and elastic properties, we manipulated nanomechanically these species by atomic force spectroscopy. The single-molecule scale identification and characterization of the fundamental unit of amyloid assemblies provide insights into early events underlying their formation and shed light on opportunities for therapeutic intervention at the early stages of aberrant protein self-assembly.

Molecular Dynamic Simulation of Space and Earth-Grown Crystal Structures of Thermostable T1 Lipase Geobacillus zalihae Revealed a Better Structure

Less sedimentation and convection in a microgravity environment has become a well-suited condition for growing high quality protein crystals. Thermostable T1 lipase derived from bacterium Geobacilluszalihae has been crystallized using the counter diffusion method under space and earth conditions. Preliminary study using YASARA molecular modeling structure program for both structures showed differences in number of hydrogen bond, ionic interaction, and conformation. The space-grown crystal structure contains more hydrogen bonds as compared with the earth-grown crystal structure. A molecular dynamics simulation study was used to provide insight on the fluctuations and conformational changes of both T1 lipase structures. The analysis of root mean square deviation (RMSD), radius of gyration, and root mean square fluctuation (RMSF) showed that space-grown structure is more stable than the earth-grown structure. Space-structure also showed more hydrogen bonds and ion interactions compared to the earth-grown structure. Further analysis also revealed that the space-grown structure has long-lived interactions, hence it is considered as the more stable structure. This study provides the conformational dynamics of T1 lipase crystal structure grown in space and earth condition.

Interaction of the N-AcAβ(13–23)NH2 segment of the beta amyloid peptide with beta-sheet-blocking peptides: site and edge specificity

The region encompassing residues 13–23 of the amyloid beta peptide (Aβ(13–23)) of Alzheimer’s disease is the self-recognition site that initiates toxic oligomerization and fibrillization and also is the site of interaction of Aβ with many other proteins. We describe herein a study by molecular dynamics of the complexes formed by R (= N-AcAβ(13–23)NH2(N-CH3C(O)HHQKLVFFAEDNH2)) with several pseudopeptides designed to form β-sheets with Aβ(1-40,42) and prevent oligomer and fibril formation. Adhesion to both edges of the R β-strand is examined by structure analysis. Umbrella sampling along a dissociation pathway reveals approximate free energies of binding in the submicromolar range. One of the three pseudopeptides binds strongly to one edge of the R β-strand and another to the opposite edge, while the third displays strong binding to both edges. It is desirable to block both edges of the self-recognition site of Aβ to prevent oligomer formation. The study reveals that this may be accomplished by a single pseudopeptide or two in combination. Thus the pseudopeptides, used singly or in pairs, may be competitive inhibitors of Aβ oligomerization at stoichiometric concentrations.

Molecular dynamics study of the monomers and dimers of N-AcAβ(13–23)NH2: on the effect of pH on the aggregation of the amyloid beta peptide of Alzheimer’s disease

The region encompassed by residues 13–23 of the amyloid beta peptide (Aβ(13–23)) of Alzheimer’s disease is the self-recognition site that initiates toxic oligomerization and fibrillization and also is the site of interaction of Aβ with many other proteins. We describe herein a study by molecular dynamics of N-AcAβ(13–23)NH2 (N-CH3C(O)HHQKLVFFAEDNH2) as a model of full-length Aβ(1–40) or Aβ(1–42) and of its dimers. The effect of pH at or below physiological (pH 7.4) is assessed by protonation of one or more of the His residues. The major conformation of the monomer of the systems is a flexible folded structure. Protonation of one or both His residues does not change the conformation in any significant way. The dimers of protonated and unprotonated systems exist almost exclusively as stable antiparallel β-sheets anchored at both ends by intermolecular salt bridges between Lys16 of one chain and the C-terminal residues Glu22 and Asp23 of the other. We also employ the technique of “umbrella sampling” whereby relative binding affinities of the complexes could be determined. In the case of unsymmetrically protonated species, each complex begins dissociation by releasing the weaker salt bridge, breaking interstrand hydrogen bonds, and losing the β-sheet character. The stronger salt bridge is the last to release and presumably is the first to form in the reverse process of aggregation. Umbrella sampling yields the free energy profiles of the dissociation as a function of the separation of the centres of mass. For each system, the dissociation profile has only a shallow maximum. By implication, the reverse process of assembly has almost no barrier. This is an example of entropy–enthalpy compensation that arises naturally during the molecular dynamics – umbrella sampling simulation.

Insights into the unfolding pathway and identification of thermally sensitive regions of phytase from Aspergillus niger by molecular dynamics simulations

Thermal stability is of great importance in the application of commercial phytases. Phytase A (PhyA) is a monomeric protein comprising twelve α-helices and ten β-sheets. Comparative molecular dynamics (MD) simulations (at 310, 350, 400, and 500 K) revealed that the thermal stability of PhyA from Aspergillus niger (A. niger) is associated with its conformational rigidity. The most thermally sensitive regions were identified as loops 8 (residues 83-106), 10 (161-174), 14 (224-230), 17 (306-331), and 24 (442-444), which are present on the surface of the protein. It was observed that solvent-exposed loops denature before or show higher flexibility than buried residues. We observed that PhyA begins to unfold at loops 8 and 14, which further extends to loop 24 at the C-terminus. The intense movement of loop 8 causes the helix H2 and beta-sheet B3 to fluctuate at high temperature. The high flexibility of the H2, H10, and H12 helices at high temperature resulted in complete denaturation. The high mobility of loop 14 easily transfers to the adjacent helices H7, H8, and H9, which fluctuate and partially unfold at high temperature (500 K). It was also observed that the salt bridges Asp110-Lys149, Asp205-Lys277, Asp335-Arg136, Asp416-Arg420, and Glu387-Arg400 are important influences on the structural stability but not the thermostability, as the lengths of these salt bridges did not increase with rising temperature. The salt bridges Glu125-Arg163, Asp299-Arg136, Asp266-Arg219, Asp339-Lys278, Asp335-Arg136, and Asp424-Arg428 are all important for thermostability, as the lengths of these bridges increased dramatically with increasing temperature. Here, for the first time, we have computationally identified the thermolabile regions of PhyA, and this information could be used to engineer novel thermostable phytases. Numerous homologous phytases of fungal as well as bacterial origin are known, and these homologs show high sequence similarity. Our findings could prove useful in attempts to increase the thermostability of homologous phytases via protein engineering.

Engineering and kinetic stabilization of the therapeutic enzyme Anabeana variabilis phenylalanine ammonia lyase

Anabeana variabilis phenylalanine ammonia lyase has just recently been discovered and introduced in clinical trials of phenylketonuria enzyme replacement therapy for its outstanding kinetic properties. In the present study, kinetic stabilization of this therapeutically important enzyme has been explored by introduction of a disulfide bond into the structure. Site-directed mutagenesis was performed with quick-change PCR method. Recombinant wild-type and mutated enzymes were expressed in Escherichia coli, and his-tagged proteins were affinity purified. Formation of disulfide bond was confirmed by Ellman's method, and then chemical unfolding, kinetic behavior, and thermal inactivation of mutated enzyme were compared with the wild type. Based on our results, the Q292C mutation resulted in a significant improvement in kinetic stability and resistance against chemical unfolding of the enzyme while kinetic parameters and pH profile of enzyme activity were remained unaffected. The results of the present study provided an insight towards designing phenylalanine ammonia lyases with higher stability.

Theoretical study on the polar hydrogen-π (Hp-π) interactions between protein side chains

BACKGROUND

In the study of biomolecular structures and interactions the polar hydrogen-π bonds (Hp-π) are an extensive molecular interaction type. In proteins 11 of 20 natural amino acids and in DNA (or RNA) all four nucleic acids are involved in this type interaction.

RESULTS

The Hp-π in proteins are studied using high level QM method CCSD/6-311 + G(d,p) + H-Bq (ghost hydrogen basis functions) in vacuum and in solutions (water, acetonitrile, and cyclohexane). Three quantum chemical methods (B3LYP, CCSD, and CCSD(T)) and three basis sets (6-311 + G(d,p), TZVP, and cc-pVTZ) are compared. The Hp-π donors include R2NH, RNH2, ROH, and C6H5OH; and the acceptors are aromatic amino acids, peptide bond unit, and small conjugate π-groups. The Hp-π interaction energies of four amino acid pairs (Ser-Phe, Lys-Phe, His-Phe, and Tyr-Phe) are quantitatively calculated.

CONCLUSIONS

Five conclusion points are abstracted from the calculation results. (1) The common DFT method B3LYP fails in describing the Hp-π interactions. On the other hand, CCSD/6-311 + G(d,p) plus ghost atom H-Bq can yield better results, very close to the state-of-the-art method CCSD(T)/cc-pVTZ. (2) The Hp-π interactions are point to π-plane interactions, possessing much more interaction conformations and broader energy range than other interaction types, such as common hydrogen bond and electrostatic interactions. (3) In proteins the Hp-π interaction energies are in the range 10 to 30 kJ/mol, comparable or even larger than common hydrogen bond interactions. (4) The bond length of Hp-π interactions are in the region from 2.30 to 3.00 Å at the perpendicular direction to the π-plane, much longer than the common hydrogen bonds (~1.9 Å). (5) Like common hydrogen bond interactions, the Hp-π interactions are less affected by solvation effects.

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.

Relationship between Escherichia coli AppA phytase's thermostability and salt bridges

In order to study on the relationship between Escherichia coli AppA phytase's thermostability and salt bridges, and indicate an effective technical route of which factor to think about and where to modify at AppA for enhancing its thermostability, a salt bridge subtraction mutant E31Q and a salt bridge addition mutant Q307D were constructed by site-directed mutagenesis. The residual activities of the wild-type AppA phytase, E31Q and Q307D were 31.42%, 17.46%, and 40.57%, respectively, after being heated at 80°C for 10 min. The salt bridge subtraction mutant E31Q showed 13.96% thermostability decreasement, and the salt bridge addition mutant Q307D showed 9.15% thermostability enhancement than the wild-type both without the pH and temperature optimum changed. It proved salt bridges play a key role in E. coli AppA phytase's thermostability and the α/β-domain of AppA may be sensitive to heat. Salt bridges and the α/β-domain of AppA should have high priority to think about to enhance AppA's thermostability for commercial application. Besides, molecular dynamics simulation was used for salt bridges analysis.

Mutations from a family-shuffling-library reveal amino acid residues responsible for the thermostability of endoglucanase CelA from Clostridium thermocellum

We constructed a library of chimeras from the major endoglucanase, CelA, of Clostridium thermocellum and a less stable endoglucanase CelB from Clostridium josui with multiple point mutations using low-fidelity family-shuffling method. Mutations that inactivated the enzyme were rapidly eliminated with high-throughput screening. The activities and thermostabilities of selected variants were evaluated, and four amino acid substitutions, K249R, P258S, S329N and E355G, were identified as having significant impact on the thermostability of CelA without affecting enzymatic activity. In the crystal structure of CelA, most of them are away from the activity cleft and are responsible for the stabilization of secondary structures.

Design of short linear peptides that show hydrogen bonding constraints in water

Using a combination of an aromatic amino acid, a homoserine side chain, and a d-amino acid, a series of linear tetrapeptides were designed that adopt an "Hse turn" in water. The conformation was stabilized by intramolecular hydrogen bonds even in the presence of surrounding water molecules. In particular, the peptide with sequence H-Abz-Homoser-Ser-d-Gln-NH(2) showed significant through-space interactions and its free energy of folding is estimated to be on the order of -4 kcal/mol. We report the design of the tetrapeptides using a novel mimicry approach and their characterization based on NMR spectroscopy and MD simulations.

Assembly of mutations for improving thermostability of Escherichia coli AppA2 phytase

We previously identified a number of mutations in Escherichia coli AppA2 phytase for enhancing its thermostability. The objective of the present study was to determine if these mutations (K46E, K65E, G103S, D112N, D144N, S209G, V227A, and G344D) could be sequentially added to further improve the thermostability of AppA2. Compared with the wild-type enzyme, two variants (D144N/V227A and D144N/V227A/G344D) out of the eight resulting mutants showed 15% enhancement in thermostability (as measured by residual activity after being heated at 80 degrees C for 10 min) and 4 to 5 degrees C increases in the melting temperatures (T (m)). Based on the structural predictions with a highly homologous AppA phytase, the substitution D144N introduces a side-chain-side-chain hydrogen bond, thereby stabilizing the loop region (Gln137-Asn144), and the V227A substitution might eliminate structural hindrance between Val222 and Val227 that face each other in the beta-hairpin structure. In addition, overall catalytic efficiency (k (cat)/K (m)) of the two mutants was also improved (P < 0.05) compared to the wild type. However, no further improvement in thermostability was observed by adding other mutations to D144N/V227A/G344D, which might result from unfavorable electrostatic interactions or structural perturbation. In conclusion, our results underscore the potential as well as difficulty of predicting synergistic effects of multiple mutations on thermostability within phytase.

Contribution of charged groups to the enthalpic stabilization of the folded states of globular proteins

In this paper we use the results from all-atom molecular dynamics (MD) simulations of proteins and peptides to assess the individual contribution of charged atomic groups to the enthalpic stability of the native state of globular proteins and investigate how the distribution of charged atomic groups in terms of solvent accessibility relates to protein enthalpic stability. The contributions of charged groups is calculated using a comparison of nonbonded interaction energy terms from equilibrium simulations of charged amino acid dipeptides in water (the "unfolded state") and charged amino acids in globular proteins (the "folded state"). Contrary to expectation, the analysis shows that many buried, charged atomic groups contribute favorably to protein enthalpic stability. The strongest enthalpic contributions favoring the folded state come from the carboxylate (COO(-)) groups of either Glu or Asp. The contributions from Arg guanidinium groups are generally somewhat stabilizing, while N(+)(3) groups from Lys contribute little toward stabilizing the folded state. The average enthalpic gain due to the transfer of a methyl group in an apolar amino acid from solution to the protein interior is described for comparison. Notably, charged groups that are less exposed to solvent contribute more favorably to protein native-state enthalpic stability than charged groups that are solvent exposed. While solvent reorganization/release has favorable contributions to folding for all charged atomic groups, the variation in folded state stability among proteins comes mainly from the change in the nonbonded interaction energy of charged groups between the unfolded and folded states. A key outcome is that the calculated enthalpic stabilization is found to be inversely proportional to the excess charge density on the surface, in support of an hypothesis proposed previously.

Strong and weak hydrogen bonds in the protein-ligand interface

The characteristics of N--H...O, O--H...O, and C--H...O hydrogen bonds and other weak intermolecular interactions are analyzed in a large and diverse group of 251 protein-ligand complexes using a new computer program that was developed in-house for this purpose. The interactions examined in the present study are those which occur in the active sites, defined here as a sphere of 10 A radius around the ligand. Notably, N--H...O and O--H...O bonds tend towards linearity. Multifurcated interactions are especially common, especially multifurcated acceptors, and the average degree of furcation is 2.6 hydrogen bonds per furcated acceptor. A significant aspect of this study is that we have been able to assess the reliability of hydrogen bond geometry as a function of crystallographic resolution. Thresholds of 2.3 and 2.0 A are established for strong and weak hydrogen bonds, below which hydrogen bond geometries may be safely considered for detailed analysis. Interactions involving water as donor or acceptor, and C--H...O bonds with Gly and Tyr as donors are ubiquitous in the active site. A similar trend was observed in an external test set of 233 protein-ligand complexes belonging to the kinase family. Weaker interactions like X--H...pi (X = C, N, O) and those involving halogen atoms as electrophiles or nucleophiles have also been studied. We conclude that the strong and weak hydrogen bonds are ubiquitous in protein-ligand recognition, and that with suitable computational tools very large numbers of strong and weak intermolecular interactions in the ligand-protein interface may be analyzed reliably. Results confirm earlier trends reported previously by us but the extended nature of the present data set mean that the observed trends are more reliable.

Molecular extensibility of mini-dystrophins and a dystrophin rod construct

Muscular dystrophies arise with various mutations in dystrophin, implicating this protein in force transmission in normal muscle. With 24 three-helix, spectrin repeats interspersed with proline-rich hinges, dystrophin's large size is an impediment to gene therapy, prompting the construction of mini-dystrophins. Results thus far in dystrophic mice suggest that at least one hinge between repeats is necessary though not sufficient for palliative effect. One such mini-dystrophin is studied here in forced extension at the single molecule level. Delta2331 consists of repeats (R) and hinges (H) H1-R1-2 approximately H3 approximately R22-24-H4 linked by native (-) and non-native (approximately) sequence. This is compared to its core fragment R2 approximately H3 approximately R22 as well as an eight-repeat rod fragment middle (RFM: R8-15). We show by atomic force microscopy that all repeats extend and unfold at forces comparable to those that a few myosin molecules can generate. The hinge regions most often extend and transmit force while limiting tandem repeat unfolding. From 23-42 degrees C, the dystrophin constructs also appear less temperature-sensitive in unfolding compared to a well-studied betaI-spectrin construct. The results thus reveal new modes of dystrophin flexibility that may prove central to functions of both dystrophin and mini-dystrophins.

Pathway shifts and thermal softening in temperature-coupled forced unfolding of spectrin domains

Pathways of unfolding a protein depend in principle on the perturbation-whether it is temperature, denaturant, or even forced extension. Widely-shared, helical-bundle spectrin repeats are known to melt at temperatures as low as 40-45 degrees C and are also known to unfold via multiple pathways as single molecules in atomic force microscopy. Given the varied roles of spectrin family proteins in cell deformability, we sought to determine the coupled effects of temperature on forced unfolding. Bimodal distributions of unfolding intervals are seen at all temperatures for the four-repeat beta(1-4) spectrin-an alpha-actinin homolog. The major unfolding length corresponds to unfolding of a single repeat, and a minor peak at twice the length corresponds to tandem repeats. Increasing temperature shows fewer tandem events but has no effect on unfolding intervals. As T approaches T(m), however, mean unfolding forces in atomic force microscopy also decrease; and circular dichroism studies demonstrate a nearly proportional decrease of helical content in solution. The results imply a thermal softening of a helical linker between repeats which otherwise propagates a helix-to-coil transition to adjacent repeats. In sum, structural changes with temperature correlate with both single-molecule unfolding forces and shifts in unfolding pathways.

Phenomenological similarities between protein denaturation and small-molecule dissolution: Insights into the mechanism driving the thermal resistance of globular proteins

This article shows that the stability profiles of thermophilic proteins are significantly displaced toward higher temperatures as compared to those of mesophilic proteins. A similar trend characterizes the aqueous transfer of N-alkyl amides. In fact, as a general feature of transfer processes, liquid dissolution profiles are centered at temperatures higher than those of solid ones. This behavior is governed by packing contributions. A partition of the unfolding thermodynamics based on the analysis of phenomenological temperatures common to dissolution and unfolding phenomena provides a clue to understanding the mechanism of thermal stabilization. In fact, the position of stability profiles along the temperature axis does not appear to depend on solvation of internal residues. Instead, it is notably affected by solidlike components, whose progressive decrease appears to drive the heat denaturation temperature increase of most thermostable proteins. As a corollary, it is shown that there are actually two limiting mechanisms of thermal stabilization.

On the dissection of the unfolding reaction by the dissolution thermodynamics of N-alkyl amides

Dissection of the unfolding thermodynamics based on small molecule dissolution led to controversial interpretations. It is proposed here that uncertainty about water transfer processes to be used in first approximation analyses of protein stability may be removed (i) separating liquid-dissolution-like effects from solid-like packing contributions; (ii) taking into account both peptide and side chain dissolution; and (iii) analysing the water-dependent part of the denaturation reaction by the dissolution thermodynamics of liquid N-alkyl amides. Based on these criteria, this paper analyses the entropy of the aqueous transfer of liquid N-alkyl amides filling a gap in a recent model of the unfolding energetics, which was limited to the enthalpy. Both enthalpic and entropic changes accompanying the liquid-dissolution-like immersion of internal amino acid residues in water during unfolding may be unambiguously described within this context. Although the model developed does not deepen our knowledge of protein unfolding, it may be of help in the analysis of whether liquid-dissolution-like effects or solid-like packing contributions play the major role in determining protein stability at elevated temperatures.