Structures of mesophilic and extremophilic citrate synthases reveal rigidity and flexibility for function

Comparative Modeling and Analysis of Extremophilic D-Ala-D-Ala Carboxypeptidases

Understanding the molecular adaptations of organisms to extreme environments requires a comparative analysis of protein structure, function, and dynamics across species found in different environmental conditions. Computational studies can be particularly useful in this pursuit, allowing exploratory studies of large numbers of proteins under different thermal and chemical conditions that would be infeasible to carry out experimentally. Here, we perform such a study of the MEROPS family S11, S12, and S13 proteases from psychophilic, mesophilic, and thermophilic bacteria. Using a combination of protein structure prediction, atomistic molecular dynamics, and trajectory analysis, we examine both conserved features and trends across thermal groups. Our findings suggest a number of hypotheses for experimental investigation.

Insights on Rigidity and Flexibility at the Global and Local Levels of Protein Structures and Their Roles in Homologous Psychrophilic, Mesophilic, and Thermophilic Proteins: A Computational Study

The rigidity and flexibility of homologous psychrophilic (P), mesophilic (M), and thermophilic (T) proteins have been investigated at the global and local levels in terms of "packing factors" and "atomic fluctuations" obtained from B-factors. For comparison of atomic fluctuations, correction of errors by considering errors in B-factors from all sources in a consolidated manner and conversion of the fluctuations to the same temperature have been suggested and validated. The results indicate no differences in the global values like the average packing factor among the three classes of protein homologues, but at local levels there are differences. A comparison of homologous protein triplets show that the average atomic fluctuations at a given temperature mainly obey the order P > M > T. Packing factors and the atomic fluctuations are anti-correlated, suggesting that altering the rigidity of the active site might be a potential strategy to make tailor-made psychrophilic or thermophilic proteins from their mesophilic homologues. The computer codes developed and used in this work are available at the link https://github.com/Munna-Sarkar/proteins-rigidity-flexibility.git.

Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism

Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/β hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7's closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7's cold adaptation rather than changes in dynamics.

Monoclonal antibody stability can be usefully monitored using the excitation-energy-dependent fluorescence edge-shift

Among the major challenges in the development of biopharmaceuticals are structural heterogeneity and aggregation. The development of a successful therapeutic monoclonal antibody (mAb) requires both a highly active and also stable molecule. Whilst a range of experimental (biophysical) approaches exist to track changes in stability of proteins, routine prediction of stability remains challenging. The fluorescence red edge excitation shift (REES) phenomenon is sensitive to a range of changes in protein structure. Based on recent work, we have found that quantifying the REES effect is extremely sensitive to changes in protein conformational state and dynamics. Given the extreme sensitivity, potentially this tool could provide a ‘fingerprint’ of the structure and stability of a protein. Such a tool would be useful in the discovery and development of biopharamceuticals and so we have explored our hypothesis with a panel of therapeutic mAbs. We demonstrate that the quantified REES data show remarkable sensitivity, being able to discern between structurally identical antibodies and showing sensitivity to unfolding and aggregation. The approach works across a broad concentration range (µg–mg/ml) and is highly consistent. We show that the approach can be applied alongside traditional characterisation testing within the context of a forced degradation study (FDS). Most importantly, we demonstrate the approach is able to predict the stability of mAbs both in the short (hours), medium (days) and long-term (months). The quantified REES data will find immediate use in the biopharmaceutical industry in quality assurance, formulation and development. The approach benefits from low technical complexity, is rapid and uses instrumentation which exists in most biochemistry laboratories without modification.

Using geometric simulation software ‘GASP’ to model conformational flexibility in a family of zinc metal–organic frameworks

Synthesis with different solvents leads the same metals and ligands to assemble into geometrically and topologically distinct frameworks, aided by intrinsic flexibility of the ligand.

Salt bridge impact on global rigidity and thermostability in thermophilic citrate synthase

