Heat-resistant cytosolic malate dehydrogenases (cMDHs) of thermophilic intertidal snails (genus Echinolittorina): protein underpinnings of tolerance to body temperatures reaching 55°C

Phylogenetic history and temperature adaptation contribute to structural and functional stability of proteins in marine mollusks

Teasing apart the influences of phylogenetic history from thermal adaptation is a focal challenge in understanding the factors driving change in protein stability. This study conducted comprehensive comparative analyses between the phylogenetic relationships and functional/structural stabilities at protein and mRNA levels of cytosolic malate dehydrogenase (cMDH) orthologs of 41 marine mollusks living at widely different environmental temperatures. At the protein level, a significant negative correlation between adaptation temperature and heat-induced movements of the cMDH backbone was found. The movement fluctuation of individual residue varied similarly among cMDH orthologs. At the mRNA level, the free energy that occurs during the formation of the ensemble of mRNA secondary structure was significantly positively correlated with adaptation temperature. The fraction of guanine and cytosine increased with adaptation temperature. The proportion of variance in adaptation temperature that can be explained by the thermal stability (R2) was decreased after phylogenetic generalized least squares but was almost significant at both protein and mRNA levels (P < 0.05). Those analyses reveal the phylogenetic influence on the thermal adaptation of species. Our findings indicated that multi-level analysis of orthologous proteins should be considered alongside phylogenetic history to permit the development of a more comprehensive understanding of protein thermal adaptation.

Computer-Aided Reconstruction and Application of Bacillus halodurans S7 Xylanase with Heat and Alkali Resistance

β-1,4-Endoxylanase is the most critical hydrolase for xylan degradation during lignocellulosic biomass utilization. However, its poor stability and activity in hot and alkaline environments hinder its widespread application. In this study, BhS7Xyl from Bacillus halodurans S7 was improved using a computer-aided design through isothermal compressibility (βT) perturbation engineering and by combining three thermostability prediction algorithms (ICPE-TPA). The best variant with remarkable improvement in specific activity, heat resistance (70 °C), and alkaline resistance (both pH 9.0 and 70 °C), R69F/E137M/E145L, exhibited a 4.9-fold increase by wild-type in specific activity (1368.6 U/mg), a 39.4-fold increase in temperature half-life (458.1 min), and a 57.6-fold increase in pH half-life (383.1 min). Furthermore, R69F/E137M/E145L was applied to the hydrolysis of agricultural waste (corncob and hardwood pulp) to efficiently obtain a higher yield of high-value xylooligosaccharides. Overall, the ICPE-TPA strategy has the potential to improve the functional performance of enzymes under extreme conditions for the high-value utilization of lignocellulosic biomass.

Roles of multi-level temperature-adaptive responses and microhabitat variation in establishing distributions of intertidal species

How intertidal species survive their harsh environment and how best to evaluate and forecast range shifts in species distribution are two important and closely related questions for intertidal ecologists and global change biologists. Adaptive variation in responses of organisms to environmental change across all levels of biological organization - from behavior to molecular systems - is of key importance in setting distribution patterns, yet studies often neglect the interactions of diverse types of biological variation (e.g. differences in thermal optima owing to genetic and acclimation-induced effects) with environmental variation, notably at the scale of microhabitats. Intertidal species have to cope with extreme and frequently changing thermal stress, and have shown high variation in thermal sensitivities and adaptive responses at different levels of biological organization. Here, I review the physiological and biochemical adaptations of intertidal species to environmental temperature on multiple spatial and temporal scales. With fine-scale datasets for the thermal limits of individuals and for environmental temperature variation at the microhabitat scale, we can map the thermal sensitivity for each individual in different microhabitats, and then scale up the thermal sensitivity analysis to the population level and, finally, to the species level by incorporating physiological traits into species distribution models. These more refined mechanistic models that include consideration of physiological variations have higher predictive power than models that neglect these variations, and they will be crucial to answering the questions posed above concerning adaptive mechanisms and the roles they play in governing distribution patterns in a rapidly changing world.

Exploring the adaptability of the secondary structure of mRNA to temperature in intertidal snails based on SHAPE experiments

