High-resolution structure prediction of a circular permutation loop

Methods for rapid and reliable design and structure prediction of linker loops would facilitate a variety of protein engineering applications. Circular permutation, in which the existing termini of a protein are linked by the polypeptide chain and new termini are created, is one such application that has been employed for decreasing proteolytic susceptibility and other functional purposes. The length and sequence of the linker can impact the expression level, solubility, structure and function of the permuted variants. Hence it is desirable to achieve atomic‐level accuracy in linker design. Here, we describe the use of RosettaRemodel for design and structure prediction of circular permutation linkers on a model protein. A crystal structure of one of the permuted variants confirmed the accuracy of the computational prediction, where the all‐atom rmsd of the linker region was 0.89 Å between the model and the crystal structure. This result suggests that RosettaRemodel may be generally useful for the design and structure prediction of protein loop regions for circular permutations or other structure‐function manipulations.     

The alkaline transition of blue copper proteins, Cucumis sativus plastocyanin and Pseudomonas aeruginosa azurin

Autoreduction of Cucumis sativus plastocyanin and Pseudomonas aeruginosa azurin took place at alkaline pHs, having been accompanied by the decrease in the intensities of the charge transfer band, Cys-S- (pi)-->Cu(II) at 597 and 626 nm, and the Cu(II)-EPR signals with small AII values of 6.5 x 10(-3) and 5.3 x 10(-3) cm(-1) for plastocyanin and azurin, respectively. Further, an extra Cu(II)-EPR signal with a large AII value of 21 x 10(-3) cm(-1) also reversibly emerged with increasing pH. Plastocyanin and azurin are in an equilibrium of the three forms at alkaline pHs.     

Alteration of the disulfide-coupled folding pathway of BPTI by circular permutation

The kinetics of disulfide-coupled folding and unfolding of four circularly permuted forms of bovine pancreatic trypsin inhibitor (BPTI) were studied and compared with previously published results for both wild-type BPTI and a cyclized form. Each of the permuted proteins was found to be less stable than either the wild-type or circular proteins, by 3-8 kcal/mole. These stability differences were used to estimate effective concentrations of the chain termini in the native proteins, which were 1 mM for the wild-type protein and 2.5 to 4000 M for the permuted forms. The circular permutations increased the rates of unfolding and caused a variety of effects on the kinetics of refolding. For two of the proteins, the rates of a direct disulfide-formation pathway were dramatically increased, making this process as fast or faster than the competing disulfide rearrangement mechanism that predominates in the folding of the wild-type protein. These two permutations break the covalent connectivity among the beta-strands of the native protein, and removal of these constraints appears to facilitate direct formation and reduction of nearby disulfides that are buried in the folded structure. The effects on folding kinetics and mechanism do not appear to be correlated with relative contact order, a measure of overall topological complexity. These observations are consistent with the results of other recent experimental and computational studies suggesting that circular permutation may generally influence folding mechanisms by favoring or disfavoring specific interactions that promote alternative pathways, rather than through effects on the overall topology of the native protein.     

Construction of hybrid gene libraries involving the circular permutation of DNA

The method of incremental truncation for the creation of hybrid enzymes (ITCHY) allows the creation of comprehensive fusion libraries between 5 and 3 fragments of two genes in a manner that is independent of DNA sequence homology. A methodology is presented for the creation of ITCHY libraries called circularly permuted ITCHY (CP-ITCHY) that allows the creation of ITCHY libraries in a manner that does not require extensive time point sampling. In addition, CP-ITCHY requires only a single vector and productively biases the library towards those fusions that are approximately the same size as the original genes. In the model system of creating fusions between fragments of the Escherichia coli and human glycinamide ribonucleotide transformylase genes, the CP-ITCHY libraries are shown to contain a diverse set of active fusions including those in regions of low-homology. In addition, a high percentage of active fusions were temperature-sensitive as they complemented an auxotrophic strain of Escherichia coli at 22 °C but not at 37 °C.     

