Relationship between stability and flexibility in the most flexible region of Photinus pyralis luciferase

Firefly luciferase is a protein with a large N-terminal and a small C-terminal domain. B-factor analysis shows that its C-terminal is much more flexible than its N-terminal. Studies on hyperthermophile proteins have been shown that the increased thermal stability of hyperthermophile proteins is due to their enhanced conformational rigidity and the relationship between flexibility, stability and function in most of proteins is on debate. Two mutations (D474K and D476N) in the most flexible region of firefly luciferase were designed. Thermostability analysis shows that D476N mutation doesn't have any significant effect but D474K mutation destabilized protein. On the other hand, flexibility analysis using dynamic quenching and limited proteolysis demonstrates that D474K mutation became much more flexible than wild type although D476N doesn't have any significant difference. Intrinsic and ANS fluorescence studies demonstrate that D476N mutation is brought about by structural changes without significant effect on thermostability and flexibility. Molecular modeling reveals that disruption of a salt bridge between D⁴⁷⁴ and K⁴⁴⁵ accompanying with some H-bond deletion may be involved in destabilization of D474K mutant.     

A comparison of structural relationships among α-crystallin, human Hsp27, γ-crystallins and βB2-crystallin

The 3D structures of α-crystallin, a major eye lens protein, and related small heat shock proteins are unresolved. It has been assumed that α-crystallin is primarily a β-sheet globular protein similar to γ-crystallin (Siezen and Argos, Biochim. Biophys. Acta, 1983, 748, 56–67) containing sequence repeats in its two domains (Wistow, FEBS Lett. 1985, 181, 1–6). Positional flexibility of amino acid residues and far UV-circular dichroism spectroscopy were used to investigate structural relationships among these proteins. The utility of flexibility plots for predicting protein structure is demonstrated by the excellent correlation of these plots with the known 3D X-ray structures of β/γ-crystallins. Similar analyses of α-crystallin subunits, αA and αB, and human heat shock protein 27 show that the C-terminal domains and connecting segments of these proteins are very similar while the N-terminal domains have significant structural differences. Unlike β/γ-crystallins, both Hsp27 and α-crystallin subunits are asymmetrical with highly flexible C-terminal domains. Flexibility is considered essential for protein functional activity. Therefore, the C-terminal region may play an active role in α-crystallin and small heat shock protein function. Differences in flexibility profiles and estimated secondary structure distribution in α-crystallin by three recent/updated algorithms from far UV-CD spectra support our predicted 3D structure and the concept that α-crystallin and members of β/γ-superfamily are structurally dissimilar.     

Role of N- and C-terminal domains and non-homologous region in co-refolding of Thermotoga maritima β-glucosidase

Thermotoga maritima β-glucosidase consists of three structural regions with 721 amino acids: the N-terminal domain, middle non-homologous region and a C-terminal domain. To investigate the role of these domains in the co-refolding of two fragments into catalytically active form, five sites coding the amino acid residue at 244, 331 in the N-terminal domain, 403 in the non-homologous region, 476 and 521 in the C-terminal domain were selected to split the gene. All the 10 resultant individual fragments were obtained as insoluble inclusion bodies and found to be catalytically inactive. However, the catalytic activity was recovered when the two fragments derived from N-terminal and C-terminal peptides were co-refolded together. It is quite interesting to find that not only the complement polypeptides such as N476/477C but also the truncated combination (N476/522C, amino acid residues from 477 to 521 is truncated) and overlapped combination (N476/245C and N476/404C, amino acid residues from 245 to 476 and from 404 to 476 are overlapped) also gave catalytically active enzymes. Our results showed that folding motifs consisted of the complete N-terminal domain play an important role in the co-refolding of the polypeptides into the catalytically active form.     

Mapping the heparin-binding domain of human hepatic lipase