RNA-based thermal regulation is an important strategy for organisms to cope with temperature changes. Inhabiting the intertidal rocky shore, a key interface of the ocean, atmosphere and terrestrial environments, intertidal species have developed variable thermal adaptation mechanisms; however, adaptions at the RNA level remain largely uninvestigated. To examine the relationship between mRNA structural stability and species distribution, in the present study, the secondary structure of cytosolic malate dehydrogenase (cMDH) mRNA of Echinolittorina malaccana, Echinolittorina radiata and Littorina brevicula was determined using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), and the change in folding free energy of formation (ΔGfold) was calculated. The results showed that ΔGfold increased as the temperature increased. The difference in ΔGfold (ΔΔGfold) between two specific temperatures (25 versus 0°C, 37 versus 0°C and 57 versus 0°C) differed among the three species, and the ΔΔGfold value of E. malaccana was significantly lower than those of E. radiata and L. brevicula. The number of stems of cMDH mRNA of the snails decreased with increasing temperature, and the breakpoint temperature of E. malaccana was the highest among these. The number of loops was also reduced with increasing temperature, while the length of the loop structure increased accordingly. Consequently, these structural changes can potentially affect the translational efficiency of mRNA. These results imply that there were interspecific differences in the thermal stability of RNA secondary structures in intertidal snails, and these differences may be related to snail distribution.

Three amino acid substitutions contributing to thermostability of phosphoglucose isomerase in the Glanville fritillary butterfly

Temperature is one of the most important environmental factors that affect organisms, especially ectotherms, due to its effects on protein stability. Understanding the general rules that govern thermostability changes in proteins to adapt high-temperature environments is crucial. Here, we report the amino acid substitutions of phosphoglucose isomerase (PGI) related to thermostability in the Glanville fritillary butterfly (Melitaea cinxia, Lepidoptera: Nymphalidae). The PGI encoded by the most common allele in M. cinxia in the Chinese population (G3-PGI), which is more thermal tolerant, is more stable under heat stress than that in the Finnish population (D1-PGI). There are 5 amino acid substitutions between G3-PGI and D1-PGI. Site-directed mutagenesis revealed that the combination of amino acid substitutions of H35Q, M49T, and I64V may increase PGI thermostability. These substitutions alter the 3D structure to increase the interaction between 2 monomers of PGI. Through molecular dynamics simulations, it was found that the amino acid at site 421 is more stable in G3-PGI, confining the motion of the α-helix 420-441 and stabilizing the interaction between 2 PGI monomers. The strategy for high-temperature adaptation through these 3 amino acid substitutions is also adopted by other butterfly species (Boloria eunomia, Aglais urticae, Colias erate, and Polycaena lua) concurrent with M. cinxia in the Tianshan Mountains of China, i.e., convergent evolution in butterflies.

A new ALK inhibitor overcomes resistance to first- and second-generation inhibitors in NSCLC

More than 60% of nonsmall cell lung cancer (NSCLC) patients show a positive response to the first ALK inhibitor, crizotinib, which has been used as the standard treatment for newly diagnosed patients with ALK rearrangement. However, most patients inevitably develop crizotinib resistance due to acquired secondary mutations in the ALK kinase domain, such as the gatekeeper mutation L1196M and the most refractory mutation, G1202R. Here, we develop XMU-MP-5 as a new-generation ALK inhibitor to overcome crizotinib resistance mutations, including L1196M and G1202R. XMU-MP-5 blocks ALK signaling pathways and inhibits the proliferation of cells harboring either wild-type or mutant EML4-ALK in vitro and suppresses tumor growth in xenograft mouse models in vivo. Structural analysis provides insights into the mode of action of XMU-MP-5. In addition, XMU-MP-5 induces significant regression of lung tumors in two genetically engineered mouse (GEM) models, further demonstrating its pharmacological efficacy and potential for clinical application. These preclinical data support XMU-MP-5 as a novel selective ALK inhibitor with high potency and selectivity. XMU-MP-5 holds great promise as a new therapeutic against clinically relevant secondary ALK mutations.

The Goldilocks Principle: A Unifying Perspective on Biochemical Adaptation to Abiotic Stressors in the Sea

The ability of marine organisms to thrive over wide ranges of environmental stressors that perturb structures of proteins, nucleic acids, and lipids illustrates the effectiveness of adaptation at the biochemical level. A critical role of these adaptations is to achieve a proper balance between structural rigidity, which is necessary for maintaining three-dimensional conformation, and flexibility, which is required to allow changes in conformation during function. The Goldilocks principle refers to this balancing act, wherein structural stability and functional properties are poised at values that are just right for the environment the organism faces. Achieving this balance involves changes in macromolecular sequence and adaptive change in the composition of the aqueous or lipid milieu in which macromolecules function. This article traces the development of the field of biochemical adaptation throughout my career and shows how comparative studies of marine animals from diverse habitats have shed light on fundamental properties of life common to all organisms.

