Qingpin Xiao1, Xiaohui Hu2, Zhiyi Wei3, Kin Yip Tam2. 1. 1. Faculty of Health Sciences, University of Macau, Taipa, Macau, China; 2. Department of Biology, Southern University of Science and Technology, Shenzhen, 518055, China. 2. 1. Faculty of Health Sciences, University of Macau, Taipa, Macau, China. 3. 2. Department of Biology, Southern University of Science and Technology, Shenzhen, 518055, China.
Abstract
Cells make use of molecular motors to transport small molecules, macromolecules and cellular organelles to target region to execute biological functions, which is utmost important for polarized cells, such as neurons. In particular, cytoskeleton motors play fundamental roles in neuron polarization, extension, shape and neurotransmission. Cytoskeleton motors comprise of myosin, kinesin and cytoplasmic dynein. F-actin filaments act as myosin track, while kinesin and cytoplasmic dynein move on microtubules. Cytoskeleton motors work together to build a highly polarized and regulated system in neuronal cells via different molecular mechanisms and functional regulations. This review discusses the structures and working mechanisms of the cytoskeleton motors in neurons.
Cells make use of molecular motors to transport small molecules, macromolecules and cellular organelles to target region to execute biological functions, which is utmost important for polarized cells, such as neurons. In particular, cytoskeleton motors play fundamental roles in neuron polarization, extension, shape and neurotransmission. Cytoskeleton motors comprise of myosin, kinesin and cytoplasmic dynein. F-actin filaments act as myosin track, while kinesin and cytoplasmic dynein move on microtubules. Cytoskeleton motors work together to build a highly polarized and regulated system in neuronal cells via different molecular mechanisms and functional regulations. This review discusses the structures and working mechanisms of the cytoskeleton motors in neurons.
A neuron cell is highly polarized, composing of a cell body, an axon and dendrites. Macromolecules, such as proteins and RNA, are synthesized in neuronal cell body, and are then transported to the axon or dendrites. Cytoskeleton motors are vital to transport cargos in the biological processes by using cytoskeleton as walking tracks and hydrolyze ATP to provide mechanical energy. Myosin is one superfamily of cytoskeletal motors, walking along the F-actin filaments. Kinesin and Dynein are both microtubule motors, walking towards microtubule's plus-end and minus-end directions respectively. The dysfunction of these motor proteins in neuron system is related to many humanneuronal diseases. In this work, we review the structural and functional characteristics of cytoskeleton motors and discuss the roles of cytoskeleton motors in neuronal polarization, development and neurotransmission processes. Finally, some possible future directions of cytoskeleton motors study will be highlighted.
Actin-based motors
According to their gene sequence homology, myosin can be classified into 35 classes in eukaryote 1. Mammals have 40 myosin genes, whose expression products fall into 13 classes 2. Figure 1 shows the selected myosin families, which have been reported and expressed in neuronal cells 3,4. Myosins compose of three common domains 5: (1) the N-terminal motor domains that interact with actin and hydrolyze ATP to generate force 6, (2) the neck domains or lever arms that transduce and amplify the force generated by motor domain 7,8 and (3) the tail domains offer specific binding sites for different cargos 9.
Figure 1
Schematic Diagrams of Myosin superfamily representations in neuron. Different myosins possess conserved N-terminal motor domains and diverse C-terminal tail regions. The neck domain binds to calmodulin or calmodulin-like light chains.
The motor domain of myosin has an actin-binding site and a nucleotide-binding site, which are mutually exclusive. Following the motor domain is lever arm, which usually contains some IQ motifs (consensus sequence IQxxxRGxxR) for calmodulin (CaM) or calmodulin-like proteins binding 10,11. The numbers of IQ motifs vary among different myosin classes. The length of lever arm determines the step size along actin. The most essential difference among myosin motor proteins lies in their C-terminal globular tail domain (GTD). This domain recognizes various cargos through direct interactions or via adaptor proteins. Usually, GTD binds to the adaptor protein, which can interact with organelle-specific Rab family GTPase protein and therefore link to the specific organelle or vesicle 12,13.
Microtubule-based motors
Kinesins (also referred to as KIFs) propel cargos along the microtubule. There are about 45 KIF genes in mammal genome 14. A standard kinesin nomenclature system had been developed for the kinesin family members based on phylogenic analysis 15. Accordingly, the kinesin proteins can be classified into 14 subfamilies, namely kinesin 1 to kinesin 14. Figure 2A gives representative kinesin of each family found in neuron 16. All kinesins contain a highly conserved globular motor domain comprising a microtubule-binding site and an ATP-binding site as well as a diverse tail domain, which is responsible for cargo recognition and binding via the adaptor or scaffolding proteins 17. Unlike myosin whose motor domain is always located in N-terminus, the motor domain of kinesins can be found in their N-terminal (N-type), middle (M-type), or C-terminal (C-type) region 18. On microtubule, N-type kinesin moves to the plus end of microtubule, C-type kinesin moves to minus end of microtubule and M-type kinesin can destabilize microtubule tracks 19. The energy supplying for kinesin processive motility comes from hydrolyzing ATP in the motor domain.
Figure 2
Schematic Diagrams of Microtubule-based motor proteins. (A) Kinesin superfamily representations in neuron. Like myosins, kinesins have conserved motor (head) domains and diverse tails. According to the positions of the motor domain, kinesins can be grouped into N-type, M-type and C-type kinesins. (B) Cytoplasmic dynein. The cytoplasmic dynein consists of C-terminal motor domain and N-terminal tail domain. The tail domains not only homodimerize the dynein heavy chains (purple color), but also recruit five other subunits (LIC, IC, LC8, LC7, TCTEX1) to accomplish dynein assembly.
Another microtubule-based motor is cytoplasmic dynein (Figure 2B). Compared to myosin and kinesin, dynein is larger in size with more complicated composition and its working mechanisms are not well understood. The motor domain of dynein is located in C-terminal region, while tail domain is situated in the N-terminal region. The motor domain composes of six ATPases associated with diverse cellular activities (AAA) domains. The six AAA domains arrange in a large doughnut-shape ring 20,21. AAA1 and AAA3 are responsible for ATP binding and hydrolysis, AAA2 and AAA4 play essential regulatory roles 22,23. The function of AAA5 and AAA6 is not clear, because they lacked the required ATP-binding residues and exhibited other functions 21. Extended from AAA4 domain is an antiparallel coiled-coil region, at the tip of which is the microtubule-binding domain. Microtubule-binding domain, coupled with ATP-hydrolysis cycle, moves the protein towards minus-end microtubule cytoskeleton. The dynein tail domain is responsible for dynein homodimerization and serves as binding sites with other subunits of dynein, including intermediate chain (IC), light intermediate chain (LIC), light chain 8 (LC8), roadblock (also known as LC7) and T-complex testis-specific protein 1 (TCTEX1). These dynein subunits can interact with adaptor proteins to form bridges 24-27. With these bridges, dynein tail domain can thus carry the target cargos.
Structural biology and cytoskeleton motors in neuron
A large number of structural studies have greatly enhanced our understanding of neurological functions of many cytoskeleton motors. We summarize a few recent cases as examples to highlight the power and uniqueness of structural biology in uncovering the molecular mechanisms of cytoskeleton motors.Many mutations in cytoskeleton motor proteins directly lead to diseases and aberrant developments. In neuron system, mutations found in myosin VIIa tail region result in the neuron disease named Usher syndrome I (USH1), which is characterized by deafness and visual impairment 28,29. Crystallographic study shows that myosin VIIa tail region MyTH4-FERM domain binds to Sans central domain (Figure 3A) 30. Sans is an adaptor protein interacting with both myosin VIIa and hormonin to form the ternary complex, which plays pivotal role in sensory neuroepithelia stereocilia formation 31. In the complex structure, numerous amino acid mutations documented in Usher syndromepatients, such as E1327K, A1340P, E1349K, are located in the MyTH4-FERM domain of myosin VIIa and affect the interaction between the motor and Sans, leading to the disruption of the hormonin targeting to the neuroepithelia stereocilia of the inner ear. Based on the current study, we can further explore how these mutations individually impact on these interactions. For example, some mutations may function to weaken the interaction, while others may involve in regulation of other partners or signal pathways.
Figure 3
Structural characteristics of cytoskeleton motors. The graphics of 3D views were drawn by PyMOL. (A) Complex structure of myosin VIIa/Sans. Cartoon representation indicates Sans-CEN1 domain binds to myosin VIIa-FERM domain which consists of F1, F2 and F3 lobes. (B) Structure of KIF13B-FHA domain bound to CENTA1. Cartoon diagram shows KIF13B-FHA domain interacts with dimerized CENTA1 PH1 and ArfGAP domains.
During axonal development, Phosphatidylinositol (3,4,5)-trisphosphate (PIP3)-containing vesicles function as the polarization factor, and need to be transported away from the soma to the tip of axon (more details will be given later). Kinesin III family KIF13B binds to PIP3 by using Centaurin-α1 (CENTA1) as adaptor protein to undertake the anterograde transport in vivo 32,33. Structure of KIF13B-CENTA1 complex had been solved in the presence of phosphatidylinositol (Figure 3B) 34. KIF13B binds to the first pleckstrin homology (PH1) domain and ArfGAP domain of CENTA1 through forkhead-associated (FHA) domain. According to the study, PH1 and PH2 domain of CENTA1 can bind simultaneously to PIP3 head group (IP4) 34. Because of homodimerization of full-length KIF13B, it is expected to form a supercomplex with CENTA1 and PIP3-containing vesicle and move along microtubule towards the axonal terminus. This complex structure not only built the first molecular model of kinesin-adaptor interaction, but confirmed the long-speculated phosphorylation-independent interaction between CENTA1 and KIF13B-FHA, which represents the only known phosphothreonine-specific recognition domain but is unique among kinesin FHA domain 35.Dynactin (dynein activator) is one of the most important dynein's adaptors, which is a very complicated molecule, comprising of 11 different subunits 36. Dynactin can bind to dynein N-terminal tail region and also microtubule 36,37. In neuron cells, the ternary complex, consisting of dynein, dynactin and microtubule, is necessary for dendritic polarization formation and axonal uniformed microtubule development 38. According to recently published Cryo-electron microscopy (EM) structures 39, dynactin consists of four parts: pointed end, filament, barbed end and shoulder (Figure 4A). Pointed end is composed of pointed-end proteins (p25, p27, p62) and actin-related protein (Arp11). Filament contains Arp1 and β-actin. Barbed end consists of actin capping proteins (CapZα and CapZβ). The shoulder contains p150Glued, p50 and p24. The pointed end interacts with cargo, while the filament binds with dynein tail domain and conserved motor adaptor Bicaudal D2 (BICD2) (Figure 4B), which can stabilize the interaction between dynein tail and dynactin by forming stable ternary complex. These complex structures by Cryo-EM showed us the assembly of dynactin and the detail of its interaction with dynein. Considering the facts that human dynein mutations are related with nervous diseases, such as sensory, motor and cognitive deficits 40, the structural elucidation of dynein and its adaptors may provide fundamental knowledge to help the rationalization of the prognosis of these diseases. Clearly, the above-mentioned work highlighted the power of cryo-EM in the elucidation of complicated macromolecule interaction involved in large cytoskeleton motors, such as dynein.
Figure 4
Cryo-EM Structure of dynactin and its complex with dynein tail. The graphics of 3D views were drawn by PyMOL. (A) Cryo-EM structure of dynactin. According to the molecular model, complicated dynactin is composed of four parts (shoulder, pointed end, filament and barbed end), each of these parts contains several subunits. (B) Complex structure of dynactin interacting with BICD2 and dynein tail. Structural architecture shows dynactin interacting with dynein tail region. The other subunits of dynactin can interact with other proteins or cargos, which facilitates the dynein-dependent cargo trafficking intracellular.
Roles in neuron
The processes of neuron development and maturation are very complicated. The polarizing factors and neurotrophin factors play essential roles in neuronal development and maturation processes. Notably, the neurotransmission is the most important function for neuron system. However, it remains ambiguous about what roles cytoskeleton motors play during these processes and what the structure-based mechanisms lie behind these roles. Here, we seek to summarize the current understanding about cytoskeleton motors trafficking in neuron function.
Polarizing-related factors transport
During the early stage of neuron development, polarization-related factors are essential for axon/dendrite determination and formation 41. As shown in Figure 5, many of these factors need to be transported to the target region and then affect the neuron polarization through a variety of signal pathways or mechanisms.
Figure 5
Diagram of cytoskeleton motors in neuron. During neuron growth and development processes, cytoskeleton molecular motors get involved in signal molecules and growth factors transport in neurons. When conducting neurotransmission, cytoskeleton motors are responsible for kinds of neurotransmitter receptors trafficking.
Transport of PIP3 by KIF13B
It has been well-established that PIP3 can stimulate self-generation by activating Rho GTPases, which can subsequently activate PI3K to produce more PIP3 42. When treated by PI3K inhibitors to disrupt this feedback circuit, the asymmetrical polarization of neuron cells was reduced dramatically. This observation established that PIP3 is critical for neuronal polarization. In cultured hippocampal neuron cells, PIP3 has been found to accumulate in the tip of the longest neurite, which developed into an axon 43. Later, KIF13B was found to bind to PIP3 by using Centaurin-α1 as adaptor protein and carried them away from the soma to the tip of axon 32,33. As mentioned before, molecular mechanism of KIF13B-mediated phosphoinositide trafficking has been revealed by crystallographic study 34. It has been reported that KIF13B can be phosphorylated by Par1, leading to a reduction in PIP3 transportation. However, Par1 can be phosphorylated and inactivated by Par3/Par6/atypical PKC (aPKC) 44. This regulation mechanism indicates that there may exist some unknown connections between cellular signaling (not limited to phosphorylation cascade) and the motor trafficking processes. These deserve further investigation.
Transport of shootin 1 by KIF20B
Shootin 1 was discovered through proteomic analysis of cultured hippocampal neurons 45. Shootin 1 showed fluctuating asymmetrical location in multiple neurites, but eventually, it anchored on a single neurite and accumulated in this neurite, which subsequently grew and developed into the axon 45. Shootin 1 overexpression disrupted its asymmetrical accumulation and resulted in the surplus axons formation. Conversely, downregulated shootin 1 expression led to the neuronal polarization inhibition. It has been demonstrated that shootin 1 could regulate the asymmetrical accumulation of PIP3, thus the PI3K activity in the growth cone, but the reverse was not true. Conclusively, shootin 1 and PI3K are involved in the cell polarity pathway, during which shootin 1 is located in the upstream of PI3K and affects asymmetrical localization of PI3K 45. During the axon determination, shootin 1 was observed to be actively transported by KIF20B from the soma cell body to the growth cone 46. Disrupting shootin 1 anterograde transport decreases shootin 1 mobilizing to the growth cone, affecting asymmetrical distribution of PI3K activity in the neurite. However, the molecular mechanism of shootin 1 interacting with KIF20B requires further investigation.
Transport of Par3/Par6/atypical PKC (aPKC) polarity complex by KIF3A
Partitioning defective 3 homolog (Par3), Partitioning defective 6 homolog (Par6) and atypical protein kinase C (aPKC) form a polarity complex in epithelial or neuronal cell development 47. In cultured hippocampal neurons, the axonal determination requires the localized activity of PI3K and aPKC 48. The abnormal expression of Par3 or Par6 disrupts the axon specification and neuronal polarity. The polarity complex uses Par6 to bind to the small GTPase Cdc42 and uses Par3 to interact with Rac-specific guanine nucleotide-exchange factors (Rac GEF). Therefore, Par3/Par6/aPKC mediates Cdc42 induced Rac activation, which can subsequently induce PI3K activity to produce PIP3 41,49. In rat hippocampal neurons, Par3 was found to accumulate at the tip of developing axon. Assays show that KIF3A can anterograde transport the polarity complex to the tip of axonal neurite through directly binding to Par3 50. The localization of the polarity complex is essential to neuron polarity. Inhibition of Par3 and KIF3A interaction or disruption of Par3/Par6/aPKC transport to the growing axon can abolish the neuronal polarity development.
Neurotrophin-related factors transport
Neurotropins are the family of proteins implicated to induce neuron survival, development and modulating neuron plasticity 51. Neurotrophin family contains four members: nerve growth factor (NGF), brain-derived neurotrophic factor (BDNF), neurotrophin-3 (NT-3) and neurotrophin-4 (NT-4) 52. After their formation, these neurotrophins interact with the plasma membrane neurotrophin receptors, Trk and p75NTR, to regulate the neuronal survival 53,54. Neurotrophin receptors bind with neurotrophin factors with different specificities 55: TrkA binds tendentiously to NGF, TrkB binds preferentially to BDNF and NT-4, and TrkC to NT-3; while p75NTR can bind to all neurotrophin factors and their precursors. The activated neurotrophin and receptor complex can promote neuron growth through regulating gene expression, only when the complex has been transported to neuron soma from the terminal of axon or presynapse (Figure 5). 56,57. It appeared that neurotrophin and its receptor are firstly endocytosed into cells and then transported to the cell nucleus periphery by interacting with cytoplasmic dynein in some unknown manner, which allows extracellular signals to be sent to the cell nucleus 58.
NGF/TrkA signaling endosome
NGF binding to TrkA extracellular domain leads to a conformation change in its cytoplasmic domain and autophosphorylation its tyrosine residues, which enhance TrkA kinase activity 59. With surrounding plasma membrane, the activated NGF/TrkA complex forms the signaling endosome, which is endocytosed or internalized into the cell and activates the downstream signaling cascades 60. NGF/TrkA signaling endosomes has been shown to require dynein motor protein for retrograde transport and for survival signaling propagation 61. Disrupting the dynein trafficking through overexpression of dynamitin (p50) suppressed the perinuclear MAPK signaling pathway, which is activated by NGF/Trk, resulting in neuron cell atrophy and death 62.
BDNF/TrkB signaling endosome
BDNF has been demonstrated to regulate dendritic arborization and outgrowth during neuronal development 63. BDNF binds to receptor TrkB, which also distributes on the axonal terminal membrane surface. Initiated by BDNF binding to TrkB, TrkB becomes phosphorylated, triggering internalization of BDNF/TrkB signaling endosome 58,64,65. The signaling endosomes then associate with dynein motor to undergo axonal retrograde transport to cell body to induce nuclear signal pathway 66. It has been shown that Snapin acts as a dynein adaptor, linking dynein and BDNF/TrkB complex and ensuring the signaling endosome axonal transport 67. Disrupting Snapin and dynein interaction leads to aberrant BDNF/Trk motility and decreased dendrites growth.
Neurotransmission-related vesicles transport
Neurotransmission is the most important function of neuron system, a process of systematic communication between two neuronal cells. During the transmission, the neurotransmitters in the presynapse are released into the synaptic cleft and bind to the specific receptors on the postsynapse membrane. In presynapse, the synaptic vesicles containing neurotransmitters need to undergo a long axonal anterograde transport to reach the presynapse 68. KIF1A and KIF1Bβ transport synaptic vesicles in axon, such as synaptotagmin, synaptophysin and Rab3A 68,69. Additionally, myosin V has also been shown to associate with synaptic vesicles. It has been demonstrated that on the surface of synaptic vesicles, the synaptobrevin-synaptophysin complex binds to myosin V through the cytoplasmic domain of synaptobrevin/synaptophysin to form a multimeric complex 70,71. In hippocampal neurons, myosin Va associates with the postsynaptic density-95 (PSD-95) through directly interacting with guanylate kinase domain-associated protein (GKAP), which binds to PSD-95 72,73. Simultaneously, myosin VI also binds to PSD-95 through its binding partner synapse-associated protein 97 (SAP97) interacting with GKAP 74,75. As a result, both myosin V and myosin VI can interact with the postsynaptic density protein PSD-95 and may work in concert to transport PSD-95-containing vesicles. In postsynapse (Figure 5), the receptors need to be transported to and from the postsynaptic membrane 76. Myosin VI has been reported to interact with alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptor (AMPAR), a glutamate receptor, through binding directly to synapse-associated protein 97 (SAP97) 77. Interestingly, KIF5 can also interact with AMPAR through adaptor glutamate receptor-interacting protein 1 (GRIP1) 78. Hence, KIF5 carries AMPAR in soma and dendrites along microtubule; while myosin VI transports AMPAR on F-actin. It implies that the motors can recognize the cargo through specific interactions with adaptor proteins. Also, the different kinds of motors may act in a coherent manner during the cargo transportation process 79,80. Besides AMPAR, both myosin VI and dynein can bind to muskelin that acts as an adaptor to link myosin VI and dynein to gamma-aminobutyric acid receptor (GABAR), and can regulate GABAR endocytosis process 80. During this process, muskelin acts as the interconnector between actin-based retrograde trafficking and microtubule-based retrograde trafficking of GABAR. Another transmitter's receptor, N-methyl-D-aspartate receptor (NMDAR), can be transported by KIF17 through bridging by the complex LIN10-LIN2-LIN7 81. But the interaction between KIF17 and LIN10 can be suppressed by septin 9 (SEPT9) which directly binds to the tail region of KIF17 and competes with LIN10 82.Some cytoskeleton motors, including myosin Ic, myosin III and myosin IX, play other roles as mechanoforce sensor, tether and signaling regulator. Myosin Ic directly binds to cell membranes through the C-terminal pleckstrin homology (PH) motif, which allows it to crosslink the membranes and actin filaments and transduce the mechanical force from tail-binding membranes to head-binding filaments in inner ear stereocilia 83. A recent structural study showed that myosin Ic undergoes CaM-mediated conformational transformation upon mechanical force transduction 84. Myosin III was also found to stimulate the elongation of stereocilia in hair cells through tethering and activating Espin1 at the tips of stereocilia, the plus ends of actin filaments 85. Gene mutations in encoding these proteins can lead to hear loss disease 86,87. The recent solved myosin III/Espin complex structure provided a mechanistic explanation of myosin III-mediated stimulation of actin bundling formation 88. Myosin IX contains a unique RhoGAP (GTPase-activating protein) domain in their tail region 89, which can switch active GTP-bound Rho to inactive GDP-bound Rho 90. The structural study of myosin IXbRhoGAP domain indicates a regulation role of myosin IXb in SLIT/ROBO pathway 91. In general, these diverse role manifests an integrated and complex function of myosin in cellular processes.
Perspectives
We have summarized and discussed the roles of cytoskeleton motors during neuronal development and transmission function. However, the structural basis for intracellular trafficking and regulations are not yet fully understood. Study on cytoskeleton motors' trafficking in neuron can help us to understand the cause of a variety of brain and nervous disorders. Although much progress has been made in past decades, structural study of neuronal-intracellular trafficking is still an attractive topic, especially in light of recent advance in macromolecule complex structure solvation technology, such as single particle cryo-EM 92.To gain further mechanistic insights, we envisage that the following three aspects of intracellular trafficking pertaining to cytoskeleton motors deserve special attentions:What are the regulations between cargo transport and trafficking network? It has been an attractive and complex question. To address this issue calls for the structural and functional studies on motors inactive/active transition, motors competitive/cooperative coordinating during trafficking processes and cargos loading/unloading regulations.How do cytoskeleton motors and the asymmetrical localization of differentiation factors play out in specific cell development? During neuron differentiation process, motors are required to determine the asymmetrical localization of many development relevant factors, which may be also indispensable for epithelial cells or other dendritic cells maturation process. Many of these motor proteins have the characteristic of tissue-specific alternative splicing (also known as isoforms), which implies that the tissue-specific motor proteins expression may be related to the specific cell type development.What's the potential therapeutic implication of cytoskeleton motors for manipulating the vesicle transport processes? It has been shown that many human pathogens, including bacterium, virus, fungus and prion, can spread into host cells by motor trafficking, and move to nucleus periphery of host cells by moving along cytoskeletons. Therefore, further studies to understand the mechanisms and regulations may offer ideas for developing novel antipathogen therapies.
Authors: H Wu; C Reissner; S Kuhlendahl; B Coblentz; S Reuver; S Kindler; E D Gundelfinger; C C Garner Journal: EMBO J Date: 2000-11-01 Impact factor: 11.598
Authors: Katrin Deinhardt; Sara Salinas; Carole Verastegui; Rose Watson; Daniel Worth; Sarah Hanrahan; Cecilia Bucci; Giampietro Schiavo Journal: Neuron Date: 2006-10-19 Impact factor: 17.173
Authors: Kalpana Mandal; Katarzyna Pogoda; Satabdi Nandi; Samuel Mathieu; Amal Kasri; Eric Klein; François Radvanyi; Bruno Goud; Paul A Janmey; Jean-Baptiste Manneville Journal: Nano Lett Date: 2019-10-07 Impact factor: 11.189
Authors: Christian Duellberg; Albert Auer; Nikola Canigova; Katrin Loibl; Martin Loose Journal: Proc Natl Acad Sci U S A Date: 2020-12-18 Impact factor: 11.205
Authors: Bárbara Matos; Stephen J Publicover; Luis Filipe C Castro; Pedro J Esteves; Margarida Fardilha Journal: Open Biol Date: 2021-06-02 Impact factor: 6.411