Human hepatic lipase (HL) is known to bind to the cell surface of hepatocytes and the sinusoidal endothelium of the liver. In each case, it appears that the enzyme remains associated with the cell surface through an ionic interaction with heparan sulfate proteoglycans. However, it remains unclear as to which residues are responsible for this critical function of the enzyme. In the present study, we have used a systematic approach to map the heparin-binding regions of human HL by utilizing peptide arrays spanning the complete sequence of the mature protein. Following probing with biotin-heparin, six peptides spanning residues 301-320 and 465-476 were identified as regions binding to heparin. Probing of an additional array containing these six parent peptides and a comprehensive series of mutant peptides identified two putative HL heparin-binding domains. The first was composed of residues R310, K312, K314, and R315 at the distal N-terminal domain and the second was composed of residues R473, K474, and R476 at the C-terminal end of the protein.     

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.     

The role of the N-domain in the ATPase activity of the mammalian AAA ATPase p97/VCP

p97/valosin-containing protein (VCP) is a type II ATPase associated with various cellular activities that forms a homohexamer with each protomer containing an N-terminal domain (N-domain); two ATPase domains, D1 and D2; and a disordered C-terminal region. Little is known about the role of the N-domain or the C-terminal region in the p97 ATPase cycle. In the p97-associated human disease inclusion body myopathy associated with Paget disease of bone and frontotemporal dementia, the majority of missense mutations are located at the N-domain D1 interface. Structure-based predictions suggest that such mutations affect the interaction of the N-domain with D1. Here we have tested ten major inclusion body myopathy associated with Paget disease of bone and frontotemporal dementia-linked mutants for ATPase activity and found that all have increased activity over the wild type, with one mutant, p97(A232E), having three times higher activity. Further mutagenesis of p97(A232E) shows that the increase in ATPase activity is mediated through D2 and requires both the N-domain and a flexible ND1 linker. A disulfide mutation that locks the N-domain to D1 in a coplanar position reversibly abrogates ATPase activity. A cryo-EM reconstruction of p97(A232E) suggests that the N-domains are flexible. Removal of the C-terminal region also reduces ATPase activity. Taken together, our data suggest that the conformation of the N-domain in relation to the D1-D2 hexamer is directly linked to ATP hydrolysis and that the C-terminal region is required for hexamer stability. This leads us to propose a model where the N-domain adopts either of two conformations: a flexible conformation compatible with ATP hydrolysis or a coplanar conformation that is inactive.     

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.     

Thermal Stability and Unfolding Pathways of Sso7d and its Mutant F31A: Insight from Molecular Dynamics Simulation

Abstract

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.     

Favorable contribution of the C-terminal residue K97 to the stability of a hyperthermophilic archaeal [P62A]Ssh10b

The role of residue K97 at the C-terminal end of archaeal [P62A]Ssh10b in the hyperthermostability of the protein is investigated using three K97-mutant variants: K97E-, K97A-, and ΔK97-mutant [P62A]Ssh10b. The thermal- and GdmHCl-induced denaturation of the three mutant variants has been monitored by circular dichroism. The results reveal that the K97E mutation leads to a stronger destabilization effect than the K97A mutation by disturbing the electrostatic interaction of the salt-bridge D63-K97 and drawing an unfavorable charge-charge repulsive interaction into the structure. However, ΔK97-[P62A]Ssh10b shows much lower stability than K97E- and K97A-mutant [P62A]Ssh10b. Analysis suggests that residue K97 at the C-terminal end makes the favorable contributions to the stability of hyperthermophilic [P62A]Ssh10b not only by the favorable electrostatic interactions with residues in close vicinity but also through maintaining the side chain packing of the surrounding residues in the C-terminal area of the protein.     

SxIP binding disrupts the constitutive homodimer interface of EB1 and stabilizes EB1 monomer

