Homology modeling and molecular dynamics study of West Nile virus NS3 protease: a molecular basis for the catalytic activity increased by the NS2B cofactor

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.     

Phosphonate inhibitors of West Nile virus NS2B/NS3 protease

West Nile virus (WNV) is a member of the flavivirus genus belonging to the Flaviviridae family. The viral serine protease NS2B/NS3 has been considered an attractive target for the development of anti-WNV agents. Although several NS2B/NS3 protease inhibitors have been described so far, most of them are reversible inhibitors. Herein, we present a series of α-aminoalkylphosphonate diphenyl esters and their peptidyl derivatives as potent inhibitors of the NS2B/NS3 protease. The most potent inhibitor identified was Cbz-Lys-Arg-(4-GuPhe)^P(OPh)₂ displaying K_i and k₂/K_i values of 0.4 µM and 28 265 M^?1s^?1, respectively, with no significant inhibition of trypsin, cathepsin G, and HAT protease.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.     

Complex of NS3 protease and NS4A peptide of BK strain hepatitis C virus: a 2.2 A resolution structure in a hexagonal crystal form.

The crystal structure of the NS3 protease of the hepatitis C virus (BK strain) has been determined in the space group P6(3)22 to a resolution of 2.2 A. This protease is bound with a 14-mer peptide representing the central region of the NS4A protein. There are two molecules of the NS3(1-180)-NS4A(21'-34') complex per asymmetric unit. Each displays a familiar chymotrypsin-like fold that includes two beta-barrel domains and four short alpha-helices. The catalytic triad (Ser-139, His-57, and Asp-81) is located in the crevice between the beta-barrel domains. The NS4A peptide forms an almost completely enclosed peptide surface association with the protease. In contrast to the reported H strain complex of NS3 protease-NS4A peptide in a trigonal crystal form (Kim JL et al., 1996, Cell 87:343-355), the N-terminus of the NS3 protease is well-ordered in both molecules in the asymmetric unit of our hexagonal crystal form. The folding of the N-terminal region of the NS3 protease is due to the formation of a three-helix bundle as a result of crystal packing. When compared with the unbound structure (Love RA et al., 1996, Cell 87:331-342), the binding of the NS4A peptide leads to the ordering of the N-terminal 28 residues of the NS3 protease into a beta-strand and an alpha-helix and also causes local rearrangements important for a catalytically favorable conformation at the active site. Our analysis provides experimental support for the proposal that binding of an NS4A-mimicking peptide, which increases catalytic rates, is necessary but not sufficient for formation of a well-ordered, compact and, hence, highly active protease molecule.     

Designing cyclic peptide inhibitor of dengue virus NS3-NS2B protease by using molecular docking approach

Peptides are preferred for designing inhibitors because of their high activity and specificity. Seven cyclopentapeptide inhibitors were designed in this study against dengue virus type 2 (DEN-2) NS3-NS2B protease: CKRRC, CGRRC, CRGRC, CRTRC, CTRRC, CKRKC and CRRKC. Docking analysis was performed to study the enzyme-inhibitor binding interactions. The free energy binding and estimated Ki values for all the inhibitors were found to be small (within micromolar range), indicating that the inhibitors bind considerably well to the binding site. The results showed that the cyclopentapeptide CKRKC was the best peptide inhibitor candidate with estimated free binding energy of -8.39 kcal/mol and Ki of 0.707 μM when compared to the standard inhibitor Bz-Nle-Lys-Arg-Arg-H that has been experimentally tested and shown to exhibit Ki value of 5.8 μM. Several modes of weak interactions were observed between the cyclopentapeptide CKRKC and the active site of DEN-2 NS3-NS2B protease. Thus, the cyclopentapeptide is proposed as a potential inhibitor to the NS3-NS2B protease activities of DEN-2. While these preliminary results are promising, further experimental investigation is necessary to validate the results.     

Activation of the West Nile virus NS3 protease: Molecular dynamics evidence for a conformational selection mechanism

