Neutron Scattering Reveals the Dynamic Basis of Protein Adaptation to Extreme Temperature

The Characterization of Ancient Methanococcales Malate Dehydrogenases Reveals That Strong Thermal Stability Prevents Unfolding Under Intense γ-Irradiation

Malate dehydrogenases (MalDHs) (EC.1.1.1.37), which are involved in the conversion of oxaloacetate to pyruvate in the tricarboxylic acid cycle, are a relevant model for the study of enzyme evolution and adaptation. Likewise, a recent study showed that Methanococcales, a major lineage of Archaea, is a good model to study the molecular processes of proteome thermoadaptation in prokaryotes. Here, we use ancestral sequence reconstruction and paleoenzymology to characterize both ancient and extant MalDHs. We observe a good correlation between inferred optimal growth temperatures and experimental optimal temperatures for activity (A-Topt). In particular, we show that the MalDH present in the ancestor of Methanococcales was hyperthermostable and had an A-Topt of 80 °C, consistent with a hyperthermophilic lifestyle. This ancestor gave rise to two lineages with different thermal constraints: one remained hyperthermophilic, while the other underwent several independent adaptations to colder environments. Surprisingly, the enzymes of the first lineage have retained a thermoresistant behavior (i.e. strong thermostability and high A-Topt), whereas the ancestor of the second lineage shows a strong thermostability, but a reduced A-Topt. Using mutants, we mimic the adaptation trajectory toward mesophily and show that it is possible to significantly reduce the A-Topt without altering the thermostability of the enzyme by introducing a few mutations. Finally, we reveal an unexpected link between thermostability and the ability to resist γ-irradiation-induced unfolding.

Dynamic Personality of Proteins and Effect of the Molecular Environment

Protein dynamics display distinct traits that are linked to their specific biological function. However, the interplay between intrinsic dynamics and the molecular environment on protein stability remains poorly understood. In this study, we investigate, by incoherent neutron scattering, the subnanosecond time scale dynamics of three model proteins: the mesophilic lysozyme, the thermophilic thermolysin, and the intrinsically disordered β-casein. Moreover, we address the influence of water, glycerol, and glucose, which create progressively more viscous matrices around the protein surface. By comparing the protein thermal fluctuations, we find that the internal dynamics of thermolysin are less affected by the environment compared to lysozyme and β-casein. We ascribe this behavior to the protein dynamic personality, i.e., to the stiffer dynamics of the thermophilic protein that contrasts the influence of the environment. Remarkably, lysozyme and thermolysin in all molecular environments reach a critical common flexibility when approaching the calorimetric melting temperature.

Decoding the Role of the Global Proteome Dynamics for Cellular Thermal Stability

Molecular mechanisms underlying the thermal response of cells remain elusive. On the basis of the recent result that the short-time diffusive dynamics of the Escherichia coli proteome is an excellent indicator of temperature-dependent bacterial metabolism and death, we used neutron scattering (NS) spectroscopy and molecular dynamics (MD) simulations to investigate the sub-nanosecond proteome mobility in psychro-, meso-, and hyperthermophilic bacteria over a wide temperature range. The magnitude of thermal fluctuations, measured by atomic mean square displacements, is similar among all studied bacteria at their respective thermal cell death. Global roto-translational motions turn out to be the main factor distinguishing the bacterial dynamical properties. We ascribe this behavior to the difference in the average proteome net charge, which becomes less negative for increasing bacterial thermal stability. We propose that the chemical-physical properties of the cytoplasm and the global dynamics of the resulting proteome are fine-tuned by evolution to uphold optimal thermal stability conditions.

Marine enzymes: Classification and application in various industries

Oceans are regarded as a plentiful and sustainable source of biological compounds. Enzymes are a group of marine biomaterials that have recently drawn more attention because they are produced in harsh environmental conditions such as high salinity, extensive pH, a wide temperature range, and high pressure. Hence, marine-derived enzymes are capable of exhibiting remarkable properties due to their unique composition. In this review, we overviewed and discussed characteristics of marine enzymes as well as the sources of marine enzymes, ranging from primitive organisms to vertebrates, and presented the importance, advantages, and challenges of using marine enzymes with a summary of their applications in a variety of industries. Current biotechnological advancements need the study of novel marine enzymes that could be applied in a variety of ways. Resources of marine enzyme can benefit greatly for biotechnological applications duo to their biocompatible, ecofriendly and high effectiveness. It is beneficial to use the unique characteristics offered by marine enzymes to either develop new processes and products or improve existing ones. As a result, marine-derived enzymes have promising potential and are an excellent candidate for a variety of biotechnology applications and a future rise in the use of marine enzymes is to be anticipated.

