Stabilization provided by neighboring strands is critical for the mechanical stability of proteins

Single-molecule force spectroscopy studies and steered molecular dynamics simulations have revealed that protein topology and pulling geometry play important roles in determining the mechanical stability of proteins. Most studies have focused on local interactions that are associated with the force-bearing β-strands. Interactions mediated by neighboring strands are often overlooked. Here we use Top7 and barstar as model systems to illustrate the critical importance of the stabilization effect provided by neighboring β-strands on the mechanical stability. Using single-molecule atomic force microscopy, we showed that Top7 and barstar, which have similar topology in their force-bearing region, exhibit vastly different mechanical-stability characteristics. Top7 is mechanically stable and unfolds at ∼150 pN, whereas barstar is mechanically labile and unfolds largely below 50 pN. Steered molecular dynamics simulations revealed that stretching force peels one force-bearing strand away from barstar to trigger unfolding, whereas Top7 unfolds via a substructure-sliding mechanism. This previously overlooked stabilization effect from neighboring β-strands is likely to be a general mechanism in protein mechanics and can serve as a guideline for the de novo design of proteins with significant mechanical stability and novel protein topology.     

Introduction of a thermophile-sourced ion pair network in the fourth beta/alpha unit of a psychophile-derived triosephosphate isomerase from Methanococcoides burtonii significantly increases its kinetic thermal stability

Hyperthermophile proteins commonly have higher numbers of surface ionic interactions than homologous proteins from other domains of life. PfuTIM, a triosephosphate isomerase (TIM) from the hyperthermophile archaeon, Pyrococcus furiosus, contains an intricate network of 4 ion pairs in its 4th beta/alpha unit, (β/α)4, whereas MbuTIM, a triosephosphate isomerase from a psychrophile archaeon, Methanococcoides burtonii, lacks this network. Notably, (β/α)4 is the first element of the structure formed during folding of certain TIM-type (beta/alpha)8 barrel proteins. Previously, we have shown that elimination of PfuTIM's ion pair network in PfuTIM significantly decreases its kinetic structural stability. Here, we describe the reciprocal experiment in which this ion pair network is introduced into MbuTIM, to produce MutMbuTIM. Recombinant MbuTIM displays multi-state unfolding with apparent Tm values of autonomous structural elements approaching, or above, 70 °C, when a temperature scanning rate of 90 °C/h is used. The protein displays significant intrinsic kinetic stability, i.e., there is a marked temperature scan rate-dependence of the Tm values associated with unfolding transitions. The Tm values drop by as much as ~ 10 °C when the temperature scanning rate is lowered to 5 °C/h. MutMbuTIM, incorporating PfuTIM's ion pair network, shows significantly higher apparent Tm values (raised by 4–6 °C over those displayed by MbuTIM). MutMbuTIM also displays significantly higher kinetic thermal stability. Thus, it appears that the thermal stability of triosephosphate isomerase can be increased, or decreased, by either enhancing, or reducing, the strength of ion pair interactions stabilizing (β/α)4, presumably through reduced cooperativity (and increased autonomy) in unfolding transitions.     

Unfolding thermodynamics of the Delta-domain in the prohead I subunit of phage HK97: determination by factor analysis of Raman spectra

An early step in the morphogenesis of the double-stranded DNA (dsDNA) bacteriophage HK97 is the assembly of a precursor shell (prohead I) from 420 copies of a 384-residue subunit (gp5). Although formation of prohead I requires direct participation of gp5 residues 2-103 (Δ-domain), this domain is eliminated by viral protease prior to subsequent shell maturation and DNA packaging. The prohead I Δ-domain is thought to resemble a phage scaffolding protein, by virtue of its highly α-helical secondary structure and a tertiary fold that projects inward from the interior surface of the shell. Here, we employ factor analysis of temperature-dependent Raman spectra to characterize the thermostability of the Δ-domain secondary structure and to quantify the thermodynamic parameters of Δ-domain unfolding. The results are compared for the Δ-domain within the prohead I architecture (in situ) and for a recombinantly expressed 111-residue peptide (in vitro). We find that the α-helicity (∼ 70%), median melting temperature (T_m = 58 °C), enthalpy (ΔH_m = 50 ± 5 kcal mol^− 1), entropy (ΔS_m = 150 ± 10 cal mol^− 1 K^− 1), and average cooperative melting unit (〈n_c〉 ∼ 3.5) of the in situ Δ-domain are altered in vitro, indicating specific interdomain interactions within prohead I. Thus, the in vitro Δ-domain, despite an enhanced helical secondary structure (∼ 90% α-helix), exhibits diminished thermostability (T_m = 40 °C; ΔH_m = 27 ± 2 kcal mol^− 1; ΔS_m = 86 ± 6 cal mol^− 1 K^− 1) and noncooperative unfolding (〈n_c〉 ∼ 1) vis-à-vis the in situ Δ-domain. Temperature-dependent Raman markers of subunit side chains, particularly those of Phe and Trp residues, also confirm different local interactions for the in situ and in vitro Δ-domains. The present results clarify the key role of the gp5 Δ-domain in prohead I architecture by providing direct evidence of domain structure stabilization and interdomain interactions within the assembled shell.     

