A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

Significance of the Disulfide Bridge in the Structure and Stability of Metalloprotein Azurin

Metalloproteins make up a class of proteins that incorporate metal ions into their structures, enabling them to perform essential functions in biological systems, such as catalysis and electron transport. Azurin is one such metalloprotein with copper cofactor, having a β-barrel structure with exceptional thermal stability. The copper metal ion is coordinated at one end of the β-barrel structure, and there is a disulfide bond at the opposite end. In this study, we explore the effect of this disulfide bond in the high thermal stability of azurin by analyzing both the native S-S bonded and S-S nonbonded (S-S open) forms using temperature replica exchange molecular dynamics (REMD). Similar to experimental observations, we find a 35 K decrease in denaturation temperature for S-S open azurin compared to that of the native holo form (420 K). As observed in the case of native holo azurin, the unfolding process of the S-S open form also started with disruptions of the α-helix. The free energy surfaces of the unfolding process revealed that the denaturation event of the S-S open form progresses through different sets of conformational ensembles. Subsequently, we compared the stabilities of individual β-sheet strands of both the S-S bonded and the S-S nonbonded forms of azurin. Further, we examined the contacts between individual residues for the central structures from the free energy surfaces of the S-S nonbonded form. The microscopic origin of the lowering in the denaturation temperature is further supplemented by thermodynamic analysis.

Role of Metal Cofactor in Enhanced Thermal Stability of Azurin

Metal cofactors are critical centers for different biochemical processes of metalloproteins, and often, this metal coordination renders additional structural stability. In this study, we explore the additional stability conferred by the copper ion on azurin by analyzing both the apo and holo forms using temperature replica exchange molecular dynamics (REMD) data. We find a 14 K decrease in denaturation temperature for apo (406 K) azurin relative to that of holo (420 K), indicating a copper ion-induced additional thermal stability for holo azurin. The unfolding of apo azurin begins with the melting of α-helix and β-sheet V, similar to that of holo form. β-Sheets IV, VII, and VIII are comparatively more stable than other β-strands and melt at higher temperatures. Similar to holo azurin, the strong hydrophobic interactions among the apolar residues in the protein core is the key factor that renders high stability to apo protein as well. We construct free energy surfaces at different temperatures to capture the major conformations along the unfolding basins of the protein. Using contact maps from different basins we show the changes in the interaction between different residues along the unfolding pathway. Furthermore, we compare the Cα root-mean-square fluctuations (Cα-RMSF) and B-factor of all residues of apo and holo forms to understand the flexibility of different regions. The concerted displacement of α-helix and β-sheets V and VI from the protein core is another distinction we observe for apo compared to the holo form, where β-sheet VI was relatively stable.

Molecular Insight into the High Thermal Stability of Metalloprotein Azurin

We investigate the events characterizing the steps of the unfolding pathway of blue copper metalloprotein azurin using replica exchange molecular dynamics (REMD). Our studies show that the unfolding of azurin begins with the melting of α-helix and β-sheets II and V. This is followed by the melting of other β-sheets and the exposure of hydrophobic protein core to the solvent, resulting in disruptions of its tertiary structure. Free energy surfaces constructed at different temperatures portray different basins that signify the stability of different melted structures in the unfolding process. The contact maps at different temperatures reveal that the strong hydrophobic interaction within the core of the protein is the vital force that renders high stability to this protein. Analysis of the individual β-sheets by looking into their amino acid sequence shows that β-sheets with charged side chains on the surface melt fast compared to others. The β-barrel of azurin is able to dynamically rearrange, and it helps the protein to preserve its hydrophobic core, holding back the native topology from melting fast. B-factor analysis shows that residues of β-sheets III, IV, and VII deviate less from their initial structure at the transition temperature.

Octarellin VI: using rosetta to design a putative artificial (β/α)8 protein

