Role of non-glycine residues in left-handed helical conformation for the conformational stability of human lysozyme

Binding site prediction between lysozyme and glucose-regulated protein 78, a hope to fight amyloidosis

Amyloidosis is an extraordinarily vigorous and heterogeneous group of disorders that causes numerous organ failures due to the precipitation of misfolded proteins. Many of these damaged proteins are discarded before causing any fatal diseases due to the contribution of the protein quality control (PQC) system and its chaperons, including glucose-regulated protein (GRP78). One of the most important enzymatic proteins inside the body is lysozyme, which is reported to have many mutated variants that may cause amyloid fibrils. This study used structural bioinformatics and molecular dynamics simulations to test and suggest binding sites for the human lysozyme protein with GRP78. Multiple sequence alignment (MSA) shows that part of the lysozyme envelope protein (C65-C81 cyclic region) has high similarities (30.77% identity) with the cyclic Pep42. Additionally, the binding between the lysozyme cyclic region (C65-C81) and GRP78 substrate binding domain (SBD) is found favorable. The number and types of interactions vary between each of the mutant isoforms of lysozyme. The more significant the conformational changes in the mutation, the greater its probability of aggregation and the formation of amyloid fibrils. Each mutation leads to different interactions and binding patterns with GRP78. The present computational study suggests a lysozyme-GRP78 binding site, thus paving the way for drug designers to construct suitable carriers that can collect misfolded lysozyme proteins and eliminate them from the body, preventing their aggregation and amyloidogenesis.Communicated by Ramaswamy H. Sarma.

Analyses of displacements resulting from a point mutation in proteins

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

Effect of heterochiral inversions on the structure of a β-hairpin peptide

We study computationally a family of β-hairpin peptides with systematically introduced chiral inversions, in explicit water, and we investigate the extent to which the backbone structure is able to fold in the presence of heterochiral perturbations. In contrast to the recently investigated case of a helical peptide, we do not find a monotonic change in secondary structure content as a function of the number of L- to D-inversions. The effects of L- to D-inversions are instead found to be highly position-specific. Additionally, in contrast to the helical peptide, some inversions increase the stability of the folded peptide: in such cases, we compute an increase in β-sheet content in the aqueous solution equilibrium ensemble. However, the tertiary structures of the stable (folded) configurations for peptides for which inversions cause an increase in β-sheet content show differences from one another, as well as from the native fold of the nonchirally perturbed β-hairpin. Our results suggest that although some chiral perturbations can increase folding stability, chirally perturbed proteins may still underperform functionally, given the relationship between structure and function.

Engineering Proteins for Thermostability with iRDP Web Server

Engineering protein molecules with desired structure and biological functions has been an elusive goal. Development of industrially viable proteins with improved properties such as stability, catalytic activity and altered specificity by modifying the structure of an existing protein has widely been targeted through rational protein engineering. Although a range of factors contributing to thermal stability have been identified and widely researched, the in silico implementation of these as strategies directed towards enhancement of protein stability has not yet been explored extensively. A wide range of structural analysis tools is currently available for in silico protein engineering. However these tools concentrate on only a limited number of factors or individual protein structures, resulting in cumbersome and time-consuming analysis. The iRDP web server presented here provides a unified platform comprising of iCAPS, iStability and iMutants modules. Each module addresses different facets of effective rational engineering of proteins aiming towards enhanced stability. While iCAPS aids in selection of target protein based on factors contributing to structural stability, iStability uniquely offers in silico implementation of known thermostabilization strategies in proteins for identification and stability prediction of potential stabilizing mutation sites. iMutants aims to assess mutants based on changes in local interaction network and degree of residue conservation at the mutation sites. Each module was validated using an extensively diverse dataset. The server is freely accessible at http://irdp.ncl.res.in and has no login requirements.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (C(n), where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C(10)(CO(δ)(i) …NH(i + 2)), C(13)(CO(δ)(i) …NH(i + 3)) and C(17)(N(δ)H(i) …CO(i - 4)) structures. In contrast, Gln predominantly forms C(16)(CO(ε)(i) …NH(i - 3)), C(12)(N(ε)H(i) …CO(i - 2)), C(15)(N(ε)H(i) …CO(i - 3)) and C(18)(N(ε)H(i) …CO(i - 4)) motifs, with only the C(18) motif being analogous to the Asn C(17) structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone β-turns and α-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered.