Changing the Mechanical Unfolding Pathway of FnIII10 by Tuning the Pulling Strength

We investigate the mechanical unfolding of the tenth type III domain from fibronectin (FnIII₁₀) both at constant force and at constant pulling velocity, by all-atom Monte Carlo simulations. We observe both apparent two-state unfolding and several unfolding pathways involving one of three major, mutually exclusive intermediate states. All three major intermediates lack two of seven native β-strands, and share a quite similar extension. The unfolding behavior is found to depend strongly on the pulling conditions. In particular, we observe large variations in the relative frequencies of occurrence for the intermediates. At low constant force or low constant velocity, all three major intermediates occur with a significant frequency. At high constant force or high constant velocity, one of them, with the N- and C-terminal β-strands detached, dominates over the other two. Using the extended Jarzynski equality, we also estimate the equilibrium free-energy landscape, calculated as a function of chain extension. The application of a constant pulling force leads to a free-energy profile with three major local minima. Two of these correspond to the native and fully unfolded states, respectively, whereas the third one can be associated with the major unfolding intermediates.     

Single-molecule force spectroscopy reveals a stepwise unfolding of Caenorhabditis elegans giant protein kinase domains

Myofibril assembly and disassembly are complex processes that regulate overall muscle mass. Titin kinase has been implicated as an initiating catalyst in signaling pathways that ultimately result in myofibril growth. In titin, the kinase domain is in an ideal position to sense mechanical strain that occurs during muscle activity. The enzyme is negatively regulated by intramolecular interactions occurring between the kinase catalytic core and autoinhibitory/regulatory region. Molecular dynamics simulations suggest that human titin kinase acts as a force sensor. However, the precise mechanism(s) resulting in the conformational changes that relieve the kinase of this autoinhibition are unknown. Here we measured the mechanical properties of the kinase domain and flanking Ig/Fn domains of the Caenorhabditis elegans titin-like proteins twitchin and TTN-1 using single-molecule atomic force microscopy. Our results show that these kinase domains have significant mechanical resistance, unfolding at forces similar to those for Ig/Fn β-sandwich domains (30-150 pN). Further, our atomic force microscopy data is consistent with molecular dynamic simulations, which show that these kinases unfold in a stepwise fashion, first an unwinding of the autoinhibitory region, followed by a two-step unfolding of the catalytic core. These data support the hypothesis that titin kinase may function as an effective force sensor.     

NMR characterization of a 264-residue hyperthermostable endo-β-1,3-glucanase

Insight into the hyperthermostable endo-β-1,3-glucanase pfLamA from Pyrococcus furiosus is obtained by using NMR spectroscopy. pfLamA functions optimally at 104 °C and recently the X-ray structure of pfLamA has been obtained at 20 °C, a temperature at which the enzyme is inactive. In this study, near-complete (>99%) NMR assignments are presented of chemical shifts of pfLamA in presence and absence of calcium at 62 °C, a temperature at which the enzyme is biologically active. The protein contains calcium and the effects of calcium on the protein are assessed. Calcium binding results in relatively small chemical shift changes in a region distant from the active site of pfLamA and thus causes only minor conformational modifications. Removal of calcium does not significantly alter the denaturation temperature of pfLamA, implying that calcium does not stabilize the enzyme against global unfolding. The data obtained form the basis for elucidation of the molecular origins involved in conformational stability and biological activity of hyperthermophilic endo-β-1,3-glucanases at extreme temperatures.     

Characterization of the protein unfolding processes induced by urea and temperature

