Solution structure and dynamics of a designed hydrophobic core variant of ubiquitin.

BACKGROUND: The recent merger of computation and protein design has resulted in a burst of success in the generation of novel proteins with native-like properties. A critical component of this coupling between theory and experiment is a detailed analysis of the structures and stabilities of designed proteins to assess and improve the accuracy of design algorithms. RESULTS: Here we report the solution structure of a hydrophobic core variant of ubiquitin, referred to as 1D7, which was designed with the core-repacking algorithm ROC. As a measure of conformational specificity, we also present amide exchange protection factors and backbone and sidechain dynamics. The results indicate that 1D7 is similar to wild-type (WT) ubiquitin in backbone structure and degree of conformational specificity. We also observe a good correlation between experimentally determined sidechain structures and those predicted by ROC. However, evaluation of the core sidechain conformations indicates that, in general, 1D7 has more sidechains in less statistically favorable conformations than WT. CONCLUSIONS: Our results provide an explanation for the lower stability of 1D7 compared to WT, and suggest modifications to design algorithms that may improve the accuracy with which structure and stability are predicted. The results also demonstrate that core packing can affect conformational flexibility in subtle ways that are likely to be important for the design of function and protein-ligand interactions.     

Hydrophobic core fluidity of homologous protein domains: relation of side-chain dynamics to core composition and packing

The side-chain dynamics of methyl groups in two structurally related proteins from the fibronectin type III (fnIII) superfamily, the third fnIII domain from human tenascin (TNfn3) and the tenth fnIII domain from human fibronectin (FNfn10), have been studied by NMR spectroscopy. Side-chain order parameters reveal that the hydrophobic cores of the two proteins have substantially different mobilities. The core of TNfn3 is very dynamic, with exceptionally low order parameters for the most deeply buried residues, while that of FNfn10 is more like those of other proteins which have been studied with this technique, having a relatively rigid core with uniformly distributed dynamics. The unusually dynamic core of TNfn3 appears to be related to its amino acid composition, which makes it more fluid-like. A further explanation for the mobility of the TNfn3 core may be found in the negative correlation between the order parameter and excess packing volume, which shows that the core of TNfn3 is less densely packed and consequently has lower methyl order parameters for its buried residues. Rotameric transitions, presumably facilitated by the lower packing density, appear to make an important contribution to lowering the order parameters, and have been probed by measuring three-bond scalar couplings. Overall, although backbone dynamics is generally similar for proteins with the same topology on a fast time scale (picoseconds to nanoseconds), this study shows that a single fold can accommodate a wide variation in the dynamic properties of its core.     

Dynamically-Driven Enhancement of the Catalytic Machinery of the SARS 3C-Like Protease by the S284-T285-I286/A Mutations on the Extra Domain

Previously we revealed that the extra domain of SARS 3CLpro mediated the catalysis via different mechanisms. While the R298A mutation completely abolished the dimerization, thus resulting in the inactive catalytic machinery, N214A inactivated the enzyme by altering its dynamics without significantly perturbing its structure. Here we studied another mutant with S284-T285-I286 replaced by Ala (STI/A) with a 3.6-fold activity increase and slightly enhanced dimerization. We determined its crystal structure, which still adopts the dimeric structure almost identical to that of the wild-type (WT), except for slightly tighter packing between two extra-domains. We then conducted 100-ns molecular dynamics (MD) simulations for both STI/A and WT, the longest reported so far for 3CLpro. In the simulations, two STI/A extra domains become further tightly packed, leading to a significant volume reduction of the nano-channel formed by residues from both catalytic and extra domains. The enhanced packing appears to slightly increase the dynamic stability of the N-finger and the first helix residues, which subsequently triggers the redistribution of dynamics over residues directly contacting them. This ultimately enhances the dynamical stability of the residues constituting the catalytic dyad and substrate-binding pockets. Further correlation analysis reveals that a global network of the correlated motions exists in the protease, whose components include all residues identified so far to be critical for the dimerization and catalysis. Most strikingly, the N214A mutation globally decouples this network while the STI/A mutation alters the correlation pattern. Together with previous results, the present study establishes that besides the classic structural allostery, the dynamic allostery also operates in the SARS 3CLpro, which is surprisingly able to relay the perturbations on the extra domain onto the catalytic machinery to manifest opposite catalytic effects. Our results thus imply a promising avenue to design specific inhibitors for 3CL proteases by disrupting their dynamic correlation network.     

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.     