The flaviviral nonstructural 3 protease (NS3pro) is essential for virus replication and is therefore a pharmaceutically relevant target to fight Dengue and West Nile virus (WNV). NS3pro is a chymotrypsin‐like serine protease which requires a polypeptide cofactor (NS2B) for activation. Recent X‐ray crystallography studies have led to the suggestion that the substrate binds to the two‐component NS2B‐NS3pro enzyme by an induced‐fit mechanism. Here, multiple explicit water molecular dynamics simulations of the WNV NS2B‐NS3pro enzyme show that the active conformation of the NS2B cofactor (in which its β‐loop is part of the substrate binding site) is stable over a 50‐ns time scale even in the absence of the inhibitor. The partial and reversible opening of the NSB2 β‐loop and its correlated motion with an adjacent NS3pro loop, both observed in the simulations started from the active conformation, are likely to facilitate substrate binding and product release. Moreover, in five of eight simulations without inhibitor (started from two X‐ray structures both with improperly formed oxyanion hole) the Thr132‐Gly133 peptide bond flips spontaneously thereby promoting the formation of the catalytically competent oxyanion hole, which then stays stable until the end of the runs. The simulation results provide evidence at atomic level of detail that the substrate binds to the NS2B‐NS3pro enzyme by conformational selection, rather than induced‐fit mechanism.     

Antiviral potential of 4-hydroxypanduratin A,secondary metabolite of Fingerroot,Boesenbergia pandurata (Schult.), towards Japanese Encephalitis virus NS2B/NS3 protease

4-hydroxypanduratin A is a secondary metabolite of Boesenbergia pandurata Schult. (Fingerroot) plant with various pharmacological activities such as neuroprotective, potent antioxidant, antibacterial and antifungal. Flaviviral NS2B/NS3 protease activity is essential for polyprotein processing and viral replication for Japanese Encephalitis Virus (JEV), a major cause of Acute Encephaltis in Asia. Inhibition of formation of this complex by arresting the binding of NS2B with NS3 would reduce the enzyme''s activity to meager proportions and hence would prevent further viral proliferation. The automated 3D structure of NS2B protein of the JEV GP78 was predicted based on the sequence-to-structure-to-function paradigm using I-TASSER and the function of NS2B protein was inferred by matching to other known proteins. The stereochemical quality of predicted structure was checked by PROCHECK. The antiviral activity of 4-hydroxypanduratin A against NS2B protein as a potential drug has been elucidated in this paper. Docking simulation analysis showed 4-hydroxypanduratin A as potential inhibitor of NS2B protein/cofactor which is necessary for NS3 protease activity. 220 derivatives of 4-hydroxypanduratin A were virtually screened with rigid criteria of Lipinski''s rule of 5 using Autodock4.2. 4-hydroxypanduratin A was found interacting with target hydrophilic domain in NS2B protein by two Hbonds (Gly80 and Asp81) with active residues, several hydrophobic interactions, Log P value of 5.6, inhibition constant (Ki) of 51.07nM and lowest binding energy of -9.95Kcal/Mol. Hence, 4-hydroxypanduratin A targeted to Site 2 will have sufficient profound effect to inhibit protease activity to abrogate viral replication. It could be a promising potential drug candidate for JEV infections using NS2B Site 2 as a Drug target.     

Characterisation of the role of zinc in the hepatitis C virus NS2/3 auto-cleavage and NS3 protease activities

Cleavage of the hepatitis C virus polyprotein between the non-structural NS2 and NS3 proteins is mediated by a poorly characterised auto-proteolytic activity that maps to the C terminus of NS2 and the N terminus of NS3, but is distinct from the NS3 protease activity responsible for downstream cleavages in the polyprotein. We have exploited the fact that the minimal precursor (residues 904-1206 of the HCV polyprotein) can be expressed as an insoluble protein in Escherichia coli and subsequently refolded into a form active for both auto-cleavage and NS3 protease activity, to further characterise the NS2/3 auto-cleavage activity. We show that both activities are zinc-dependent and show an absolute requirement for cysteine residues 1123, 1125 and 1171 within NS3. In contrast cysteine 922 (within NS2) is only required for NS2/3 auto-cleavage activity and histidine 1175 is only required for NS3 activity. Although the complete NS3 protease domain (including the C-terminal alpha-helix) is required for NS2/3 auto-cleavage, the activity of the NS3 protease is not essential. Lastly we show that the NS2/3 auto-cleavage activity is more sensitive to zinc chelation by 1,10-phenanthroline than the NS3 protease activity. This observation is consistent with different conformations of the precursor competent for either NS2/3 auto-cleavage or NS3 protease activity; these two conformations can be distinguished by their relative strength and geometry of zinc coordination.     

Mutagenesis of D80-82 and G83 residues in West Nile Virus NS2B: Effects on NS2B-NS3 activity and viral replication

