The role of calcium in the conformational dynamics and thermal stability of the D-galactose/D-glucose-binding protein from Escherichia coli

We have characterized stability and conformational dynamics of the calcium depleted D-galactose/D-glucose-binding protein (GGBP) from Escherichia coli. The structural stability of the protein was investigated by steady state and time resolved fluorescence, and far-UV circular dichroism in the temperature range from 20 degrees C to 70 degrees C. We have found that the absence of the Ca(2+) ion results in a significant destabilization of the C-terminal domain of the protein. In particular, the melting temperature decreases by about 10 degrees C with the simultaneous loss of the melting cooperativity. Time resolved fluorescence quenching revealed significant loosening of the protein when highly shielded Trp residue(s) became accessible to acrylamide at higher temperatures. We have documented a significant stabilizing effect of glucose that mostly reverts the effect of calcium, that is, the thermal stability of the protein increases by about 10 degrees C and the melting cooperativity is restored. Moreover, the protein structure remains compact with low amplitude of the segmental mobility up to high temperatures. We have used molecular dynamics to identify the structural feature responsible for changes in the temperature stability. Disintegration of the Ca(2+)-binding loop seems to be responsible for the loss of the stability in the absence of calcium. The new insights on the structural properties and temperature stability of the calcium depleted GGBP contribute to better understanding of the protein function and constitute important information for the development of new biotechnological applications of this class of proteins.     

Time-resolved fluorescence spectroscopy and molecular dynamics simulations point out the effects of pressure on the stability and dynamics of the porcine odorant-binding protein

The effects of hydrostatic pressure on the structure and stability of porcine odorant-binding protein (pOBP) in the presence and absence of the odorant molecule 2-isobutyl-3-methoxypyrazine (IBMP) were studied by steady-state and time-resolved fluorescence spectroscopy as well as by molecular dynamics simulation. The authors found that the application of moderate values of hydrostatic pressure to pOBP solutions perturbed the microenvironment of Trp(16) and disrupted its highly quenched complex with Met(39). In addition, compared with the protein in the absence of IBMP, the MD simulations experiments carried out at different pressures highlighted the role of this ligand in stabilizing the Trp(16)/Met(39) interaction even at 2000 bar. The obtained results will assist for the tailoring of this protein as specific sensing element in a new class of fluorescence-based biosensors for the detection of explosives.     

D-trehalose/D-maltose-binding protein from the hyperthermophilic archaeon Thermococcus litoralis: the binding of trehalose and maltose results in different protein conformational states

In this work, we used fluorescence spectroscopy, molecular dynamics simulation, and Fourier transform infrared spectroscopy for investigating the effect of trehalose binding and maltose binding on the structural properties and the physical parameters of the recombinant D-trehalose/D-maltose binding protein (TMBP) from the hyperthermophilic archaeon Thermococcus litoralis. The binding of the two sugars to TMBP was studied in the temperature range 20 degrees-100 degrees C. The results show that TMBP possesses remarkable temperature stability and its secondary structure does not melt up to 90 degrees C. Although both the secondary structure itself and the sequence of melting events were not significantly affected by the sugar binding, the protein assumes different conformations with different physical properties depending whether maltose or trehalose is bound to the protein. At low and moderate temperatures, TMBP possesses a structure that is highly compact both in the absence and in the presence of two sugars. At about 90 degrees C, the structure of the unliganded TMBP partially relaxes whereas both the TMBP/maltose and the TMBP/trehalose complexes remain in the compact state. In addition, Fourier transform infrared results show that the population of alpha-helices exposed to the solvent was smaller in the absence than in the presence of the two sugars. The spectroscopic results are supported by molecular dynamics simulations. Our data on dynamics and stability of TMBP can contribute to a better understanding of transport-related functions of TMBP and constitute ground for targeted modifications of this protein for potential biotechnological applications.     

