Enhancing extracellular protein production in Escherichia coli by deleting the d‐alanyl‐d‐alanine carboxypeptidase gene dacC

d ‐Alanyl‐d ‐alanine carboxypeptidase DacC is important for synthesis and stabilization of the peptidoglycan layer of Escherichia coli. In this work, dacC of E. coli BL21 (DE3) was successfully deleted, and the effects of this deletion on extracellular protein production in E. coli were investigated. The extracellular activities and fluorescence value of recombinant amylase, green fluorescent protein, and α‐galactosidase of the deletion mutants were increased by 82.3, 29.1, and 37.7%, respectively, compared with that of control cells. The outer membrane permeability and intracellular soluble peptidoglycan accumulation of deletion mutant were also enhanced compared with those of control cells, respectively. Based on fluorescence‐assisted cell sorting analyses, we found that the morphology of the E. coli deletion mutant cells was altered compared with that of control cells. Local transparent bulges in the poles of the E. coli mutant with deletion of the dacC gene were found by transmission electron microscopy analysis. These bulges in the poles could explain the improvement in the production of extracellular protein by the E. coli mutant with deletion of the dacC gene. These findings provide important insights into the extracellular production of proteins using E. coli as microbial cell factories.     

pH‐sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short‐interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH‐dependent residues that affect CIRV p19–siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue. At high pH, the deprotonation of Lys60, Lys67, Lys71, and Cys134 has the largest effect on the binding stability. Similarly, deprotonation of several acidic residues (Asp9, Glu12, Asp20, Glu35, and/or Glu41) at low pH results in a decrease in binding stability. At neutral pH, residues Glu17 and His132 provide a small increase in the binding stability and we find that the optimal pH range for siRNA binding is between 7.0 and 10.0. Overall, our findings further inform recent experiments and are in excellent agreement with data on the pH‐dependent binding profile.     

Role of AP‐endonuclease (Ape1) active site residues in stabilization of the reactant enzyme‐DNA complex

Apurinic/apyrimidinic endonuclease 1 (Ape1) is an important metal‐dependent enzyme in the base excision repair mechanism, responsible for the backbone cleavage of abasic DNA through a phosphate hydrolysis reaction. Molecular dynamics simulations of Ape1 complexed to its substrate DNA performed for models containing 1 or 2 Mg²⁺‐ions as cofactor located at different positions show a complex with 1 metal ion bound on the leaving group site of the scissile phosphate to be the most likely reaction‐competent conformation. Active‐site residue His309 is found to be protonated based on pKa calculations and the higher conformational stability of the Ape1‐DNA substrate complex compared to scenarios with neutral His309. Simulations of the D210N mutant further support the prevalence of protonated His309 and strongly suggest Asp210 as the general base for proton acceptance by a nucleophilic water molecule.     

The role of the 6 lysines and the terminal amine of Escherichia coli single-strand binding protein in its binding of single-stranded DNA

Differential chemical modification of the lysines and amino-terminus of Escherichia coli single-strand binding (SSB) protein was used to determine their roles in the binding of SSB to single-stranded DNA (ssDNA). A combination of isotope labeling and mass spectrometry was used to determine the rates at which SSB was acetylated by acetic anhydride. First, SSB was labeled by deuterated acetic anhydride for given lengths of time in the presence or absence of single-stranded ssDNA. Then, the protein was denatured and completely acetylated by nondeuterated acetic anhydride. Enzymatic digests of the completely acetylated, isotopically labeled SSB were analyzed by electrospray ionization mass spectrometry. The intensities of the deuterated and nondeuterated forms of acetylated peptides provided accurate quantification of the reactivity of the amines in native SSB, either free or bound to ssDNA. Acetylation rate constants were determined from time course measurements. In the absence of ssDNA, the terminal alpha-amine of SSB was 10-fold more reactive than Lys residues at positions 43, 62, 73, and 87. The reactivities of Lys 7 and 49 were much lower yet, suggesting that they have very limited access to solution under any condition. In the presence of ssDNA, the reactivities of the amino-terminus and Lys residues 43, 62, 73, and 87 were reduced by factors of 3.7-25, indicating that the environments around all of these amines is substantially altered by binding of SSB to ssDNA. Three of these residues are located near putative ssDNA binding sites, whereas Lys 87 is located at the monomer-monomer interface.     

Molecular dynamics simulations of apo,holo, and inactivator bound GABA‐at reveal the role of active site residues in PLP dependent enzymes

