Solution structure of the Cys74 to Ala74 mutant of the recombinant catalytic domain of Zoocin A

Zoocin A is a Zn‐metallopeptidase secreted by Streptococcus zooepidemicus strain 4881. Its catalytic domain is responsible for cleaving the D‐alanyl‐L‐alanine peptide bond in streptococcal peptidoglycan. The solution NMR structure of the Cys74 to Ala74 mutant of the recombinant catalytic domain (rCAT C74A) has been determined. With a previous structure determination for the recombinant target recognition domain (rTRD), this completes the 3D structure of zoocin A. While the structure of rCAT C74A resembles those of the catalytic domains of lysostaphin and LytM, the substrate binding groove is wider and no tyrosine residue was observed in the active site. Proteins 2016; 85:177–181. © 2016 Wiley Periodicals, Inc.     

Proposed docking interface between peptidoglycan and the target recognition domain of zoocin A

A docking model is proposed for the target recognition domain of the lytic exoenzyme zoocin A with the peptidoglycan on the outer cell surface of sensitive bacterial strains. Solubilized fragments from such peptidoglycans perturb specific backbone and side chain amide resonances in the recombinant form of the domain designated rTRD as detected in two-dimensional ¹H–¹⁵N correlation NMR spectra. The affected residues comprise a shallow surface cleft on the protein surface near W115, N53, N117, and Q105 among others, which interacts with the peptide portion of the peptidoglycan. Calculations with AutoDock Vina provide models of the docking interface. There is approximate homology between the rTDR-peptidoglycan docking site and the antigen binding site of Fab antibodies with the immunoglobin fold. EDTA was also found to bind to rTRD, but at a site distinct from the proposed peptidoglycan docking site.     

The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

The metal binding site of zoocin A

Direct metal analysis of the bacteriolytic exoenzyme zoocin A failed to unequivocally identify a putative metal cofactor; hence, indirect experiments utilizing NMR were undertaken to settle this question. Cd(2+) as a surrogate metal ion was reconstituted into EDTA-treated, metal-free recombinant zoocin, and (113)Cd-NMR was employed to explore binding in the protein for this ion. The Cd-substituted enzyme was found to have 80-85% of native streptococcolytic activity. A major (113)Cd resonance at 113.6 ppm was observed which with time split into resonances at 113.6 and 107.2 ppm. A minor (113)Cd resonance at 87.3 ppm was observed which increased in intensity with time. These Cd chemical shifts are indicative of two N atoms and two O atoms ligating directly to the metal site.On the basis of conserved amino acid residues in a homologous protein of known structure, LytM, the ligands in zoocin are tentatively assigned to H45, D49, H133, and some combination of water or buffer ions as the fourth oxygen donor in zoocin A. Comparison of the combined intensities for (113)Cd-substituted zoocin with a known quantity of another Cd-substituted protein gave Cd binding as approximately stoichiometric (1.2 +/- 0.2) with protein. Additional metal-removal and reconstitution experiments on the recombinant catalytic domain of zoocin implicate Zn(2+) as the metal cofactor. Therefore, the evidence supports zoocin as a single Zn(2+) ion binding metalloenzyme.     

NMR evidence for independent domain structures in zoocin A, an antibacterial exoenzyme

NMR was used to obtain spectroscopic evidence supporting a two domain model for zoocin A in which an N-terminal catalytic domain is linked by a threonine-proline rich linker to a target recognition domain responsible for recognizing the cell wall of bacteria susceptible to the bacteriolytic action of the enzyme. When cloned and separately expressed, each domain retains the folding found in the whole enzyme. Additionally, spectroscopy suggests that the target recognition domain has a conformation typical of a soluble globular protein, while the catalytic domain aggregates at low millimolar concentrations.     

Temperature‐induced partially unfolded state of hUBF HMG Box‐5: Conformational and dynamic investigations of the Box‐5 thermal intermediate ensemble

