Tryptophan phosphorescence study of enzyme flexibility and unfolding in laboratory-evolved thermostable esterases

Directed evolution of p-nitrobenzyl esterase (pNB E) has yielded eight generations of increasingly thermostable variants. The most stable esterase, 8G8, has 13 amino acid substitutions, a melting temperature 17 degrees C higher than the wild-type enzyme, and increased hydrolytic activity toward p-nitrophenyl acetate (pNPA), the substrate used for evolution, at all temperatures. Room-temperature activities of the evolved thermostable variants range from 3.5 times greater to 4.0 times less than wild type. The relationships between enzyme stability, catalytic activity, and flexibility for the esterases were investigated using tryptophan phosphorescence. We observed no correlation between catalytic activity and enzyme flexibility in the vicinity of the tryptophan (Trp) residues. Increases in stability, however, are often accompanied by decreases in flexibility, as measured by Trp phosphorescence. Phosphorescence data also suggest that the N- and C-terminal regions of pNB E unfold independently. The N-terminal region appears more thermolabile, yet most of the thermostabilizing mutations are located in the C-terminal region. Mutational studies show that the effects of the N-terminal mutations depend on one or more mutations in the C-terminal region. Thus, the pNB E mutants are stabilized by long-range, cooperative interactions between distant parts of the enzyme.     

Time resolved fluorescence and phosphorescence properties of the individual tryptophan residues of barnase: evidence for protein-protein interactions.

Steady-state and time-resolved fluorescence, as well as phosphorescence measurements, were used to resolve the luminescence properties of the three individual tryptophan residues of barnase. Assignment of the fluorescence properties was performed using single-tryptophan-containing mutants and the results were compared with the information available from the study of wild-type and two-tryptophan-containing mutants (Willaert, Lowenthal, Sancho, Froeyen, Fersht, Engelborghs, Biochemistry 1992;31:711-716). The fluorescence and the phosphorescence emission of wild-type barnase is dominated by Trp35, although Trp71 has the strongest intrinsic fluorescence when present alone. Fluorescence emission of these two tryptophan residues is blue-shifted and pH-independent. The fluorescence decay parameters of Trp94 are pH-dependent, and an intramolecular collision frequency of 2 to 5 x 10(9) s(-1) between Trp94 and His18 is calculated. Fluorescence emission of Trp94 is red-shifted. Fluorescence anisotropy decay reveals the local mobility of the individual tryptophan residues and this result correlates well with their phosphorescence properties. Trp35 and Trp71 display a single phosphorescence lifetime, which reflects the rigidity of their environment. Surface Trp94 does not exhibit detectable phosphorescence emission. The existence of energy transfer between Trp71 and Trp94, as previously detected by fluorescence measurements, is also observed in the phosphorescence emission of barnase. Dynamic quenching causes the phosphorescence intensity to be protein-concentration dependent. In addition, fluorescence anisotropy shows concentration dependency, and this can be described by the formation of trimers in solution.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Tryptophan phosphorescence study of the internal dynamics of human erythrocyte membrane proteins upon spectrin modification

Room-temperature tryptophan phosphorescence has been used analyze the slow (millisecond) internal dynamics of proteins in isolated native human erythrocyte membranes, after removal of 95% of spectrin, and after thermal denaturation of spectrin or medium acidification to pH 6.0–4.0, as well as the internal dynamics of spectrin extracted from the membrane in solution. The integral membrane proteins prove to differ sharply from spectrin in their structural and dynamic state. The millisecond movements of structural elements in integral proteins are considerably hindered as compared with spectrin. Removal of the bulk of spectrin from membranes leads to amplification of slow fluctuations in the structure of integral proteins. This suggests involvement of spectrin in the control of the structural and dynamic state of the erythrocyte membrane proteins. The acidification of the medium to pH 6.0–4.0 decreases the internal dynamics of native membrane proteins, which is explained by the pH-induced aggregation of spectrin. After thermal denaturation of spectrin, there is no pH-induced increase in the rigidity of the structure of membrane proteins.     

Disease-associated mutations impacting BC-loop flexibility trigger long-range transthyretin tetramer destabilization and aggregation