The pyridoxal 5‐phosphate (PLP) cofactor is a significant organic molecule in medicinal chemistry. It is often found covalently bound to lysine residues in proteins to form PLP dependent enzymes. An example of this family of PLP dependent enzymes is γ‐aminobutyric acid aminotransferase (GABA‐AT) which is responsible for the degradation of the neurotransmitter GABA. Its inhibition or inactivation can be used to prevent the reduction of GABA concentration in brain which is the source of several neurological disorders. As a test case for PLP dependent enzymes, we have performed molecular dynamics simulations of GABA‐AT to reveal the roles of the protein residues and its cofactor. Three different states have been considered: the apoenzyme, the holoenzyme, and the inactive state obtained after the suicide inhibition by vigabatrin. Different protonation states have also been considered for PLP and two key active site residues: Asp298 and His190. Together, 24 independent molecular dynamics trajectories have been simulated for a cumulative total of 2.88 µs. Our results indicate that, unlike in aqueous solution, the PLP pyridine moiety is protonated in GABA‐AT. This is a consequence of a pK_a shift triggered by a strong charge–charge interaction with an ionic “diad” formed by Asp298 and His190 that would help the activation of the first half‐reaction of the catalytic mechanism in GABA‐AT: the conversion of PLP to free pyridoxamine phosphate (PMP). In addition, our MD simulations exhibit additional strong hydrogen bond networks between the protein and PLP: the phosphate group is held in place by the donation of at least three hydrogen bonds while the carbonyl oxygen of the pyridine ring interacts with Gln301; Phe181 forms a π–π stacking interaction with the pyridine ring and works as a gate keeper with the assistance of Val300. All these interactions are hypothesized to help maintain free PMP in place inside the protein active site to facilitate the second half‐reaction in GABA‐AT: the regeneration of PLP‐bound GABA‐AT (i.e., the holoenzyme). Proteins 2016; 84:875–891. © 2016 Wiley Periodicals, Inc.     

Outside‐binding site mutations modify the active site's shapes in neuraminidase from influenza A H1N1

The recent occurrence of 2009 influenza A (H1N1) pandemic as well as others has raised concern of a far more dangerous outcome should this virus becomes resistant to current drug therapies. The number of clinical cases that are resistant to oseltamivir (Tamiflu®) is larger than the limited number of neuraminidase (NA) mutations (H275Y, N295S, and I223R) that have been identified at the active site and that are associated to oseltamivir resistance. In this study, we have performed a comparative analysis between a set of NAs that have the most representative mutations located outside the active site. The recently crystallized NA‐oseltamivir complex (PDB ID: 3NSS) was used as a wild‐type structure. After selecting the target NA sequences, their three‐dimensional (3D) structure was built using 3NSS as a template by homology modeling. The 3D NA models were refined by molecular dynamics (MD) simulations. The refined models were used to perform a docking study, using oseltamivir as a ligand. Furthermore, the docking results were refined by free‐energy analysis using the MM‐PBSA method. The analysis of the MD simulation results showed that the NA models reached convergence during the first 10 ns. Visual inspection and structural measures showed that the mutated NA active sites show structural variations. The docking and MM‐PBSA results from the complexes showed different binding modes and free energy values. These results suggest that distant mutations located outside the active site of NA affect its structure and could be considered to be a new source of resistance to oseltamivir, which agrees with reports in the clinical literature. © 2012 Wiley Periodicals, Inc.     

Role of active‐site residues Tyr55 and Tyr114 in catalysis and substrate specificity of Corynebacterium diphtheriae C‐S lyase

In recent years, there has been increased interest in bacterial methionine biosynthesis enzymes as antimicrobial targets because of their pivotal role in cell metabolism. C‐S lyase from Corynebacterium diphtheriae is a pyridoxal 5′‐phosphate‐dependent enzyme in the transsulfuration pathway that catalyzes the α,β‐elimination of sulfur‐containing amino acids, such as l ‐cystathionine, to generate ammonia, pyruvate, and homocysteine, the immediate precursor of L ‐methionine. In order to gain deeper insight into the functional and dynamic properties of the enzyme, mutants of two highly conserved active‐site residues, Y55F and Y114F, were characterized by UV‐visible absorbance, fluorescence, and CD spectroscopy in the absence and presence of substrates and substrate analogs, as well as by steady‐state kinetic studies. Substitution of Tyr55 with Phe apparently causes a 130‐fold decrease in at pH 8.5 providing evidence that Tyr55 plays a role in cofactor binding. Moreover, spectral data show that the mutant accumulates the external aldimine intermediate suggesting that the absence of interaction between the hydroxyl moiety and PLP‐binding residue Lys222 causes a decrease in the rate of substrate deprotonation. Mutation of Tyr114 with Phe slightly influences hydrolysis of l ‐cystathionine, and causes a change in substrate specificity towards l ‐serine and O‐acetyl‐l ‐serine compared to the wild type enzyme. These findings, together with computational data, provide useful insights in the substrate specificity of C‐S lyase, which seems to be regulated by active‐site architecture and by the specific conformation in which substrates are bound, and will aid in development of inhibitors. Proteins 2015; 83:78–90. © 2014 Wiley Periodicals, Inc.     