The differences in the microenvironment of the two tryptophan residues of the glutamine-binding protein from Escherichia coli shed light on the binding properties and the structural dynamics of the protein

Glutamine-binding protein (GlnBP) from Escherichia coli is a monomer (26 kDa) that is responsible for the first step in the active transport of L-glutamine across the cytoplasmic membrane. GlnBP consists of two domains (termed large and small) linked by two antiparallel beta-strands. The large domain is similar to the small domain but it contains two additional alpha-helices and three more short antiparallel beta-strands. The deep cleft formed between the two domains contains the ligand-binding site. The binding of L-glutamine leads to cleft closing and a significant structural change with the formation of the so-called "closed form" structure. The protein contains two tryptophan residues (W32 and W220) and 10 tyrosine residues. We used phosphorescence spectroscopy measurements to characterize the role of the two tryptophan residues in the protein structure in the absence and the presence of glutamine. Our results pointed out that the phosphorescence of GlnBP is easily detected in fluid solutions where the emission of the two tryptophan residues is readily discriminated by the drastic difference in the phosphorescence lifetime allowing the assignments of the short lifetime to W220 and the long lifetime to W32. In addition, our results showed that the triplet lifetime of the superficial W220 is unusually short because of intramolecular quenching by the proximal Y163. On the contrary, the lifetime of W32 is several hundred milliseconds long, implicating a well-ordered, compact fold of the surrounding polypeptide. The spectroscopic data were analyzed and discussed together with a detailed inspection of the 3D structure of GlnBP.     

Crystal structure of the dipeptide binding protein from Escherichia coli involved in active transport and chemotaxis.

The Escherichia coli periplasmic dipeptide binding protein functions in both peptide transport and taxis toward peptides. The structure of the dipeptide binding protein in complex with Gly-Leu (glycyl-L-leucine) has been determined at 3.2 A resolution. The binding site for dipeptides is designed to recognize the ligand's backbone while providing space to accommodate a variety of side chains. Some repositioning of protein side chains lining the binding site must occur when the dipeptide's second residue is larger than leucine. The protein's fold is very similar to that of the Salmonella typhimurium oligopeptide binding protein, and a comparison of the structures reveals the structural basis for the dipeptide binding protein's preference for shorter peptides.     

Ligand binding specificity of the Escherichia coli periplasmic histidine binding protein,HisJ

The HisJ protein from Escherichia coli and related Gram negative bacteria is the periplasmic component of a bacterial ATP‐cassette (ABC) transporter system. Together these proteins form a transmembrane complex that can take up L‐histidine from the environment and translocate it into the cytosol. We have studied the specificity of HisJ for binding L‐His and many related naturally occurring compounds. Our data confirm that L‐His is the preferred ligand, but that 1‐methyl‐L‐His and 3‐methyl‐L‐His can also bind, while the dipeptide carnosine binds weakly and D‐histidine and the histidine degradation products, histamine, urocanic acid and imidazole do not bind. L‐Arg, homo‐L‐Arg, and post‐translationally modified methylated Arg‐analogs also bind with reasonable avidity, with the exception of symmetric dimethylated‐L‐Arg. In contrast, L‐Lys and L‐Orn have considerably weaker interactions with HisJ and methylated and acetylated Lys variants show relatively poor binding. It was also observed that the carboxylate group of these amino acids and their variants was very important for proper recognition of the ligand. Taken together our results are a key step towards designing HisJ as a specific protein‐based reagentless biosensor.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

Conformational stability of the full‐atom hexameric model of the ClpB chaperone from Escherichia coli