Thermodynamic and structural consequences of flexible loop deletion by circular permutation in the streptavidin-biotin system.

A circularly permuted streptavidin (CP51/46) has been designed to remove the flexible polypeptide loop that undergoes an open to closed conformational change when biotin is bound. The original termini have been joined by a tetrapeptide linker, and four loop residues have been removed, resulting in the creation of new N- and C-termini. Isothermal titration calorimetric studies show that the association constant has been reduced approximately six orders of magnitude below that of wild-type streptavidin to 10(7) M(-1). The deltaH degrees of biotin association for CP51/46 is reduced by 11.1 kcal/mol. Crystal structures of CP51/46 and its biotin complex show no significant alterations in the binding site upon removal of the loop. A hydrogen bond between Ser45 and Ser52 found in the absence of biotin is broken in the closed conformation as the side-chain hydroxyl of Ser45 moves to hydrogen bond to a ureido nitrogen of biotin. This is true in both the wild-type and CP51/46 forms of the protein, and the hydrogen bonding interaction might thus help nucleate closure of the loop. The reduced entropic cost of binding biotin to CP51/46 is consistent with the removal of this loop and a reduction in entropic costs associated with loop closure and immobilization. The reduced enthalpic contribution to the free energy of binding is not readily explainable in terms of the molecular structure, as the binding contacts are nearly entirely conserved, and only small differences in solvent accessible surfaces are observed relative to wild-type streptavidin.     

Vibrational and electronic circular dichroism monitoring of copper(II) coordination with a chiral ligand

Novel copper(II) coordination compounds with chiral macrocyclic imine ligands derived from R-/S-camphor were asymmetrically synthesized and characterized with the aid of chiroptical spectroscopies. Crystal structures of both enantiomers were determined by single crystal X-ray diffraction analysis. Circular dichroism (CD) spectra were analyzed using a simplified exciton model as well as quantum chemical computations. The absolute configuration of the copper(II) coordination compounds determined from CD was found consistent with the crystal data. The copper(II) complexes were further investigated by vibrational CD (VCD) measurement combined with density functional theory calculation. The complex formation was evidenced by spectral shifts of the characteristic bands in the CD and VCD spectra.     

Anti-cooperativity in the two electron oxidation of the S118C disulfide dimer of azurin

We have constructed a disulfide dimer of S118C azurin, in which two copper centers are coupled through a relatively short covalent pathway, and studied its electron transfer properties. The dimer exhibits intriguing mechanistic properties. Due to the strain in the molecule, caused by the limited accessibility of Cys118, anti-cooperativity occurs in the two step oxidation of the dimer with a difference in redox potential between the two half reactions of 33 mV. Upon oxidation, the dimer favours the semi-reduced over the fully oxidized state, as the Cu(I) site in the semi-reduced dimer is able to stabilize the strained dimer complex. The internal electron transfer is surprisingly slow, which could be partially due to an increase in reorganization energy.     

Importance of polarization effect in the study of metalloproteins: Application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein‐specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild‐type azurin (ΔE_cal) were calculated and compared with experimental values. The ΔE_cal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein‐specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations. Proteins 2014; 82:2209–2219. © 2014 Wiley Periodicals, Inc.     

Effect of circular permutations on transient partial unfolding in proteins

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three‐dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N‐termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native‐state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild‐type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild‐type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild‐type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.     

Spectral magnitude effects on the analyses of secondary structure from circular dichroism spectroscopic data

The effects of spectral magnitude on the calculated secondary structures derived from circular dichroism (CD) spectra were examined for a number of the most commonly used algorithms and reference databases. Proteins with different secondary structures, ranging from mostly helical to mostly beta-sheet, but which were not components of existing reference databases, were used as test systems. These proteins had known crystal structures, so it was possible to ascertain the effects of magnitude on both the accuracy of determining the secondary structure and the goodness-of-fit of the calculated structures to the experimental data. It was found that most algorithms are highly sensitive to spectral magnitude, and that the goodness-of-fit parameter may be a useful tool in assessing the correct scaling of the data. This means that parameters that affect magnitude, including calibration of the instrument, the spectral cell pathlength, and the protein concentration, must be accurately determined to obtain correct secondary structural analyses of proteins from CD data using empirical methods.     