Human upstream binding factor (hUBF) HMG Box‐5 is a highly conserved protein domain, containing 84 amino acids and belonging to the family of the nonspecific DNA‐binding HMG boxes. Its native structure adopts a twisted L shape, which consists of three α‐helices and two hydrophobic cores: the major wing and the minor wing. In this article, we report a reversible three‐state thermal unfolding equilibrium of hUBF HMG Box‐5, which is investigated by differential scanning calorimetry (DSC), circular dichroism spectroscopy, fluorescence spectroscopy, and NMR spectroscopy. DSC data show that Box‐5 unfolds reversibly in two separate stages. Spectroscopic analyses suggest that different structural elements exhibit noncooperative transitions during the unfolding process and that the major form of the Box‐5 thermal intermediate ensemble at 55°C shows partially unfolded characteristics. Compared with previous thermal stability studies of other boxes, it appears that Box‐5 possesses a more stable major wing and two well separated subdomains. NMR chemical shift index and sequential ¹H^N_i‐¹H^N_i+1 NOE analyses indicate that helices 1 and 2 are native‐like in the thermal intermediate ensemble, while helix 3 is partially unfolded. Detailed NMR relaxation dynamics are compared between the native state and the intermediate ensemble. Our results implicate a fluid helix‐turn‐helix folding model of Box‐5, where helices 1 and 2 potentially form the helix 1‐turn‐helix 2 motif in the intermediate, while helix 3 is consolidated only as two hydrophobic cores form to stabilize the native structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Noncollagenous region of the streptococcal collagen‐like protein is a trimerization domain that supports refolding of adjacent homologous and heterologous collagenous domains

Proper folding of the (Gly‐Xaa‐Yaa)_n sequence of animal collagens requires adjacent N‐ or C‐terminal noncollagenous trimerization domains which often contain coiled‐coil or beta sheet structure. Collagen‐like proteins have been found recently in a number of bacteria, but little is known about their folding mechanism. The Scl2 collagen‐like protein from Streptococcus pyogenes has an N‐terminal globular domain, designated V_sp, adjacent to its triple‐helix domain. The V_sp domain is required for proper refolding of the Scl2 protein in vitro. Here, recombinant V_sp domain alone is shown to form trimers with a significant α‐helix content and to have a thermal stability of T_m = 45°C. Examination of a new construct shows that the V_sp domain facilitates efficient in vitro refolding only when it is located N‐terminal to the triple‐helix domain but not when C‐terminal to the triple‐helix domain. Fusion of the V_sp domain N‐terminal to a heterologous (Gly‐Xaa‐Yaa)_n sequence from Clostridium perfringens led to correct folding and refolding of this triple‐helix, which was unable to fold into a triple‐helical, soluble protein on its own. These results suggest that placement of a functional trimerization module adjacent to a heterologous Gly‐Xaa‐Yaa repeating sequence can lead to proper folding in some cases but also shows specificity in the relative location of the trimerization and triple‐helix domains. This information about their modular nature can be used in the production of novel types of bacterial collagen for biomaterial applications.     

Crystal structure of BinB: A receptor binding component of the binary toxin from Lysinibacillus sphaericus

The binary toxin (Bin), produced by Lysinibacillus sphaericus, is composed of BinA (42 kDa) and BinB (51 kDa) proteins, which are both required for full toxicity against Culex and Anopheles mosquito larvae. Specificity of Bin toxin is determined by the binding of BinB component to a receptor present on the midgut epithelial membranes, while BinA is proposed to be a toxic component. Here, we determined the first crystal structure of the active form of BinB at a resolution of 1.75 Å. BinB possesses two distinct structural domains in its N‐ and C‐termini. The globular N‐terminal domain has a β‐trefoil scaffold which is a highly conserved architecture of some sugar binding proteins or lectins, suggesting a role of this domain in receptor‐binding. The BinB β‐rich C‐terminal domain shares similar three‐dimensional folding with aerolysin type β‐pore forming toxins, despite a low sequence identity. The BinB structure, therefore, is a new member of the aerolysin‐like toxin family, with probably similarities in the cytolytic mechanism that takes place via pore formation. Proteins 2014; 82:2703–2712. © 2014 Wiley Periodicals, Inc.     

Chemical synthesis and evaluation of a backbone‐cyclized minimized 2‐helix Z‐domain