Temperature-resistant and solvent-tolerant lipases as industrial biocatalysts: Biotechnological approaches and applications

The development of new biocatalytic systems to replace the chemical catalysts, with suitable characteristics in terms of efficiency, stability under high temperature reactions and in the presence of organic solvents, reusability, and eco-friendliness is considered a very important step to move towards the green processes. From this basis, the use of lipase as a catalyst is highly desired for many industrial applications because it offers the reactions in which could be used, stability in harsh conditions, reusability and a greener process. Therefore, the introduction of temperature-resistant and solvent-tolerant lipases have become essential and ideal for industrial applications. Temperature-resistant and solvent-tolerant lipases have been involved in many large-scale applications including biodiesel, detergent, food, pharmaceutical, organic synthesis, biosensing, pulp and paper, textile, animal feed, cosmetics, and leather industry. So, the present review provides a comprehensive overview of the industrial use of lipase. Moreover, special interest in biotechnological and biochemical techniques for enhancing temperature-resistance and solvent-tolerance of lipases to be suitable for the industrial uses.

Improving the kinetic stability of a hyperthermostable β-mannanase by a rationally combined strategy

Feasible and easily accessible methods for the rational design of enzyme engineering strategies remain to be established. Thus, a new rationally combined strategy based on disulfide bond engineering and HotSpot Wizard 3.0 was proposed and experimentally demonstrated to be effective using a hyperthermostable β-mannanase. Ten of 42 mutants showed prominent enhancement of kinetic stability with 26.4%-39.9% increases in t_1/2 (75 °C) compared with the parent enzyme ManAKH. The best mutant, D273-V308, showed apparent increases in both optimal temperature (5 °C) and T₅₀ (6.8 °C), as well as advanced catalytic efficiency. The low rate of inactive mutants and the high rate of positive mutants indicated that newly introduced screening factors (distance from catalytic residues, Gibbs free energy term, molecular simulation, and visual inspections) greatly enhance the design of thermostable β-mannanase. Moreover, these findings further advance the industrial application of β-mannanase (ManAK) in food and food-related applications.

Structural changes induced by ligand binding drastically increase the thermostability of the Ser/Thr protein kinase TpkD from Thermus thermophilus HB8

Thermophilic proteins maintain their structure at high temperatures through a combination of various factors. Here, we report the ligand-induced stabilization of a thermophilic Ser/Thr protein kinase. Thermus thermophilus TpkD unfolds completely at 55 °C despite the optimum growth temperature of 75 °C. Unexpectedly, we found that the TpkD structure is drastically stabilized by its natural ligands ATP and ADP, as evidenced by the increase in the melting temperature to 80 °C. Such a striking effect of a substrate on thermostability has not been reported for other protein kinases. Conformational changes upon ATP binding were observed in fluorescence quenching and limited proteolysis experiments. Urea denaturation of Trp mutants suggested that ATP binding affects not only the ATP-binding site, but also the remote regions. Our findings shed light on thermoadaptation of thermophilic proteins.

Revealing the Structural Contributions to Thermal Adaptation of the TATA-Box Binding Protein: Molecular Dynamics and QSPR Analyses

The TATA-box binding protein (TBP) is an important element of the transcription machinery in archaea and eukaryotic organisms. TBP is expressed in organisms adapted to different temperatures, indicating a robust structure, and experimental studies have shown that the mid-unfolding temperature (Tm) of TBP is directly correlated with the optimal growth temperature (OGT) of the organism. To understand which are the relevant structural requirements for its stability, we present the first structural and dynamic computational study of TBPs, combining molecular dynamics (MD) simulations and a quantitative structure-property relationship (QSPR) over a set of TBPs of organisms adapted to different temperatures. We found that the main structural properties of TBP used to adapt to high temperatures are an increase in the ease of desolvation of charged residues at the surface, an increase in the local resiliency, the presence of Leu clusters in the protein core, and an increase in the loss of hydrophobic packing in the N-terminal subdomain. In view of our results, we consider that TBP is a good model to study thermal adaptation, and our analysis opens the possibility of performing protein engineering on TBPs to study transcription at high or low temperatures.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Properties and Applications of Extremozymes from Deep-Sea Extremophilic Microorganisms: A Mini Review