Design and structural analysis of an engineered thermostable chicken lysozyme.

A hyperstable (hs) variant of chicken egg-white lysozyme with enhanced thermal (delta Tm approximately +10.5 degrees C) and chemical (delta Cm for guanidine hydrochloride denaturation = +1.3 M) stabilities relative to wild-type (WT) was constructed by combining several individual stabilizing substitutions. The free energy difference between the native and denatured states of the hs variant is 3.1 (GdnHCl, 25 degrees C) to 4.0 (differential scanning calorimetry, 74 degrees C) kcal mol-1 greater than that of WT. The specific activity of the hs variant is 2.5-fold greater than that of WT. The choice of mutations came from diverse sources: (1) The I55L/S91T core construct with delta Tm = 3.3 degrees C from WT was available from the accompanying study (Shih P, Holland DR, Kirsch JF, 1995, Protein Sci 4:2050-2062). (2) The A31V mutation was suggested by the better atomic packing in the human lysozyme structure where the Ala 31 equivalent is Leu. (3) The H15L and R114H substitutions were selected on the basis of sequence comparisons with pheasant lysozymes that are more stable than the chicken enzyme. (4) The D101S variant was identified from a screen of mutants previously prepared in this laboratory. The effects of the individual mutations on stability are cumulative and nearly additive.     

Hydrophobic sliding: a possible mechanism for drug resistance in human immunodeficiency virus type 1 protease

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.     

The role played by environmental residues on sidechain torsional angles within homologous families of proteins: A new method of sidechain modeling

We investigated the conservation of sidechain conformation for each residue within a homologous family of proteins in the Protein Data Bank (PDB) and performed sidechain modeling using this information. The information was represented by the probability of conserved sidechain torsional angles obtained from many families of proteins, and these were calculated for a pair of residues at topologically equivalent positions as a result of structural alignment. Probabilities were obtained for a pair of same amino acids and for a pair of different amino acids. The correlation between environmental residues and the fluctuation of probability was examined for the pair of same amino acid residues, and the simple probability was calculated for the pair of different amino acids. From the results on the same amino acid pairs, 17 amino acids, except for Ala, Gly, and Pro, were divided into two types: those that were influenced and those that were not influenced by the environmental residues. From results on different amino acid pairs, a replacement between large residues, such as Trp, Phe, and Tyr, was performed assuming conservation of their torsional angles within a homologous family of proteins. We performed sidechain modeling for 11 known proteins from their native and modeled backbones, respectively. With the native backbones, the percentage of the χ₁ angle correct within 30° was found to be 67% and 80% for all and core residues, respectively. With the modeled backbones, the percentage of the correct χ₁ angle was found to be 60% and 72% for all and core residues, respectively. To estimate an upper limit on the accuracy for predicting sidechain conformations, we investigated the probability of conserved sidechain torsional angles for highly similar proteins having > 90% sequence identity and <2.5-Å X-ray resolution. In those proteins, 83% of the sidechain conformations were conserved for the χ₁ angle. Proteins 31:355–369, 1998. © 1998 Wiley-Liss, Inc.     

The primary structure of the Laburnum alpinum seed lectin

Complete proton and carbon sidechain assignments are reported for 22 lysine and 11 leucine residues in staphylococcal nuclease, an enzyme with 149 residues. These assignments are readily obtained in a direct manner from the correlations observed in the 3D HCCH-COSY and HCCH-TOCSY spectra and the known protein backbone assignments. These assignments open the way to detailed studies of the sidechain structure and dynamics at the active site, in the hydrophobic core and on the surface of the protein.     

Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils

Membrane-embedded protein domains frequently exist as α-helical bundles, as exemplified by photosynthetic reaction centers, bacteriorhodopsin, and cytochrome C oxidase. The sidechain packing between their transmembrane helices was investigated by a nearest-neighbor analysis which identified sets of interfacial residues for each analyzed helix–helix interface. For the left-handed helix–helix pairs, the interfacial residues almost exclusively occupy positions a, d, e, or g within a heptad motif (abcdefg) which is repeated two to three times for each interacting helical surface. The connectivity between the interfacial residues of adjacent helices conforms to the knobs-into-holes type of sidechain packing known from soluble coiled coils. These results demonstrate on a quantitative basis that the geometry of sidechain packing is similar for left-handed helix–helix pairs embedded in membranes and coiled coils of soluble proteins. The transmembrane helix–helix interfaces studied are somewhat less compact and regular as compared to soluble coiled coils and tolerate all hydrophobic amino acid types to similar degrees. The results are discussed with respect to previous experimental findings which demonstrate that specific interactions between transmembrane helices are important for membrane protein folding and/or oligomerization. Proteins 31:150–159, 1998. © 1998 Wiley-Liss, Inc.     

Conformational Dynamics Allows Sampling of an “Active-like” State by Oncogenic K-Ras-GDP

Mutations in K-Ras GTPase replacing Gly12 with either Asp or Val are common in cancer. These mutations decelerate intrinsic and catalyzed GTP hydrolysis, leading to accumulation of K-Ras-GTP in cells. Signaling cascades initiated by K-Ras-GTP promote cell proliferation, survival, and invasion. Despite functional differences between the most frequent G12D mutation and the most aggressive and chemotherapy resistant G12V mutation, their long-suspected distinct structural features remain elusive. Using NMR, X-ray structures, and computational methods, we found that oncogenic mutants of K-Ras4B, the predominant splice variant of K-Ras, exhibit distinct conformational dynamics when GDP-bound, visiting the “active-like” conformational state similar to the one observed in GTP-bound K-Ras. This behavior distinguishes G12V from wild type and G12D K-Ras4B-GDP. The likely reason is interactions between the aliphatic sidechain of V12 and the Switch II region of K-Ras4B^G12V-GDP, which are distinct in K-Ras4B^G12D-GDP. In the X-ray structures, crystal contacts reduce the dynamics of the sidechain at position 12 by stabilizing the Switch I region of the protein. This explains why structural differences between G12V and G12D K-Ras have yet not been reported. Together, our results suggest a previously unknown mechanism of K-Ras activation. This mechanism relies on conformational dynamics caused by specific oncogenic mutations in the GDP-bound state. Our findings also imply that the therapeutic strategies decreasing the level of K-Ras-GTP by interfering with nucleotide exchange or by expediting GTP hydrolysis may work differently in different oncogenic mutants.     

De novo design of the hydrophobic core of ubiquitin.

We have previously reported the development and evaluation of a computational program to assist in the design of hydrophobic cores of proteins. In an effort to investigate the role of core packing in protein structure, we have used this program, referred to as Repacking of Cores (ROC), to design several variants of the protein ubiquitin. Nine ubiquitin variants containing from three to eight hydrophobic core mutations were constructed, purified, and characterized in terms of their stability and their ability to adopt a uniquely folded native-like conformation. In general, designed ubiquitin variants are more stable than control variants in which the hydrophobic core was chosen randomly. However, in contrast to previous results with 434 cro, all designs are destabilized relative to the wild-type (WT) protein. This raises the possibility that beta-sheet structures have more stringent packing requirements than alpha-helical proteins. A more striking observation is that all variants, including random controls, adopt fairly well-defined conformations, regardless of their stability. This result supports conclusions from the cro studies that non-core residues contribute significantly to the conformational uniqueness of these proteins while core packing largely affects protein stability and has less impact on the nature or uniqueness of the fold. Concurrent with the above work, we used stability data on the nine ubiquitin variants to evaluate and improve the predictive ability of our core packing algorithm. Additional versions of the program were generated that differ in potential function parameters and sampling of side chain conformers. Reasonable correlations between experimental and predicted stabilities suggest the program will be useful in future studies to design variants with stabilities closer to that of the native protein. Taken together, the present study provides further clarification of the role of specific packing interactions in protein structure and stability, and demonstrates the benefit of using systematic computational methods to predict core packing arrangements for the design of proteins.     

Protein sequence optimization with a pairwise decomposable penalty for buried unsatisfied hydrogen bonds

