Thermostability in rubredoxin and its relationship to mechanical rigidity

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Uncoupling of Voltage- and Ligand-Induced Activation in HCN2 Channels by Glycine Inserts

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are tetramers that generate electrical rhythmicity in special brain neurons and cardiomyocytes. The channels are activated by membrane hyperpolarization. The binding of cAMP to the four available cyclic nucleotide-binding domains (CNBD) enhances channel activation. We analyzed in the present study the mechanism of how the effect of cAMP binding is transmitted to the pore domain. Our strategy was to uncouple the C-linker (CL) from the channel core by inserting one to five glycine residues between the S6 gate and the A′-helix (constructs 1G to 5G). We quantified in full-length HCN2 channels the resulting functional effects of the inserted glycines by current activation as well as the structural dynamics and statics using molecular dynamics simulations and Constraint Network Analysis. We show functionally that already in 1G the cAMP effect on activation is lost and that with the exception of 3G and 5G the concentration-activation relationships are shifted to depolarized voltages with respect to HCN2. The strongest effect was found for 4G. Accordingly, the activation kinetics were accelerated by all constructs, again with the strongest effect in 4G. The simulations reveal that the average residue mobility of the CL and CNBD domains is increased in all constructs and that the junction between the S6 and A′-helix is turned into a flexible hinge, resulting in a destabilized gate in all constructs. Moreover, for 3G and 4G, there is a stronger downward displacement of the CL-CNBD than in HCN2 and the other constructs, resulting in an increased kink angle between S6 and A′-helix, which in turn loosens contacts between the S4-helix and the CL. This is suggested to promote a downward movement of the S4-helix, similar to the effect of hyperpolarization. In addition, exclusively in 4G, the selectivity filter in the upper pore region and parts of the S4-helix are destabilized. The results provide new insights into the intricate activation of HCN2 channels.

Allosteric signaling in C-linker and cyclic nucleotide-binding domain of HCN2 channels

Opening of hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels is controlled by membrane hyperpolarization and binding of cyclic nucleotides to the tetrameric cyclic nucleotide-binding domain (CNBD), attached to the C-linker (CL) disk. Confocal patch-clamp fluorometry revealed pronounced cooperativity of ligand binding among protomers. However, by which pathways allosteric signal transmission occurs remained elusive. Here, we investigate how changes in the structural dynamics of the CL-CNBD of mouse HCN2 upon cAMP binding relate to inter- and intrasubunit signal transmission. Applying a rigidity-theory-based approach, we identify two intersubunit and one intrasubunit pathways that differ in allosteric coupling strength between cAMP-binding sites or toward the CL. These predictions agree with results from electrophysiological and patch-clamp fluorometry experiments. Our results map out distinct routes within the CL-CNBD that modulate different cAMP-binding responses in HCN2 channels. They signify that functionally relevant submodules may exist within and across structurally discernable subunits in HCN channels.

Systematically Scrutinizing the Impact of Substitution Sites on Thermostability and Detergent Tolerance for Bacillus subtilis Lipase A

Improving an enzyme's (thermo-)stability or tolerance against solvents and detergents is highly relevant in protein engineering and biotechnology. Recent developments have tended toward data-driven approaches, where available knowledge about the protein is used to identify substitution sites with high potential to yield protein variants with improved stability, and subsequently, substitutions are engineered by site-directed or site-saturation (SSM) mutagenesis. However, the development and validation of algorithms for data-driven approaches have been hampered by the lack of availability of large-scale data measured in a uniform way and being unbiased with respect to substitution types and locations. Here, we extend our knowledge on guidelines for protein engineering following a data-driven approach by scrutinizing the impact of substitution sites on thermostability or/and detergent tolerance for Bacillus subtilis lipase A (BsLipA) at very large scale. We systematically analyze a complete experimental SSM library of BsLipA containing all 3439 possible single variants, which was evaluated as to thermostability and tolerances against four detergents under respectively uniform conditions. Our results provide systematic and unbiased reference data at unprecedented scale for a biotechnologically important protein, identify consistently defined hot spot types for evaluating the performance of data-driven protein-engineering approaches, and show that the rigidity theory and ensemble-based approach Constraint Network Analysis yields hot spot predictions with an up to ninefold gain in precision over random classification.

Basal Histamine H4 Receptor Activation: Agonist Mimicry by the Diphenylalanine Motif