Hereditary transthyretin amyloidosis (ATTR) is an autosomal dominant disease characterized by the extracellular deposition of the transport protein transthyretin (TTR) as amyloid fibrils. Despite the progress achieved in recent years, understanding why different TTR residue substitutions lead to different clinical manifestations remains elusive. Here, we studied the molecular basis of disease-causing missense mutations affecting residues R34 and K35. R34G and K35T variants cause vitreous amyloidosis, whereas R34T and K35N mutations result in amyloid polyneuropathy and restrictive cardiomyopathy. All variants are more sensitive to pH-induced dissociation and amyloid formation than the wild-type (WT)-TTR counterpart, specifically in the variants deposited in the eyes amyloid formation occurs close to physiological pHs. Chemical denaturation experiments indicate that all the mutants are less stable than WT-TTR, with the vitreous amyloidosis variants, R34G and K35T, being highly destabilized. Sequence-induced stabilization of the dimer–dimer interface with T119M rendered tetramers containing R34G or K35T mutations resistant to pH-induced aggregation. Because R34 and K35 are among the residues more distant to the TTR interface, their impact in this region is therefore theorized to occur at long range. The crystal structures of double mutants, R34G/T119M and K35T/T119M, together with molecular dynamics simulations indicate that their strong destabilizing effect is initiated locally at the BC loop, increasing its flexibility in a mutation-dependent manner. Overall, the present findings help us to understand the sequence-dynamic-structural mechanistic details of TTR amyloid aggregation triggered by R34 and K35 variants and to link the degree of mutation-induced conformational flexibility to protein aggregation propensity.     

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all‐atom molecular dynamics simulations to study how gradual shortening a very long or solvent‐exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability. Proteins 2015; 83:2137–2146. © 2015 Wiley Periodicals, Inc.     

Molecular dynamics simulations of Hydrogenobacter thermophilus cytochrome c552: comparisons of the wild-type protein, a b-type variant, and the apo state

Molecular dynamic simulations have been performed for wild-type Hydrogenobacter thermophilus cytochrome c(552), a b-type variant of the protein, and the apo state with the heme prosthetic group removed. In the b-type variant, Cys 10 and Cys 13 were mutated to alanine residues, and so the heme group was no longer covalently bound to the protein. Two 8-ns simulations have been performed for each system at 298 and 360 K. The simulations of the wild-type protein at 298 K show a very close agreement with experimental NMR data. A fluxional process involving the side chain of Met 59, which coordinates to the heme iron, is observed in accord with proposals from NMR studies. Overall, the structure and dynamical behavior of the protein during the simulations of the b-type variant is closely similar to that of the wild-type protein. However, side chains in the heme-binding site show larger fluctuations in the b-type variant simulation at 360 K. In addition, structural changes are seen for a number of residues close to the heme group, particularly Gly 22 and Ser 51. The simulations of the apo state show significant conformational changes for residues 50-59. These residues form a loop region, which packs over the heme group in the wild-type protein and hydrogen bonds to the heme propionate groups. In the absence of heme, in the apo state simulations, these residues form short but persistent regions of beta-sheet secondary structure. These could provide nucleation sites for the conversion to amyloid fibrils.     

Early events during folding of wild-type staphylococcal nuclease and a single-tryptophan variant studied by ultrarapid mixing