Thermal adaptation of mRNA secondary structure: stability versus lability

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species' normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability-lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔG_fold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔG_fold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.

Physiological determinants of biogeography: The importance of metabolic depression to heat tolerance

A quantitative understanding of physiological thermal responses is vital for forecasting species distributional shifts in response to climate change. Many studies have focused on metabolic rate as a global metric for analyzing the sublethal effects of changing environments on physiology. Thermal performance curves (TPCs) have been suggested as a viable analytical framework, but standard TPCs may not fully capture physiological responses, due in part to failure to consider the process of metabolic depression. We derived a model based on the nonlinear regression of biological temperature-dependent rate processes and built a heart rate data set for 26 species of intertidal molluscs distributed from 33°S to ~40°N. We then calculated physiological thermal performance limits with continuous heating using T 1 / 2 H , the temperature at which heart rate is decreased to 50% of the maximal rate, as a more realistic measure of upper thermal limits. Results indicate that heat-induced metabolic depression of cardiac performance is a common adaptive response that allows tolerance of harsh environments. Furthermore, our model accounted for the high inter-individual variability in the shape of cardiac TPCs. We then used these TPCs to calculate physiological thermal safety margins (pTSM), the difference between the maximal operative temperature (95th percentile of field temperatures) and T 1 / 2 H of each individual. Using pTSMs, we developed a physiological species distribution model (pSDM) to forecast future geographic distributions. pSDM results indicate that climate-induced species range shifts are potentially less severe than predicted by a simple correlative SDM. Species with metabolic depression below the optimum temperature will be more thermal resistant at their warm trailing edges. High intraspecific variability further suggests that models based on species-level vulnerability to environmental change may be problematic. This multi-scale, mechanistic understanding that incorporates metabolic depression and inter-individual variability in thermal response enables better predictions about the relationship between thermal stress and species distributions.

Behavioral repertoire of high-shore littorinid snails reveals novel adaptations to an extreme environment

Species that inhabit high-shore environments on rocky shores survive prolonged periods of emersion and thermal stress. Using two Hong Kong high-shore littorinids (Echinolittorina malaccana and E. radiata) as models, we examined their behavioral repertoire to survive these variable and extreme conditions. Environmental temperatures ranged from 4°C in the cool season to 55.5°C in the hot season, with strong seasonal and daily fluctuations. In the hot season, both species allocated >35% of their activity budgets to stress-mitigating thermoregulatory behaviors (e.g. standing, towering) and relatively small proportions to foraging (<20%) and reproduction (<10%). In the assumedly benign cool season, greater proportions (>70%) of activity budgets were allocated to stress mitigation behaviors (crevice occupation, aggregation formation). Both species exhibited multifunctional behaviors that optimized time use during their tidally-constrained activity window in the hot season. Females mated while foraging when awash by the rising tide, and some males crawled on top of females prior to ceasing movement to form 'towers', which have both thermoregulatory benefits and reduce searching time for mates during subsequent activity. The function of such behaviors varies in a state-dependent manner, for example, the function of trail following changes over an activity cycle from mate searching on rising tides, to stress mitigation on falling tides (aiding aggregation formation), and to both functions through tower formation just before movement stops. Many of these behavioral responses are, therefore, multifunctional and can vary according to local conditions, allowing snails in this family to successfully colonize the extreme high-shore environment.

Temperature adaptations of the thermophilic snail Echinolittorina malaccana: insights from metabolomic analysis