The Escherichia coli heat shock protein ClpB, a member of the Hsp100 family, plays a crucial role in cellular thermotolerance. In co‐operation with the Hsp70 chaperone system, it is able to solubilize proteins aggregated by heat shock conditions and refold them into the native state in an ATP‐dependent way. It was established that the mechanism of ClpB action depends on the formation of a ring‐shaped hexameric structure and the translocation of a protein substrate through an axial channel. The structural aspects of this process are not fully known. By means of homology modeling and protein–protein docking, we obtained a model of the hexameric arrangement of the full‐length ClpB protein complexed with ATP. A molecular dynamics simulation of this model was performed to assess its flexibility and conformational stability. The high mobility of the “linker” M‐domain, essential for the renaturing activity of ClpB, was demonstrated, and the size and shape of central channel were analyzed. In this model, we propose the coordinates for a loop between b4 and B6 structural elements, not defined in previous structural research, which faces the inside of the channel and may therefore play a role in substrate translocation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 47–60, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Binding of glutamine to glutamine-binding protein from Escherichia coli induces changes in protein structure and increases protein stability

Glutamine-binding protein (GlnBP) from Escherichia coli is a monomeric protein localized in the periplasmic space of the bacterium. It is responsible for the first step in the active transport of L-glutamine across the cytoplasmic membrane. The protein consists of two similar globular domains linked by two peptide hinges, and X-ray crystallographic data indicate that the two domains undergo large movements upon ligand binding. Fourier transform infrared spectroscopy (FTIR) was used to analyze the structure and thermal stability of the protein in detail. The data indicate that glutamine binding induces small changes in the secondary structure of the protein and that it renders the structure more thermostable and less flexible. Detailed analyses of IR spectra show a lower thermal sensitivity of alpha-helices than beta-sheets in the protein both in the absence and in the presence of glutamine. Generalized two-dimensional (2D) analyses of IR spectra reveal the same sequence of unfolding events in the protein in the absence and in the presence of glutamine, indicating that the amino acid does not affect the unfolding pathway of the protein. The data give new insight into the structural characteristics of GlnBP that are useful for both basic knowledge and biotechnological applications.     