A continuous-flow mixing device with a dead time of 100 micros coupled with intrinsic tryptophan and 1-anilinonaphthalene-8-sulfonate (ANS) fluorescence was used to monitor structure formation during early stages of the folding of staphylococcal nuclease (SNase). A variant with a unique tryptophan fluorophore in the N-terminal beta-barrel domain (Trp76 SNase) was obtained by replacing the single Trp140 in wild-type SNase with His in combination with Trp substitution of Phe76. A common background of P47G, P117G and H124L mutations was chosen in order to stabilize the protein and prevent accumulation of cis proline isomers under native conditions. In contrast to WT(*) SNase, which shows no changes in tryptophan fluorescence prior to the rate-limiting folding step ( approximately 100 ms), the F76W/W140H variant shows additional changes (enhancement) during an early folding phase with a time constant of 75 micros. Both proteins exhibit a major increase in ANS fluorescence and identical rates for this early folding event. These findings are consistent with the rapid accumulation of an ensemble of states containing a loosely packed hydrophobic core involving primarily the beta-barrel domain while the specific interactions in the alpha-helical domain involving Trp140 are formed only during the final stages of folding. The fact that both variants exhibit the same number of kinetic phases with very similar rates confirms that the folding mechanism is not perturbed by the F76W/W140H mutations. However, the Trp at position 76 reports on the rapid formation of a hydrophobic cluster in the N-terminal beta-sheet region while the wild-type Trp140 is silent during this early stage of folding. Quantitative modeling of the (un)folding kinetics and thermodynamics of these two proteins versus urea concentration revealed that the F76W/W140H mutation selectively destabilizes the native state relative to WT(*) SNase while the stability of transient intermediates remains unchanged, leading to accumulation of intermediates under equilibrium conditions at moderate denaturant concentrations.     

Effect of single point mutations in a form of systemic amyloidosis

Amyloid deposits of light-chain proteins are associated with the most common form of systemic amyloidosis. We have studied the effects of single point mutations on amyloid formation of these proteins using explicit solvent model molecular dynamics simulations. For this purpose, we compare the stability of the wild-type immunoglobulin light-chain protein REI in its native and amyloid forms with that of four mutants: R61N, G68D, D82I, and A84T. We argue that the experimentally observed differences in the propensity for amyloid formation result from two effects. First, the mutant dimers have a lower stability than the wild-type dimer due to increase exposure of certain hydrophobic residues. The second effect is a shift in equilibrium between monomers with amyloid-like structure and such with native structures. Hence, when developing drugs against light-chain associated systemic amyloidosis, one should look for components that either stabilize the dimer by binding to the dimer interface or reduce for the monomers the probability of the amyloid form.     

A statistical approach to the interpretation of molecular dynamics simulations of calmodulin equilibrium dynamics

A sample of 35 independent molecular dynamics (MD) simulations of calmodulin (CaM) equilibrium dynamics was prepared from different but equally plausible initial conditions (20 simulations of the wild-type protein and 15 simulations of the D129N mutant). CaM's radius of gyration and backbone mean-square fluctuations were analyzed for the effect of the D129N mutation, and simulations were compared with experiments. Statistical tests were employed for quantitative comparisons at the desired error level. The computational model predicted statistically significant compaction of CaM relative to the crystal structure, consistent with the results of small-angle X-ray scattering (SAXS) experiments. This effect was not observed in several previously reported studies of (Ca2+)(4)-CaM, which relied on a single MD run. In contrast to radius of gyration, backbone mean-square fluctuations showed a distinctly non-normal and positively skewed distribution for nearly all residues. Furthermore, the D129N mutation affected the backbone dynamics in a complex manner and reduced the mobility of Glu123, Met124, Ile125, Arg126, and Glu127 located in the adjacent alpha-helix G. The implications of these observations for the comparisons of MD simulations with experiments are discussed. The proposed approach may be useful in studies of protein equilibrium dynamics where MD simulations fall short of properly sampling the conformational space, and when the comparison with experiments is affected by the reproducibility of the computational model.     

Key residue-dominated protein folding dynamics

A “key-residue” hypothesis that a few residues’ characteristics contain the essential dynamics of the whole protein is proposed for the study of side-chain relaxation near native states. Molecular dynamics simulation is performed on the folding of Trp-cage, and four key residues are discovered and shown to be highly sensitive to the change of state of the protein away from the native state. Order parameters that characterize the geometrical properties of key residues are shown to form valuable phase plane on which one distinguishes different reaction pathways. Furthermore, one of the key residues, Trp6, is observed to display two reconfiguration processes, in which one is induced by an unconstrained torsion of the side-chain of Trp6, with a rate faster by almost an order of magnitude than the other one described by Kussell’s model. The faster process seems to occur more frequently in our simulation and thus represent a significant mechanism in folding dynamics.     

Effect of heavy water on protein flexibility