Histamine H4 receptor (H4 R) orthologues are G-protein-coupled receptors (GPCRs) that exhibit species-dependent basal activity. In contrast to the basally inactive mouse H4 R (mH4 R), human H4 R (hH4 R) shows a high degree of basal activity. We have performed long-timescale molecular dynamics simulations and rigidity analyses on wild-type hH4 R, the experimentally characterized hH4 R variants S179M, F169V, F169V+S179M, F168A, and on mH4 R to investigate the molecular nature of the differential basal activity. H4 R variant-dependent differences between essential motifs of GPCR activation and structural stabilities correlate with experimentally determined basal activities and provide a molecular explanation for the differences in basal activation. Strikingly, during the MD simulations, F16945.55 dips into the orthosteric binding pocket only in the case of hH4 R, thus adopting the role of an agonist and contributing to the stabilization of the active state. The results shed new light on the molecular mechanism of basal H4 R activation that are of importance for other GPCRs.

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.

Improvement of BsAPA Aspartic Protease Thermostability via Autocatalysis-Resistant Mutation

An aspartic protease gene (Bsapa) was cloned from Bispora sp. MEY-1 and expressed in Pichia pastoris. The recombinant BsAPA showed maximal activity at pH 3.0 and 75 °C and remained stable at 70 °C and below, indicating the thermostable nature of BsAPA. However, heat inactivation still limits the application of BsAPA. To further improve its thermostability, an autocatalysis site (L205-F206) in BsAPA was identified and three mutants (F193W, K204P, and A371V) were generated based on the analysis of the structure neighboring the autocatalysis site. These mutants have improved thermostability, and their half-life at 75 °C increased by 0.5-, 0.2-, and 0.3-fold, respectively. A triple-site mutant (F193W/K204P/A371V) was generated, with 1.5-fold increased half-life at 80 and a 10.7 °C increased T_m, compared with those of the wild-type. These results indicate that autocatalysis of aspartic protease reduces enzyme thermostability. Furthermore, site-directed mutagenesis at regions near the autocatalysis site is an efficient approach to improve aspartic protease thermostability.

Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration

We previously reported that inactivation of the transmembrane taurine transporter (TauT or solute carrier 6a6) causes early retinal degeneration in mice. Compatible with taurine's indispensability for cell volume homeostasis, protein stabilization, cytoprotection, antioxidation, and immuno- and neuromodulation, mice develop multisystemic dysfunctions (hearing loss; liver fibrosis; and behavioral, heart, and skeletal muscle abnormalities) later on. Here, by genetic, cell biologic, in vivo ¹H-magnetic resonance spectroscopy and molecular dynamics simulation studies, we conducted in-depth characterization of a novel disorder: human TAUT deficiency. Loss of TAUT function due to a homozygous missense mutation caused panretinal degeneration in 2 brothers. TAUT_p.A78E still localized in the plasma membrane but is predicted to impact structural stabilization. ³H-taurine uptake by peripheral blood mononuclear cells was reduced by 95%, and taurine levels were severely reduced in plasma, skeletal muscle, and brain. Extraocular dysfunctions were not yet detected, but significantly increased urinary excretion of 8-oxo-7,8-dihydroguanosine indicated generally enhanced (yet clinically unapparent) oxidative stress and RNA oxidation, warranting continuous broad surveillance.-Preising, M. N., Görg, B., Friedburg, C., Qvartskhava, N., Budde, B. S., Bonus, M., Toliat, M. R., Pfleger, C., Altmüller, J., Herebian, D., Beyer, M., Zöllner, H. J., Wittsack, H.-J., Schaper, J., Klee, D., Zechner, U., Nürnberg, P., Schipper, J., Schnitzler, A., Gohlke, H., Lorenz, B., Häussinger, D., Bolz, H. J. Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration.

Improvements of thermophilic enzymes: From genetic modifications to applications

Thermozymes (from thermophiles or hyperthermophiles) offer obvious advantages due to their excellent thermostability, broad pH adaptation, and hydrolysis ability, resulting in diverse industrial applications including food, paper, and textile processing, biofuel production. However, natural thermozymes with low yield and poor adaptability severely hinder their large-scale applications. Extensive studies demonstrated that using genetic modifications such as directed evolution, semi-rational design, and rational design, expression regulations and chemical modifications effectively improved enzyme's yield, thermostability and catalytic efficiency. However, mechanism-based techniques for thermozymes improvements and applications need more attention. In this review, stabilizing mechanisms of thermozymes are summarized for thermozymes improvements, and these improved thermozymes eventually have large-scale industrial applications.