It has been suggested that structural rigidity is connected to thermostability, e.g. in enzymes from thermophilic microorganisms. We examine the importance of correctly handling salt bridges, and interactions which we term 'strong polars', when constructing the constraint network for global rigidity analysis in these systems. Through a comparison of rigidity in citrate synthases, we clarify the relationship between rigidity and thermostability. In particular, with our corrected handling of strong polar interactions, the difference in rigidity between mesophilic and thermophilic structures is detected more clearly than in previous studies. The increase in rigidity did not detract from the functional flexibility of the active site in all systems once their respective temperature range had been reached. We then examine the distribution of salt bridges in thermophiles that were previously unaccounted for in flexibility studies. We show that in hyperthermophiles these have stabilising roles in the active site; occuring in close proximity to key residues involved in catalysis and binding of the protein.

Structure and function of L-threonine-3-dehydrogenase from the parasitic protozoan Trypanosoma brucei revealed by X-ray crystallography and geometric simulations

Two of the world's most neglected tropical diseases, human African trypanosomiasis (HAT) and Chagas disease, are caused by protozoan parasites of the genus Trypanosoma. These organisms possess specialized metabolic pathways, frequently distinct from those in humans, which have potential to be exploited as novel drug targets. This study elucidates the structure and function of L-threonine-3-dehydrogenase (TDH) from T. brucei, the causative pathogen of HAT. TDH is a key enzyme in the metabolism of L-threonine, and an inhibitor of TDH has been shown to have trypanocidal activity in the procyclic form of T. brucei. TDH is a nonfunctional pseudogene in humans, suggesting that it may be possible to rationally design safe and specific therapies for trypanosomiasis by targeting this parasite enzyme. As an initial step, the TDH gene from T. brucei was expressed and the three-dimensional structure of the enzyme was solved by X-ray crystallography. In multiple crystallographic structures, T. brucei TDH is revealed to be a dimeric short-chain dehydrogenase that displays a considerable degree of conformational variation in its ligand-binding regions. Geometric simulations of the structure have provided insight into the dynamic behaviour of this enzyme. Furthermore, structures of TDH bound to its natural substrates and known inhibitors have been determined, giving an indication of the mechanism of catalysis of the enzyme. Collectively, these results provide vital details for future drug design to target TDH or related enzymes.

In Silico and Structural Analyses Demonstrate That Intrinsic Protein Motions Guide T Cell Receptor Complementarity Determining Region Loop Flexibility

T-cell immunity is controlled by T cell receptor (TCR) binding to peptide major histocompatibility complexes (pMHCs). The nature of the interaction between these two proteins has been the subject of many investigations because of its central role in immunity against pathogens, cancer, in autoimmunity, and during organ transplant rejection. Crystal structures comparing unbound and pMHC-bound TCRs have revealed flexibility at the interaction interface, particularly from the perspective of the TCR. However, crystal structures represent only a snapshot of protein conformation that could be influenced through biologically irrelevant crystal lattice contacts and other factors. Here, we solved the structures of three unbound TCRs from multiple crystals. Superposition of identical TCR structures from different crystals revealed some conformation differences of up to 5 Å in individual complementarity determining region (CDR) loops that are similar to those that have previously been attributed to antigen engagement. We then used a combination of rigidity analysis and simulations of protein motion to reveal the theoretical potential of TCR CDR loop flexibility in unbound state. These simulations of protein motion support the notion that crystal structures may only offer an artifactual indication of TCR flexibility, influenced by crystallization conditions and crystal packing that is inconsistent with the theoretical potential of intrinsic TCR motions.

A complete thermodynamic analysis of enzyme turnover links the free energy landscape to enzyme catalysis