Changes in secondary structure follow the dissociation of human insulin hexamers: a circular dichroism study

Vacuum UV circular dichroism spectra measured down to 178 nm for hexameric 2-zinc human insulin, zinc-free human insulin, and the two engineered and biologically active monomeric mutants, [B/S9D] and [B/S9D,T27E] human insulin, show significant differences. The secondary structure analysis of the 2-zinc human insulin (T6) in neutral solution was determined: 57% helix, 1% beta-strand, 18% turn, and 24% random coil. This is very close to the corresponding crystal structure showing that the solution and solid structures are similar. The secondary structure of the monomer shows a 10-15% increase in antiparallel beta-structure and a corresponding reduction in random coil structure. These structural changes are consistent with an independent analysis of the corresponding difference spectra. The advantage of secondary structure analyses of difference spectra is that the contribution of odd spectral features stemming mainly from side chain chromophores is minimized and the sensitivity of the analyses improved. Analysis of the CD spectra of T6 2-zinc, zinc-free human insulin and monomeric mutant insulin by singular value decomposition indicates that the secondary structure changes following the dissociation of hexamers into dimers and monomers are two-state processes.     

The metal site of Pseudomonas aeruginosa azurin,revealed by a crystal structure determination of the co(II) derivative and co-EPR spectroscopy

The crystal structure of cobalt-substituted azurin from Pseudomonas aeruginosa has been determined to final crystallographic R value of 0.175 at 1.9 Å resolution. There are four molecules in the asymmetric unit in the structure, and these four molecules are packed as a dimer of dimers. The dimer packing is very similar to that of the wild-type Pseudomonas aeruginosa azurin dimer. Replacement of the native copper by the cobalt ion has only small effects on the metal binding site presumably because of the existence of an extensive network of hydrogen bonds in its immediate neighborhood. Some differences are obvious, however. In wild-type azurin the copper atom occupies a distorted trigonal bipyramidal site, while cobalt similar to zinc and nickel occupy a distorted tetrahedral site, in which the distance to the Met121,S^δ atom is increased to 3.3–3.5 Å and the distance to the carbonyl oxygen of Gly45 has decreased to 2.1–2.4 Å. The X-band EPR spectrum of the high-spin Co(II) in azurin is well resolved (apparent g values g_x′ = 5.23; g_y′ = 3.83; g_z′ = 1.995, and hyperfine splittings A_x′ = 31; A_y′ = 20–30; A_z′ = 53 G) and indicates that the ligand field is close to axial. Proteins 27:385–394, 1997. © 1997 Wiley-Liss, Inc.     

Structural, dynamical and functional aspects of the inner motions in the blue copper protein azurin

Molecular dynamics was applied to dissect out the internal motions of azurin, a copper protein performing electron transfer. Simulations of 16.5 ns were analyzed in search of coordinated displacements of amino acid residues that are important for the protein function. A region with high conformational instability was found in the 'southern' end of the molecule, far away from the copper site and the binding sites for the redox partners of azurin. By excluding the 'southern' region from the subsequent analysis, correlated motions were identified in the hydrophobic patch that surrounds the protein active site. The simulation results are in excellent agreement with recent NMR data on azurin in solution [A. V. Zhuravleva, D. M. Korzhnev, E. Kupce, A. S. Arseniev, M. Billeter, V. Y. Orekhov, Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures, J. Mol. Biol. 342 (2004) 1599-1611] and suggest a rationale for cooperative displacements of protein residues that are thought to be critical for the electron transfer process. A number of other structural and dynamic features of azurin are discussed in the context of the blue copper protein family and an explanation is proposed to account for the variability/conservation of some regions in the cupredoxins.     