Kinematic Flexibility Analysis: Hydrogen Bonding Patterns Impart a Spatial Hierarchy of Protein Motion

Elastic network models (ENMs) and constraint-based, topological rigidity analysis are two distinct, coarse-grained approaches to study conformational flexibility of macromolecules. In the two decades since their introduction, both have contributed significantly to insights into protein molecular mechanisms and function. However, despite a shared purpose of these approaches, the topological nature of rigidity analysis, and thereby the absence of motion modes, has impeded a direct comparison. Here, we present an alternative, kinematic approach to rigidity analysis, which circumvents these drawbacks. We introduce a novel protein hydrogen bond network spectral decomposition, which provides an orthonormal basis for collective motions modulated by noncovalent interactions, analogous to the eigenspectrum of normal modes. The zero modes decompose proteins into rigid clusters identical to those from topological rigidity, while nonzero modes rank protein motions by their hydrogen bond collective energy penalty. Our kinematic flexibility analysis bridges topological rigidity theory and ENM, enabling a detailed analysis of motion modes obtained from both approaches. Analysis of a large, structurally diverse data set revealed that collectivity of protein motions, reported by the Shannon entropy, is significantly reduced for rigidity theory compared to normal mode approaches. Strikingly, kinematic flexibility analysis suggests that the hydrogen bonding network encodes a protein-fold specific, spatial hierarchy of motions, which goes nearly undetected in ENM. This hierarchy reveals distinct motion regimes that rationalize experimental and simulated protein stiffness variations. Kinematic motion modes highly correlate with reported crystallographic B factors and molecular dynamics simulations of adenylate kinase. A formal expression for changes in free energy derived from the spectral decomposition indicates that motions across nearly 40% of modes obey enthalpy-entropy compensation. Taken together, our results suggest that hydrogen bond networks have evolved to modulate protein structure and dynamics, which can be efficiently probed by kinematic flexibility analysis.

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.

Improving the thermal stability of azoreductase from Halomonas elongata by introducing a disulfide bond via site-directed mutagenesis

Azoreductases mainly reduce azo dyes, the largest class of colorants, to colorless aromatic amines. AzoH, a new azoreductase from the halophilic bacterium, Halomonas elongata, has been recently cloned and expressed in Escherichia coli. The aim of this study was to improve thermal stability of this enzyme by introducing new disulfide bonds. Since X-ray crystallography was not available, homology modeling and molecular dynamics was used to construct the enzyme three-dimensional structure. Potential disulfide bonds for increasing thermal stability were found using DIScover online software. Appropriate mutations (L49C/D108C) to form a disulfide bond were introduced by the Quik-Change method. Mutant protein expressed in E. coli showed increased thermal stability at 50 °C (increased half-life from 12.6 Min in AzoH to 26.66 Min in a mutated enzyme). The mutated enzyme could also tolerate 5% (w/v) NaCl and retained 30% of original activity after 24 H incubation, whereas the wild-type enzyme was completely inactivated. According to circular dichroism studies, the secondary structure was not altered by this mutation; however, a blue shift in intrinsic florescent graph revealed changes in the tertiary structure. This is the first study to improve thermal stability and salt tolerance of a halophilic azoreductase.

Engineering the Enantioselectivity and Thermostability of a (+)-γ-Lactamase from Microbacterium hydrocarbonoxydans for Kinetic Resolution of Vince Lactam (2-Azabicyclo[2.2.1]hept-5-en-3-one)