Our understanding of how enzymes work is coloured by static structure depictions where the enzyme scaffold is presented as either immobile, or in equilibrium between well-defined static conformations. Proteins, however, exhibit a large degree of motion over a broad range of timescales and magnitudes and this is defined thermodynamically by the enzyme free energy landscape (FEL). The role and importance of enzyme motion is extremely contentious. Much of the challenge is in the experimental detection of so called 'conformational sampling' involved in enzyme turnover. Herein we apply combined pressure and temperature kinetics studies to elucidate the full suite of thermodynamic parameters defining an enzyme FEL as it relates to enzyme turnover. We find that the key thermodynamic parameters governing vibrational modes related to enzyme turnover are the isobaric expansivity term and the change in heat capacity for enzyme catalysis. Variation in the enzyme FEL affects these terms. Our analysis is supported by a range of biophysical and computational approaches that specifically capture information on protein vibrational modes and the FEL (all atom flexibility calculations, red edge excitation shift spectroscopy and viscosity studies) that provide independent evidence for our findings. Our data suggest that restricting the enzyme FEL may be a powerful strategy when attempting to rationally engineer enzymes, particularly to alter thermal activity. Moreover, we demonstrate how rational predictions can be made with a rapid computational approach.

Evidence of Correlated Static Disorder in the Fenna-Matthews-Olson Complex

Observation of excitonic quantum beats in photosynthetic antennae has prompted wide debate regarding the function of excitonic coherence in pigment-protein complexes. Much of this work focuses on the interactions of excitons with the femto-to-picosecond dynamical fluctuations of their environment. However, in experiments these effects can be masked by static disorder of the excited-state energies across ensembles, whose microscopic origins are challenging to predict. Here the excited-state properties of ∼2000 atom clusters of the Fenna-Matthews-Olson complex are simulated using a unique combination of linear-scaling density functional theory and constrained geometric dynamics. While slow, large amplitude protein motion leads to large variations in the Q_y transitions of two pigments, we identify pigment-protein correlations that greatly reduce variations in the energy gap across the ensemble, which is consistent with experimental observations of suppressed inhomogeneous dephasing of quantum beats.

Structure of the family B DNA polymerase from the hyperthermophilic archaeon Pyrobaculum calidifontis

The family B DNA polymerase from Pyrobaculum calidifontis (Pc-polymerase) consists of 783 amino acids and is magnesium-ion dependent. It has an optimal pH of 8.5, an optimal temperature of 75°C and a half-life of 4.5 h at 95°C, giving it greater thermostability than the widely used Taq DNA polymerase. The enzyme is also capable of PCR-amplifying larger DNA fragments of up to 7.5 kb in length. It was shown to have functional, error-correcting 3'-5' exonuclease activity, as do the related high-fidelity DNA polymerases from Pyrococcus furiosus, Thermococcus kodakarensis KOD1 and Thermococcus gorgonarius, which have extensive commercial applications. Pc-polymerase has a quite low sequence identity of approximately 37% to these enzymes, which, in contrast, have very high sequence identity to each other, suggesting that the P. calidifontis enzyme is distinct. Here, the structure determination of Pc-polymerase is reported, which has been refined to an R factor of 24.47% and an R_free of 28.81% at 2.80 Å resolution. The domains of the enzyme are arranged in a circular fashion to form a disc with a narrow central channel. One face of the disc has a number of connected crevices in it, which allow the protein to bind duplex and single-stranded DNA. The central channel is thought to allow incoming nucleoside triphosphates to access the active site. The enzyme has a number of unique structural features which distinguish it from other archaeal DNA polymerases and may account for its high processivity. A model of the complex with the primer-template duplex of DNA indicates that the largest conformational change that occurs upon DNA binding is the movement of the thumb domain, which rotates by 7.6° and moves by 10.0 Å. The surface potential of the enzyme is dominated by acidic groups in the central region of the molecule, where catalytic magnesium ions bind at the polymerase and exonuclease active sites. The outer regions are richer in basic amino acids that presumably interact with the sugar-phosphate backbone of DNA. The large number of salt bridges may contribute to the high thermal stability of this enzyme.

Structural studies of substrate and product complexes of 5-aminolaevulinic acid dehydratase from humans,Escherichia coliand the hyperthermophilePyrobaculum calidifontis