Pressure effect on the stability and the conformational dynamics of the D-Galactose/D-Glucose-binding protein from Escherichia coli

The effect of the pressure on the structure and stability of the D-Galactose/D-Glucose binding protein (GGBP) from Escherichia coli was studied by steady-state and time-resolved fluorescence spectroscopy, and the ability of glucose ligand to stabilize the GGBP structure was also investigated. Steady-state fluorescence experiments showed a marked quenching of fluorescence emission of GGBP in the absence of glucose. Instead, the presence of glucose seems to stabilize the structure of GGBP at low and moderate pressure values. Time-resolved fluorescence measurements showed that the GGBP taumean in the absence of glucose varies significantly up to 600 bar, while in the presence of the ligand it is almost unaffected by pressure increase up to 600 bar. The effect of the pressure on GGBP was also studied by molecular dynamics simulations. The simulation data support the spectroscopic results and confirm that the presence of glucose is able to contrast the negative effects of pressure on the protein structure. Taken together, the spectroscopic and computer simulation studies suggest that at pressure values up to 2000 bar the structure of GGBP in the absence of glucose remains folded, but a significant perturbation of the protein secondary structures can be detected. The binding of glucose reduces the negative effect of pressure on protein structure and confers protection from perturbation especially at moderate pressure values.     

Molecular recognition of Fc‐specific ligands binding onto the consensus binding site of IgG: insights from molecular simulation

Immunoglobulin G (IgG) plays an important role in clinical diagnosis and therapeutics. Meanwhile, the consensus binding site (CBS) on the Fc domain of IgG is responsible for ligand recognition, especially for Fc‐specific ligands. In this study, molecular simulation methods were used to investigate molecular interactions between the CBS of the Fc domain and seven natural Fc‐specific ligands. The analysis on the binding energy of the Fc–ligand complex indicated that hydrophobic interactions provide the main driving force for the Fc–ligand binding processes. The hot spots on the ligands and Fc were identified with the computational alanine scanning approach. It was found that the residues of tryptophan and tyrosine on the ligands have significant contributions for the Fc–ligand binding, while Met252, Ile253, Asn434, His435, and Tyr436 are the key residues of Fc. Moreover, two binding modes based on tryptophan or tyrosine were summarized and constructed according to the pairwise interaction analysis. Guidelines for the rational design of CBS‐specific ligands with high affinity and specificity were proposed. Copyright © 2014 John Wiley & Sons, Ltd.     

Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate

We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.     

Exploring the early stages of the pH‐induced conformational change of influenza hemagglutinin

Hemagglutinin (HA) mediates the membrane fusion process of influenza virus through its pH‐induced conformational change. However, it remains challenging to study its structure reorganization pathways in atomic details. Here, we first applied continuous constant pH molecular dynamics approach to predict the pK_a values of titratable residues in H2 subtype HA. The calculated net‐charges in HA1 globular heads increase from 0e (pH 7.5) to +14e (pH 4.5), indicating that the charge repulsion drives the detrimerization of HA globular domains. In HA2 stem regions, critical pH sensors, such as Glu103₂, His18₁, and Glu89₁, are identified to facilitate the essential structural reorganizations in the fusing pathways, including fusion peptide release and interhelical loop transition. To probe the contribution of identified pH sensors and unveil the early steps of pH‐induced conformational change, we carried out conventional molecular dynamics simulations in explicit water with determined protonation state for each titratable residue in different environmental pH conditions. Particularly, energy barriers involving previously uncharacterized hydrogen bonds and hydrophobic interactions are identified in the fusion peptide release pathway. Nevertheless, comprehensive comparisons across HA family members indicate that different HA subtypes might employ diverse pH sensor groups along with different fusion pathways. Finally, we explored the fusion inhibition mechanism of antibody CR6261 and small molecular inhibitor TBHQ, and discovered a novel druggable pocket in H2 and H5 subtypes. Our results provide the underlying mechanism for the pH‐driven conformational changes and also novel insight for anti‐flu drug development. Proteins 2014; 82:2412–2428. © 2014 Wiley Periodicals, Inc.     