To produce promising biocatalysts, natural enzymes often need to be engineered to increase their catalytic performance. In this study, the enantioselectivity and thermostability of a (+)-γ-lactamase from Microbacterium hydrocarbonoxydans as the catalyst in the kinetic resolution of Vince lactam (2-azabicyclo[2.2.1]hept-5-en-3-one) were improved. Enantiomerically pure (-)-Vince lactam is the key synthon in the synthesis of antiviral drugs, such as carbovir and abacavir, which are used to fight against HIV and hepatitis B virus. The work was initialized by using the combinatorial active-site saturation test strategy to engineer the enantioselectivity of the enzyme. The approach resulted in two mutants, Val54Ser and Val54Leu, which catalyzed the hydrolysis of Vince lactam to give (-)-Vince lactam, with 99.2% (enantiomeric ratio [E] > 200) enantiomeric excess (ee) and 99.5% ee (E > 200), respectively. To improve the thermostability of the enzyme, 11 residues with high temperature factors (B-factors) calculated by B-FITTER or high root mean square fluctuation (RMSF) values from the molecular dynamics simulation were selected. Six mutants with increased thermostability were obtained. Finally, the mutants generated with improved enantioselectivity and mutants evolved for enhanced thermostability were combined. Several variants showing (+)-selectivity (E value > 200) and improved thermostability were observed. These engineered enzymes are good candidates to serve as enantioselective catalysts for the preparation of enantiomerically pure Vince lactam.IMPORTANCE Enzymatic kinetic resolution of the racemic Vince lactam using (+)-γ-lactamase is the most often utilized means of resolving the enantiomers for the preparation of carbocyclic nucleoside compounds. The efficiency of the native enzymes could be improved by using protein engineering methods, such as directed evolution and rational design. In our study, two properties (enantioselectivity and thermostability) of a γ-lactamase identified from Microbacterium hydrocarbonoxydans were tackled using a semirational design. The protein engineering was initialized by combinatorial active-site saturation test to improve the enantioselectivity. At the same time, two strategies were applied to identify mutation candidates to enhance the thermostability based on calculations from both a static (B-FITTER based on the crystal structure) and a dynamic (root mean square fluctuation [RMSF] values based on molecular dynamics simulations) way. After combining the mutants, we successfully obtained the final mutants showing better properties in both properties. The engineered (+)-lactamase could be a candidate for the preparation of (-)-Vince lactam.

Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery

Allostery describes the functional coupling between sites in biomolecules. Recently, the role of changes in protein dynamics for allosteric communication has been highlighted. A quantitative and predictive description of allostery is fundamental for understanding biological processes. Here, we integrate an ensemble-based perturbation approach with the analysis of biomolecular rigidity and flexibility to construct a model of dynamic allostery. Our model, by definition, excludes the possibility of conformational changes, evaluates static, not dynamic, properties of molecular systems, and describes allosteric effects due to ligand binding in terms of a novel free-energy measure. We validated our model on three distinct biomolecular systems: eglin c, protein tyrosine phosphatase 1B, and the lymphocyte function-associated antigen 1 domain. In all cases, it successfully identified key residues for signal transmission in very good agreement with the experiment. It correctly and quantitatively discriminated between positively or negatively cooperative effects for one of the systems. Our model should be a promising tool for the rational discovery of novel allosteric drugs.

Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch

Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.

Matching Multiple Rigid Domain Decompositions of Proteins

We describe efficient methods for consistently coloring and visualizing collections of rigid cluster decompositions obtained from variations of a protein structure, and lay the foundation for more complex setups, that may involve different computational and experimental methods. The focus here is on three biological applications: the conceptually simpler problems of visualizing results of dilution and mutation analyses, and the more complex task of matching decompositions of multiple Nucleic Magnetic Resonance (NMR) models of the same protein. Implemented into the KINematics And RIgidity (KINARI) web server application, the improved visualization techniques give useful information about protein folding cores, help examining the effect of mutations on protein flexibility and function, and provide insights into the structural motions of Protein Data Bank proteins solved with solution NMR. These tools have been developed with the goal of improving and validating rigidity analysis as a credible coarse-grained model capturing essential information about a protein's slow motions near the native state.

Rational Design of Disulfide Bonds Increases Thermostability of a Mesophilic 1,3-1,4-β-Glucanase from Bacillus terquilensis