The deep sea, which is defined as sea water below a depth of 1000 m, is one of the largest biomes on the Earth, and is recognised as an extreme environment due to its range of challenging physical parameters, such as pressure, salinity, temperature, chemicals and metals (such as hydrogen sulphide, copper and arsenic). For surviving in such extreme conditions, deep-sea extremophilic microorganisms employ a variety of adaptive strategies, such as the production of extremozymes, which exhibit outstanding thermal or cold adaptability, salt tolerance and/or pressure tolerance. Owing to their great stability, deep-sea extremozymes have numerous potential applications in a wide range of industries, such as the agricultural, food, chemical, pharmaceutical and biotechnological sectors. This enormous economic potential combined with recent advances in sampling and molecular and omics technologies has led to the emergence of research regarding deep-sea extremozymes and their primary applications in recent decades. In the present review, we introduced recent advances in research regarding deep-sea extremophiles and the enzymes they produce and discussed their potential industrial applications, with special emphasis on thermophilic, psychrophilic, halophilic and piezophilic enzymes.

Ternary Complex Formation and Photoactivation of a Photoenzyme Results in Altered Protein Dynamics

The interplay between protein dynamics and catalysis remains a fundamental question in enzymology. We here investigate the ns-timescale dynamics of a light-dependent NADPH:protochlorophyllide oxidoreductase (LPOR), a photoenzyme crucial for chlorophyll synthesis. LPORs catalyze the light-triggered trans addition of a hydride and a proton across the C17═C18 double bond of the chlorophyll precursor protochlorophyllide (Pchlide). Because of the lack of an LPOR structure, the global structural and dynamic consequences of LPOR/Pchlide/NADPH ternary complex formation remain elusive. Moreover, photoactivation of LPORs by low-light preillumination is controversially discussed as unequivocal proof for this phenomenon is lacking. By employing quasielastic neutron spectroscopy (QENS), we show that the formation of the ternary holoprotein complex as well as photoactivation lead to progressive rigidification of the protein. These findings are supported by thermostability measurements, which reveal different melting behavior and thermostabilities for the apo- and holoprotein ternary complexes. Molecular dynamics simulations in good agreement with the experimental QENS results suggest that the increased flexibility observed for the apoprotein stems from structural fluctuations of the NADPH and Pchlide substrate binding sites of the enzyme. On the basis of our results, in conjunction with activity and stability measurements, we provide independent proof for LPOR photoactivation, defined as a process that modifies the protein structure and dynamics, resulting in an increased substrate turnover. Our findings advance the structural and dynamic understanding of LPORs and provide a first link between protein dynamics and catalysis for this enzyme class.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Pulse EPR of Triarylmethyl Probes: A New Approach for the Investigation of Molecular Motions in Soft Matter

Triarylmethyl (TAM) radicals have become widely used free radicals in the past few years. Their electron spins have long relaxation times and narrow electron paramagnetic resonance (EPR) lines, which make them an important class of probes and tags in biological applications and materials science. In this work, we propose a new approach to characterize librations by means of TAM radicals. The temperature dependence of motional parameter ⟨α²⟩τ_c, where ⟨α²⟩ is the mean-squared amplitude of librations and τ_c is their characteristic time, is obtained by comparison of the 1/ T_m phase-relaxation rates at X- and Q-band EPR frequencies. We study three soft matrixes, viz., glassy trehalose and two ionic liquids, using TAMs with optimized relaxation properties OX063D and a dodeca- n-butyl homologue of Finland trityl (DBT). The motional parameters ⟨α²⟩τ_c obtained using TAMs are in excellent agreement with those obtained by means of nitroxide radicals. At the same time, the new TAM-based approach has (1) greater sensitivity due to the narrower EPR spectrum and (2) greater measuring accuracy and broader temperature range due to longer relaxation times. The developed approach may be fruitfully implemented to probe low-temperature molecular motions of TAM-labeled biopolymers, membrane systems, polymers, molecules in glassy media, and ionic liquids.

Trehalose Mediated Inhibition of Lactate Dehydrogenase from Rabbit Muscle. The Application of Kramers' Theory in Enzyme Catalysis

Lactate dehydrogenase (LDH) catalyzes the reduction of pyruvate to lactate by using NADH. LDH kinetics has been proposed to be dependent on the dynamics of a loop over the active site. Kramers' theory has been useful in the study of enzyme catalysis dependent on large structural dynamics. In this work, LDH kinetics was studied in the presence of trehalose and at different temperatures. In the absence of trehalose, temperature increase raised exponentially the LDH V_max and revealed a sigmoid transition of K_m toward a low-affinity state similar to protein unfolding. Notably, LDH V_max diminished when in the presence of trehalose, while pyruvate affinity increased and the temperature-mediated binding site transition was hindered. The effect of trehalose on k_cat was viscosity dependent as described by Kramers' theory since V_max correlated inversely with the viscosity of the medium. As a result, activation energy ( E_a) for pyruvate reduction was dramatically increased by trehalose presence. This work provides experimental evidence that the dynamics of a structural component in LDH is essential for catalysis, i.e., the closing of the loop on the active site. While the trehalose mediated-increased of pyruvate affinity is proposed to be due to the compaction and/or increase of structural order at the binding site.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