Tautomeric states of the active-site histidines of phosphorylated and unphosphorylated IIIGlc, a signal-transducing protein from Escherichia coli, using two-dimensional heteronuclear NMR techniques.

IIIGlc is an 18.1-kDa signal-transducing phosphocarrier protein of the phosphoenolpyruvate:glycose phosphotransferase system from Escherichia coli. The 1H, 15N, and 13C histidine ring NMR signals of both the phosphorylated and unphosphorylated forms of IIIGlc have been assigned using two-dimensional 1H-15N and 1H-13C heteronuclear multiple-quantum coherence (HMQC) experiments and a two-dimensional 13C-13C-1H correlation spectroscopy via JCC coupling experiment. The data were acquired on uniformly 15N-labeled and uniformly 15N/13C-labeled protein samples. The experiments rely on one-bond and two-bond J couplings that allowed for assignment of the signals without the need for the analysis of through-space (nuclear Overhauser effect spectroscopy) correlations. The 15N and 13C chemical shifts were used to determine that His-75 exists predominantly in the N epsilon 2-H tautomeric state in both the phosphorylated and unphosphorylated forms of IIIGlc, and that His-90 exists primarily in the N delta 1-H state in the unphosphorylated protein. Upon phosphorylation of the N epsilon 2 nitrogen of His-90, the N delta 1 nitrogen remains protonated, resulting in the formation of a charged phospho-His-90 moiety. The 1H, 15N, and 13C signals of the phosphorylated and unphosphorylated proteins showed only minor shifts in the pH range from 6.0 to 9.0. These data indicate that the pK alpha values for both His-75 and His-90 in IIIGlc and His-75 in phospho-IIIGlc are less than 5.0, and that the pK alpha value for phospho-His-90 is greater than 10. The results are presented in relation to previously obtained structural data on IIIGlc, and implications for proposed mechanisms of phosphoryl transfer are discussed.     

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     

Titration_DB: Storage and analysis of NMR‐monitored protein pH titration curves

NMR‐monitored pH titration experiments are routinely used to measure site‐specific protein pKa values. Accurate experimental pKa values are essential in dissecting enzyme catalysis, in studying the pH‐dependence of protein stability and ligand binding, in benchmarking pKa prediction algorithms, and ultimately in understanding electrostatic effects in proteins. However, due to the complex ways in which pH‐dependent electrostatic and structural changes manifest themselves in NMR spectra, reported apparent pKa values are often dependent on the way that NMR pH‐titration curves are analyzed. It is therefore important to retain the raw NMR spectroscopic data to allow for documentation and possible re‐interpretation. We have constructed a database of primary NMR pH‐titration data, which is accessible via a web interface. Here, we report statistics of the database contents and analyze the data with a global perspective to provide guidelines on best practice for fitting NMR titration curves. Titration_DB is available at http://enzyme.ucd.ie/Titration_DB . Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The concerted action of a positive charge and hydrogen bonds dynamically regulates the pKa of the nucleophilic cysteine in the NrdH‐redoxin family

NrdH‐redoxins shuffle electrons from the NADPH pool in the cell to Class Ib ribonucleotide reductases, which in turn provide the precursors for DNA replication and repair. NrdH‐redoxins have a CVQC active site motif and belong to the thioredoxin‐fold protein family. As for other thioredoxin‐fold proteins, the pK_a of the nucleophilic cysteine of NrdH‐redoxins is of particular interest since it affects the catalytic reaction rate of the enzymes. Recently, the pK_a value of this cysteine in Corynebacterium glutamicum and Mycobacterium tuberculosis NrdH‐redoxins were determined, but structural insights explaining the relatively low pK_a remained elusive. We subjected C. glutamicum NrdH‐redoxin to an extensive molecular dynamics simulation to expose the factors regulating the pK_a of the nucleophilic cysteine. We found that the nucleophilic cysteine receives three hydrogen bonds from residues within the CVQC active site motif. Additionally, a fourth hydrogen bond with a lysine located N‐terminal of the active site further lowers the cysteine pK_a. However, site‐directed mutagenesis data show that the major contribution to the lowering of the cysteine pK_a comes from the positive charge of the lysine and not from the additional Lys‐Cys hydrogen bond. In 12% of the NrdH‐redoxin family, this lysine is replaced by an arginine that also lowers the cysteine pK_a. All together, the four hydrogen bonds and the electrostatic effect of a lysine or an arginine located N‐terminally of the active site dynamically regulate the pK_a of the nucleophilic cysteine in NrdH‐redoxins.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

A multi‐resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand‐binding site

In multi‐resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse‐grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi‐resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902–1913. © 2016 Wiley Periodicals, Inc.     

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.