The periwinkle snail Echinolittorina malaccana, for which the upper lethal temperature is near 55°C, is one of the most heat-tolerant eukaryotes known. We conducted a multi-level investigation - including cardiac physiology, enzyme activity, and targeted and untargeted metabolomic analyses - that elucidated a spectrum of adaptations to extreme heat in this organism. All systems examined showed heat intensity-dependent responses. Under moderate heat stress (37-45°C), the snail depressed cardiac activity and entered a state of metabolic depression. The global metabolomic and enzymatic analyses revealed production of metabolites characteristic of oxygen-independent pathways of ATP generation (lactate and succinate) in the depressed metabolic state, which suggests that anaerobic metabolism was the main energy supply pathway under heat stress (37-52°C). The metabolomic analyses also revealed alterations in glycerophospholipid metabolism under extreme heat stress (52°C), which likely reflected adaptive changes to maintain membrane structure. Small-molecular-mass organic osmolytes (glycine betaine, choline and carnitine) showed complex changes in concentration that were consistent with a role of these protein-stabilizing solutes in protection of the proteome under heat stress. This thermophilic species can thus deploy a wide array of adaptive strategies to acclimatize to extremely high temperatures.

Impact of short-term elevated temperature stress on winter-acclimated individuals of the marine gastropod Crepidula fornicata

The intertidal zone is an especially stressful thermal habitat, typically exposing residents to air temperatures for up to 6 h at a time, twice daily. Tolerance to elevated temperatures has been particularly well-studied for a variety of intertidal species, especially with regard to upper thermal limits during summers. However, in recent years, as climates have been changing around the world, temperate zone intertidal organisms have sometimes been exposed to periods of unusually high air temperatures during the winter. The present study sought to examine the impact of elevated temperatures on survival and clearance rates of winter-acclimated intertidal individuals of the sedentary marine suspension-feeding gastropod Crepidula fornicata. Individuals were collected intertidally from Nahant, Massachusetts from late January to early April each year for 5 years, maintained in the laboratory at the acclimation temperature of 6 °C, and exposed in the laboratory for 3 h to temperatures as high as 37 °C in seawater either once or twice, 24 h apart. Although mean clearance rates were substantially reduced for at least the next 12-24 h after individuals were returned to the 6 °C control condition following exposures to elevated temperatures as low as 21-26 °C, we saw little mortality even following two 3 h exposures to 35 °C, or single exposures to 37 °C. Mortality was substantial, however, in one experiment following a double exposure to 37 °C. Smaller individuals (~5-12 mm in shell length) were somewhat more sensitive to the thermal stress than adults were. Intertidal members of C. fornicata in Massachusetts seem well-prepared to deal with the increasing range of winter air temperatures associated with the global climate confusion predicted for future years. Additional studies will be required to understand the physiological and biochemical mechanisms used by winter-acclimated individuals of this species to tolerate such periodic substantial temperature increases of 29-31 °C.

Thermal adaptation in the honeybee (Apis mellifera) via changes to the structure of malate dehydrogenase

In honeybees there are three alleles of cytosolic malate dehydrogenase gene: F, M and S. Allele frequencies are correlated with environmental temperature, suggesting that the alleles have temperature-dependent fitness benefits. We determined the enzyme activity of each allele across a range of temperatures in vitro The F and S alleles have higher activity and are less sensitive to high temperatures than the M allele, which loses activity after incubation at temperatures found in the thorax of foraging bees in hot climates. Next, we predicted the protein structure of each allele and used molecular dynamics simulations to investigate their molecular flexibility. The M allozyme is more flexible than the S and F allozymes at 50°C, suggesting a plausible explanation for its loss of activity at high temperatures, and has the greatest structural flexibility at 15°C, suggesting that it can retain some enzyme activity at cooler temperatures. MM bees recovered from 2 h of cold narcosis significantly better than all other genotypes. Combined, these results explain clinal variation in malate dehydrogenase allele frequencies in the honeybee at the molecular level.

High-throughput quantification of protein structural change reveals potential mechanisms of temperature adaptation in Mytilus mussels

Background

Temperature exerts a strong influence on protein evolution: species living in thermally distinct environments often exhibit adaptive differences in protein structure and function. However, previous research on protein temperature adaptation has focused on small numbers of proteins and on proteins adapted to extreme temperatures. Consequently, less is known about the types and quantity of evolutionary change that occurs to proteins when organisms adapt to small shifts in environmental temperature. In this study, these uncertainties were addressed by developing software that enabled comparison of structural changes associated with temperature adaptation (hydrogen bonding, salt bridge formation, and amino acid use) among large numbers of proteins from warm- and cold-adapted species of marine mussels, Mytilus galloprovincialis and Mytilus trossulus, respectively.

Results