The Z‐molecule is a small, engineered IgG‐binding affinity protein derived from the immunoglobulin‐binding domain B of Staphylococcus aureus protein A. The Z‐domain consists of 58 amino acids forming a well‐defined antiparallel three‐helix structure. Two of the three helices are involved in ligand binding, whereas the third helix provides structural support to the three‐helix bundle. The small size and the stable three‐helix structure are two attractive properties comprised in the Z‐domain, but a further reduction in size of the protein is valuable for several reasons. Reduction in size facilitates synthetic production of any protein‐based molecule, which is beneficial from an economical viewpoint. In addition, a smaller protein is easier to manipulate through chemical modifications. By omitting the third stabilizing helix from the Z‐domain and joining the N‐ and C‐termini by a native peptide bond, the affinity protein obtains the advantageous properties of a smaller scaffold and in addition becomes resistant to exoproteases. We here demonstrate the synthesis and evaluation of a novel cyclic two‐helix Z‐domain. The molecule has retained affinity for its target protein, is resistant to heat treatment, and lacks both N‐ and C‐termini. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Crystallization and X-ray analysis of a single fab binding domain from protein L of Peptostreptococcus magnus

Protein L is a multi domain cell wall constituent of certain strains of Peptostreptococcus magnus which binds to the variable domain of immunoglobulin κ-light chains. A single immunoglobulin-binding domain of M_r = 9000 from this protein has been isolated and crystallized. The crystals are of space group P4₂2₁2, with cell dimensions a = b = 66.9 Å, c = 68.3 Å, and diffract to at least 2.2 Å resolution. The asymmetric unit of the crystal contains two molecules of the protein L domain, related by a noncrystallographic 2-fold axis, as revealed by a self-rotation function calculated with native diffraction data. © 1995 Wiley-Liss, Inc.     

Occurrence and characterization of PEND proteins in angiosperms

The PEND protein is a DNA-binding protein in the inner envelope membrane of the developing chloroplast. It consists of a short pre-sequence, an N-terminal DNA-binding domain (cbZIP), a central repeat domain, and a C-terminal transmembrane domain. PEND homologs have been detected in various angiosperms, including Arabidopsis thaliana, Brassica napus, Medicago truncatula, cucumber and cherry. Monocot homologs have also been detected in barley and rice, but sequence conservation was low in monocots. PEND-related sequences have not been detected in non-flowering plants and algae. Green fluorescent protein fusions consisting of the N-terminal as well as full-length PEND homologs in A. thaliana and B. napus were targeted to chloroplasts, and localized to nucleoids and chloroplast periphery, respectively. Immunoblot analysis suggested that crucifer homologs were present in chloroplasts probably as a dimer, as in the case of pea. These results suggest that PEND protein is present in angiosperms, and the homologs in crucifers are functionally analogous to the PEND protein in pea.     

Identification of mycobacterial lectins from genomic data

Sixty‐four sequences containing lectin domains with homologs of known three‐dimensional structure were identified through a search of mycobacterial genomes. They appear to belong to the β‐prism II, the C‐type, the Microcystis virdis (MV), and the β‐trefoil lectin folds. The first three always occur in conjunction with the LysM, the PI‐PLC, and the β‐grasp domains, respectively while mycobacterial β‐trefoil lectins are unaccompanied by any other domain. Thirty heparin binding hemagglutinins (HBHA), already annotated, have also been included in the study although they have no homologs of known three‐dimensional structure. The biological role of HBHA has been well characterized. A comparison between the sequences of the lectin from pathogenic and nonpathogenic mycobacteria provides insights into the carbohydrate binding region of the molecule, but the structure of the molecule is yet to be determined. A reasonable picture of the structural features of other mycobacterial proteins containing one or the other of the four lectin domains can be gleaned through the examination of homologs proteins, although the structure of none of them is available. Their biological role is also yet to be elucidated. The work presented here is among the first steps towards exploring the almost unexplored area of the structural biology of mycobacterial lectins. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long‐range hydrophobic interactions