1,3-1,4-β-glucanase is an important biocatalyst in brewing industry and animal feed industry, while its low thermostability often reduces its application performance. In this study, the thermostability of a mesophilic β-glucanase from Bacillus terquilensis was enhanced by rational design and engineering of disulfide bonds in the protein structure. Protein spatial configuration was analyzed to pre-exclude the residues pairs which negatively conflicted with the protein structure and ensure the contact of catalytic center. The changes in protein overall and local flexibility among the wild-type enzyme and the designated mutants were predicted to select the potential disulfide bonds for enhancement of thermostability. Two residue pairs (N31C-T187C and P102C-N125C) were chosen as engineering targets and both of them were proved to significantly enhance the protein thermostability. After combinational mutagenesis, the double mutant N31C-T187C/P102C-N125C showed a 48.3% increase in half-life value at 60°C and a 4.1°C rise in melting temperature (Tm) compared to wild-type enzyme. The catalytic property of N31C-T187C/P102C-N125C mutant was similar to that of wild-type enzyme. Interestingly, the optimal pH of double mutant was shifted from pH6.5 to pH6.0, which could also increase its industrial application. By comparison with mutants with single-Cys substitutions, the introduction of disulfide bonds and the induced new hydrogen bonds were proved to result in both local and overall rigidification and should be responsible for the improved thermostability. Therefore, the introduction of disulfide bonds for thermostability improvement could be rationally and highly-effectively designed by combination with spatial configuration analysis and molecular dynamics simulation.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

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.

Does deamidation cause protein unfolding? A top-down tandem mass spectrometry study

Deamidation is a nonenzymatic post-translational modification of asparagine to aspartic acid or glutamine to glutamic acid, converting an uncharged amino acid to a negatively charged residue. It is plausible that deamidation of asparagine and glutamine residues would result in disruption of a proteins' hydrogen bonding network and thus lead to protein unfolding. To test this hypothesis Calmodulin and B2M were deamidated and analyzed using tandem mass spectrometry on a Fourier transform ion cyclotron resonance mass spectrometer (FTICR-MS). The gas phase hydrogen bonding networks of deamidated and nondeamidated protein isoforms were probed by varying the infra-red multi-photon dissociation laser power in a linear fashion and plotting the resulting electron capture dissociation fragment intensities as a melting curve at each amino acid residue. Analysis of the unfolding maps highlighted increased fragmentation at lower laser powers localized around heavily deamidated regions of the proteins. In addition fragment intensities were decreased across the rest of the proteins which we propose is because of the formation of salt-bridges strengthening the intramolecular interactions of the central regions. These results were supported by a computational flexibility analysis of the mutant and unmodified proteins, which would suggest that deamidation can affect the global structure of a protein via modification of the hydrogen bonding network near the deamidation site and that top down FTICR-MS is an appropriate technique for studying protein folding.

Protein thermal stability enhancement by designing salt bridges: a combined computational and experimental study

Protein thermal stability is an important factor considered in medical and industrial applications. Many structural characteristics related to protein thermal stability have been elucidated, and increasing salt bridges is considered as one of the most efficient strategies to increase protein thermal stability. However, the accurate simulation of salt bridges remains difficult. In this study, a novel method for salt-bridge design was proposed based on the statistical analysis of 10,556 surface salt bridges on 6,493 X-ray protein structures. These salt bridges were first categorized based on pairing residues, secondary structure locations, and Cα–Cα distances. Pairing preferences generalized from statistical analysis were used to construct a salt-bridge pair index and utilized in a weighted electrostatic attraction model to find the effective pairings for designing salt bridges. The model was also coupled with B-factor, weighted contact number, relative solvent accessibility, and conservation prescreening to determine the residues appropriate for the thermal adaptive design of salt bridges. According to our method, eight putative salt-bridges were designed on a mesophilic β-glucosidase and 24 variants were constructed to verify the predictions. Six putative salt-bridges leaded to the increase of the enzyme thermal stability. A significant increase in melting temperature of 8.8, 4.8, 3.7, 1.3, 1.2, and 0.7°C of the putative salt-bridges N437K–D49, E96R–D28, E96K–D28, S440K–E70, T231K–D388, and Q277E–D282 was detected, respectively. Reversing the polarity of T231K–D388 to T231D–D388K resulted in a further increase in melting temperatures by 3.6°C, which may be caused by the transformation of an intra-subunit electrostatic interaction into an inter-subunit one depending on the local environment. The combination of the thermostable variants (N437K, E96R, T231D and D388K) generated a melting temperature increase of 15.7°C. Thus, this study demonstrated a novel method for the thermal adaptive design of salt bridges through inference of suitable positions and substitutions.

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