Small differences in habitat temperature that characterize the evolutionary history of Mytilus mussels were sufficient to cause protein structural changes consistent with temperature adaptation. Hydrogen bonds and salt bridges that increase stability and protect against heat-induced denaturation were more abundant in proteins from warm-adapted M. galloprovincialis compared with proteins from cold-adapted M. trossulus. These structural changes were related to deviations in the use of polar and charged amino acids that facilitate formation of hydrogen bonds and salt bridges within proteins, respectively. Enzymes, in particular those within antioxidant and cell death pathways, were over-represented among proteins with the most hydrogen bonds and salt bridges in warm-adapted M. galloprovincialis. Unlike extremophile proteins, temperature adaptation in Mytilus proteins did not involve substantial changes in the number of hydrophobic or large volume amino acids, nor in the content of glycine or proline.

Conclusions

Small shifts in organism temperature tolerance, such as that needed to cope with climate warming, may result from structural and functional changes to a small percentage of the proteome. Proteins in which function is dependent on large conformational change, notably enzymes, may be particularly sensitive to temperature perturbation and represent foci for natural selection. Protein temperature adaptation can occur through different types and frequencies of structural change, and adaptive mechanisms used to cope with small shifts in habitat temperature appear different from mechanisms used to retain protein function at temperature extremes.

The Place for Enzymes and Biologically Active Peptides from Marine Organisms for Application in Industrial and Pharmaceutical Biotechnology

Since the beginning of written history, diverse texts have reported the use of enzymatic preparations in food processing and have described the medicinal properties of crude and fractionated venoms to treat various diseases and injuries. With the biochemical characterization of enzymes from distinct sources and bioactive polypeptides from animal venoms, the last sixty years have testified the advent of industrial enzymology and protein therapeutics, which are currently applicable in a wide variety of industrial processes, household products, and pharmaceuticals. Bioprospecting of novel biocatalysts and bioactive peptides is propelled by their unsurpassed properties that are applicable for current and future green industrial processes, biotechnology, and biomedicine. The demand for both novel enzymes with desired characteristics and novel peptides that lead to drug development, has experienced a steady increase in response to the expanding global market for industrial enzymes and peptidebased drugs. Moreover, although largely unexplored, oceans and marine realms, with their unique ecosystems inhabited by a large variety of species, including a considerable number of venomous animals, are recognized as untapped reservoirs of molecules and macromolecules (enzymes and bioactive venom-derived peptides) that can potentially be converted into highly valuable biopharmaceutical products. In this review, we have focused on enzymes and animal venom (poly)peptides that are presently in biotechnological use, and considering the state of prospection of marine resources, on the discovery of useful industrial biocatalysts and drug leads with novel structures exhibiting selectivity and improved performance.

Non-reversible and Reversible Heat Tolerance Plasticity in Tropical Intertidal Animals: Responding to Habitat Temperature Heterogeneity

The theory for thermal plasticity of tropical ectotherms has centered on terrestrial and open-water marine animals which experience reduced variation in diurnal and seasonal temperatures, conditions constraining plasticity selection. Tropical marine intertidal animals, however, experience complex habitat thermal heterogeneity, circumstances encouraging thermal plasticity selection. Using the tropical rocky-intertidal gastropod, Echinolittorina malaccana, we investigated heat tolerance plasticity in terms of laboratory acclimation and natural acclimatization of populations from thermally-dissimilar nearby shorelines. Laboratory treatments yielded similar capacities of snails from either population to acclimate their lethal thermal limit (LT50 variation was ∼2°C). However, the populations differed in the temperature range over which acclimatory adjustments could be made; LT50 plasticity occurred over a higher temperature range in the warm-shore snails compared to the cool-shore snails, giving an overall acclimation capacity for the populations combined of 2.9°C. In addition to confirming significant heat tolerance plasticity in tropical intertidal animals, these findings reveal two plasticity forms, reversible (laboratory acclimation) and non-reversible (population or shoreline specific) plasticity. The plasticity forms should account for different spatiotemporal scales of the environmental temperature variation; reversible plasticity for daily and tidal variations in microhabitat temperature and non-reversible plasticity for lifelong, shoreline temperature conditions. Non-reversible heat tolerance plasticity, likely established after larvae settle on the shore, should be energetically beneficial in preventing heat shock protein overexpression, but also should facilitate widespread colonization of coasts that support thermally-diverse shorelines. This first demonstration of different plasticity forms in benthic intertidal animals supports the hypothesis that habitat heterogeneity (irrespective of latitude) drives thermal plasticity selection. It further suggests that studies not making reference to different spatial scales of thermal heterogeneity, nor seeking how these may drive different thermal plasticity forms, risk misinterpreting ectothermic responses to environmental warming.