Molecular adaptation and salt stress response of Halobacterium salinarum cells revealed by neutron spectroscopy

Halobacterium salinarum is an extreme halophile archaeon with an absolute requirement for a multimolar salt environment. It accumulates molar concentrations of KCl in the cytosol to counterbalance the external osmotic pressure imposed by the molar NaCl. As a consequence, cytosolic proteins are permanently exposed to low water activity and highly ionic conditions. In non-adapted systems, such conditions would promote protein aggregation, precipitation, and denaturation. In contrast, in vitro studies showed that proteins from extreme halophilic cells are themselves obligate halophiles. In this paper, adaptation via dynamics to low-salt stress in H. salinarum cells was measured by neutron scattering experiments coupled with microbiological characterization. The molecular dynamic properties of a proteome represent a good indicator for environmental adaptation and the neutron/microbiology approach has been shown to be well tailored to characterize these modifications. In their natural setting, halophilic organisms often have to face important variations in environmental salt concentration. The results showed deleterious effects already occur in the H. salinarum proteome, even when the external salt concentration is still relatively high, suggesting the onset of survival mechanisms quite early when the environmental salt concentration decreases.

Adaptation of Extremophilic Proteins with Temperature and Pressure: Evidence from Initiation Factor 6

In this work, we study dynamical properties of an extremophilic protein, Initiation Factor 6 (IF6), produced by the archeabacterium Methanocaldococcus jannascii, which thrives close to deep-sea hydrothermal vents where temperatures reach 80 °C and the pressure is up to 750 bar. Molecular dynamics simulations (MD) and quasi-elastic neutron scattering (QENS) measurements give new insights into the dynamical properties of this protein with respect to its eukaryotic and mesophilic homologue. Results obtained by MD are supported by QENS data and are interpreted within the framework of a fractional Brownian dynamics model for the characterization of protein relaxation dynamics. IF6 from M. jannaschii at high temperature and pressure shares similar flexibility with its eukaryotic homologue from S. cerevisieae under ambient conditions. This work shows for the first time, to our knowledge, that the very common pattern of corresponding states for thermophilic protein adaptation can be extended to thermo-barophilic proteins. A detailed analysis of dynamic properties and of local structural fluctuations reveals a complex pattern for "corresponding" structural flexibilities. In particular, in the case of IF6, the latter seems to be strongly related to the entropic contribution given by an additional, C-terminal, 20 amino-acid tail which is evolutionary conserved in all mesophilic IF6s.

Marine extremophiles: a source of hydrolases for biotechnological applications

The marine environment covers almost three quarters of the planet and is where evolution took its first steps. Extremophile microorganisms are found in several extreme marine environments, such as hydrothermal vents, hot springs, salty lakes and deep-sea floors. The ability of these microorganisms to support extremes of temperature, salinity and pressure demonstrates their great potential for biotechnological processes. Hydrolases including amylases, cellulases, peptidases and lipases from hyperthermophiles, psychrophiles, halophiles and piezophiles have been investigated for these reasons. Extremozymes are adapted to work in harsh physical-chemical conditions and their use in various industrial applications such as the biofuel, pharmaceutical, fine chemicals and food industries has increased. The understanding of the specific factors that confer the ability to withstand extreme habitats on such enzymes has become a priority for their biotechnological use. The most studied marine extremophiles are prokaryotes and in this review, we present the most studied archaea and bacteria extremophiles and their hydrolases, and discuss their use for industrial applications.

Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations

Improving the thermostability of industrial enzymes is an important protein engineering challenge. Point mutations, induced to increase thermostability, affect the structure and dynamics of the target protein in several ways and thus can also affect its activity. There appears to be no general rules for improving the thermostabilty of enzymes without adversely affecting their enzymatic activity. We report MD simulations, of wild type Bacillus subtilis lipase (WT) and its six progressively thermostable mutants (2M, 3M, 4M, 6M, 9M, and 12M), performed at different temperatures, to address this issue. Less thermostable mutants (LTMs), 2M to 6M, show WT-like dynamics at all simulation temperatures. However, the two more thermostable mutants (MTMs) show the required flexibility at appropriate temperature ranges and maintain conformational stability at high temperature. They show a deep and rugged free-energy landscape, confining them within a near-native conformational space by conserving noncovalent interactions, and thus protecting them from possible aggregation. In contrast, the LTMs having marginally higher thermostabilities than WT show greater probabilities of accessing non-native conformations, which, due to aggregation, have reduced possibilities of reverting to their respective native states under refolding conditions. Our analysis indicates the possibility of nonadditive effects of point mutations on the conformational stability of LTMs.