The structural role of the copper-coordinating and surface-exposed histidine residue in the blue copper protein azurin

Copper K-edge extended X-ray absorption fine structure (EXAFS) spectroscopy and (15)N NMR relaxation studies were performed on samples of a variant azurin in which the surface-exposed histidine ligand of the copper atom (His117) has been replaced by glycine. The experiments were performed to probe the structure of the active site and the protein dynamics. The cavity in the protein structure created by the His-->Gly replacement could be filled by external ligands, which can either restore the spectroscopic properties of the original type-1 copper site or create a new type-2 copper site. The binding of external ligands occurs only when the copper atom is in its oxidised state. In the reduced form, the binding is abolished. From the EXAFS experiments, it is concluded that for the oxidised type-1 copper sites the protein plus external ligand (L) provide an NSS*L donor set deriving from His46, Cys112, Met121 and the external ligand. The type-2 copper site features an S(N/O)(3) donor set in which the S-donor derives from Cys112, one N-donor from His46 and the remaining two N or O donors from one or more external ligands. Upon reduction of the type-1 as well as the type-2 site, the external ligand drops out of the copper site and the coordination reduces to 3-fold with an SS*N donor set deriving from His46, Cys112 and Met121. The Cu-S(delta)(Met) distance is reduced from about 3.2 to 2.3 A. Analysis of the NMR data shows that the hydrophobic patch around His117 has gained fluxionality when compared to wild-type azurin, which may explain why the His117Gly variant is able to accommodate a variety of external ligands of different sizes and with different chelating properties. On the other hand, the structure and dynamics of the beta-sandwich, which comprises the main body of the protein, is only slightly affected by the mutation. The unusually high reduction potential of the His117Gly azurin is discussed in light of the present results.     

The circular dichroism of tumor necrosis factor-α: Measurement into the vacuum UV and analysis for secondary structure

The circular dichroism (CD) spectrum of tumor necrosis factor-α has been measured into the vacuum UV to 168 nm. Analysis of the CD for secondary structure is in good agreement with X-ray diffraction results, but the analysis is somewhat unstable. Adding the CD of this protein together with its X-ray determined secondary structure to the basis set should improve subsequent analyses of CD spectra for other all-β proteins.     

Complete sequential 1H and 15N nuclear magnetic resonance assignments and solution secondary structure of the blue copper protein azurin from Pseudomonas aeruginosa.

Complete sequential 1H and 15N resonance assignments for the reduced Cu(I) form of the blue copper protein azurin (M(r) 14,000, 128 residues) from Pseudomonas aeruginosa have been obtained at pH 5.5 and 40 degrees C by using homo- and heteronuclear two-dimensional (2D) and three-dimensional (3D) nuclear magnetic resonance spectroscopic experiments. Combined analysis of a 3D homonuclear 1H Hartmann-Hahn nuclear Overhauser (3D 1H HOHAHA-NOESY) spectrum and a 3D heteronuclear 1H nuclear Overhauser 1H[15N] single-quantum coherence (3D 1H[15N] NOESY-HSQC) spectrum proved especially useful. The latter spectrum was recorded without irradiation of the water signal and provided for differential main chain amide (NH) exchange rates. NMR data were used to determine the secondary structure of azurin in solution. Comparison with the secondary structure of azurin obtained from X-ray analysis shows a virtually complete resemblance; the two beta-sheets and a 3(10)-alpha-3(10) helix are preserved at 40 degrees C, and most loops contain well-defined turns. Special findings are the unexpectedly slow exchange of the Asn-47 and Phe-114 NH's and the observation of His-46 and His-117 N epsilon 2H resonances. The implications of these observations for the assignment of azurin resonance Raman spectra, the rigidity of the blue copper site, and the electron transfer mechanism of azurin are discussed.     

Effect of lysine modification on the secondary structure of ovalbumin