Citrate synthase (CS) catalyses the entry of carbon into the citric acid cycle and is highly-conserved structurally across the tree of life. Crystal structures of dimeric CSs are known in both "open" and "closed" forms, which differ by a substantial domain motion that closes the substrate-binding clefts. We explore both the static rigidity and the dynamic flexibility of CS structures from mesophilic and extremophilic organisms from all three evolutionary domains. The computational expense of this wide-ranging exploration is kept to a minimum by the use of rigidity analysis and rapid all-atom simulations of flexible motion, combining geometric simulation and elastic network modeling. CS structures from thermophiles display increased structural rigidity compared with the mesophilic enzyme. A CS structure from a psychrophile, stabilized by strong ionic interactions, appears to display likewise increased rigidity in conventional rigidity analysis; however, a novel modified analysis, taking into account the weakening of the hydrophobic effect at low temperatures, shows a more appropriate decreased rigidity. These rigidity variations do not, however, affect the character of the flexible dynamics, which are well conserved across all the structures studied. Simulation trajectories not only duplicate the crystallographically observed symmetric open-to-closed transitions, but also identify motions describing a previously unidentified antisymmetric functional motion. This antisymmetric motion would not be directly observed in crystallography but is revealed as an intrinsic property of the CS structure by modeling of flexible motion. This suggests that the functional motion closing the binding clefts in CS may be independent rather than symmetric and cooperative.

Redistribution of flexibility in stabilizing antibody fragment mutants follows Le Châtelier's principle

Le Châtelier's principle is the cornerstone of our understanding of chemical equilibria. When a system at equilibrium undergoes a change in concentration or thermodynamic state (i.e., temperature, pressure, etc.), La Châtelier's principle states that an equilibrium shift will occur to offset the perturbation and a new equilibrium is established. We demonstrate that the effects of stabilizing mutations on the rigidity ⇔ flexibility equilibrium within the native state ensemble manifest themselves through enthalpy-entropy compensation as the protein structure adjusts to restore the global balance between the two. Specifically, we characterize the effects of mutation to single chain fragments of the anti-lymphotoxin-β receptor antibody using a computational Distance Constraint Model. Statistically significant changes in the distribution of both rigidity and flexibility within the molecular structure is typically observed, where the local perturbations often lead to distal shifts in flexibility and rigidity profiles. Nevertheless, the net gain or loss in flexibility of individual mutants can be skewed. Despite all mutants being exclusively stabilizing in this dataset, increased flexibility is slightly more common than increased rigidity. Mechanistically the redistribution of flexibility is largely controlled by changes in the H-bond network. For example, a stabilizing mutation can induce an increase in rigidity locally due to the formation of new H-bonds, and simultaneously break H-bonds elsewhere leading to increased flexibility distant from the mutation site via Le Châtelier. Increased flexibility within the VH β4/β5 loop is a noteworthy illustration of this long-range effect.

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

Efficient and robust analysis of biomacromolecular flexibility using ensembles of network topologies based on fuzzy noncovalent constraints

We describe an approach (ENT(FNC)) for performing rigidity analyses of biomacromolecules on ensembles of network topologies (ENT) generated from a single input structure. The ENT is based on fuzzy noncovalent constraints, which considers thermal fluctuations of biomacromolecules without actually sampling conformations. Definitions for fuzzy noncovalent constraints were derived from persistency data from molecular dynamics (MD) simulations. A very good agreement between local flexibility and rigidity characteristics from ENT(FNC) and MD simulations-generated ensembles is found. Regarding global characteristics, convincing results were obtained when relative thermostabilities of citrate synthase and lipase A structures were computed. The ENT(FNC) approach significantly improves the robustness of rigidity analyses, is highly efficient, and does not require a protein-specific parameterization. Its low computational demand makes it especially valuable for the analysis of large data sets, e.g., for data-driven protein engineering.

Constraint Network Analysis (CNA): a Python software package for efficiently linking biomacromolecular structure, flexibility, (thermo-)stability, and function