Structural relationships in the lysozyme superfamily: significant evidence for glycoside hydrolase signature motifs

BACKGROUND

Chitin is a polysaccharide that forms the hard, outer shell of arthropods and the cell walls of fungi and some algae. Peptidoglycan is a polymer of sugars and amino acids constituting the cell walls of most bacteria. Enzymes that are able to hydrolyze these cell membrane polymers generally play important roles for protecting plants and animals against infection with insects and pathogens. A particular group of such glycoside hydrolase enzymes share some common features in their three-dimensional structure and in their molecular mechanism, forming the lysozyme superfamily.

RESULTS

Besides having a similar fold, all known catalytic domains of glycoside hydrolase proteins of lysozyme superfamily (families and subfamilies GH19, GH22, GH23, GH24 and GH46) share in common two structural elements: the central helix of the all-α domain, which invariably contains the catalytic glutamate residue acting as general-acid catalyst, and a β-hairpin pointed towards the substrate binding cleft. The invariant β-hairpin structure is interestingly found to display the highest amino acid conservation in aligned sequences of a given family, thereby allowing to define signature motifs for each GH family. Most of such signature motifs are found to have promising performances for searching sequence databases. Our structural analysis further indicates that the GH motifs participate in enzymatic catalysis essentially by containing the catalytic water positioning residue of inverting mechanism.

CONCLUSIONS

The seven families and subfamilies of the lysozyme superfamily all have in common a β-hairpin structure which displays a family-specific sequence motif. These GH β-hairpin motifs contain potentially important residues for the catalytic activity, thereby suggesting the participation of the GH motif to catalysis and also revealing a common catalytic scheme utilized by enzymes of the lysozyme superfamily.

Functional characterization of the re-face loop spanning residues 536-541 and its interactions with the cofactor in the flavin mononucleotide-binding domain of flavocytochrome P450 from Bacillus megaterium

Flavocytochrome P450BM-3, a bacterial monooxygenase, contains a flavin mononucleotide-binding domain bearing a strong structural homology to the bacterial flavodoxin. The flavin mononucleotide (FMN) serves as the one-electron donor to the heme iron, but in contrast to the electron transfer mechanism of mammalian cytochrome P450 reductase, the FMN semiquinone state is not thermodynamically stable and appears transiently as the anionic rather than the neutral form. A unique loop region comprised of residues (536)Y-N-G-H-P-P(541), which forms a type I' reverse turn and provides several interactions with the FMN isoalloxazine ring, was targeted in this study. Nuclear magnetic resonance studies support the presence of a strong hydrogen bond between the backbone amide of Asn537 and FMN N5, the anionic ionization state of the hydroquinone, and for a change in the hybridization state of the N5 upon reduction. Replacement of Tyr536, which flanks the flavin ring, with a basic residue (histidine or arginine) did not significantly influence the redox properties of the FMN or the accumulation of the anionic semiquinone. The central residues of the type I' turn (Asn-Gly) were replaced with various combinations of glycine and alanine as a means of altering the turn and its interactions. Gly538 was found to be crucial in maintaining the type I' turn conformation of the loop and the strong H-bonding interaction at N5. The functional role of the tandem Pro-Pro sequence which anchors and possible "rigidifies" the loop was investigated through alanine replacements. Despite changes in the stabilities of the oxidized and hydroquinone redox states of the FMN, none of the replacements studied significantly altered the two-electron midpoint potentials. Pro541 does contribute to some degree to the strength of the N5 interaction and the formation of the anionic semiquinone. Unlike that of the flavodoxin, it would appear that the conformation of the FMN rather than the loop changes in response to reduction in this flavoprotein.

Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a robust monomeric protein under kinetic control, which possesses some proline residues at the N-terminal of alpha-helices. Proline residue at the N-terminal of an alpha-helix is thought to stabilize a protein. In this work, the thermostability and folding kinetics of Tk-RNase HII were measured for mutant proteins in which a proline residue is introduced (Xaa to Pro) or removed (Pro to Ala) at the N-terminal of alpha-helices. In the folding experiments, the mutant proteins examined exhibit little influence on the remarkably slow unfolding of Tk-RNase HII. In contrast, E111P and K199P exhibit some thermostabilization, whereas P46A, P70A and P174A have some thermodestabilization. E111P/K199P and P46A/P70A double mutations cause cumulative changes in stability. We conclude that the proline effect on protein thermostability is observed in a hyperthermophilic protein, but each proline residue at the N-terminal of an alpha-helix slightly contributes to the thermostability. The present results also mean that even a natural hyperthermophilic protein can acquire improved thermostability.