Interface matters: the stiffness route to stability of a thermophilic tetrameric malate dehydrogenase

In this work we investigate by computational means the behavior of two orthologous bacterial proteins, a mesophilic and a thermophilic tetrameric malate dehydrogenase (MalDH), at different temperatures. Namely, we quantify how protein mechanical rigidity at different length- and time-scales correlates to protein thermophilicity as commonly believed. In particular by using a clustering analysis strategy to explore the conformational space of the folded proteins, we show that at ambient conditions and at the molecular length-scale the thermophilic variant is indeed more rigid that the mesophilic one. This rigidification is the result of more efficient inter-domain interactions, the strength of which is further quantified via ad hoc free energy calculations. When considered isolated, the thermophilic domain is indeed more flexible than the respective mesophilic one. Upon oligomerization, the induced stiffening of the thermophilic protein propagates from the interface to the active site where the loop, controlling the access to the catalytic pocket, anchors down via an extended network of ion-pairs. On the contrary in the mesophilic tetramer the loop is highly mobile. Simulations at high temperature, could not re-activate the mobility of the loop in the thermophile. This finding opens questions on the similarities of the binding processes for these two homologues at their optimal working temperature and suggests for the thermophilic variant a possible cooperative role of cofactor/substrate.

Correlation between supercoiling and conformational motions of the bacterial flagellar filament

The bacterial flagellar filament is a very large macromolecular assembly of a single protein, flagellin. Various supercoiled states of the filament exist, which are formed by two structurally different conformations of flagellin in different ratios. We investigated the correlation between supercoiling of the protofilaments and molecular dynamics in the flagellar filament using quasielastic and elastic incoherent neutron scattering on the picosecond and nanosecond timescales. Thermal fluctuations in the straight L- and R-type filaments were measured and compared to the resting state of the wild-type filament. Amplitudes of motion on the picosecond timescale were found to be similar in the different conformational states. Mean-square displacements and protein resilience on the 0.1 ns timescale demonstrate that the L-type state is more flexible and less resilient than the R-type, whereas the wild-type state lies in between. Our results provide strong support that supercoiling of the protofilaments in the flagellar filament is determined by the strength of molecular forces in and between the flagellin subunits.

How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain

Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.

Neutron scattering techniques and applications in structural biology

Neutron scattering is exquisitely sensitive to the position, concentration, and dynamics of hydrogen atoms in materials and is a powerful tool for the characterization of structure-function and interfacial relationships in biological systems. Modern neutron scattering facilities offer access to a sophisticated, nondestructive suite of instruments for biophysical characterization that provides spatial and dynamic information spanning from Ångstroms to microns and from picoseconds to microseconds, respectively. Applications in structural biology range from the atomic-resolution analysis of individual hydrogen atoms in enzymes through to meso- and macro-scale analysis of complex biological structures, membranes, and assemblies. The large difference in neutron scattering length between hydrogen and deuterium allows contrast variation experiments to be performed and enables H/D isotopic labeling to be used for selective and systematic analysis of the local structure, dynamics, and interactions of multi-component systems. This overview describes the available techniques and summarizes their practical application to the study of biomolecular systems.

Perspectives in biological physics: the nDDB project for a neutron Dynamics Data Bank for biological macromolecules

Neutron spectroscopy provides experimental data on time-dependent trajectories, which can be directly compared to molecular dynamics simulations. Its importance in helping us to understand biological macromolecules at a molecular level is demonstrated by the results of a literature survey over the last two to three decades. Around 300 articles in refereed journals relate to neutron scattering studies of biological macromolecular dynamics, and the results of the survey are presented here. The scope of the publications ranges from the general physics of protein and solvent dynamics, to the biologically relevant dynamics-function relationships in live cells. As a result of the survey we are currently setting up a neutron Dynamics Data Bank (nDDB) with the aim to make the neutron data on biological systems widely available. This will benefit, in particular, the MD simulation community to validate and improve their force fields. The aim of the database is to expose and give easy access to a body of experimental data to the scientific community. The database will be populated with as much of the existing data as possible. In the future it will give value, as part of a bigger whole, to high throughput data, as well as more detailed studies. A range and volume of experimental data will be of interest in determining how quantitatively MD simulations can reproduce trends across a range of systems and to what extent such trends may depend on sample preparation and data reduction and analysis methods. In this context, we strongly encourage researchers in the field to deposit their data in the nDDB.

Enhancing the promiscuous phosphotriesterase activity of a thermostable lactonase (GkaP) for the efficient degradation of organophosphate pesticides