SxIP is a microtubule tip localizing signal found in many +TIP proteins that bind to the hydrophobic cavity of the C-terminal domain of end binding protein 1 (EB1) and then positively regulate the microtubule plus-end tracking of EBs. However, the exact mechanism of microtubule activation of EBs in the presence of SxIP signaling motif is not known. Here, we studied the effect of SxIP peptide on the native conformation of EB1 in solution. Using various NMR experiments, we found that SxIP peptide promoted the dissociation of natively formed EB1 dimer. We also discovered that I224A mutation of EB1 resulted in an unfolded C-terminal domain, which upon binding with the SxIP motif folded to its native structure. Molecular dynamics simulations also confirmed the relative structural stability of EB1 monomer in the SxIP bound state. Residual dipolar couplings and heteronuclear NOE analysis suggested that the binding of SxIP peptide at the C-terminal domain of EB1 decreased the dynamics and conformational flexibility of the N-terminal domain involved in EB1-microtubule interaction. The SxIP-induced disruption of the dimeric interactions in EB1, coupled with the reduction in conformational flexibility of the N-terminal domain of EB1, might facilitate the microtubule association of EB1.     

Proteolytic cleavage of MLL generates a complex of N- and C-terminal fragments that confers protein stability and subnuclear localization

The mixed-lineage leukemia gene (MLL, ALL1, HRX) encodes a 3,969-amino-acid nuclear protein homologous to Drosophila trithorax and is required to maintain proper Hox gene expression. Chromosome translocations in human leukemia disrupt MLL (11q23), generating chimeric proteins between the N terminus of MLL and multiple translocation partners. Here we report that MLL is normally cleaved at two conserved sites (D/GADD and D/GVDD) and that mutation of these sites abolishes the proteolysis. MLL cleavage generates N-terminal p320 (N320) and C-terminal p180 (C180) fragments, which form a stable complex that localizes to a subnuclear compartment. The FYRN domain of N320 directly interacts with the FYRC and SET domains of C180. Disrupting the interaction between N320 and C180 leads to a marked decrease in the level of N320 and a redistribution of C180 to a diffuse nuclear pattern. These data suggest a model in which a dynamic post-cleavage association confers stability to N320 and correct nuclear sublocalization of the complex, to control the availability of N320 for target genes. This predicts that MLL fusion proteins of leukemia which would lose the ability to complex with C180 have their stability conferred instead by the fusion partners, thus providing one mechanism for altered target gene expression.     

Increase in bioluminescence intensity of firefly luciferase using genetic modification

Firefly luciferase is widely used for enzymatic measurement of ATP, and its gene is used as a reporter for gene expression experiments. From our mutant library, we selected novel mutations in Photinus pyralis luciferase with higher luminescence intensity. These included mutations at Ile423, Asp436, and Leu530. Luciferase is structurally composed of a large N-terminal active site domain (residues 1-436), a flexible linker (residues 436-440) peptide, and a small C-terminal domain (residues 440-550) facing the N domain. Thus, the mutations are located at the junction of the N-terminal domain and the flexible linker, in the flexible linker peptide, and in the tip of the C-terminal domain, respectively. Substitution of Asp436 with a nonbulky amino acid such as Gly remarkably increased the substrate affinity for ATP and d-luciferin. Substitution of Ile423 with a hydrophobic amino acid such as Leu and that of Leu530 with a positively charged amino acid such as Arg increased the substrate affinity and the turnover rate. Combining these mutations, we obtained luciferases that generate more than 10-fold higher luminescence intensity than the wild-type enzyme.     

Introducing Wilson disease mutations into the zinc-transporting P-type ATPase of Escherichia coli. The mutation P634L in the 'hinge' motif (GDGXNDXP) perturbs the formation of the E2P state.