Correct folding is critical for the biological activities of proteins. As a contribution to a better understanding of the protein (un)folding problem, we studied the effect of temperature and of urea on peptostreptococcal Protein L destructuration. We performed standard molecular dynamics simulations at 300 K, 350 K, 400 K, and 480 K, both in 10 M urea and in water. Protein L followed at least two alternative unfolding pathways. Urea caused the loss of secondary structure acting preferentially on the β-sheets, while leaving the α-helices almost intact; on the contrary, high temperature preserved the β-sheets and led to a complete loss of the α-helices. These data suggest that urea and high temperature act through different unfolding mechanisms, and protein secondary motives reveal a differential sensitivity to various denaturant treatments. As further validation of our results, replica-exchange molecular dynamics simulations of the temperature-induced unfolding process in the presence of urea were performed. This set of simulations allowed us to compute the thermodynamical parameters of the process and confirmed that, in the configurational space of Protein L unfolding, both of the above pathways are accessible, although to a different relative extent.     

Cloning and characterization of a new β-mannosidase from Streptomyces sp. S27

A new β-mannosidase gene, designated as man2S27, was cloned from Streptomyces sp. S27 using the colony PCR method and expressed in Escherichia coli BL21 (DE3). The full-length gene consists of 2499 bp and encodes 832 amino acids with a calculated molecular mass of 92.6 kDa. The amino acid sequence shares highest identity of 62.6% with the mannosidase Man2A from Cellulomonas fimi which belongs to the glycoside hydrolase family 2. Purified recombinant Man2S27 showed optimal activity at pH 7.0 and 50 °C. The specific activity, K_m, and k_cat values for p-nitrophenyl-β-d-mannopyranoside (p-NP-β-MP) were 35.3 U mg^-1, 0.23 mM, and 305 s^-1, respectively. Low transglycosylation activity was observed when Man2S27 was incubated with p-NP-β-MP (glycosyl donor) and methyl-α-d-mannopyranoside (p-NP-α-MP) (acceptor) at 50 °C and pH 7.0, and a small amount of methylmannobioside was synthesized. Using locust bean gum as the substrate, more reducing sugars were liberated by the synergistic action of Man2S27 and β-mannanase (Man5S27), and the synergy degree in sequential reactions with Man5S27 firstly and Man2S27 secondly was higher than that in the simultaneous reactions.     

Conformational stability of α-amylase from malted sorghum (Sorghum bicolor): Reversible unfolding by denaturants

α-Amylase from Sorghum bicolor, is reversibly unfolded by chemical denaturants at pH 7.0 in 50 mM Hepes containing 13.6 mM calcium and 15 mM DTT. The isothermal equilibrium unfolding at 27 °C is characterized by two state transition with ΔG (H₂O) of 16.5 kJ mol⁻¹ and 22 kJ mol⁻¹, respectively, at pH 4.8 and pH 7.0 for GuHCl and ΔG (H₂O) of 25.2 kJ mol⁻¹ at pH 4.8 for urea. The conformational stability indicators such as the change in excess heat capacity (ΔC_p), the unfolding enthalpy (H_g) and the temperature at ΔG = 0 (T_g) are 17.9 ± 0.7 kJ mol⁻¹ K⁻¹, 501.2 ± 18.2 kJ mol⁻¹ and 337.3 ± 6.9 K at pH 4.8 and 14.3 ± 0.5 kJ mol⁻¹ K⁻¹, 509.3 ± 21.7 kJ mol⁻¹ and 345.4 ± 4.8 K at pH 7.0, respectively. The reactivity of the conserved cysteine residues, during unfolding, indicates that unfolding starts from the ‘B’ domain of the enzyme. The oxidation of cysteine residues, during unfolding, can be prevented by the addition of DTT. The conserved cysteine residues are essential for enzyme activity but not for the secondary and tertiary fold acquired during refolding of the denatured enzyme. The pH dependent stability described by ΔG (H₂O) and the effect of salt on urea induced unfolding confirm the role of electrostatic interactions in enzyme stability.     

Ca2+ Binding Enhanced Mechanical Stability of an Archaeal Crystallin