The phosphotriesterase-like lactonase (PLL) enzymes in the amidohydrolase superfamily hydrolyze various lactones and exhibit latent phosphotriesterase activities. These enzymes serve as attractive templates for in vitro evolution of neurotoxic organophosphates (OPs) with hydrolytic capabilities that can be used as bioremediation tools. Here, a thermostable PLL from Geobacillus kaustophilus HTA426 (GkaP) was targeted for joint laboratory evolution with the aim of enhancing its catalytic efficiency against OP pesticides. By a combination of site saturation mutagenesis and whole-gene error-prone PCR approaches, several improved variants were isolated. The most active variant, 26A8C, accumulated eight amino acid substitutions and demonstrated a 232-fold improvement over the wild-type enzyme in reactivity (k(cat)/K(m)) for the OP pesticide ethyl-paraoxon. Concomitantly, this variant showed a 767-fold decrease in lactonase activity with δ-decanolactone, imparting a specificity switch of 1.8 × 10(5)-fold. 26A8C also exhibited high hydrolytic activities (19- to 497-fold) for several OP pesticides, including parathion, diazinon, and chlorpyrifos. Analysis of the mutagenesis sites on the GkaP structure revealed that most mutations are located in loop 8, which determines substrate specificity in the amidohydrolase superfamily. Molecular dynamics simulation shed light on why 26A8C lost its native lactonase activity and improved the promiscuous phosphotriesterase activity. These results permit us to obtain further insights into the divergent evolution of promiscuous enzymes and suggest that laboratory evolution of GkaP may lead to potential biological solutions for the efficient decontamination of neurotoxic OP compounds.

Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus porosus) were investigated using circular dichroism, incoherent elastic neutron scattering and coarse-grained Brownian dynamics simulations. The experimental results from Hb solutions revealed a direct correlation between protein resilience, melting temperature and average body temperature of the different species on the 0.1 ns time scale. Molecular forces appeared to be adapted to permit conformational fluctuations with a root mean square displacement close to 1.2 Å at the corresponding average body temperature of the endotherms. Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit over the entire temperature range in which the animal lives. In fully hydrated powder samples of human and chicken, Hb mean square displacements and effective force constants on the 1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, in contrast to the solution case. A complementary result of the study, therefore, is that one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the pico- to nanosecond time scale which might be relevant for biological function. Coarse-grained Brownian dynamics simulations permitted to explore residue-specific effects. They indicated that temperature sensing of human and chicken Hb occurs mainly at residues lining internal cavities in the β-subunits.

Dynamics-stability relationships in apo- and holomyoglobin: a combined neutron scattering and molecular dynamics simulations study

The removal of the heme group from myoglobin (Mb) results in a destabilization of the protein structure. The dynamic basis of the destabilization was followed by comparative measurements on holo- (holo-Mb) and apomyoglobin (apo-Mb). Mean-squared displacements (MSD) and protein resilience on the picosecond-to-nanosecond timescale were measured by elastic incoherent neutron scattering. Differences in thermodynamic parameters, MSD, and resilience were observed for both proteins. The resilience of holo-Mb was significantly lower than that of apo-Mb, indicating entropic stabilization by a higher degree of conformational sampling in the heme-bound folded protein. Molecular dynamics simulations provided site-specific information. Averaged over the whole structure, the molecular dynamics simulations yielded similar MSD and resilience values for the two proteins. The mobility of residues around the heme group in holo-Mb showed a smaller MSD and higher resilience compared to the same residue group in apo-Mb. It is of interest that in holo-Mb, higher MSD values are observed for the residues outside the heme pocket, indicating an entropic contribution to protein stabilization by heme binding, which is in agreement with experimental results.

Temperature adaptations in psychrophilic, mesophilic and thermophilic chloride-dependent alpha-amylases

The functional and structural adaptations to temperature have been addressed in homologous chloride-dependent α-amylases from a psychrophilic Antarctic bacterium, the ectothermic fruit fly, the homeothermic pig and from a thermophilic actinomycete. This series covers nearly all temperatures encountered by living organisms. We report a striking continuum in the functional properties of these enzymes coupled to their structural stability and related to the thermal regime of the source organism. In particular, thermal stability recorded by intrinsic fluorescence, circular dichroism and differential scanning calorimetry appears to be a compromise between the requirement for a stable native state and the proper structural dynamics to sustain the function at the environmental/physiological temperatures. The thermodependence of activity, the kinetic parameters, the activations parameters and fluorescence quenching support these activity-stability relationships in the investigated α-amylases.

Dynamics of thermodynamically stable, kinetically trapped, and inhibitor-bound states of pepsin