ZntA, a bacterial zinc-transporting P-type ATPase, is homologous to two human ATPases mutated in Menkes and Wilson diseases. To explore the roles of the bacterial ATPase residues homologous to those involved in the human diseases, we have introduced several point mutations into ZntA. The mutants P401L, D628A and P634L correspond to the Wilson disease mutations P992L, D1267A and P1273L, respectively. The mutations D628A and P634L are located in the C-terminal part of the phosphorylation domain in the so-called hinge motif conserved in all P-type ATPases. P401L resides near the N-terminal portion of the phosphorylation domain whereas the mutations H475Q and P476L affect the heavy metal ATPase-specific HP motif in the nucleotide binding domain. All mutants show reduced ATPase activity corresponding 0-37% of the wild-type activity. The mutants P401L, H475Q and P476L are poorly phosphorylated by both ATP and P(i). Their dephosphorylation rates are slow. The D628A mutant is inactive and cannot be phosphorylated at all. In contrast, the mutant P634L six residues apart in the same domain shows normal phosphorylation by ATP. However, phosphorylation by P(i) is almost absent. In the absence of added ADP the P634L mutant dephosphorylates much more slowly than the wild-type, whereas in the presence of ADP the dephosphorylation rate is faster than that of the wild-type. We conclude that the mutation P634L affects the conversion between the states E1P and E2P so that the mutant favors the E1 or E1P state.     

Molecular Dynamics Studies of the Inhibitor C34 Binding to the Wild-Type and Mutant HIV-1 gp41: Inhibitory and Drug Resistant Mechanism

Mutations on NHR (N-terminal heptad repeat) associated with resistance to fusion inhibitor were observed. In addition, mutations on CHR (C-terminal heptad repeat) accompanied NHR mutations of gp41 are noted in many cases, like N43D/S138A double mutation. In this work, we explored the drug resistant mechanism of N43D mutation and the role of S138A second mutation in drug resistance. The binding modes of the wild type gp41 and the two mutants, N43D and N43D/S138A, with the HIV-1 fusion inhibitor C34, a 34-residue peptide mimicking CHR of gp41, were carried out by using molecular dynamics simulations. Based on the MD simulations, N43D mutation affects not only the stability of C34 binding, but also the binding energy of the inhibitor C34. Because N43D mutation may also affect the stable conformation of 6-HB, we introduced S138A second mutation into CHR of gp41 and determined the impact of this mutation. Through the comparative analysis of MD results of the N43D mutant and the N43D/S138A mutant, we found that CHR with S138A mutation shown more favorable affinity to NHR. Compelling differences in structures have been observed for these two mutants, particularly in the binding modes and in the hydrophobic interactions of the CHR (C34) located near the hydrophobic groove of the NHR. Because the conformational stability of 6-HB is important to HIV-1 infection, we suggested a hypothetical mechanism for the drug resistance: N43D single mutation not only impact the binding of inhibitor, but also affect the affinity between NHR and CHR of gp41, thus may reduce the rate of membrane fusion; compensatory mutation S138A would induce greater hydrophobic interactions between NHR and CHR, and render the CHR more compatible to NHR than inhibitors.     

Site-site communication in the EF-hand Ca2+-binding protein calbindin D9k

The cooperative binding of Ca2+ ions is an essential functional property of the EF-hand family of Ca2+-binding proteins. To understand how these proteins function, it is essential to characterize intermediate binding states in addition to the apo- and holo-proteins. The three-dimensional solution structure and fast time scale internal motional dynamics of the backbone have been determined for the half-saturated state of the N56A mutant of calbindin D9k with Ca2+ bound only in the N-terminal site. The extent of conformational reorganization and a loss of flexibility in the C-terminal EF-hand upon binding of an ion in the N-terminal EF-hand provide clear evidence of the importance of site-site interactions in this family of proteins, and demonstrates the strength of long-range effects in the cooperative EF-hand Ca2+-binding domain.     

Crystal structures of highly constrained substrate and hydrolysis products bound to HIV-1 protease. Implications for the catalytic mechanism