Comparing mutagenesis and simulations as tools for identifying functionally important sequence changes for protein thermal adaptation

Comparative studies of orthologous proteins of species evolved at different temperatures have revealed consistent patterns of temperature-related variation in thermal stabilities of structure and function. However, the precise mechanisms by which interspecific variations in sequence foster these adaptive changes remain largely unknown. Here, we compare orthologs of cytosolic malate dehydrogenase (cMDH) from marine molluscs adapted to temperatures ranging from -1.9 °C (Antarctica) to ∼55 °C (South China coast) and show how amino acid usage in different regions of the enzyme (surface, intermediate depth, and protein core) varies with adaptation temperature. This eukaryotic enzyme follows some but not all of the rules established in comparisons of archaeal and bacterial proteins. To link the effects of specific amino acid substitutions with adaptive variations in enzyme thermal stability, we combined site-directed mutagenesis (SDM) and in vitro protein experimentation with in silico mutagenesis using molecular dynamics simulation (MDS) techniques. SDM and MDS methods generally but not invariably yielded common effects on protein stability. MDS analysis is shown to provide insights into how specific amino acid substitutions affect the conformational flexibilities of mobile regions (MRs) of the enzyme that are essential for binding and catalysis. Whereas these substitutions invariably lie outside of the MRs, they effectively transmit their flexibility-modulating effects to the MRs through linked interactions among surface residues. This discovery illustrates that regions of the protein surface lying outside of the site of catalysis can help establish an enzyme's thermal responses and foster evolutionary adaptation of function.

Substantial heat tolerance acclimation capacity in tropical thermophilic snails, but to what benefit?

The theory for thermal acclimation of ectotherms suggests that (1) heat tolerance is traded off for thermal acclimation in thermophilic species and that (2) plasticity is constrained in tropically distributed ectotherms, which commonly experience relatively thermally stable environments. We observed substantial heat tolerance plasticity in a test of this theory using tropical, thermophilic marine intertidal snails that inhabit extremely hot and highly variable thermal environments. The implication of these results is that plasticity selection is largely driven by habitat temperature conditions irrespective of basal heat tolerance or latitude. However, heat tolerance of field-fresh snails was comparable with that of laboratory warm-acclimated snails, suggesting that snails in the field may often be unable to improve heat hardening with further environmental warming. These findings suggest that field referencing is crucial to using laboratory-measured acclimation capacity when inferring climate warming vulnerability in ectotherms, and overall they question how well current thermal biology theory predicts the outcomes of global change in intertidal environments.

Structural flexibility and protein adaptation to temperature: Molecular dynamics analysis of malate dehydrogenases of marine molluscs

Orthologous proteins of species adapted to different temperatures exhibit differences in stability and function that are interpreted to reflect adaptive variation in structural "flexibility." However, quantifying flexibility and comparing flexibility across proteins has remained a challenge. To address this issue, we examined temperature effects on cytosolic malate dehydrogenase (cMDH) orthologs from differently thermally adapted congeners of five genera of marine molluscs whose field body temperatures span a range of ∼60 °C. We describe consistent patterns of convergent evolution in adaptation of function [temperature effects on K_M of cofactor (NADH)] and structural stability (rate of heat denaturation of activity). To determine how these differences depend on flexibilities of overall structure and of regions known to be important in binding and catalysis, we performed molecular dynamics simulation (MDS) analyses. MDS analyses revealed a significant negative correlation between adaptation temperature and heat-induced increase of backbone atom movements [root mean square deviation (rmsd) of main-chain atoms]. Root mean square fluctuations (RMSFs) of movement by individual amino acid residues varied across the sequence in a qualitatively similar pattern among orthologs. Regions of sequence involved in ligand binding and catalysis-termed mobile regions 1 and 2 (MR1 and MR2), respectively-showed the largest values for RMSF. Heat-induced changes in RMSF values across the sequence and, importantly, in MR1 and MR2 were greatest in cold-adapted species. MDS methods are shown to provide powerful tools for examining adaptation of enzymes by providing a quantitative index of protein flexibility and identifying sequence regions where adaptive change in flexibility occurs.