A number of X-ray analyses of an enzyme involved in a key early stage of tetrapyrrole biosynthesis are reported. Two structures of human 5-aminolaevulinate dehydratase (ALAD), native and recombinant, have been determined at 2.8 Å resolution, showing that the enzyme adopts an octameric quaternary structure in accord with previously published analyses of the enzyme from a range of other species. However, this is in contrast to the finding that a disease-related F12L mutant of the human enzyme uniquely forms hexamers [Breiniget al.(2003),Nature Struct. Biol.10, 757–763]. Monomers of all ALADs adopt the TIM-barrel fold; the subunit conformation that assembles into the octamer includes the N-terminal tail of one monomer curled around the (α/β)8barrel of a neighbouring monomer. Both crystal forms of the human enzyme possess two monomers per asymmetric unit, termedAandB. In the native enzyme there are a number of distinct structural differences between theAandBmonomers, with the latter exhibiting greater disorder in a number of loop regions and in the active site. In contrast, the second monomer of the recombinant enzyme appears to be better defined and the active site of both monomers clearly possesses a zinc ion which is bound by three conserved cysteine residues. In native human ALAD, theAmonomer also has a ligand resembling the substrate ALA which is covalently bound by a Schiff base to one of the active-site lysines (Lys252) and is held in place by an ordered active-site loop. In contrast, these features of the active-site structure are disordered or absent in theBsubunit of the native human enzyme. The octameric structure of the zinc-dependent ALAD from the hyperthermophilePyrobaculum calidifontisis also reported at a somewhat lower resolution of 3.5 Å. Finally, the details are presented of a high-resolution structure of theEscherichia coliALAD enzyme co-crystallized with a noncovalently bound moiety of the product, porphobilinogen (PBG). This structure reveals that the pyrrole side-chain amino group is datively bound to the active-site zinc ion and that the PBG carboxylates interact with the enzymeviahydrogen bonds and salt bridges with invariant residues. A number of hydrogen-bond interactions that were previously observed in the structure of yeast ALAD with a cyclic intermediate resembling the product PBG appear to be weaker in the new structure, suggesting that these interactions are only optimal in the transition state.

The flexibility and dynamics of protein disulfide isomerase

We have studied the mobility of the multidomain folding catalyst, protein disulfide isomerase (PDI), by a coarse‐graining approach based on flexibility. We analyze our simulations of yeast PDI (yPDI) using measures of backbone movement, relative positions and orientations of domains, and distances between functional sites. We find that there is interdomain flexibility at every interdomain junction but these show very different characteristics. The extent of interdomain flexibility is such that yPDI's two active sites can approach much more closely than is found in crystal structures—and indeed hinge motion to bring these sites into proximity is the lowest energy normal mode of motion of the protein. The flexibility predicted for yPDI (based on one structure) includes the other known conformation of yPDI and is consistent with (i) the mobility observed experimentally for mammalian PDI and (ii) molecular dynamics. We also observe intradomain flexibility and clear differences between the domains in their propensity for internal motion. Our results suggest that PDI flexibility enables it to interact with many different partner molecules of widely different sizes and shapes, and highlights considerable similarities of yPDI and mammalian PDI. Proteins 2016; 84:1776–1785. © 2016 Wiley Periodicals, Inc.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal β-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

Unraveling the B. pseudomallei Heptokinase WcbL: From Structure to Drug Discovery

Gram-negative bacteria utilize heptoses as part of their repertoire of extracellular polysaccharide virulence determinants. Disruption of heptose biosynthesis offers an attractive target for novel antimicrobials. A critical step in the synthesis of heptoses is their 1-O phosphorylation, mediated by kinases such as HldE or WcbL. Here, we present the structure of WcbL from Burkholderia pseudomallei. We report that WcbL operates through a sequential ordered Bi-Bi mechanism, loading the heptose first and then ATP. We show that dimeric WcbL binds ATP anti-cooperatively in the absence of heptose, and cooperatively in its presence. Modeling of WcbL suggests that heptose binding causes an elegant switch in the hydrogen-bonding network, facilitating the binding of a second ATP molecule. Finally, we screened a library of drug-like fragments, identifying hits that potently inhibit WcbL. Our results provide a novel mechanism for control of substrate binding and emphasize WcbL as an attractive anti-microbial target for Gram-negative bacteria.