Structural topology plays an important role in protein mechanical stability. Proteins with β-sandwich topology consisting of Greek key structural motifs, for example, I27 of muscle titin and ¹⁰FNIII of fibronectin, are mechanically resistant as shown by single-molecule force spectroscopy (SMFS). In proteins with β-sandwich topology, if the terminal strands are directly connected by backbone H-bonding then this geometry can serve as a “mechanical clamp”. Proteins with this geometry are shown to have very high unfolding forces. Here, we set out to explore the mechanical properties of a protein, M-crystallin, which belongs to β-sandwich topology consisting of Greek key motifs but its overall structure lacks the “mechanical clamp” geometry at the termini. M-crystallin is a Ca²⁺ binding protein from Methanosarcina acetivorans that is evolutionarily related to the vertebrate eye lens β and γ-crystallins. We constructed an octamer of crystallin, (M-crystallin)₈, and using SMFS, we show that M-crystallin unfolds in a two-state manner with an unfolding force ∼90 pN (at a pulling speed of 1000 nm/sec), which is much lower than that of I27. Our study highlights that the β-sandwich topology proteins with a different strand-connectivity than that of I27 and ¹⁰FNIII, as well as lacking “mechanical clamp” geometry, can be mechanically resistant. Furthermore, Ca²⁺ binding not only stabilizes M-crystallin by 11.4 kcal/mol but also increases its unfolding force by ∼35 pN at the same pulling speed. The differences in the mechanical properties of apo and holo M-crystallins are further characterized using pulling speed dependent measurements and they show that Ca²⁺ binding reduces the unfolding potential width from 0.55 nm to 0.38 nm. These results are explained using a simple two-state unfolding energy landscape.     

Identification of a Mechanical Rheostat in the Hydrophobic Core of Protein L

The ability of proteins and their complexes to withstand or respond to mechanical stimuli is vital for cells to maintain their structural organisation, to relay external signals and to facilitate unfolding and remodelling. Force spectroscopy using the atomic force microscope allows the behaviour of single protein molecules under an applied extension to be investigated and their mechanical strength to be quantified. protein L, a simple model protein, displays moderate mechanical strength and is thought to unfold by the shearing of two mechanical sub-domains. Here, we investigate the importance of side-chain packing for the mechanical strength of protein L by measuring the mechanical strength of a series of protein L variants containing single conservative hydrophobic volume deletion mutants. Of the five thermodynamically destabilised variants characterised, only one residue (I60V) close to the interface between two mechanical sub-domains was found to differ in mechanical properties to wild type (ΔF_I60V-WT = − 36 pN at 447 nm s^− 1, Δx_uI60V-WT = 0.2 nm). Φ-value analysis of the unfolding data revealed a highly native transition state. To test whether the number of hydrophobic contacts across the mechanical interface does affect the mechanical strength of protein L, we measured the mechanical properties of two further variants. protein L L10F, which increases core packing but does not enhance interfacial contacts, increased mechanical strength by 13 ± 11 pN at 447 nm s^− 1. By contrast, protein L I60F, which increases both core and cross-interface contacts, increased mechanical strength by 72 ± 13 pN at 447 nm s^− 1. These data suggest a method by which nature can evolve a varied mechanical response from a limited number of topologies and demonstrate a generic but facile method by which the mechanical strength of proteins can be rationally modified.     

Cloning,characterization, hypoxia and heat shock response of hypoxia inducible factor-1 (HIF-1) from the small abalone Haliotis diversicolor

In this study, hypoxia inducible factor-1α (HIF-1α) and hypoxia inducible factor-1β (HIF-1β) from small abalone Haliotis diversicolor were cloned. The cDNA of H. diversicolor HIF-1α (HdHIF-1α) is 2833 bp encoding a protein of 711aa and H. diversicolor HIF-1β (HdHIF-1β) is 1919 bp encoding a protein of 590aa. Similar to other species' HIF-1, HdHIF-1 has one basic helix–loop–helix (bHLH) domain and two Per-Arnt-Sim (PAS) domains, and HdHIF-1α has a oxygen-dependent degradation domain (ODDD) with two proline hydroxylation motifs and a C-terminal transactivation domain (C-TAD) with an asparagine hydroxylation motif. Under normoxic conditions, HdHIF-1α and HdHIF-1β mRNAs were constitutively present in all examined tissues. Under hypoxia (2.0 mg/L DO at 25 °C) stress, HdHIF-1α expression was up-regulated in gills at 4 h, 24 h and 96 h, and in hemocytes at 24 h and 96 h, while HdHIF-1β remained relatively constant. Under thermal stress (31 °C), HdHIF-1α expression was significantly increased in gills at 4 h, and hemocytes at 0 h and 4 h, while HdHIF-1β expression still remained relatively constant. These results suggested that HIF-1α may play an important role in adaption to poor environment in H. diversicolor.     

Enzymatic hydrolysis of cellulose by the cellobiohydrolase domain of CelB from the hyperthermophilic bacterium Caldicellulosiruptor saccharolyticus