The pepsin folding mechanism involves a prosegment (PS) domain that catalyzes folding, which is then removed, resulting in a kinetically trapped native state. Although native pepsin (Np) is kinetically stable, it is irreversibly denatured due to a large folding barrier, and in the absence of the PS it folds to a more thermodynamically stable denatured state, termed refolded pepsin (Rp). This system serves as a model to understand the nature of kinetic barriers and folding transitions between compact states. Quasielastic neutron scattering (QENS) was used to characterize and compare the flexibility of Np, as a kinetically trapped state, with that of Rp, as a thermodynamically stable fold. Additionally, the dynamics of Np were compared with those of a partially unfolded form and a thermally stabilized, inhibitor-bound form. QENS revealed length-scale-dependent differences between Np and Rp on a picosecond timescale and indicated greater flexibility in Np, leading to the conclusion that kinetic stabilization likely does not correspond to reduced internal dynamics. Furthermore, large differences were observed upon inhibition, indicating that QENS of proteins in solution may prove useful for examining the role of conformational entropy changes in ligand binding.

Crowding induces differences in the diffusion of thermophilic and mesophilic proteins: a new look at neutron scattering results

The dynamical basis underlying the increased thermal stability of thermophilic proteins remains uncertain. Here, we challenge the new paradigm established by neutron scattering experiments in solution, in which the adaptation of thermophilic proteins to high temperatures lies in the lower sensitivity of their flexibility to temperature changes. By means of a combination of molecular dynamics and Brownian dynamics simulations, we report a reinterpretation of those experiments and show evidence that under crowding conditions, such as in vivo, thermophilic and homolog mesophilic proteins have diffusional properties with different thermal behavior.

Thermophilic proteins: insight and perspective from in silico experiments

Proteins from thermophilic and hyperthermophilic organisms are stable and function at high temperatures (50-100 °C). The importance of understanding the microscopic mechanisms underlying this thermal resistance is twofold: it is key for acquiring general clues on how proteins maintain their fold stable and for targeting those medical and industrial applications that aim at designing enzymes that can work under harsh conditions. In this tutorial review we first provide the general background of protein thermostability by specifically focusing on the structural and thermodynamic peculiarities; next, we discuss how computational studies based on Molecular Dynamics simulations can broaden and refine our knowledge on such special class of proteins.

Improving the promiscuous nerve agent hydrolase activity of a thermostable archaeal lactonase

The thermostable Phosphotriesterase-Like Lactonase from Sulfolobus solfataricus (SsoPox) hydrolyzes lactones and, at a lower rate, neurotoxic organophosphorus compounds. The persistent demand of detoxification tools in the field of agricultural wastes and restoring of conditions after terrorist acts prompted us to exploit SsoPox as a "starter" to evolve its ancillary nerve agents hydrolytic capability. A directed evolution strategy yielded, among several variants, the single mutant W263F with k(cat) and specificity constant against paraoxon 16- and 6-fold enhanced, respectively, compared to the wild type. Furthermore, a phenomenon of enzyme activation by SDS has been observed, which allowed to increase those values 150- and 28-fold, respectively. The activity of SsoPox against the deadly nerve gas Cyclosarin has been reported for the first time and proved to be substantially unaffected for variant W263F. Finally, outperforming efficiency of W263F was demonstrated, under severe stressing conditions, with respect to the best known phosphotriesterase PTE from Brevundimonas diminuta.

Thermus thermophilus proteins that are differentially expressed in response to growth temperature and their implication in thermoadaptation

As a kind of important extremophiles to realize the adaptation of life at high temperatures, thermophiles have attracted extensive studies. However, the pathways of thermophile proteins related to thermoadaptation remain to be addressed. Our study showed that there existed two types of protein profiles for the thermophile Thermus thermophilus wl in response to temperature change. One of them came from cultures growing below 65 degrees C, which was close to the optimal growth temperature, and another from cultures at or above 65 degrees C. These protein profiles were confirmed by Northern blots. On the basis of the proteomic and computational analyses, it was found that the thermophile proteins related to thermoadaptation might be involved in metabolic pathways as well as the stabilities and modifications of DNA and proteins. Interestingly, for the basic metabolism glycolysis, the phosphoglucomutase was up-regulated at below-optimum temperature, while the glyceraldehyde-3-phosphate dehydrogenase was up-regulated at above-optimum temperature, suggesting that different regulations might be used for basic metabolism at different temperatures. To characterize the proteins in response to high temperatures, superoxide dismutase (SOD), an important enzyme in organism to remove free radical produced in stress environment such as high temperature, was selected as a target protein for this investigation. SOD was inactivated to construct a SOD mutant. The results showed that the SOD protein was essential in thermoadaptation of T. thermophilus. Our study, therefore, presented the thermophile proteins required for thermoadaptation and their possible pathways in thermoadaptation.

Neutron spin-echo studies of hemoglobin and myoglobin: multiscale internal dynamics