In aqueous solution, polar groups make hydrogen bonds with water, and hence burial of such groups in the interior of a protein is unfavorable unless the loss of hydrogen bonds with water is compensated by formation of new ones with other protein groups. For this reason, buried “unsatisfied” polar groups making no hydrogen bonds are very rare in proteins. Efficiently representing the energetic cost of unsatisfied hydrogen bonds with a pairwise-decomposable energy term during protein design is challenging since whether or not a group is satisfied depends on all of its neighbors. Here we describe a method for assigning a pairwise-decomposable energy to sidechain rotamers such that following combinatorial sidechain packing, buried unsaturated polar atoms are penalized. The penalty can be any quadratic function of the number of unsatisfied polar groups, and can be computed very rapidly. We show that inclusion of this term in Rosetta sidechain packing calculations substantially reduces the number of buried unsatisfied polar groups.     

Effects of the mur1 mutation on xyloglucans produced by suspension-cultured Arabidopsis thaliana cells

Mutation of the Arabidopsis thaliana (L.) Heynh. gene MUR1, which encodes an isoform of GDP-D-mannose-4,6-dehydratase, affects the biosynthetic conversion of GDP-mannose to GDP-fucose. Cell walls in the aerial tissues of mur1 plants are almost devoid of alpha-L-fucosyl residues, which are partially replaced by closely related alpha-L-galactosyl residues. A line of suspension-cultured A. thaliana cells was generated from leaves of mur1 plants and the structure of the xyloglucan in the walls of these cells was structurally characterized. Xyloglucan fractions were prepared from the walls of both wild-type (WT) and mur1 cells by sequential extraction with a xyloglucan-specific endoglucanase (XEG) and aqueous KOH. Structural analysis of these fractions revealed that xyloglucan produced by cultured mur1 cells is similar, but not identical to that isolated from leaves of mur1 plants. As previously reported for mur1 leaves, the xyloglucan from cultured mur1 cells contains less than 5% of the fucose present in the xyloglucan from WT cells. Fucosylation of the xyloglucan is substantially restored when mur1 cells are grown in medium supplemented with L-fucose. Xyloglucan isolated from leaves contains more oligosaccharide subunits in which the central sidechain is terminated with a beta-D-galactosyl residue than does xyloglucan prepared from cultured cells. This was observed for both mur1 and WT plants, indicating that this correlation is independent of the mur1 mutation and that it is possible to distinguish changes due to genetic mutation from those due to the physiological state of the cells in culture. Suspension-cultured cells thus provide a convenient source of genetically altered cell wall material, facilitating the biochemical characterization of mutations that affect cell wall structure.     

Kinetic stabilization of Bacillus licheniformis alpha-amylase through introduction of hydrophobic residues at the surface

It is generally assumed that in proteins hydrophobic residues are not favorable at solvent-exposed sites, and that amino acid substitutions on the surface have little effect on protein thermostability. Contrary to these assumptions, we have identified hyperthermostable variants of Bacillus licheniformis alpha-amylase (BLA) that result from the incorporation of hydrophobic residues at the surface. Under highly destabilizing conditions, a variant combining five stabilizing mutations unfolds 32 times more slowly and at a temperature 13 degrees C higher than the wild-type. Crystal structure analysis at 1.7 A resolution suggests that stabilization is achieved through (a) extension of the concept of increased hydrophobic packing, usually applied to cavities, to surface indentations, (b) introduction of favorable aromatic-aromatic interactions on the surface, (c) specific stabilization of intrinsic metal binding sites, and (d) stabilization of a beta-sheet by introducing a residue with high beta-sheet forming propensity. All mutated residues are involved in forming complex, cooperative interaction networks that extend from the interior of the protein to its surface and which may therefore constitute "weak points" where BLA unfolding is initiated. This might explain the unexpectedly large effect induced by some of the substitutions on the kinetic stability of BLA. Our study shows that substantial protein stabilization can be achieved by stabilizing surface positions that participate in underlying cooperatively formed substructures. At such positions, even the apparently thermodynamically unfavorable introduction of hydrophobic residues should be explored.     

Molecular determinants of ATP-sensitive potassium channel MgATPase activity: diabetes risk variants and diazoxide sensitivity