In order to determine the effect of chemical modification of the -amino groups on the secondary structure of ovalbumin, we prepared six acetylated (17, 36, 54, 70, 82, and 98%) and four succinylated derivatives (25, 50, 72, and 97%) of the protein. Native ovalbumin and the acylated derivatives were homogeneous as revealed by the electrophoretic pattern. The UV-absorption and fluorescence spectra changed progressively with the extent of modification. However, circular dichroic (CD) studies indicated that acylation of 15 of the 20 lysine residues had little effect on the secondary structure of ovalbumin. Acylation of the remaining five lysine residues resulted in a fairly severe change in the secondary structure. The -helical content decreased from about 31% in the native state to 16.5% in the 97% succinylated ovalbumin and to 21.5% in the 98% acetylated derivative. A comparison of these data with the spectral and hydrodynamic data of Qasim and Salahuddin (1978) suggested that the secondary structure of ovalbumin is more resistant to acylation than is the tertiary structure and, thus, the tertiary and the secondary structures are, to some extent, mutually independent. Raising thepH to 11.2 did not alter the secondary structure of ovalbumin and increasing the ionic strength by more than 20-fold did not reverse the loss of helical structure in 97% succinylated protein. These two observations suggest that the change in secondary structure upon maximal acylation may not only involve electrostatic effects, but also certain other factors, such as steric hindrance due to the entering bulky groups.     

Structure of heavy and light chain subunits of type A botulinum neurotoxin analyzed by circular dichroism and fluorescence measurements

Summary The secondary and tertiary structural features of botulinum neurotoxin (NT) serotype A, a dichain protein (M_r 145 000), and its two subunits, the heavy (H) and light (L) chains (M_r 97 000 and 53 000, respectively) were examined using circular dichroism and fluorescence spectorscopy. Nearly 70% of the amino acid residues in each of the three polypeptide preparations were found in ordered structure (sum of helix, sheet and turns). Also, the helix, sheet, turns and random coil contents of the dichain NT were nearly equal to the weighted mean of each of these secondary structure parameters of the L and H chains; e.g., sum of helix of L chain (22%) and H chain (18.7%), as weighted mean, 19.8% was similar to that of NT (20%). These agreements suggested that the secondary structures of the subunits of the dichain NT do not significantly change when they are separated as isolated L and H chains. Fluorescence emission maximum of L chain, 4 nm less (blue shift) than that of H chain, suggested relatively more hydrophobic environment of fluorescent tryptophan residue(s) of L chain. Tryptophan fluorescence quantum yields of L chain, H chain and the NT, 0.072, 0.174 and 0.197, respectively, suggested that a) an alteration in the micro-environment of the tryptophan residues was possibly caused by interactions of L and H chain subunits of the NT and b) quantum yields for L and H chains were altered when they are together as subunits of the NT. Possible implications of structural features of the L and H chains, their interactions and the molecular mechanism of action of botulinum NT are assessed.     

Design and initial characterization of a circular permuted variant of the potent HIV-inactivating protein cyanovirin-N.

A circular permuted variant of the potent human immunodeficiency virus (HIV)-inactivating protein cyanovirin-N (CV-N) was constructed. New N- and C-termini were introduced into an exposed helical loop, and the original termini were linked using residues of the original loop. Since the three-dimensional structure of wild-type cyanovirin-N is a pseudodimer, the mutant essentially exhibits a swap between the two pseudo-symmetrically related halves. The expressed protein, which accumulates in the insoluble fraction, was purified, and conditions for in vitro refolding were established. During refolding, a transient dimeric species is also formed that converts to a monomer. Similar to the wild-type CV-N, the monomeric circular permuted protein exhibits reversible thermal unfolding and urea denaturation. The mutant is moderately less stable than the wild-type protein, but it displays significantly reduced anti-HIV activity. Using nuclear magnetic resonance spectroscopy, we demonstrate that this circular permuted monomeric molecule adopts the same fold as the wild-type protein. Characterization of these two architecturally very similar molecules allows us to embark, for the first time, on a structure guided focused mutational study, aimed at delineating crucial features for the extraordinary difference in the activity of these molecules.