A 20‐residue peptide, IG(42–61), derived from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three‐state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a β‐hairpin shape over a wide range of temperatures (283–313 K). Our results show that IG (42–61) possesses a well‐organized three‐dimensional structure stabilized by long‐range hydrophobic interactions (Tyr50 ··· Phe57 and Trp48 ··· Val59) at T = 283 K and (Trp48 ··· Val59) at 305 and 313 K. The mechanism of β‐hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a C_α structure‐based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well‐separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non‐native contact interactions in different folding scenarios. These findings strongly correlate with the protein free‐energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.Proteins 2013; 81:1727–1737. © 2013 Wiley Periodicals, Inc.     

Structural characterization of the heme‐based oxygen sensor,AfGcHK, its interactions with the cognate response regulator,and their combined mechanism of action in a bacterial two‐component signaling system

The oxygen sensor histidine kinase AfGcHK from the bacterium Anaeromyxobacter sp. Fw 109‐5 forms a two‐component signal transduction system together with its cognate response regulator (RR). The binding of oxygen to the heme iron of its N‐terminal sensor domain causes the C‐terminal kinase domain of AfGcHK to autophosphorylate at His183 and then transfer this phosphate to Asp52 or Asp169 of the RR protein. Analytical ultracentrifugation revealed that AfGcHK and the RR protein form a complex with 2:1 stoichiometry. Hydrogen‐deuterium exchange coupled to mass spectrometry (HDX‐MS) suggested that the most flexible part of the whole AfGcHK protein is a loop that connects the two domains and that the heme distal side of AfGcHK, which is responsible for oxygen binding, is the only flexible part of the sensor domain. HDX‐MS studies on the AfGcHK:RR complex also showed that the N‐side of the H9 helix in the dimerization domain of the AfGcHK kinase domain interacts with the helix H1 and the β‐strand B2 area of the RR protein's Rec1 domain, and that the C‐side of the H8 helix region in the dimerization domain of the AfGcHK protein interacts mostly with the helix H5 and β‐strand B6 area of the Rec1 domain. The Rec1 domain containing the phosphorylable Asp52 of the RR protein probably has a significantly higher affinity for AfGcHK than the Rec2 domain. We speculate that phosphorylation at Asp52 changes the overall structure of RR such that the Rec2 area containing the second phosphorylation site (Asp169) can also interact with AfGcHK. Proteins 2016; 84:1375–1389. © 2016 Wiley Periodicals, Inc.     

A tobacco protein kinase, NPK2, has a domain homologous to a domain found in activators of mitogen-activated protein kinases (MAPKKs)

A cDNA (cNPK2) that encodes a protein of 518 amino acids was isolated from a library prepared from poly(A)⁺ RNAs of tobacco cells in suspension culture. The N-terminal half of the predicted NPK2 protein is similar in amino acid sequence to the catalytic domains of kinases that activate mitogen-activated protein kinases (designated here MAPKKs) from various animals and to those of yeast homologs of MAPKKs. The N-terminal domain of NPK2 was produced as a fusion protein in Escherichia coli, and the purified fusion protein was found to be capable of autophosphorylation of threonine and serine residues. These results indicate that the N-terminal domain of NPK2 has activity of a serine/threonine protein kinase. Southern blot analysis showed that genomic DNAs from various plant species, including Arabidopsis thaliana and sweet potato, hybridized strongly with cNPK2, indicating that these plants also have genes that are closely related to the gene for NPK2. The structural similarity between the catalytic domain of NPK2 and those of MAPKKs and their homologs suggests that tobacco NPK2 corresponds to MAPKKs of other organisms. Given the existence of plant homologs of an MAP kinase and tobacco NPK1, which is structurally and functionally homologous to one of the activator kinases of yeast homologs of MAPKK (MAPKKKs), it seems likely that a signal transduction pathway mediated by a protein kinase cascade that is analogous to the MAP kinase cascades proposed in yeasts and animals, is also conserved in plants.     

Co‐evolutionary analysis implies auxiliary functions of HSP110 in Plasmodium falciparum