Flaviviral NS2B is a required cofactor for NS3 serine protease activity and plays an important role in promoting functional NS2B-NS3 protease configuration and maintaining critical interactions with protease catalysis substrates. The residues D⁸⁰DDG in West Nile virus (WNV) NS2B are important for protease activity. To investigate the effects of D⁸⁰DDG in NS2B on protease activity and viral replication, the negatively charged region D⁸⁰DD and the conserved residue G83 of NS2B were mutated (D⁸⁰DD/E⁸⁰EE, D⁸⁰DD/K⁸⁰KK, D⁸⁰DD/A⁸⁰AA, G83F, G83S, G83D, G83K, and G83A), and NS3 D75A was designated as the negative control. The effects of the mutations on NS2B-NS3 activity, viral translation, and viral RNA replication were analyzed using kinetic analysis of site-directed enzymes and a transient replicon assay. All substitutions resulted in significantly decreased enzyme activity and blocked RNA replication. The negative charge of D⁸⁰DD is not important for maintaining NS2B function, but side chain changes in G83 have dramatic effects on protease activity and RNA replication. These results demonstrate that NS2B is important for viral replication and that D⁸⁰DD and G83 substitutions prevent replication; they will be useful for understanding the relationship between NS2B and NS3.     

Structural evidence for regulation and specificity of flaviviral proteases and evolution of the Flaviviridae fold

Pathogenic members of the flavivirus family, including West Nile Virus (WNV) and Dengue Virus (DV), are growing global threats for which there are no specific treatments. The two-component flaviviral enzyme NS2B-NS3 cleaves the viral polyprotein precursor within the host cell, a process that is required for viral replication. Here, we report the crystal structure of WNV NS2B-NS3pro both in a substrate-free form and in complex with the trypsin inhibitor aprotinin/BPTI. We show that aprotinin binds in a substrate-mimetic fashion in which the productive conformation of the protease is fully formed, providing evidence for an "induced fit" mechanism of catalysis and allowing us to rationalize the distinct substrate specificities of WNV and DV proteases. We also show that the NS2B cofactor of WNV can adopt two very distinct conformations and that this is likely to be a general feature of flaviviral proteases, providing further opportunities for regulation. Finally, by comparing the flaviviral proteases with the more distantly related Hepatitis C virus, we provide insights into the evolution of the Flaviviridae fold. Our work should expedite the design of protease inhibitors to treat a range of flaviviral infections.     

Identification of novel thiadiazoloacrylamide analogues as inhibitors of dengue-2 virus NS2B/NS3 protease

Dengue virus is endemic throughout tropical and subtropical regions, and cause severe epidemic diseases. The NS2B/NS3 protease is a promising drug target for dengue virus. Herein, we report the discovery and modification of a novel class of thiadiazoloacrylamide derivatives with potent inhibitory activity against the NS2B/NS3 protease. Thiadiazolopyrimidinone 1 was firstly determined as a new chemical structure against NS2B/NS3 from a commercial compound library. Then, we sought to identify similar compounds with the thiadiazoloacrylamide core that would exhibit better activity. A series of analogues were synthesized and fourteen of them were identified with strong inhibitory activities, in which the nitrile group in the linker part was discovered as an essential group for the inhibitory activity. The best of these (8b) demonstrated an IC₅₀ at 2.24 μM based on in vitro DENV2 NS2B-NS3pro assays.     

Effects of NS2B-NS3 protease and furin inhibition on West Nile and Dengue virus replication

West Nile virus (WNV) and Dengue virus (DENV) replication depends on the viral NS2B-NS3 protease and the host enzyme furin, which emerged as potential drug targets. Modification of our previously described WNV protease inhibitors by basic phenylalanine analogs provided compounds with reduced potency against the WNV and DENV protease. In a second series, their decarboxylated P1-trans-(4-guanidino)cyclohexylamide was replaced by an arginyl-amide moiety. Compound 4-(guanidinomethyl)-phenylacetyl-Lys-Lys-Arg-NH₂ inhibits the NS2B-NS3 protease of WNV with an inhibition constant of 0.11?µM. Due to the similarity in substrate specificity, we have also tested the potency of our previously described multibasic furin inhibitors. Their further modification provided chimeric inhibitors with additional potency against the WNV and DENV proteases. A strong inhibition of WNV and DENV replication in cell culture was observed for the specific furin inhibitors, which reduced virus titers up to 10,000-fold. These studies reveal that potent inhibitors of furin can block the replication of DENV and WNV.     

Ligand-bound structures of the dengue virus protease reveal the active conformation