The celB gene of Caldicellulosiruptor saccharolyticus was cloned and expressed in Escherichia coli to create a recombinant biocatalyst for hydrolyzing lignocellulosic biomass at high temperature. The GH5 domain of CelB hydrolyzed 4-nitrophenyl-β-d-cellobioside and carboxymethyl cellulose with optimum activity at pH 4.7-5.5 and 80 °C. The recombinant GH5 and CBM3-GH5 constructs were both stable at 80 °C with half-lives of 23 h and 39 h, respectively, and retained >94% activity after 48 h at 70 °C. Enzymatic hydrolysis of corn stover and cellulose pretreated with the ionic liquid 1-ethyl-3-methylimidazolium acetate showed that GH5 and CBM3-GH5 primarily produce cellobiose, with product yields for CBM3-GH5 being 1.2- to 2-fold higher than those for GH5. Confocal microscopy of bound protein on cellulose confirmed tighter binding of CBM3-GH5 to cellulose than GH5, indicating that the enhancement of enzymatic activity on solid substrates may be due to the substrate binding activity of CBM3 domain.     

Biochemical characterization of a 27 kDa 1,3-β-d-glucanase from Trichoderma asperellum induced by cell wall of Rhizoctonia solani

Trichoderma asperellum produces two extracellular 1,3-β-d-glucanase upon induction with cell walls from Rhizoctonia solani. A minor 1,3-β-d-glucanase was purified to homogeneity by ion exchange chromatography on Q-Sepharose and gel filtration on Sephacryl S-100. A typical procedure provided 13.8-fold purification with 70% yield. SDS-PAGE of the purified enzyme showed a single protein band of molecular weight 27 kDa. The enzyme exhibited optimum catalytic activity at pH 3.6 and 45 °C. It was thermostable at 40 °C, and retained 75% activity after 60 min at 45 °C. The K_m and V_max values for 1,3-β-d-glucanase, using laminarin as substrate, were 0.323 mg ml⁻¹ and 0.315 U min⁻¹, respectively. The enzyme was strongly inhibited by Hg²⁺ and SDS. The enzyme was only active toward glucans containing β-1,3-linkages. Peptide sequences showed similarity with two endo-1,3(4)-β-d-glucanases from Aspergillus fumigatus Af293when compared against GenBank non-redundant database.     

The cellulose-binding domain of cellobiohydrolase Cel7A from Trichoderma reesei is also a thermostabilizing domain

The thermostability of cellobiohydrolase I Cel7A from Trichoderma reesei was investigated using dynamic light scattering. While the whole enzyme displayed a melting point of 59 °C, the catalytic domain obtained via papain-catalyzed proteolysis was shown to denature at 51 °C and the cellulose-binding domain (with linker attached) melted at 65-66 °C. This variation in individual melting temperatures is proposed to account for the full retention of binding capacity of Cel7A at 50 °C, along with a loss of catalytic activity observed for the catalytic domain alone. Thus, the cellulose-binding domain of Cel7A acts as a thermostabilizing domain for the enzyme. The effect of reducing agents on the protein melting behavior was also investigated.     

Order Statistics Theory of Unfolding of Multimeric Proteins

Dynamic force spectroscopy has become indispensable for the exploration of the mechanical properties of proteins. In force-ramp experiments, performed by utilizing a time-dependent pulling force, the peak forces for unfolding transitions in a multimeric protein (D)_N are used to map the free energy landscape for unfolding for a protein domain D. We show that theoretical modeling of unfolding transitions based on combining the observed first (f₁), second (f₂), …, N^th (f_N) unfolding forces for a protein tandem of fixed length N, and pooling the force data for tandems of different length, n₁ < n₂ < … < N, leads to an inaccurate estimation of the distribution of unfolding forces for the protein D, ψ_D(f). This problem can be overcome by using Order statistics theory, which, in conjunction with analytically tractable models, can be used to resolve the molecular characteristics that determine the unfolding micromechanics. We present a simple method of estimation of the parent distribution, ψ_D(f), based on analyzing the force data for a tandem (D)_n of arbitrary length n. Order statistics theory is exemplified through a detailed analysis and modeling of the unfolding forces obtained from pulling simulations of the monomer and oligomers of the all-β-sheet WW domain.     

Phenotypic variation in hatchling Chinese ratsnakes (Zaocys dhumnades) from eggs incubated at constant temperatures