Ionic Peptide Aggregation: Exploration of Conformational Dynamics in Aqueous Solution by Computational Techniques

The effects of end groups on KEK peptide conformational characteristics and self-assembling properties in water solution are investigated by using long lasting all-atom molecular dynamics simulations. The analysis of the structural macroscopic and microscopic properties and the examination of intra- and intermolecular interactions suggest, in agreement with experimental observations, the role played by side chains and terminal regions in determining the characteristic features of the assemblages. Competition between intra- and interchain interactions greatly affects the diffusivity of peptide molecules and the conformational space that they can sample, ultimately controlling the shape, size, and distribution of the aggregate configurations. Different peptide end groups influence peptide flexibility and seem to play a crucial role in determining the aggregates' supramolecular architectures.

Calorimetric Analysis of Thermodynamic Stability and Aggregation for Apo and Holo Amyotrophic Lateral Sclerosis-associated Gly-93 Mutants of Superoxide Dismutase

Differential scanning calorimetry was used to measure changes in thermodynamic stability and aggregation for glycine 93 mutants of human copper, zinc-superoxide dismutase (SOD). Glycine 93 is a conserved residue at position i + 3 of a tight turn and has been found to be a mutational hot spot in familial amyotrophic lateral sclerosis (fALS). The fALS-associated mutations, G93A, G93S, G93R, G93D, and G93V, were made in a pseudo wild-type background containing no free cysteines, which prevented the formation of aberrant disulfide bonds upon thermal unfolding, and enabled quantitative thermodynamic analysis of the effects of the mutations. Thermal unfolding was highly reversible for all the SODs in both the fully metallated (holo) and metal-free (apo) forms. The data for all the holo-SODs and for the apo-pseudo-wild-type SOD were well fit by a 2-state unfolding model for native dimer (N2) to two unfolded monomers (2U), N2 <--> 2U. The holo- and apo-forms of the mutants are significantly destabilized (by 1.5-3.5 kcal mol(-1) monomer) relative to the corresponding forms of pseudo wild-type, with the relative stabilities being correlated with statistical preferences for amino acids in this structural context. Although van't Hoff (DeltaHvH) to calorimetric (DeltaHcal) enthalpy ratios are close to unity for all the holo-SODs and for apo-pseudo-wild-type, consistent with a 2-state transition, DeltaHvH is considerably larger than DeltaHcal for all the apo-mutants. This suggests that the mutations cause apo-SOD to have an increased propensity to misfold or aggregate, which may be linked to increased toxic mutant SOD aggregation in fALS.

Structural features of thermozymes

Enzymes synthesized by thermophiles and hyperthermophiles are known as thermozymes. These enzymes are typically thermostable, or resistant to irreversible inactivation at high temperatures, and thermophilic, i.e. optimally active at elevated temperatures between 60 and 125 degrees C. Enzyme thermostability encompasses thermodynamic stability and kinetic stability. Thermodynamic stability is defined by the enzyme's free energy of stabilization (deltaG(stab)) and by its melting temperature (Tm). An enzyme's kinetic stability is often expressed as its halflife (t1/2) at defined temperature. DeltaG(stab) of thermophilic proteins is 5-20 kcal/mol higher than that of mesophilic proteins. The thermostability mechanisms for thermozymes are varied and depend on the enzyme; nevertheless, some common features can be identified as contributing to stability. These features include more interactions (i.e. hydrogen bonds, electrostatic interactions, hydrophobic interactions, disulfide bonds, metal binding) than in less stable enzymes and superior conformational structure (i.e. more rigid, higher packing efficiency, reduced entropy of unfolding, conformational strain release and stability of alpha-helix). Understanding of the stabilizing features will greatly facilitate reengineering of some of the mesozymes to more stable thermozymes.