Dengue is a mosquito-borne viral hemorrhagic disease that is a major threat to human health in tropical and subtropical regions. Here we report crystal structures of a peptide covalently bound to dengue virus serotype 3 (DENV-3) protease as well as the serine-protease inhibitor aprotinin bound to the same enzyme. These structures reveal, for the first time, a catalytically active, closed conformation of the DENV protease. In the presence of the peptide, the DENV-3 protease forms the closed conformation in which the hydrophilic β-hairpin region of NS2B wraps around the NS3 protease core, in a manner analogous to the structure of West Nile virus (WNV) protease. Our results confirm that flavivirus proteases form the closed conformation during proteolysis, as previously proposed for WNV. The current DENV-3 protease structures reveal the detailed interactions at the P4' to P3 sites of the substrate. The new structural information explains the sequence preference, particularly for long basic residues in the nonprime side, as well as the difference in substrate specificity between the WNV and DENV proteases at the prime side. Structural analysis of the DENV-3 protease-peptide complex revealed a pocket that is formed by residues from NS2B and NS3; this pocket also exists in the WNV NS2B/NS3 protease structure and could be targeted for potential antivirus development. The structural information presented in the current study is invaluable for the design of specific inhibitors of DENV protease.     

A small compound targeting the interaction between nonstructural proteins 2B and 3 inhibits dengue virus replication

The non-structural protein NS2B/NS3 serine-protease complex of the dengue virus (DENV) is required for the maturation of the viral polyprotein. Dissociation of the NS2B cofactor from NS3 diminishes the enzymatic activity of the complex. In this study, we identified a small molecule inhibitor that interferes with the interaction between NS2B and NS3 using structure-based screening and a cell-based viral replication assay. A library containing 661,417 small compounds derived from the Molecular Operating Environment lead-like database was docked to the NS2B/NS3 structural model. Thirty-nine compounds with high scores were tested in a secondary screening using a cell-based viral replication assay. SK-12 was found to inhibit replication of all DENV serotypes (EC₅₀ = 0.74–4.92 μM). In silico studies predicted that SK-12 pre-occupies the NS2B-binding site of NS3. Steady-state kinetics using a fluorogenic short peptide substrate demonstrated that SK-12 is a noncompetitive inhibitor against the NS2B/NS3 protease. Inhibition to Japanese encephalitis virus by SK-12 was relatively weak (EC₅₀ = 29.81 μM), and this lower sensitivity was due to difference in amino acid at position 27 of NS3. SK-12 is the promising small-molecule inhibitor that targets the interaction between NS2B and NS3.     

An inhibition model of BPTI to unlinked dengue virus NS2B-NS3 protease

One approach to treating the dengue virus infection is to inhibit its NS2B-NS3 protease that plays a vital role in virus maturation. However, the lack of structural information on the active conformation of the protease hindered related drug design. With a co-expression system, we obtained the active two-component protease in its unlinked form. BPTI shows strong competitive inhibitory activity (K_i = 6.5 nM) against this unlinked protease, which adopts a closed conformation. Based on the biochemical and NMR perturbation information, an inhibition model of BPTI to NS2B-NS3 protease is proposed.     

Mechanistic role of NS4A and substrate in the activation of HCV NS3 protease

Hepatitis C virus (HCV) NS3 protease is the key enzyme for its maturation. Three hypotheses have been advanced in the literature to demonstrate the mechanism of the activation of the HCV NS3 protease. A virus-encoded protein NS4A and substrate are proposed to be involved in the activation of the HCV NS3 protease. However, the three hypotheses are not completely consistent with one another. Multiple molecular dynamics simulations were performed on various NS3 protease systems: free NS3 protease, NS3/4A, NS3/inhibitor, and NS3/4A/inhibitor complexes, to further unravel the mechanism of the activation of the NS3 protease. Simulation results suggest that the binding of NS4A induces a classic serine protease conformation of the catalytic triad of the NS3 protease. NS4A rearranges the secondary structure of both the N-terminus and catalytic site of the NS3 protease, reduces the mobility of the global structure of the NS3 protease, especially the catalytic site, and provides a rigid and tight structure, except for the S1 pocket, for the binding and hydrolysis of substrates. The binding of substrate also contributes to the activation of the NS3 protease by an induced-fit of the classic serine protease catalytic triad. However, the global structure of the NS3 protease is still loose and highly flexible without stable secondary structural elements, such as helix α0 at the N-terminus and helix α1 and β-sheet E1-F1 at the catalytic site. The structure of the NS3 protease without NS4A is not suitable for the binding and hydrolysis of substrates.     

Switching the substrate specificity of the two-component NS2B-NS3 flavivirus proteinase by structure-based mutagenesis