For deriving maximal advantage from information on biomacromolecular flexibility and rigidity, results from rigidity analyses must be linked to biologically relevant characteristics of a structure. Here, we describe the Python-based software package Constraint Network Analysis (CNA) developed for this task. CNA functions as a front- and backend to the graph-based rigidity analysis software FIRST. CNA goes beyond the mere identification of flexible and rigid regions in a biomacromolecule in that it (I) provides a refined modeling of thermal unfolding simulations that also considers the temperature-dependence of hydrophobic tethers, (II) allows performing rigidity analyses on ensembles of network topologies, either generated from structural ensembles or by using the concept of fuzzy noncovalent constraints, and (III) computes a set of global and local indices for quantifying biomacromolecular stability. This leads to more robust results from rigidity analyses and extends the application domain of rigidity analyses in that phase transition points ("melting points") and unfolding nuclei ("structural weak spots") are determined automatically. Furthermore, CNA robustly handles small-molecule ligands in general. Such advancements are important for applying rigidity analysis to data-driven protein engineering and for estimating the influence of ligand molecules on biomacromolecular stability. CNA maintains the efficiency of FIRST such that the analysis of a single protein structure takes a few seconds for systems of several hundred residues on a single core. These features make CNA an interesting tool for linking biomacromolecular structure, flexibility, (thermo-)stability, and function. CNA is available from http://cpclab.uni-duesseldorf.de/software for nonprofit organizations.

Global and local indices for characterizing biomolecular flexibility and rigidity

Understanding flexibility and rigidity characteristics of biomolecules is a prerequisite for understanding biomolecular structural stability and function. Computational methods have been implemented that directly characterize biomolecular flexibility and rigidity by constraint network analysis. For deriving maximal advantage from these analyses, their results need to be linked to biologically relevant characteristics of a structure. Such links are provided by global and local measures ("indices") of biomolecular flexibility and rigidity. To date, more than 14 indices are available with sometimes overlapping or only vague definitions. We present concise definitions of these indices, analyze the relation between, and the scope and limitations of them, and compare their informative value. For this, we probe the structural stability of the calcium binding protein α-lactalbumin as a showcase, both in the "ground state" and after perturbing the system by changing the network topology. In addition, we introduce three indices for the first time that extend the application domain of flexibility and rigidity analyses. The results allow us to provide guidelines for future studies suggesting which of these indices could best be used for analyzing, understanding, and quantifying structural features that are important for biomolecular stability and function. Finally, we make suggestions for proper index notations in future studies to prevent the misinterpretation and to facilitate the comparison of results obtained from flexibility and rigidity analyses.

A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35 V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35 V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.

Inhibition of HIV-1 protease: the rigidity perspective

MOTIVATION

HIV-1 protease is a key drug target due to its role in the life cycle of the HIV-1 virus. Rigidity analysis using the software First is a computationally inexpensive method for inferring functional information from protein crystal structures. We evaluate the rigidity of 206 high-resolution (2 Å or better) X-ray crystal structures of HIV-1 protease and compare the effects of different inhibitors binding to the enzyme.

RESULTS

Inhibitor binding has little effect on the overall rigidity of the protein homodimer, including the rigidity of the active site. The principal effect of inhibitor binding on rigidity is to constrain the flexibility of the β-hairpin flaps, which move to allow access to the active site of the enzyme. We show that commercially available antiviral drugs which target HIV-1 protease can be divided into two classes, those which significantly affect flap rigidity and those which do not. The non-peptidic inhibitor tipranavir is distinctive in its consistently strong effect on flap rigidity.

CONTACT

jack.heal@warwick.ac.uk; r.roemer@warwick.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Changes in Lysozyme Flexibility upon Mutation Are Frequent, Large and Long-Ranged

We investigate changes in human c-type lysozyme flexibility upon mutation via a Distance Constraint Model, which gives a statistical mechanical treatment of network rigidity. Specifically, two dynamical metrics are tracked. Changes in flexibility index quantify differences within backbone flexibility, whereas changes in the cooperativity correlation quantify differences within pairwise mechanical couplings. Regardless of metric, the same general conclusions are drawn. That is, small structural perturbations introduced by single point mutations have a frequent and pronounced affect on lysozyme flexibility that can extend over long distances. Specifically, an appreciable change occurs in backbone flexibility for 48% of the residues, and a change in cooperativity occurs in 42% of residue pairs. The average distance from mutation to a site with a change in flexibility is 17–20 Å. Interestingly, the frequency and scale of the changes within single point mutant structures are generally larger than those observed in the hen egg white lysozyme (HEWL) ortholog, which shares 61% sequence identity with human lysozyme. For example, point mutations often lead to substantial flexibility increases within the β-subdomain, which is consistent with experimental results indicating that it is the nucleation site for amyloid formation. However, β-subdomain flexibility within the human and HEWL orthologs is more similar despite the lowered sequence identity. These results suggest compensating mutations in HEWL reestablish desired properties.