A Survey of Left-handed Helices in Protein Structures

All naturally occurring amino acids with the exception of glycine contain one or more chiral carbon atoms and can therefore occur in two different configurations, L (levo, left-handed) and D (dextro, right-handed). Proteins are almost exclusively built from L-amino acids. The stereochemical bias of nature is further reflected at the secondary structure level where right-handed helices are strongly preferred over left-handed helices. The handedness of helices has not received much attention in the past and is often overlooked during the analysis, description and deposition of experimentally solved protein structures. Therefore, an extensive survey of left-handed helices in the Protein Data Bank (PDB) was undertaken to analyse their frequency of occurrence, length, amino acid composition, conservation and possible structural or functional role. All left-handed helices (of four or more residues) in a non-redundant subset of the PDB, were identified using hydrogen-bonding analysis, comparison of related structures, and experimental electron density assessment to filter out likely spurious and artefactual hits. This analysis yielded 31 verified left-handed helices in a set of 7284 proteins. The phi angles of the residues in the left-handed helices lie between 30 degrees and 130 degrees and the psi angles lie between -50 degrees and 100 degrees . Most of the helices are short (four residues) and for 87% of them, it was possible to determine that they are important for the stability of the protein, for ligand binding, or as part of the active site. This suggests that, even though left-handed helices are rare, when they do occur, they are structurally or functionally significant. Four secondary structure assignment programs were tested for their ability to identify the handedness of the helices. Of these programs, only DSSP correctly assigns the handedness.

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.

Identification of a key structural element for protein folding within beta-hairpin turns

Specific residues in a polypeptide may be key contributors to the stability and foldability of the unique native structure. Identification and prediction of such residues is, therefore, an important area of investigation in solving the protein folding problem. Atypical main-chain conformations can help identify strains within a folded protein, and by inference, positions where unique amino acids may have a naturally high frequency of occurrence due to favorable contributions to stability and folding. Non-Gly residues located near the left-handed alpha-helical region (L-alpha) of the Ramachandran plot are a potential indicator of structural strain. Although many investigators have studied mutations at such positions, no consistent energetic or kinetic contributions to stability or folding have been elucidated. Here we report a study of the effects of Gly, Ala and Asn substitutions found within the L-alpha region at a characteristic position in defined beta-hairpin turns within human acidic fibroblast growth factor, and demonstrate consistent effects upon stability and folding kinetics. The thermodynamic and kinetic data are compared to available data for similar mutations in other proteins, with excellent agreement. The results have identified that Gly at the i+3 position within a subset of beta-hairpin turns is a key contributor towards increasing the rate of folding to the native state of the polypeptide while leaving the rate of unfolding largely unchanged.

Role of amino acid residues in left-handed helical conformation for the conformational stability of a protein

Our previous study of six non-Gly to Gly/Ala mutant human lysozymes in a left-handed helical region showed that only one non-Gly residue at a rigid site had unfavorable strain energy as compared with Gly at the same position (Takano et al., Proteins 2001; 44:233-243). To further examine the role of left-handed residues in the conformational stability of a protein, we constructed ten Gly to Ala mutant human lysozymes. Most Gly residues in human lysozyme are located in the left-handed helix region. The thermodynamic parameters for denaturation and crystal structures were determined by differential scanning calorimetry and X-ray analysis, respectively. The difference in denaturation Gibbs energy (DeltaDeltaG) for the ten Gly to Ala mutants ranged from + 1.9 to -7.5 kJ/mol, indicating that the effect of the mutation depends on the environment of the residue. We confirm that Gly in a left-handed region is more favorable at rigid sites than non-Gly, but there is little difference in energetic cost between Gly and non-Gly at flexible sites. The present results indicate that dihedral angles in the backbone conformation and also the flexibility at the position should be considered for analyses of protein stability, and protein structural determination, prediction, and design.