The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (eta(xy)), and steady state [1H]-15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 degrees C for 89-100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/eta(xy) was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to approximately 0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be approximately 9 degrees C lower (78 degrees C rather than 87 degrees C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Experimental evaluation of single‐domain antibodies predicted by molecular dynamics simulations to have elevated thermal stability

Recently Bekker et al. [Bekker G‐J et al. Protein Sci. 2019;28:429–438.] described a computational strategy of applying molecular‐dynamics simulations to estimate the relative stabilities of single‐domain antibodies, and utilized their method to design changes with the aim of increasing the stability of a single‐domain antibody with a known crystal structure. The structure from which they generated potentially stabilizing mutations is an anti‐cholera toxin single domain antibody selected from a naïve library which has relatively low thermal stability, reflected by a melting point of 48°C. Their work was purely theoretical, so to examine their predictions, we prepared the parental and predicted stabilizing mutant single domain antibodies and examined their thermal stability, ability to refold and affinity. We found that the mutation that improved stability the most (~7°C) was one which changed an amino acid in CDR1 from an asparagine to an aspartic acid. This change unfortunately was also accompanied by a reduction in affinity. Thus, while their modeling did appear to successfully predict stabilizing mutations, introducing mutations in the binding regions is problematic. Of further interest, the mutations selected via their high temperature simulations, did improve refolding, suggesting that they were successful in stabilizing the structure at high temperatures and thereby decrease aggregation. Our result should permit them to reassess and refine their model and may one day lead to a usefulin silico approach to protein stabilization.     

Role of membrane potential in protein folding and domain formation during secretion in Escherichia coli

The synthesis and processing of the periplasmic components of the leucine transport system of E coli have been studied to determine the role played by transmembrane potential in protein secretion. Both the leucine-isoleucine-valine binding protein and the leucine-specific binding protein are synthesized as precursors with 23 amino acid N-terminal leader sequences. The processing of these precursors is sensitive to the transmembrane potential. Since the amino acid sequence and the crystal structure have been determined for the leucine-isoleucine-valine binding protein, it and the closely related leucine-specific binding protein represent convenient models in which to examine the mechanism of protein secretion in E coli. A model for secretion has been proposed, suggesting a role for transmembrane potential. In this model, the N-terminal amino acid sequence of the precursor is assumed to form a hairpin of two helices. The membrane potential may orient this structure to make it accessible to processing. In addition, the model suggests that a negatively charged, folded domain of the secretory protein may electrophorese toward the trans-positive side of the membrane, thus providing an additional role for the transmembrane potential.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

A recombinant glutamine-binding protein from Escherichia coli: effect of ligand-binding on protein conformational dynamics

We have investigated the effect of the binding of glutamine on the conformational dynamics of the recombinant glutamine binding protein (GlnBP) from Escherichia coli by steady-state and time-resolved fluorescence techniques. The structural stability of the protein was also studied by far-UV circular dichroism spectroscopy in the range of temperature between 25 and 80 degrees C. The results showed that the interaction of the protein with the ligand resulted in a marked change of the structural and conformational dynamics features of the protein. In particular, the fluorescence and circular dichroism data showed that the presence of glutamine resulted in a dramatic increase of the protein thermal stability of about 10 degrees C. In addition, the fluorescence time-resolved data pointed out that both in the absence and in the presence of glutamine the protein structure was highly rigid with small amplitude of segmental motion up to 65 degrees C and a low accessibility of the protein tryptophan residues to acrylamide. The obtained results on the structural properties of the recombinant glutamine-binding protein in the absence and in the presence of glutamine can contribute to a better understanding of the transport-related functions of the protein and structurally similar periplasmic transport proteins, as well as to the design and development of new biotechnological applications of this class of proteins.     

Stability and interactions of the amino-terminal domain of ClpB from Escherichia coli

ClpB is a member of a multichaperone system in Escherichia coli (with DnaK, DnaJ, and GrpE) that reactivates aggregated proteins. The sequence of ClpB contains two ATP-binding regions that are enclosed between the N- and C-terminal extensions. Whereas it has been found that the N-terminal region of ClpB is essential for the chaperone activity, the structure of this region is not known, and its biochemical properties have not been studied. We expressed and purified the N-terminal fragment of ClpB (residues 1-147). Circular dichroism of the isolated N-terminal region showed a high content of alpha-helical structure. Differential scanning calorimetry showed that the N-terminal region of ClpB is thermodynamically stable and contains a single folding domain. The N-terminal domain is monomeric, as determined by gel-filtration chromatography, and the elution profile of the N-terminal domain does not change in the presence of the N-terminally truncated ClpB (ClpBDeltaN). This indicates that the N-terminal domain does not form strong contacts with ClpBDeltaN. Consistently, addition of the separated N-terminal domain does not reverse an inhibition of ATPase activity of ClpBDeltaN in the presence of casein. As shown by ELISA measurements, full-length ClpB and ClpBDeltaN bind protein substrates (casein, inactivated luciferase) with similar affinity. We also found that the isolated N-terminal domain of ClpB interacts with heat-inactivated luciferase. Taken together, our results indicate that the N-terminal fragment of ClpB forms a distinct domain that is not strongly associated with the ClpB core and is not required for ClpB interactions with other proteins, but may be involved in recognition of protein substrates.     

Molecular dynamics study on the stability of wild-type and the R220K mutant of human prion protein

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases related to the structure transition of α-helix into β-sheet. In order to gain more direct insight into the molecular basis of the disease, the stability of the wild-type human prion protein (hPrPc) and the R220K mutant (m-hPrPc) was studied by molecular dynamics (MD) and flow MD simulation. Both the thermodynamic stability and the mechanical properties of hPrPc were investigated in this work. It was found that β-sheet was more readily to be unfolded in m-hPrPc. In the case of hPrPc, less content of helix was preserved after water turbulence. The H-bond network formed by the mutation-related residue 220 was found to play a key role in the stability of hPrPc.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.