The effects of heavy water (D(2)O) on internal dynamics of proteins were assessed by both the intrinsic phosphorescence lifetime of deeply buried Trp residues, which reports on the local structure about the triplet probe, and the bimolecular acrylamide phosphorescence quenching rate constant that is a measure of the average acrylamide diffusion coefficient through the macromolecule. The results obtained with several protein systems (ribonuclease T1, superoxide dismutase, beta-lactoglobulin, liver alcohol dehydrogenase, alkaline phosphatase, and apo- and Cd-azurin) demonstrate that in most cases D(2)O does significantly increase the rigidity the native structure. With the exception of alkaline phosphatase, the kinetics of the structure tightening effect of deuteration are rapid compared with the rate of H/D exchange of internal protons, which would then assign the dampening of structural fluctuations in D(2)O to a solvent effect, rather than to stronger intramolecular D bonding. Structure tightening by heavy water is generally amplified at higher temperatures, supporting a mostly hydrophobic nature of the underlying interaction, and under conditions that destabilize the globular fold.     

Wild type and mutants of the HET‐s(218–289) prion show different flexibility at fibrillar ends: A simulation study

The C‐terminal segment (residues 218–289) of the HET‐s protein of the filamentous fungus Podosporina anserina is a prion‐forming domain. The structural model of the HET‐s(218–289) amyloid fibril based on solid‐state nuclear magnetic resonance (NMR) restraints shows a β solenoid topology which is comprised of a β‐sheet core and interconnecting loops. For the single‐point mutants Phe286Ala and Trp287Ala, slower aggregation rates in vitro and loss of prionic infectivity have been reported recently. Here we have used molecular dynamics to compare the flexibility of the mutants and wild type. The simulations, initiated from a trimeric aggregate extracted from the NMR structural model, show structural stability on a 100‐ns time scale for wild type and mutants. Analysis of the fluctuations along the simulations reveals that the mutants are less flexible than the wild type in the C‐terminal segment at only one of the two external monomers. Analysis of interaction energy and buried accessible surface indicates that residue Phe286 in particular is stabilized in the Trp287Ala mutant. The simulation results provide an atomistic explanation of the suggestion (based on indirect experimental evidence) that flexibility at the protofibril end(s) is required for fibril elongation. Moreover, they provide further evidence that the growth of the HET‐s amyloid fibril is directional. Proteins 2014; 82:399–404. © 2013 Wiley Periodicals, Inc.     

Comparative docking studies of CYP1b1 and its PCG-associated mutant forms

Molecular docking has been used to compare and contrast the binding modes of oestradiol with the wild-type and some disease-associated mutant forms of the human CYP1b1 protein. The receptor structures used for docking were derived from molecular dynamics simulations of homology-modelled structures. Earlier studies involving molecular dynamics and principal component analysis indicated that mutations could have a disruptive effect on function, by destabilizing the native properties of the functionally important regions, especially those of the haem-binding and substrate-binding regions, which constitute the site of catalytic activity of the enzyme. In order to gain more insights into the possible differences in substrate-binding and catalysis between the wild-type and mutant proteins, molecular docking studies were carried out. Mutants showed altered protein-ligand interactions compared with the wild-type as a consequence of changes in the geometry of the substrate-binding region and in the position of haem relative to the active site. An important difference in ligand-protein interactions between the wild-type and mutants is the presence of stacking interaction with phenyl residues in the wild-type, which is either completely absent or considerably weaker in mutants. The present study revealed essential differences in the interactions between ligand and protein in wild-type and disease mutants, and helped in understanding the deleterious nature of disease mutations at the level of molecular function.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Activity and spectroscopic properties of bacterial D-amino acid transaminase after multiple site-directed mutagenesis of a single tryptophan residue