The computational protein design protocol Rosetta has been applied successfully to a wide variety of protein engineering problems. Here the aim was to test its ability to design de novo a protein adopting the TIM-barrel fold, whose formation requires about twice as many residues as in the largest proteins successfully designed de novo to date. The designed protein, Octarellin VI, contains 216 residues. Its amino acid composition is similar to that of natural TIM-barrel proteins. When produced and purified, it showed a far-UV circular dichroism spectrum characteristic of folded proteins, with α-helical and β-sheet secondary structure. Its stable tertiary structure was confirmed by both tryptophan fluorescence and circular dichroism in the near UV. It proved heat stable up to 70°C. Dynamic light scattering experiments revealed a unique population of particles averaging 4 nm in diameter, in good agreement with our model. Although these data suggest the successful creation of an artificial α/β protein of more than 200 amino acids, Octarellin VI shows an apparent noncooperative chemical unfolding and low solubility.

Molecular dynamics simulation of phosphorylated KID post-translational modification

Background

Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood.

Methodology

To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding.

Conclusions

Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix α_B, folding of helix α_A, completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices α_A and α_B are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Φ-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification.

Identifying critical residues in protein folding: Insights from phi-value and P(fold) analysis

We apply a simulational proxy of the phi-value analysis and perform extensive mutagenesis experiments to identify the nucleating residues in the folding "reactions" of two small lattice Go polymers with different native geometries. Our findings show that for the more complex native fold (i.e., the one that is rich in nonlocal, long-range bonds), mutation of the residues that form the folding nucleus leads to a considerably larger increase in the folding time than the corresponding mutations in the geometry that is predominantly local. These results are compared to data obtained from an accurate analysis based on the reaction coordinate folding probability P(fold) and on structural clustering methods. Our study reveals a complex picture of the transition state ensemble. For both protein models, the transition state ensemble is rather heterogeneous and splits up into structurally different populations. For the more complex geometry the identified subpopulations are actually structurally disjoint. For the less complex native geometry we found a broad transition state with microscopic heterogeneity. These findings suggest that the existence of multiple transition state structures may be linked to the geometric complexity of the native fold. For both geometries, the identification of the folding nucleus via the P(fold) analysis agrees with the identification of the folding nucleus carried out with the phi-value analysis. For the most complex geometry, however, the applied methodologies give more consistent results than for the more local geometry. The study of the transition state structure reveals that the nucleus residues are not necessarily fully native in the transition state. Indeed, it is only for the more complex geometry that two of the five critical residues show a considerably high probability of having all its native bonds formed in the transition state. Therefore, one concludes that, in general, the phi-value correlates with the acceleration/deceleration of folding induced by mutation, rather than with the degree of nativeness of the transition state, and that the "traditional" interpretation of phi-values may provide a more realistic picture of the structure of the transition state only for more complex native geometries.

Engineering the thermostability of a TIM-barrel enzyme by rational family shuffling

A possible approach to generate enzymes with an engineered temperature optimum is to create chimeras of homologous enzymes with different temperature optima. We tested this approach using two family-10 xylanases from Thermotoga maritima: the thermophilic xylanase A catalytic domain (TmxAcat, T(opt)=68 degrees C), and the hyperthermophilic xylanase B (TmxB, T(opt)=102 degrees C). Twenty-one different chimeric constructs were created by mimicking family shuffling in a rational manner. The measured temperature optima of the 16 enzymatically active chimeras do not monotonically increase with the percentage of residues coming from TmxB. Only four chimeras had a higher temperature optimum than TmxAcat, the most stable variant (T(opt)=80 degrees C) being the one in which both terminal segments came from TmxB. Further analysis suggests that the interaction between the N- and C-terminal segments has a disproportionately high contribution to the overall thermostability. The results may be generalizable to other enzymes where the N- and C-termini are in contact.

Identification and Characterization of Key Substructures Involved in the Early Folding Events of a (β/α)8-barrel Protein as Studied by Experimental and Computational Methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.