HIV-1 protease is a key target in treating HIV infection and AIDS, with 10 inhibitors used clinically. Here we used an unusual hexapeptide substrate, containing two macrocyclic tripeptides constrained to mimic a beta strand conformation, linked by a scissile peptide bond, to probe the structural mechanism of proteolysis. The substrate has been cocrystallized with catalytically active synthetic HIV-1 protease and an inactive isosteric (D25N) mutant, and three-dimensional structures were determined (1.60 A). The structure of the inactive HIVPR(D25N)/substrate complex shows an intact substrate molecule in a single orientation that perfectly mimics the binding of conventional peptide ligands of HIVPR. The structure of the active HIVPR/product complex shows two monocyclic hydrolysis products trapped in the active site, revealing two molecules of the N-terminal monocyclic product bound adjacent to one another, one molecule occupying the nonprime site, as expected, and the other monocycle binding in the prime site in the reverse orientation. The results suggest that both hydrolysis products are released from the active site upon cleavage and then rebind to the enzyme. These structures reveal that N-terminal binding of ligands is preferred, that the C-terminal site is more flexible, and that HIVPR can recognize substrate shape rather than just sequence alone. The product complex reveals three carboxylic acids in an almost planar orientation, indicating an unusual hexagonal homodromic complex between three carboxylic acids. The data presented herein regarding orientation of catalytic aspartates support the cleavage mechanism proposed by Northrop. The results imply strategies for design of inhibitors targeting the N-terminal side of the cleavage site or taking advantage of the flexibility in the protease domain that accommodates substrate/inhibitor segments C-terminal to the cleavage site.     

Probing the subtle conformational state of N138ND2-Q106O hydrogen bonding deletion mutant (Asn138Asp) of staphylococcal nuclease using time of flight mass spectrometry with limited proteolysis

Recent studies indicate that the N138ND2-Q106O hydrogen bonding deletion in staphylococcal nuclease significantly alters the conformational integrity and stability of the nuclease. To find out the structural basis of the changes, mass spectrometry and limited proteolysis methods were combined to probe the subtle conformational changes in the SNaseN138D mutant and SNaseN138D-Ca2+-pdTp complex. The results reveal that the N138ND2-Q106O hydrogen bonding deletion makes the C-terminal part of alpha-helix 1 and alpha-helix 2 in the C-terminal subdomain of SNaseN138D unfold to some extent, but does not have much effect on the N-terminal part of alpha-helix 1, alpha-helix 3, and the N-terminal beta-barrel subdomain of SNaseN138D. Binding of ligands makes the alpha-helices 1 and 2 more resistant to protease Glu-C attack and converts the partially unfolded state to a native-like state. This study also demonstrates how mass spectrometry can be combined with limited proteolysis to observe conformational changes induced by ligand binding.     

Relationship between the stability of hen egg-white lysozymes mutated at sites designed to interact with alpha-helix dipoles and their secretion amounts in yeast

The positively charged lysine at the C-terminals of three long alpha-helices (5-15, 25-35, and 88-99) was replaced with alanine (K13A, K33A, K97A) or aspartic acid (K13D, K33D, K97D) in hen lysozyme by genetic engineering. The denaturation transition point (Tm) and Gibbs energy change Delta G of the mutant lysozymes decreased remarkably, suggesting that the positive charge at the C-terminals of helices is involved in the stabilization of the helix dipole. On the other hand, the non-charged asparagine at the N-terminal of the long alpha-helices (25-35 and 88-99) was replaced with negatively charged aspartic acid (N27D and N93D). The Tm and Delta G of N27D increased, suggesting that the dipole moment of the N-terminal of the helices is diminished by replacement with negatively charged amino acid strengthening the stability of the helices. The stabilities of those hen egg white lysozymes mutated at the N- or C-terminal sites of the three long alpha-helices were related with their secretion amounts in yeast (Pichia pastoris). The secretion amounts of these mutant lysozymes in yeast were closely correlated with their stability.     

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.