We incubated eggs of the Chinese ratsnake Zaocys dhumnades at four constant temperatures (24, 27, 30 and 30 °C) to examine the effects of incubation temperature on hatching success and hatchling phenotypes. Incubation length increased nonlinearly as temperature decreased, with the mean incubation length being 76.7 d at 24 °C, 57.4 d at 27 °C, 47.3 d at 30 °C, and 44.1 d at 33 °C. Hatching successes were lower at the two extreme temperatures (69% at 24 °C, and 44% at 33 °C) than at the other two moderate temperatures (96% at 27 °C, and 93% at 30 °C). Incubation temperature affected nearly all hatchling traits examined in this study. Incubation of Z. dhumnades eggs at 33 °C resulted in production of smaller hatchlings that characteristically had less-developed carcasses but contained more unutilized yolks. Hatchlings from eggs incubated at 27 and 30 °C did not differ in any examined traits. Taking the rate of embryonic development, hatching success and hatchling phenotypes into account, we conclude that the temperature range optimal for incubation of Z. dhumnades eggs is narrower than the range of 24−33 °C but should be wider than the range of 27−30 °C.     

Construction and characterization of different fusion proteins between cellulases and β-glucosidase to improve glucose production and thermostability

A β-glucosidase from Clostridium cellulovorans (CcBG) was fused with one of three different types of cellulases from Clostridium thermocellum, including a cellulosomal endoglucanase CelD (CtCD), a cellulosomal exoglucanase CBHA (CtCA) and a non-cellulosomal endoglucanase Cel9I (CtC9I). Six bifunctional enzymes were constructed with either β-glucosidase or cellulase in the upstream. CtCD-CcBG showed the favorable specific activities on phosphoric acid swollen cellulose (PASC), an amorphous cellulose, with more glucose production (2 folds) and less cellobiose accumulation (3 folds) when compared with mixture of the single enzymes. Moreover, CtCD-CcBG had significantly improved thermal stability with a melting temperature (T_m) of 10.9 °C higher than that of CcBG (54.5 °C) based on the CD unfolding experiments. This bifunctional enzyme is thus useful in industrial application to convert cellulose to glucose.     

Myoglobin with chlorophyllous chromophores: Influence on protein stability

The stabilities of myoglobin, apo-myoglobin, and of two myoglobins with chlorophyllous chromophores (Zn-pheophorbide a and Zn-bacteriopheophorbide a), have been studied by thermal and chemical denaturation. With guanidinium chloride, the stability order is myoglobin > Zn-pheophorbide-myoglobin > Zn-bacteriopheophorbide-myoglobin ∼ apo-myoglobin. The thermal behavior is more complex. The transition temperature of thermal unfolding of the apoprotein (62.4 °C) is increased by Zn-pheophorbide a (83.9 °C) and Zn-bacteriopheophorbide a (82.6 °C) to a similar degree as by the native chromophore, heme (83.5 °C). The recovery with Zn-pheophorbide (92-98%) is even higher than with heme (74-76%), while with Zn-bacteriopheophorbide (40%) it is as low as with the apoprotein (42%). Recovery also depends on the rates of heating, and in particular the time spent at high temperatures. It is concluded that irreversibility of unfolding is related to loss of the chromophores, which are required for proper re-folding.     

Response of juvenile diamond-backed terrapins (Malaclemys terrapin) to an aquatic thermal gradient

In many ectotherms, selection of environmental thermal niches may positively affect growth, nutrient assimilation rates, immune system function, and ultimately survival. Temperature preference in some turtle species may be influenced by environmental conditions, including acclimation temperature. We tested for effects of acclimation temperature (22 °C, 27 °C) on the selected temperature and movement patterns of 14 juvenile Malaclemys terrapin (Reptilia: Emydidae) in an aquatic thermal gradient of 14–34 °C and in single-temperature (22 °C, 27 °C) control tests. Among 8–10 month old terrapins, acclimation temperature influenced activity and movement patterns but did not affect temperature selection. In thermal gradient and single-temperature control tests, turtles acclimated to 27 °C used more tank chambers and relocated between chambers significantly more frequently than individuals acclimated to 22 °C. However, acclimation temperature did not affect temperature selection: both 22- and 27 °C-acclimated turtles selected the warmest temperature (34 °C), and avoided the other temperatures available, during thermal gradient tests. These results suggest that young M. terrapin are capable of detecting small temperature increments and prefer warm temperatures that may positively influence growth and metabolism.