One of the three tryptophan residues per subunit of thermostable D-amino acid transaminase, Trp-139, is close to the active-site Lys-145 in the sequence of the protein. This tryptophan has been changed to several other types of residues by site-directed mutagenesis. The only mutant protein that was sufficiently active and stable for study had Phe substituted for Trp (W139F). The spectroscopic properties of this mutant enzyme differed from those of the wild-type transaminase. For example, denatured W139F showed the expected decrease in fluorescence emission intensity at 350 nm due to the deletion of one Trp residue, but the fluorescence emission of the wild-type and W139F enzymes in the native state did not differ in intensity. This result suggests that the fluorescence of Trp-139 in the native, wild-type enzyme is not manifested perhaps due to its proximity to the coenzyme, pyridoxal phosphate. Results of energy-transfer studies at several wavelengths could also be interpreted as due to the proximity of Trp-139 and the coenzyme. Circular dichroism studies indicated that the negative Cotton effect at 420 nm due to the coenzyme was still present in W139F. However, the 280-nm optically active band present in the wild-type enzyme was greatly diminished in W139F. The mutant protein with Asp at position 139 (W139D) could not be isolated presumably because it was degraded. The other mutant enzymes, W139P, W139A, and W139H, were isolated with partial activities (15-35%) that were slowly lost upon storage at 4 degrees C. Overall, these results indicate the importance of Trp-139 in the thermostable D-amino acid transaminase.     

Enhanced thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 by rational engineering of a glycine to proline mutation

Protein thermostability can be increased by some glycine to proline mutations in a target protein. However, not all glycine to proline mutations can improve protein thermostability, and this method is suitable only at carefully selected mutation sites that can accommodate structural stabilization. In this study, homology modeling and molecular dynamics simulations were used to select appropriate glycine to proline mutations to improve protein thermostability, and the effect of the selected mutations was proved by the experiments. The structure of methyl parathion hydrolase (MPH) from Ochrobactrum sp. M231 (Ochr-MPH) was constructed by homology modeling, and molecular dynamics simulations were performed on the modeled structure. A profile of the root mean square fluctuations of Ochr-MPH was calculated at the nanosecond timescale, and an eight-amino acid loop region (residues 186-193) was identified as having high conformational fluctuation. The two glycines nearest to this region were selected as mutation targets that might affect protein flexibility in the vicinity. The structures and conformational fluctuations of two single mutants (G194P and G198P) and one double mutant (G194P/G198P) were modeled and analyzed using molecular dynamics simulations. The results predicted that the mutant G194P had the decreased conformational fluctuation in the loop region and might increase the thermostability of Ochr-MPH. The thermostability and kinetic behavior of the wild-type and three mutant enzymes were measured. The results were consistent with the computational predictions, and the mutant G194P was found to have higher thermostability than the wild-type enzyme.     

Phage P22 tailspike protein: removal of head-binding domain unmasks effects of folding mutations on native-state thermal stability.

A shortened, recombinant protein comprising residues 109-666 of the tailspike endorhamnosidase of Salmonella phage P22 was purified from Escherichia coli and crystallized. Like the full-length tailspike, the protein lacking the amino-terminal head-binding domain is an SDS-resistant, thermostable trimer. Its fluorescence and circular dichroism spectra indicate native structure. Oligosaccharide binding and endoglycosidase activities of both proteins are identical. A number of tailspike folding mutants have been obtained previously in a genetic approach to protein folding. Two temperature-sensitive-folding (tsf) mutations and the four known global second-site suppressor (su) mutations were introduced into the shortened protein and found to reduce or increase folding yields at high temperature. The mutational effects on folding yields and subunit folding kinetics parallel those observed with the full-length protein. They mirror the in vivo phenotypes and are consistent with the substitutions altering the stability of thermolabile folding intermediates. Because full-length and shortened tailspikes aggregate upon thermal denaturation, and their denaturant-induced unfolding displays hysteresis, kinetics of thermal unfolding were measured to assess the stability of the native proteins. Unfolding of the shortened wild-type protein in the presence of 2% SDS at 71 degrees C occurs at a rate of 9.2 x 10(-4) s(-1). It reflects the second kinetic phase of unfolding of the full-length protein. All six mutations were found to affect the thermal stability of the native protein. Both tsf mutations accelerate thermal unfolding about 10-fold. Two of the su mutations retard thermal unfolding up to 5-fold, while the remaining two mutations accelerate unfolding up to 5-fold. The mutational effects can be rationalized on the background of the recently determined crystal structure of the protein.