Exploration of Structural and Functional Variations Owing to Point Mutations in α-NAGA

Schindler disease is a lysosomal storage disorder caused due to deficiency or defective activity of alpha-N-acetylgalactosaminidase (α-NAGA). Mutations in gene encoding α-NAGA cause wide range of diseases, characterized with mild to severe clinical features. Molecular effects of these mutations are yet to be explored in detail. Therefore, this study was focused on four missense mutations of α-NAGA namely, S160C, E325K, R329Q and R329W. Native and mutant structures of α-NAGA were analysed to determine geometrical deviations such as the contours of root mean square deviation, root mean square fluctuation, percentage of residues in allowed regions of Ramachandran plot and solvent accessible surface area, using conformational sampling technique. Additionally, global energy-minimized structures of native and mutants were further analysed to compute their intra-molecular interactions, hydrogen bond dilution and distribution of secondary structure. In addition, docking studies were also performed to determine variations in binding energies between native and mutants. The deleterious effects of mutants were evident due to variations in their active site residues pertaining to spatial conformation and flexibility, comparatively. Hence, variations exhibited by mutants, namely S160C, E325K, R329Q and R329W to that of native, consequently, lead to the detrimental effects causing Schindler disease. This study computationally explains the underlying reasons for the pathogenesis of the disease, thereby aiding future researchers in drug development and disease management.

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.

Temperature dependent mistranslation in a hyperthermophile adapts proteins to lower temperatures

All organisms universally encode, synthesize and utilize proteins that function optimally within a subset of growth conditions. While healthy cells are thought to maintain high translational fidelity within their natural habitats, natural environments can easily fluctuate outside the optimal functional range of genetically encoded proteins. The hyperthermophilic archaeon Aeropyrum pernix (A. pernix) can grow throughout temperature variations ranging from 70 to 100°C, although the specific factors facilitating such adaptability are unknown. Here, we show that A. pernix undergoes constitutive leucine to methionine mistranslation at low growth temperatures. Low-temperature mistranslation is facilitated by the misacylation of tRNA^Leu with methionine by the methionyl-tRNA synthetase (MetRS). At low growth temperatures, the A. pernix MetRS undergoes a temperature dependent shift in tRNA charging fidelity, allowing the enzyme to conditionally charge tRNA^Leu with methionine. We demonstrate enhanced low-temperature activity for A. pernix citrate synthase that is synthesized during leucine to methionine mistranslation at low-temperature growth compared to its high-fidelity counterpart synthesized at high-temperature. Our results show that conditional leucine to methionine mistranslation can make protein adjustments capable of improving the low-temperature activity of hyperthermophilic proteins, likely by facilitating the increasing flexibility required for greater protein function at lower physiological temperatures.

Structural comparison between the open and closed forms of citrate synthase from Thermus thermophilus HB8

The crystal structures of citrate synthase from the thermophilic eubacteria Thermus thermophilus HB8 (TtCS) were determined for an open form at 1.5 Å resolution and for closed form at 2.3 Å resolution, respectively. In the absence of ligands TtCS in the open form was crystalized into a tetragonal form with a single subunit in the asymmetric unit. TtCS was also co-crystallized with citrate and coenzyme-A to form an orthorhombic crystal with two homodimers in the asymmetric unit. Citrate and CoA are found in the active site situated between the large domain and the small domain in all subunit whereas the complex shows two distinct closed conformations, the fully closed form and partially closed form. Structural comparisons are performed to describe conformational changes associated with binding of products of TtCS. Upon binding of citrate, basic residues in the active site move toward citrate and make a hydrogen bond network in the active site, inducing a large-scale rotation of the small domain relative to the large domain. CoA is sandwiched between the small and large domains and then the cysteamine tail is inserted into the active site with a cooperative rotation around mainchain dihedrals in the hinge region connecting helices M and N. According to this rotation these helices are extended to close the active site completely. The considerable flexibility and structural rearrangements in the hinge region are crucial for an ordered bibi reaction in catalysis for microbial CSs.

Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions

Large-scale conformational change is a common feature in the catalytic cycles of enzymes. Many enzymes function as homodimers with active sites that contain elements from both chains. Symmetric and anti-symmetric cooperative motions in homodimers can potentially lead to correlated active site opening and/or closure, likely to be important for ligand binding and release. Here, we examine such motions in two different domain-swapped homodimeric enzymes: the DcpS scavenger decapping enzyme and citrate synthase. We use and compare two types of all-atom simulations: conventional molecular dynamics simulations to identify physically meaningful conformational ensembles, and rapid geometric simulations of flexible motion, biased along normal mode directions, to identify relevant motions encoded in the protein structure. The results indicate that the opening/closure motions are intrinsic features of both unliganded enzymes. In DcpS, conformational change is dominated by an anti-symmetric cooperative motion, causing one active site to close as the other opens; however a symmetric motion is also significant. In CS, we identify that both symmetric (suggested by crystallography) and asymmetric motions are features of the protein structure, and as a result the behaviour in solution is largely non-cooperative. The agreement between two modelling approaches using very different levels of theory indicates that the behaviours are indeed intrinsic to the protein structures. Geometric simulations correctly identify and explore large amplitudes of motion, while molecular dynamics simulations indicate the ranges of motion that are energetically feasible. Together, the simulation approaches are able to reveal unexpected functionally relevant motions, and highlight differences between enzymes.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

X-ray, spectroscopic and normal-mode dynamics of calexcitin: structure-function studies of a neuronal calcium-signalling protein

The protein calexcitin was originally identified in molluscan photoreceptor neurons as a 20 kDa molecule which was up-regulated and phosphorylated following a Pavlovian conditioning protocol. Subsequent studies showed that calexcitin regulates the voltage-dependent potassium channel and the calcium-dependent potassium channel as well as causing the release of calcium ions from the endoplasmic reticulum (ER) by binding to the ryanodine receptor. A crystal structure of calexcitin from the squid Loligo pealei showed that the fold is similar to that of another signalling protein, calmodulin, the N- and C-terminal domains of which are known to separate upon calcium binding, allowing interactions with the target protein. Phosphorylation of calexcitin causes it to translocate to the cell membrane, where its effects on membrane excitability are exerted and, accordingly, L. pealei calexcitin contains two protein kinase C phosphorylation sites (Thr61 and Thr188). Thr-to-Asp mutations which mimic phosphorylation of the protein were introduced and crystal structures of the corresponding single and double mutants were determined, which suggest that the C-terminal phosphorylation site (Thr188) exerts the greatest effects on the protein structure. Extensive NMR studies were also conducted, which demonstrate that the wild-type protein predominantly adopts a more open conformation in solution than the crystallographic studies have indicated and, accordingly, normal-mode dynamic simulations suggest that it has considerably greater capacity for flexible motion than the X-ray studies had suggested. Like calmodulin, calexcitin consists of four EF-hand motifs, although only the first three EF-hands of calexcitin are involved in binding calcium ions; the C-terminal EF-hand lacks the appropriate amino acids. Hence, calexcitin possesses two functional EF-hands in close proximity in its N-terminal domain and one functional calcium site in its C-terminal domain. There is evidence that the protein has two markedly different affinities for calcium ions, the weaker of which is most likely to be associated with binding of calcium ions to the protein during neuronal excitation. In the current study, site-directed mutagenesis has been used to abolish each of the three calcium-binding sites of calexcitin, and these experiments suggest that it is the single calcium-binding site in the C-terminal domain of the protein which is likely to have a sensory role in the neuron.

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.