ATP-sensitive K⁺ (K_ATP) channels play an important role in insulin secretion. K_ATP channels possess intrinsic MgATPase activity that is important in regulating channel activity in response to metabolic changes, although the precise structural determinants are not clearly understood. Furthermore, the sulfonylurea receptor 1 (SUR1) S1369A diabetes risk variant increases MgATPase activity, but the molecular mechanisms remain to be determined. Therefore, we hypothesized that residue–residue interactions between 1369 and 1372, predicted from in silico modelling, influence MgATPase activity, as well as sensitivity to the clinically used drug diazoxide that is known to increase MgATPase activity. We employed a point mutagenic approach with patch-clamp and direct biochemical assays to determine interaction between residues 1369 and 1372. Mutations in residues 1369 and 1372 predicted to decrease the residue interaction elicited a significant increase in MgATPase activity, whereas mutations predicted to possess similar residue interactions to wild-type (WT) channels elicited no alterations in MgATPase activity. In contrast, mutations that were predicted to increase residue interactions resulted in significant decreases in MgATPase activity. We also determined that a single S1369K substitution in SUR1 caused MgATPase activity and diazoxide pharmacological profiles to resemble those of channels containing the SUR2A subunit isoform. Our results provide evidence, at the single residue level, for a molecular mechanism that may underlie the association of the S1369A variant with type 2 diabetes. We also show a single amino acid difference can account for the markedly different diazoxide sensitivities between channels containing either the SUR1 or SUR2A subunit isoforms.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

In vivo and in vitro analyses of single-amino acid variants of the Salmonella enterica phosphotransacetylase enzyme provide insights into the function of its N-terminal domain

The function of the N-terminal domain ( approximately 350 residues) of the Pta (phosphotransacetylase) enzyme of Salmonella enterica is unclear. Results from in vivo genetic and in vitro studies suggest that the N-terminal domain of Pta is a sensor for NADH and pyruvate. We isolated 10 single-amino acid variants of Pta that, unlike the wild-type protein, supported growth of a strain of S. enterica devoid of Acs (acetyl-CoA synthetase; AMP-forming) activity on 10 mm acetate. All mutations were mapped within the N-terminal domain of the protein. Kinetic analyses of the wild type and three variant Pta proteins showed that two of the variant proteins were faster enzymes (k(cat) 2.5-3-fold > k(cat) Pta(WT). Results from sedimentation equilibrium experiments are consistent with Pta(WT) being a trimer. Pta variants formed more hexamer than the Pta(WT) protein. NADH inhibited Pta(WT) activity by inducing a conformational change detectable by limited trypsin proteolysis; NADH did not inhibit variant protein Pta(R252H). Pyruvate stimulated Pta(WT) activity, and its effect was potentiated in the variants, being most pronounced on Pta(R252H).     

Solution structure of switch Arc,a mutant with 3(10) helices replacing a wild-type beta-ribbon

Adjacent N11L and L12N mutations in the antiparallel beta-ribbon of Arc repressor result in dramatic changes in local structure in which each beta-strand is replaced by a right-handed helix. The full solution structure of this "switch" Arc mutant shows that irregular 3(10) helices compose the new secondary structure. This structural metamorphosis conserves the number of main-chain and side-chain to main-chain hydrogen bonds and the number of fully buried core residues. Apart from a slight widening of the interhelical angle between alpha-helices A and B and changes in side-chain conformation of a few core residues in Arc, no large-scale structural adjustments in the remainder of the protein are necessary to accommodate the ribbon-to-helix change. Nevertheless, some changes in hydrogen-exchange rates are observed, even in regions that have very similar structures in the two proteins. The surface of switch Arc is packed poorly compared to wild-type, leading to approximately 1000A(2) of additional solvent-accessible surface area, and the N termini of the 3(10) helices make unfavorable head-to-head electrostatic interactions. These structural features account for the positive m value and salt dependence of the ribbon-to-helix transition in Arc-N11L, a variant that can adopt either the mutant or wild-type structures. The tertiary fold is capped in different ways in switch and wild-type Arc, showing how stepwise evolutionary transformations can arise through small changes in amino acid sequence.