The flavivirus NS2B-NS3(pro)teinase is an essential element in the proteolytic processing of the viral precursor polyprotein and therefore a potential drug target. Recently, crystal structures and substrate preferences of NS2B-NS3pro from Dengue and West Nile viruses (DV and WNV) were determined. We established that the presence of Gly-Gly at the P1'-P2' positions is optimal for cleavage by WNV NS3pro, whereas DV NS3pro tolerates well the presence of bulky residues at either P1' or P2'. Structure-based modeling suggests that Arg(76) and Pro(131)-Thr(132) limit the P1'-P2' subsites and restrict the cleavage preferences of the WNV enzyme. In turn, Leu(76) and Lys(131)-Pro(132) widen the specificity of DV NS3pro. Guided by these structural models, we expressed and purified mutant WNV NS2B-NS3pro and evaluated cleavage preferences by using positional scanning of the substrate peptides in which the P4-P1 and the P3'-P4' positions were fixed and the P1' and P2' positions were each randomized. We established that WNV R76L and P131K-T132P mutants acquired DV-like cleavage preferences, whereas T52V had no significant effect. Our work is the first instance of engineering a viral proteinase with switched cleavage preferences and should provide valuable data for the design of optimized substrates and substrate-based selective inhibitors of flaviviral proteinases.     

A FRET-based assay for the discovery of West Nile Virus NS2B-NS3 protease inhibitors

The West Nile Virus (WNV) has been a worldwide epidemic since the early 1990s. Currently there are no therapeutic treatments for WNV infections. One particular avenue of treatment is inhibition of the NS2B-NS3 protease, an enzyme that is crucial for WNV replication. In our effort to increase the number of NS2B-NS3 protease inhibitors, we report a novel FRET-based high throughput assay for the discovery of WNV NS2B-NS3 protease inhibitors. For this assay, a FRET-based peptide substrate was synthesized and kinetically characterized with the NS2B-NS3 protease. The new substrate exhibits a K_m of 3.35 ± 0.31 μM, a k_cat of 0.0717 ± 0.0016 s^?1 and a k_cat/K_m of 21,400 ± 2000 M^?1 s^?1.     

Yellow fever virus NS2B/NS3 protease: Hydrolytic Properties and Substrate Specificity

Here we report the hydrolytic behavior of recombinant YFV NS2B/NS3 protease against FRET substrates mimicking the prime and non-prime region of the natural polyprotein cleavage sites. While the P2-P′1 motif is the main factor associated with the catalytic efficiency of Dengue (DV) and West Nile Virus (WNV) protease, we show that the k_cat/K_m of YFV NS2B/NS3 varied by more than two orders of magnitude, despite the presence of the same motif in all natural substrates. The catalytic significance of this homogeneity – a unique feature among worldwide prominent flavivirus – was kinetically analyzed using FRET peptides containing all possible combinations of two and three basic amino acids in tandem, and Arg and Lys residues produced distinct effects on k_cat/K_m. The parallel of our data with those obtained in vivo by Chambers et al. (1991) restrains the idea that these sites co-evolved with the NS2B/NS3 protease to promote highly efficient hydrolysis and supports the notion that secondary substrate interaction distant from cleavage sites are the main factor associated with the different hydrolytic rates on YFV NS2B-NS3pro natural substrates.     

Mechanistic and kinetic characterization of hepatitis C virus NS3 protein interactions with NS4A and protease inhibitors

The mechanism and kinetics of the interactions between ligands and immobilized full‐length hepatitis C virus (HCV) genotype 1a NS3 have been characterized by SPR biosensor technology. The NS3 interactions for a series of NS3 protease inhibitors as well as for the NS4A cofactor, represented by a peptide corresponding to the sequence interacting with the enzyme, were found to be heterogeneous. It may represent interactions with two stable conformations of the protein. The NS3–NS4A interaction consisted of a high‐affinity (K_D = 50 nM) and a low‐affinity (K_D = 2 µM) interaction, contributing equally to the overall binding. By immobilizing NS3 alone or together with NS4A it was shown that all inhibitors had a higher affinity for NS3 in the presence of NS4A. NS4A thus has a direct effect on the binding of inhibitors to NS3 and not only on catalysis. As predicted, the mechanism‐based inhibitor VX 950 exhibited a time‐dependent interaction with a slow formation of a stable complex. BILN 2061 or ITMN‐191 showed no signs of time‐dependent interactions, but ITMN‐191 had the highest affinity of the tested compounds, with both the slowest dissociation (k_off) and fastest association rate, closely followed by BILN 2061. The k_off for the inhibitors correlated strongly with their NS3 protease inhibitory effect as well as with their effect on replication of viral proteins in replicon cell cultures, confirming the relevance of the kinetic data. This approach for obtaining kinetic and mechanistic data for NS3 protease inhibitor and cofactor interactions is expected to be of importance for understanding the characteristics of HCV NS3 functionality as well as for anti‐HCV lead discovery and optimization. Copyright © 2010 John Wiley & Sons, Ltd.