Neutron spin-echo spectroscopy was used to study structural fluctuations that occur in hemoglobin (Hb) and myoglobin (Mb) in solution. Using neutron spin-echo data up to a very high momentum transfer q ( approximately 0.62 A(-)(1)), we characterized the internal dynamics of these proteins at the levels of dynamic pair correlation function and self-correlation function in the time range of several picoseconds to a few nanoseconds. In the same protein solution, data transition from pair correlation motion to self-correlation motion as the momentum transfer q increases. At low q, coherent scattering dominates; at high q, observations are largely due to incoherent scattering. The low q data were interpreted in terms of an effective diffusion coefficient; on the other hand, the high q data were interpreted in terms of mean square displacements. Comparison of data from the two homologous proteins collected at different temperatures and protein concentrations was used to assess the contributions made by translational and rotational diffusion and internal modes of motion to the data. The temperature dependence of decay times can be attributed to changes in the viscosity and temperature of the solvent, as predicted by the Stokes-Einstein relationship. This is true for contributions from both diffusive and internal modes of motion, indicating an intimate relationship between the internal dynamics of the proteins and the viscosity of the solvent. Viscosity change associated with protein concentration can account for changes in diffusion observed at different concentrations, but is apparently not the only factor involved in the changes in internal dynamics observed with change in protein concentration. Data collected at high q indicate that internal modes in Mb are generally faster than those in Hb, perhaps due to the greater surface-to-volume ratio of Mb and the fact that surface groups tend to exhibit faster motion than buried groups. Comparison of data from Hb and data from Mb at low q indicates an unexpectedly rapid motion of Hb alphabeta dimers relative to one another. Dynamic motion of subunits is increasingly perceived as important to the allosteric behavior of Hb. Our data demonstrate that this motion is highly sensitive to protein concentration, temperature, and solvent viscosity, indicating that great care needs to be exercised in interpreting its effect on protein function.

An examination of the relationship between active site loop size and thermodynamic activation parameters for orotidine 5'-monophosphate decarboxylase from mesophilic and thermophilic organisms

Closure of the active site phosphate gripper loop of orotidine 5'-monophosphate decarboxylase from Saccharomyces cerevisiae (ScOMPDC) over the bound substrate orotidine 5'-monophosphate (OMP) activates the bound substrate for decarboxylation by at least 10(4)-fold [Amyes, T. L., Richard, J. P., and Tait, J. J. (2005) J. Am. Chem. Soc. 127, 15708-15709]. The 19-residue phosphate gripper loop of the mesophilic ScOMPDC is much larger than the nine-residue loop at the ortholog from the thermophile Methanothermobacter thermautotrophicus (MtOMPDC). This difference in loop size results in a small decrease in the total intrinsic phosphate binding energy of the phosphodianion group of OMP from 11.9 to 11.6 kcal/mol, along with a modest decrease in the extent of activation by phosphite dianion of decarboxylation of the truncated substrate 1-(beta-D-erythrofuranosyl)orotic acid. The activation parameters DeltaH(double dagger) and DeltaS(double dagger) for k(cat) for decarboxylation of OMP are 3.6 kcal/mol and 10 cal K(-1) mol(-1) more positive, respectively, for MtOMPDC than for ScOMPDC. We suggest that these differences are related to the difference in the size of the active site loops at the mesophilic ScOMPDC and the thermophilic MtOMPDC. The greater enthalpic transition state stabilization available from the more extensive loop-substrate interactions for the ScOMPDC-catalyzed reaction is largely balanced by a larger entropic requirement for immobilization of the larger loop at this enzyme.

Protein flexibility and metal coordination changes in DHAP-dependent aldolases

The mobility of rhamnulose-1-phosphate aldolase (RhuA) was analysed with a normal mode description and high level calculations on models of the active site. We report the connection between the mobility and the chemical properties of the active site, and compare them to a closely related enzyme, fuculose-1-phosphate aldolase (FucA). Calculations show that the different coordination number for the zinc ion, reported in the crystal structures of RhuA and FucA, was due to a different spatial arrangement of the residues, not to their different chemical nature. Moreover, the metal coordination change is correlated with activity. The domain mobility of the enzyme can reshape the active site of RhuA into the arrangement found in the FucA structure, and vice-versa. This has a direct influence on the energy barrier for the aldol reaction catalyzed by these enzymes, thus showing a coupling of the domain movements and the catalytic effects. Hence domain movements and the coordination chemistry of the active site metal suggest an explanation of why these enzymes have similar experimental turnover rates.

New sources and instrumentation for neutrons in biology

Neutron radiation offers significant advantages for the study of biological molecular structure and dynamics. A broad and significant effort towards instrumental and methodological development to facilitate biology experiments at neutron sources worldwide is reviewed.