The functional importance of protein dynamics is universally accepted, making the study of dynamical similarities and differences among proteins of the same function an intriguing problem. While some metrics are likely to be conserved across family, differences are also very common. In previous works we have used a Distance Constraint Model to quantify flexibility differences across sets of orthologous proteins, which reproduce this diversity. In the same manner, this work investigates changes occurring upon individual point mutations. Somewhat surprisingly, the small structural perturbations caused by mutation lead to changes throughout the protein. These changes can be quite large, actually surpassing the scale for differences between ortholog pairs. Moreover, changes in flexibility frequently occur at sites far from the mutation site. These results underscore the sensitivity of protein dynamics in connection with allostery, and help explain why differences across protein families are so common.

Flexibility analysis of biomacromolecules with application to computer-aided drug design

Flexibility characteristics of biomacromolecules can be efficiently determined down to the atomic level by a graph-theoretical technique as implemented in the FIRST (Floppy Inclusion and Rigid Substructure Topology) and ProFlex software packages. The method has been successfully applied to a series of protein and nucleic acid structures. Here, we describe practical guidelines for setting up and performing a flexibility analysis, discuss current bottlenecks of the approach, and provide sample applications as to how this technique can support computer-aided drug design approaches.

Double mutation in photosystem II reaction centers and elevated CO2 grant thermotolerance to mesophilic cyanobacterium

Photosynthetic biomass production rapidly declines in mesophilic cyanobacteria grown above their physiological temperatures largely due to the imbalance between degradation and repair of the D1 protein subunit of the heat susceptible Photosystem II reaction centers (PSIIRC). Here we show that simultaneous replacement of two conserved residues in the D1 protein of the mesophilic Synechocystis sp. PCC 6803, by the analogue residues present in the thermophilic Thermosynechococcus elongatus, enables photosynthetic growth, extensive biomass production and markedly enhanced stability and repair rate of PSIIRC for seven days even at 43 °C but only at elevated CO(2) (1%). Under the same conditions, the Synechocystis control strain initially presented very slow growth followed by a decline after 3 days. Change in the thylakoid membrane lipids, namely the saturation of the fatty acids is observed upon incubation for the different strains, but only the double mutant shows a concomitant major change of the enthalpy and entropy for the light activated Q(A)(-)→Q(B) electron transfer, rendering them similar to those of the thermophilic strain. Following these findings, computational chemistry and protein dynamics simulations we propose that the D1 double mutation increases the folding stability of the PSIIRC at elevated temperatures. This, together with the decreased impairment of D1 protein repair under increased CO(2) concentrations result in the observed photothermal tolerance of the photosynthetic machinery in the double mutant.

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.

Solution-state NMR structure and biophysical characterization of zinc-substituted rubredoxin B (Rv3250c) from Mycobacterium tuberculosis

Owing to the evolution of multi-drug-resistant and extremely drug-resistant Mycobacterium tuberculosis strains, there is an urgent need to develop new antituberculosis strategies to prevent TB epidemics in the industrial world. Among the potential new drug targets are two small nonheme iron-binding proteins, rubredoxin A (Rv3251c) and rubredoxin B (Rv3250c), which are believed to play a role in electron-transfer processes. Here, the solution structure and biophysical properties of one of these two proteins, rubredoxin B (Mt-RubB), determined in the zinc-substituted form are reported. The zinc-substituted protein was prepared by expressing Mt-RubB in minimal medium containing excess zinc acetate. Size-exclusion chromatography and NMR spectroscopy indicated that Mt-RubB was a monomer in solution. The structure (PDB entry 2kn9) was generally similar to those of other rubredoxins, containing a three-stranded antiparallel β-sheet (β2-β1-β3) and a metal tetrahedrally coordinated to the S atoms of four cysteine residues (Cys9, Cys12, Cys42 and Cys45). The first pair of cysteine residues is at the C-terminal end of the first β-strand and the second pair of cysteine residues is towards the C-terminal end of the loop between β2 and β3. The structure shows the metal buried deeply within the protein, an observation that is supported by the inability to remove the metal with excess EDTA at room temperature. Circular dichroism spectroscopy shows that this stability extends to high temperature, with essentially no change being observed in the CD spectrum of Mt-RubB upon heating to 353 K.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.