research communications
Cwp84, a Clostridium difficile cysteine protease, exhibits conformational flexibility in the absence of its propeptide
aDepartment of Biology and Biochemistry, University of Bath, Claverton Down, Bath BA2 7AY, England, and bPublic Health England, Porton Down, Salisbury SP4 0JG, England
*Correspondence e-mail: bsskra@bath.ac.uk
In recent decades, the global healthcare problems caused by Clostridium difficile have increased at an alarming rate. A greater understanding of this antibiotic-resistant bacterium, particularly with respect to how it interacts with the host, is required for the development of novel strategies for fighting C. difficile infections. The (S-layer) of C. difficile is likely to be of significant importance to host–pathogen interactions. The mature S-layer is formed by a proteinaceous array consisting of multiple copies of a high-molecular-weight and a low-molecular-weight S-layer protein. These components result from the cleavage of SlpA by Cwp84, a cysteine protease. The structure of a truncated Cwp84 active-site mutant has recently been reported and the key features have been identified, providing the first structural insights into the role of Cwp84 in the formation of the S-layer. Here, two structures of Cwp84 after propeptide cleavage are presented and the three conformational changes that are observed are discussed. These changes result in a reconfiguration of the active site and exposure of the hydrophobic pocket.
Keywords: Clostridium difficile; surface layer; Cwp84; host–pathogen interactions.
1. Introduction
The mainly nosocomially acquired bacterium Clostridium difficile is the primary aetiological agent of antibiotic-associated diarrhoea, pseudomembranous colitis and toxic megacolon (Kachrimanidou & Malisiovas, 2011; Rupnik et al., 2009). C. difficile causes tens of thousands of deaths annually and is a major global economic burden (Bouza, 2012; Dubberke & Olsen, 2012). This results in a need for further understanding of this Gram-positive bacterium, particularly towards the development of novel treatment strategies.
C. difficile presents a paracrystalline protein array as its outermost structure, known as a (S-layer; Fagan & Fairweather, 2014; Kawata et al., 1984). The S-layer is primarily comprised of two proteins: the high-molecular-weight S-layer protein (HMW SLP) and the low-molecular-weight S-layer protein (LMW SLP), which result from the cleavage of the S-layer precursor protein SlpA (Calabi et al., 2001; Cerquetti et al., 2000; Karjalainen et al., 2001). This cleavage is performed by Cwp84, a cysteine protease (de la Riva et al., 2011; Janoir et al., 2004; Kirby et al., 2009) that is one of up to 28 SlpA paralogues coded for by the C. difficile genome (Calabi et al., 2001; Fagan et al., 2011). All SlpA paralogues, or Cwps (cell-wall proteins), contain three Pfam 04122 cell-wall binding domains, as seen in HMW SLP, while many also contain a `functional' domain. As well as cleaving SlpA, Cwp84 has also been shown to be able to degrade the extracellular matrix proteins fibronectin, vitronectin and laminin, but is unable to break down gelatine (Janoir et al., 2007).
It has previously been demonstrated that knocking out or inhibiting Cwp84 activity results in severe growth defects in C. difficile but is not lethal (Dang et al., 2010; de la Riva et al., 2011; Kirby et al., 2009). This effect is likely to be due to the resultant perturbation of the S-layer structure causing stress to the bacterium, which, if it could be replicated in vivo, may significantly reduce the effects of C. difficile infection. It has been shown that a Cwp84 knockout was still able to cause symptoms in hamsters (Kirby et al., 2009), but it has been speculated that perturbation of the S-layer may result in the bacterium being more susceptible to antibiotics (Dang et al., 2010). A detailed understanding of the structure of the S-layer at the atomic level will be vital to causing a similar perturbation in vivo, preventing or reducing the severity of infections.
Recently, the et al., 2014). Key features of the structure were identified, including the propeptide, the cysteine protease domain and the previously uncharacterized `lectin-like' domain. The cysteine protease domain assumes a cathepsin L-like fold with notable differences in the occluding loop region and the proregion binding loop. The lectin-like domain, which bears significant secondary-structural and tertiary-structural similarity to carbohydrate-binding proteins, consists of a twisted β-sandwich and is stabilized by a calcium ion. The C-terminal half of the propeptide assumes a typical papain protease fold, occupying the active-site groove in the opposite direction to that of the substrate, while the N-terminal half wraps around the lectin-like domain. A method of propeptide cleavage has also been discussed (Bradshaw et al., 2014).
of a truncated Cwp84 active-site mutant has been reported at a resolution of 1.4 Å (BradshawHere, two Cwp84 structures after propeptide cleavage are presented at 1.6 and 1.84 Å resolution (Fig. 1). There are at present relatively few mature cysteine protease structures in the Protein Data Bank with no ligands bound. The structures presented here therefore provide a rare opportunity to discuss the conformational changes observed in a cysteine protease upon activation. None of the conformational changes discussed appear to be shared by papain or cathepsins, the archetypes of the papain protease family.
2. Materials and methods
2.1. Protein expression and purification
Cwp84 was expressed and purified and the propeptide was cleaved as described previously (Bradshaw et al., 2014). Briefly, a construct coding for a GST-tagged truncated Cwp84 active-site mutant, Cwp8433–497C116A, was expressed in LB and purified on a GST column. The tag and propeptide were removed by incubation with trypsin at a molar ratio of approximately 10:1 at 210 K for 45 min. The resulting protein, Cwp8492–497C116A, was further purified on a size-exclusion column, eluting in 10 mM Tris–HCl pH 8.0. Successful cleavage of the propeptide was confirmed by and N-terminal sequencing.
2.2. X-ray crystallographic studies
Crystallization conditions were screened using a Phoenix crystallization robot (Art Robbins Instruments). Crystals were observed in PACT condition D7 (0.2 M NaCl, 0.1 M Tris pH 8.0, 20% PEG 6000) and MIDAS condition F7 [20%(v/v) dimethyl sulfoxide, 20%(v/v) Jeffamine M-2070] (Molecular Dimensions). The former condition was optimized by addition of 10% Silver Bullets condition E9 [0.2%(w/v) 1,4-diaminobutane, 0.2%(w/v) cystamine dihydrochloride, 0.2%(w/v) diloxanide furoate, 0.2%(w/v) sarcosine, 0.2%(w/v) spermine, 20 mM sodium HEPES pH 6.8; Hampton Research]. Both crystals were cryoprotected by addition of 1 µl 50% reservoir solution and 25% glycerol.
X-ray diffraction data were collected on beamlines I04 (at 0.9795 Å wavelength) for the crystal obtained from PACT condition D7 and I04-1 (0.9200 Å wavelength) for the crystal obtained from MIDAS condition F7 at Diamond Light Source, Didcot, England. 360° of data were collected for both crystals with 0.1° oscillations. The data were autoprocessed with XDS (Kabsch, 2010) and xia2 (Winter et al., 2013) and rescaled with AIMLESS (Evans & Murshudov, 2013), at which point 143.5° of poorer images were excluded from the MIDAS F7 data set and 60° of images were excluded from the PACT D7 data set. The structures were solved by with Phaser (McCoy et al., 2007) using the previously determined Cwp84 structure (PDB entry 4ci7) with the propeptide removed as a starting model. This was followed by with Coot (Emsley & Cowtan, 2004) and REFMAC (Murshudov et al., 2011) and validation with MolProbity (Chen et al., 2010).
3. Results
The structure of a truncated Cwp84 active-site mutant has been determined without the propeptide in two different sets of conditions (Fig. 2). Both sets of conditions resulted in crystallization in the triclinic P1 with two molecules in the but with different unit-cell parameters and significantly different packing of the two molecules. The structures were solved at 1.84 and 1.6 Å resolution. The lower resolution structure (from the crytal obtained using MIDAS condition F7) is referred to as `structure 1', while the higher resolution structure (from the crytal obtained using PACT condition D7) is referred to as `structure 2'. Structure 1 contained two calcium ions, two Jeffamine molecules and 449 water molecules, while structure 2 contained four calcium ions, eight PEG molecules, four glycerol molecules and 791 water molecules. Crystallographic statistics are summarized in Table 1.
|
As in the previous structure with the propeptide, both structures contained a calcium ion coordinated by Leu339, Glu448, Lys460 and Asn487 (Bradshaw et al., 2014). The two extra calcium ions in structure 2, which are coordinated by Asp138, Leu139, Glu141 and Glu181, displace Lys108 (Fig. 3). The four coordinating residues are moderately conserved in other papain (Bradshaw et al., 2014; Coulombe et al., 1996; Kamphuis et al., 1984). Their positioning around a positively charged moiety is likely to play an important role in stabilizing the fold of the cysteine protease domain. However, the identity of this moiety does not appear to be crucial: the displacement of Lys108 by a calcium ion simply results in a different side-chain conformation. It does not appear to have a significant structural effect. All calcium-ion assignments were performed based on coordination (Zheng et al., 2008) and the ability of the ion to fit the density; they were then confirmed with the CheckMyMetal server (Zheng et al., 2014).
The two structures without the propeptide are largely similar when compared to that with the propeptide (Fig. 2). The six chains (each structure contains two chains) superpose on each other with root-mean-square values of between 0.12 and 0.70 Å. There are, however, three loops that undergo notable conformational changes: Met160–Ser164 and Leu315–Asp320, which both form part of the active-site groove (Fig. 4), and Thr479–Pro485, which is located on the surface of the lectin-like domain (Fig. 5).
In the previously determined structure with the propeptide, the loop formed by Met160–Ser164, which forms the S1 pocket, is involved in a hydrogen-bond network with part of the propeptide (Table 2), which is supported by a high number of van der Waals interactions (Table 3; Bradshaw et al., 2014). This array of noncovalent interactions is lost upon propeptide cleavage, allowing the backbone of Met160–Ser164 to rotate, including a rotation of approximately 160° of the peptide bond between Ser161 and Gly162. This rotation is accompanied by a movement of the loop towards the active-site residues, with a 3 Å movement of the Cα atom of Gly162 (Fig. 4).
‡The bond length is 2.92 Å in both chains. |
Secondly, in the presence of the propeptide, a loop formed by Leu315–Asp320 closely interacts with the propeptide in a similar way (Tables 2 and 3; Bradshaw et al., 2014). When the propeptide is cleaved, this stabilization is lost and the loop moves away from its previous position closer to the active-site groove, including a 4 Å movement of Asp318 (Fig. 4).
The cleavage of the propeptide and the conformational changes observed result in chain A of structure 1 exhibiting an active-site groove of 1253 Å3 and chain B of structure 2 possessing an active-site groove of 993 Å3, as measured by Swiss-PdbViewer (Guex & Peitsch, 1997). The differences in the volumes are due to slight changes in side-chain conformations. The groove narrows considerably between some of the active-site pockets, resulting in the active site being recognized as multiple separate grooves in chain B of structure 1 and chain A of structure 2.
Thirdly, a loop formed by Thr479–Pro485 is involved in the formation of the hydrophobic pocket on the surface of the lectin-like domain. With the propeptide bound, Leu36 and Val39 insert into the hydrophobic pocket, stabilizing it (Bradshaw et al., 2014); this interaction is also aided by a large number of noncovalent interactions (Tables 2 and 3). After propeptide cleavage, the pocket is left somewhat more solvent-accessible, allowing the loop to become more flexible; this results in two different conformations in the structures presented here. In structure 1, apart from a slight movement away from the position of the propeptide, the conformation is largely unchanged from that seen with the propeptide, while in structure 2 the loop assumes a markedly different conformation, with the Cα atom of Glu482 over 9.5 Å from its previous position (Fig. 5).
4. Discussion
It is known that in the C1A cysteine protease family the propeptides sit in the active-site groove in the opposite direction to the substrate (Schaschke et al., 1998), allowing them to act as inhibitors (Fox et al., 1992; Wiederanders, 2003). Cleavage of the propeptide reveals the active-site groove, allowing the substrate to bind. The two structures presented here complement the previously reported structure of Cwp84 with the propeptide (Bradshaw et al., 2014). Together, these structures allow the identification and characterization of the conformational changes that occur upon activation of Cwp84.
Conformational changes have been described for three loops: Met160–Ser164, Leu315–Asp320 and Thr479–Pro485. If the conformational changes observed were in reality samplings of conformational space, there would be an increase in the B factor of the loops and different conformations may be seen in each chain. It can be concluded from this that the first two conformational changes, which show minor variations in B factors and no significant conformational difference between any of the four chains presented here, are conformational changes that occur upon propeptide cleavage. Conversely, owing to more significant increases in B factors and different conformations being observed for the third loop in the two structures, this is more likely to represent a loss of conformational stability caused by the loss of the stabilizing effect of the propeptide on the hydrophobic pocket.
The first two conformational changes described, those of Met160–Ser164 and Leu315–Asp320, result in a reconfiguration of the active site. The former is likely to be involved in formation of the S1 pocket, while the latter is likely to form the S2 pocket. There are no symmetry-related molecules in close proximity to either of these loops in either chain of any of the three structures, so they are unlikely to be crystallographic artefacts.
The rotation of Met160–Ser164 appears to result in a fairly minor change to the surface of the S1 pocket (Fig. 4). The shape is very similar; however, the previously occluded Cα atom of Gly162 produces a small hydrophobic patch and the side chain of Ser163 points into the active-site groove rather than away from it. These changes, particularly the movement of Ser163, may result in the stabilization of the conserved P1 serine of SlpA (Qazi et al., 2009). Comparing this with the archetypal C1A cysteine the loop does not appear to exhibit any conformational change that is likely to be significant upon propeptide cleavage in papain (Kamphuis et al., 1984; Roy et al., 2012), cathepsin L (Adams-Cioaba et al., 2011; Coulombe et al., 1996) or cathepsin B (Musil et al., 1991; Podobnik et al., 1997; Turk et al., 1996).
The conformation of the S2 pocket observed in the present structures is likely to be unfavourable when the propeptide is bound, as Asp318 would lie 2.5–3.0 Å from Pro68, resulting in hydrophobic/hydrophilic repulsion. The result of this is that Leu315–Asp320 are further away from the active-site residues when the propeptide is bound. This conformation is then stabilized by a range of noncovalent interactions, producing the conformation seen in the previously reported structure with the propeptide intact (Bradshaw et al., 2014).
The stabilization of this conformation will be significantly weakened upon cleavage of the propeptide, allowing the loop to move closer to where Pro68 would have been, forming the mature S2 pocket seen in the present structures. The movement of this loop results in the occlusion of the majority of the previously described S2 pocket (Fig. 4; Bradshaw et al., 2014), leaving Asp318 and Asp320 forming a negatively charged patch to which the P2 lysine of SlpA can bind.
This loop does not show any movement upon propeptide cleavage in papain (Kamphuis et al., 1984; Roy et al., 2012) or cathepsin B (Musil et al., 1991; Podobnik et al., 1997; Turk et al., 1996), while in cathepsin L there is some movement, particularly for Ser213, but the S2 selectivity residue, Ala214, does not move (Adams-Cioaba et al., 2011; Coulombe et al., 1996).
Thr479–Pro485, the residues involved in the third conformational change, form part of the previously described hydrophobic pocket on the surface of the lectin-like domain, into which Leu36 and Val39 from the propeptide insert (Bradshaw et al., 2014). Upon removal of the propeptide, all residues forming the pocket become somewhat more accessible, potentially reducing the ability of the hydrophobic effect to drive the conformation of the loop.
In structure 1, this loop assumes a similar closed conformation to when the propeptide is present, but is likely to be more flexible. This is shown by the slightly weaker electron density and resultant higher B factors. The loop adopts a vastly different, much more open conformation in structure 2 (Fig. 5), exposing the hydrophobic pocket. There are also a small number of interactions involving two symmetry-related molecules, which are likely to have an influence on the side-chain orientations of Tyr480 and Phe483 and potentially stabilize the open conformation. However, it is likely that this conformation is only permitted by the loss of the closed-conformation stabilization of the propeptide.
5. Conclusions
Two structures of Cwp84 without its propeptide have been determined, complementing the previously determined structure with the propeptide intact. Together, these structures allow the identification and discussion of structural changes that occur upon propeptide cleavage. Propeptide cleavage causes two loops which form the S1 and S2 pockets of the active-site groove to undergo conformational changes, resulting in a reconfiguration of the pockets. The new conformations seen in the present structures are likely to facilitate binding of the P1 and P2 residues of SlpA: serine and lysine, respectively. A third loop in structure 2, found on the surface of the lectin-like domain, also exhibits a conformational change. This exposes the hydrophobic pocket that was previously occluded by the propeptide. This conformation appears to be somewhat stabilized by crystal contacts. Nonetheless, the same loop in structure 1 does assume a slightly more open conformation than when the propeptide is bound. The structural data presented here provide a detailed molecular basis for the role of the propeptide of Cwp84.
Acknowledgements
We thank Diamond Light Source for access to beamlines I04 and I04-1 (MX8922), which contributed to the results presented here. This work was supported by a postgraduate studentship from Public Health England (PHE, Porton Down, England) and the University of Bath to WJB and a Medical Research Council UK project grant (MK/K027123/1) to KRA and CCS. Author contributions are as follows. WJB performed protein expression, purification and structural biology experiments, analysed the structures and wrote the manuscript. AKR and CCS edited the manuscript. KRA supervised the study, analysed the data and edited the manuscript. All authors reviewed the manuscript.
References
Adams-Cioaba, M. A., Krupa, J. C., Xu, C., Mort, J. S. & Min, J. (2011). Nature Commun. 2, 197. Google Scholar
Bouza, E. (2012). Clin. Microbiol. Infect. 18, Suppl. 6, 5–12. Google Scholar
Bradshaw, W. J., Kirby, J. M., Thiyagarajan, N., Chambers, C. J., Davies, A. H., Roberts, A. K., Shone, C. C. & Acharya, K. R. (2014). Acta Cryst. D70, 1983–1993. Web of Science CrossRef IUCr Journals Google Scholar
Calabi, E., Ward, S., Wren, B., Paxton, T., Panico, M., Morris, H., Dell, A., Dougan, G. & Fairweather, N. (2001). Mol. Microbiol. 40, 1187–1199. Web of Science CrossRef PubMed CAS Google Scholar
Cerquetti, M., Molinari, A., Sebastianelli, A., Diociaiuti, M., Petruzzelli, R., Capo, C. & Mastrantonio, P. (2000). Microb. Pathog. 28, 363–372. Web of Science CrossRef PubMed CAS Google Scholar
Chen, V. B., Arendall, W. B., Headd, J. J., Keedy, D. A., Immormino, R. M., Kapral, G. J., Murray, L. W., Richardson, J. S. & Richardson, D. C. (2010). Acta Cryst. D66, 12–21. Web of Science CrossRef CAS IUCr Journals Google Scholar
Coulombe, R., Grochulski, P., Sivaraman, J., Ménard, R., Mort, J. S. & Cygler, M. (1996). EMBO J. 15, 5492–5503. CAS PubMed Web of Science Google Scholar
Dang, T. H., de la Riva, L., Fagan, R. P., Storck, E. M., Heal, W. P., Janoir, C., Fairweather, N. F. & Tate, E. W. (2010). ACS Chem. Biol. 5, 279–285. Web of Science CrossRef CAS PubMed Google Scholar
Dubberke, E. R. & Olsen, M. A. (2012). Clin. Infect. Dis. 55, S88–S92. Web of Science CrossRef PubMed Google Scholar
Emsley, P. & Cowtan, K. (2004). Acta Cryst. D60, 2126–2132. Web of Science CrossRef CAS IUCr Journals Google Scholar
Evans, P. R. & Murshudov, G. N. (2013). Acta Cryst. D69, 1204–1214. Web of Science CrossRef CAS IUCr Journals Google Scholar
Fagan, R. P. & Fairweather, N. F. (2014). Nature Rev. Microbiol. 12, 211–222. Web of Science CrossRef CAS Google Scholar
Fagan, R. P., Janoir, C., Collignon, A., Mastrantonio, P., Poxton, I. R. & Fairweather, N. F. (2011). J. Med. Microbiol. 60, 1225–1228. Web of Science CrossRef CAS PubMed Google Scholar
Fox, T., de Miguel, E., Mort, J. S. & Storer, A. C. (1992). Biochemistry, 31, 12571–12576. CrossRef PubMed CAS Web of Science Google Scholar
Guex, N. & Peitsch, M. C. (1997). Electrophoresis, 18, 2714–2723. Web of Science CrossRef CAS PubMed Google Scholar
Janoir, C., Grénery, J., Savariau-Lacomme, M.-P. & Collignon, A. (2004). Pathol. Biol. 52, 444–449. Web of Science CrossRef PubMed CAS Google Scholar
Janoir, C., Péchiné, S., Grosdidier, C. & Collignon, A. (2007). J. Bacteriol. 189, 7174–7180. Web of Science CrossRef PubMed CAS Google Scholar
Kabsch, W. (2010). Acta Cryst. D66, 125–132. Web of Science CrossRef CAS IUCr Journals Google Scholar
Kachrimanidou, M. & Malisiovas, N. (2011). Crit. Rev. Microbiol. 37, 178–187. Web of Science CrossRef CAS PubMed Google Scholar
Kamphuis, I. G., Kalk, K. H., Swarte, M. B. & Drenth, J. (1984). J. Mol. Biol. 179, 233–256. CrossRef CAS PubMed Web of Science Google Scholar
Karjalainen, T., Waligora-Dupriet, A. J., Cerquetti, M., Spigaglia, P., Maggioni, A., Mauri, P. & Mastrantonio, P. (2001). Infect. Immun. 69, 3442–3446. Web of Science CrossRef PubMed CAS Google Scholar
Kawata, T., Takeoka, A., Takumi, K. & Masuda, K. (1984). FEMS Microbiol. Lett. 24, 323–328. CrossRef CAS Web of Science Google Scholar
Kirby, J. M., Ahern, H., Roberts, A. K., Kumar, V., Freeman, Z., Acharya, K. R. & Shone, C. C. (2009). J. Biol. Chem. 284, 34666–34673. Web of Science CrossRef PubMed CAS Google Scholar
McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C. & Read, R. J. (2007). J. Appl. Cryst. 40, 658–674. Web of Science CrossRef CAS IUCr Journals Google Scholar
Murshudov, G. N., Skubák, P., Lebedev, A. A., Pannu, N. S., Steiner, R. A., Nicholls, R. A., Winn, M. D., Long, F. & Vagin, A. A. (2011). Acta Cryst. D67, 355–367. Web of Science CrossRef CAS IUCr Journals Google Scholar
Musil, D., Zucic, D., Turk, D., Engh, R. A., Mayr, I., Huber, R., Popovic, T., Turk, V., Towatari, T., Katunuma, N. & Bode, W. (1991). EMBO J. 10, 2321–2330. PubMed CAS Web of Science Google Scholar
Podobnik, M., Kuhelj, R., Turk, V. & Turk, D. (1997). J. Mol. Biol. 271, 774–788. CrossRef CAS PubMed Web of Science Google Scholar
Qazi, O., Hitchen, P., Tissot, B., Panico, M., Morris, H. R., Dell, A. & Fairweather, N. (2009). J. Mass Spectrom. 44, 368–374. Web of Science CrossRef PubMed CAS Google Scholar
Riva, L. de la, Willing, S. E., Tate, E. W. & Fairweather, N. F. (2011). J. Bacteriol. 193, 3276–3285. Web of Science PubMed Google Scholar
Roy, S., Choudhury, D., Aich, P., Dattagupta, J. K. & Biswas, S. (2012). Acta Cryst. D68, 1591–1603. Web of Science CrossRef IUCr Journals Google Scholar
Rupnik, M., Wilcox, M. H. & Gerding, D. N. (2009). Nature Rev. Microbiol. 7, 526–536. Web of Science CrossRef CAS Google Scholar
Schaschke, N., Assfalg-Machleidt, I., Machleidt, W. & Moroder, L. (1998). FEBS Lett. 421, 80–82. Web of Science CrossRef CAS PubMed Google Scholar
Turk, D., Podobnik, M., Kuhelj, R., Dolinar, M. & Turk, V. (1996). FEBS Lett. 384, 211–214. CrossRef CAS PubMed Web of Science Google Scholar
Wiederanders, B. (2003). Acta Biochim. Pol. 50, 691–713. Web of Science PubMed CAS Google Scholar
Winter, G., Lobley, C. M. C. & Prince, S. M. (2013). Acta Cryst. D69, 1260–1273. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zheng, H., Chordia, M. D., Cooper, D. R., Chruszcz, M., Müller, P., Sheldrick, G. M. & Minor, W. (2014). Nature Protoc. 9, 156–170. Web of Science CrossRef CAS Google Scholar
Zheng, H., Chruszcz, M., Lasota, P., Lebioda, L. & Minor, W. (2008). J. Inorg. Biochem. 102, 1765–1776. Web of Science CrossRef PubMed CAS Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.