Plasmodium falciparum encounters frequent environmental challenges during its life cycle which makes productive protein folding immensely challenging for its metastable proteome. To identify the important components of protein folding machinery involved in maintaining P. falciparum proteome, we performed a proteome‐wide phylogenetic profiling across various species. We found that except HSP110, the parasite lost all other cytosolic nucleotide exchange factors essential for regulating HSP70 which is the centrum of the protein folding network. Evolutionary and structural analysis shows that besides its canonical interaction with HSP70, PfHSP110 has acquired sequence insertions for additional dynamic interactions. Molecular co‐evolution profile depicts that the co‐evolving proteins of PfHSP110 belong to distinct pathways like genetic variation, DNA repair, fatty acid biosynthesis, protein modification/trafficking, molecular motions, and apoptosis. These proteins exhibit unique physiochemical properties like large size, high iso‐electric point, low solubility, and antigenicity, hence require PfHSP110 chaperoning to attain functional state. Co‐evolving protein interaction network suggests that PfHSP110 serves as an important hub to coordinate protein quality control, survival, and immune evasion pathways in the parasite. Overall, our findings highlight potential accessory roles of PfHSP110 that may provide survival advantage to the parasite during its lifecycle and febrile conditions. The data also open avenues for experimental validation of auxiliary functions of PfHSP110 and their exploration for design of better antimalarial strategies. Proteins 2015; 83:1513–1525. © 2015 Wiley Periodicals, Inc.     

AtIRE1A/AtIRE1B and AGB1 independently control two essential unfolded protein response pathways in Arabidopsis

The endoplasmic reticulum (ER) has the ability to maintain the balance between demand for and synthesis of secretory proteins. To ensure protein‐folding homeostasis in the ER, cells invoke signaling pathways known as the unfolded protein response (UPR). To initiate UPR, yeasts largely rely on a conserved sensor, IRE1. In metazoans, there are at least three independent UPR signalling pathways. Some UPR transducers have been identified in plants, but no genetic interaction among them has yet been examined. The Arabidopsis genome encodes two IRE1 sequence homologs, AtIRE1A and AtIRE1B. Here we provide evidence that AtIRE1A and AtIRE1B have overlapping functions that are essential for the plant UPR. A double mutant of AtIRE1A and AtIRE1B, atire1a atire1b, showed reduced ER stress tolerance and a compromised UPR activation phenotype. We have also established that Arabidopsis AGB1, a subunit of the ubiquitous heterotrimeric GTP‐binding protein family, and AtIRE1A/AtIRE1B independently control two essential plant UPR pathways. By demonstrating that atire1a atire1b has a short root phenotype that is enhanced by an agb1 loss‐of‐function mutation, we have identified a role for UPR transducers in organ growth regulation.     

The stability of Taq DNA polymerase results from a reduced entropic folding penalty; identification of other thermophilic proteins with similar folding thermodynamics

The thermal stability of Taq DNA polymerase is well known, and is the basis for its use in PCR. A comparative thermodynamic characterization of the large fragment domains of Taq (Klentaq) and E. coli (Klenow) DNA polymerases has been performed by obtaining full Gibbs‐Helmholtz stability curves of the free energy of folding (ΔG) versus temperature. This analysis provides the temperature dependencies of the folding enthalpy and entropy (ΔH and ΔS), and the heat capacity (ΔC_p) of folding. If increased or enhanced non‐covalent bonding in the native state is responsible for enhanced thermal stabilization of a protein, as is often proposed, then an enhanced favourable folding enthalpy should, in general, be observed for thermophilic proteins. However, for the Klenow–Klentaq homologous pair, the folding enthalpy (ΔH_fold) of Klentaq is considerably less favorable than that of Klenow at all temperatures. In contrast, it is found that Klentaq's extreme free energy of folding (ΔG_fold) originates from a significantly reduced entropic penalty of folding (ΔS_fold). Furthermore, the heat capacity changes upon folding are similar for Klenow and Klentaq. Along with this new data, comparable extended analysis of available thermodynamic data for 17 other mesophilic–thermophilic protein pairs (where enough applicable thermodynamic data exists) shows a similar pattern in seven of the 18 total systems. When analyzed with this approach, the more familiar “reduced ΔC_p mechanism” for protein thermal stabilization (observed in a different six of the 18 systems) frequently manifests as a temperature dependent shift from enthalpy driven stabilization to a reduced‐entropic‐penalty model. Proteins 2014; 82:785–793. © 2013 Wiley Periodicals, Inc.