key: cord-0938698-nayr807b authors: Martí, Didac; Martín-Martínez, Eduard; Torras, Juan; Bertran, Oscar; Turon, Pau; Alemán, Carlos title: In silico antibody engineering for SARS-CoV-2 detection date: 2021-10-07 journal: Comput Struct Biotechnol J DOI: 10.1016/j.csbj.2021.10.010 sha: 108a20a8299e12d037ede6d645863b58bc77c532 doc_id: 938698 cord_uid: nayr807b Engineered immunoglobulin-G molecules (IgGs) are of wide interest for the development of detection elements in protein-based biosensors with clinical applications. The strategy usually employed for the de novo design of such engineered IgGs consists on merging fragments of the three-dimensional structure of a native IgG, which is immobilized on the biosensor surface, and of an antibody with an exquisite target specificity and affinity. In this work conventional and accelerated classical molecular dynamics (cMD and aMD, respectively) simulations have been used to propose two IgG-like antibodies for COVID-19 detection. More specifically, the crystal structure of the IgG1 B12 antibody, which inactivates the human immunodeficiency virus-1, has been merged with the structure of the antibody CR3022 Fab tightly bounded to SARS-CoV-2 receptor-binding domain (RBD) and the structure of the S309 antibody Fab fragment complexed with SARS-CoV-2 RBD. The two constructed antibodies, named IgG1-CR3022 and IgG1-S309, respectively, have been immobilized on a stable gold surface through a linker. Analyses of the influence of both the merging strategy and the substrate on the stability of the two constructs indicate that the IgG1-S309 antibody better preserves the neutralizing structure than the IgG1-CR3022 one. Overall, results indicate that the IgG1-S309 is appropriated for the generation of antibody based sensors for COVID-19 diagnosis. Immunoglobulin G (IgG) antibodies, which play a key role regulating the human immune system [1] , are amongst the most exquisitely designed and engineered molecules in Nature. Because of their exceptional bio-recognition elements, which exhibit high specificity and affinity for their cognate antigen, IgG antibodies are widely used in serological sensor devices (immunosensors) for detection of pathogens and toxins [2] [3] [4] [5] [6] [7] [8] . Within this context, recombinant technology, which easily allows antibodies genetic manipulation, is a valuable and robust tool for the fabrication of immunosensors. IgG antibodies present a monomeric "H2L2" structure, consisting of two identical heavy chains (H) and two identical light (L) chains. The molecular weight of H chains (50 kDa) is approximately twice that of L chains (25 kDa), the former ones being linked among them and to a L chain each by disulfide bonds. The quaternary structure is characterized by two identical halves that joint forming a Y-like shape (Scheme 1). Each H chain contains a variable domain (V H ) and three constant domains (CH 1 , CH 2 , CH 3 ), with an additional "hinge region" between CH 1 and CH 2 . Besides, L chains involve variable domain (V L ) and a constant domain (C L ). The light chain associates with the V H and CH 1 domains to form the fragment antigen binding (Fab) arm ("Fab" = fragment antigen binding). The lower hinge region and the CH 2 /CH 3 domains form the fragment crystalline (Fc). The fragment composed of Fab regions, joined by the hinge region, is known as F(ab') 2 . The performance of manufactured antibody-based sensors depends, among others, on the procedure used to immobilize antibody while maintaining its natural activity. Thus, although antibodies can be attached to the solid support in many different orientations, immobilization through the Fc region is crucial to leave the Fab exposed for recognition of the antigen, optimizing the sensitivity and the limit of detection of the sensor [9] . Within this context, different strategies have been reported to avoid random orientations, promoting the immobilization of the Fc-specific orientation (e.g. photon-assisted methods based on ultrashort UV laser pulses [10] , controlled electric fields [11] and surface functionalization [9, 12] ). Scheme 1. Parts of the Y-like shape IgG antibodies. Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has triggered a global health crisis with high social impact through the COVID-19 disease [13] . The spike glycoprotein protein, which consists on a trimer with three monomers exhibiting identical primary structure, plays a key role in viral infection and pathogenesis [14] . SARS-CoV-2 infection undergoes a series of processes. The binding of the receptor-binding domain (RBD) to its receptor, angiotensin converting enzyme 2 (ACE2), to form an RBD/ACE2 complex is the first [15, 16] . It triggers conformational changes in the spike protein, leading to membrane fusion mediated via others part of the spike [17, 18] . This process culminates in viral entry into target cells. Within this context, the development of rapid and efficient immunosensors for early detection in the diagnosis of SARS-CoV-2 is highly desirable. The design of neutralizing IgG-like antibodies that binds the RBD of the SARS-CoV-2 protein has been found to be a promising approach against COVID-19 [19] [20] [21] [22] . Indeed, a huge amount of SARS-CoV-2 antibodies, obtained from COVID-19 patients, have been disclosed since the beginning of the pandemic [23] . In this work we focus on two of them, CR3022 and S309. CR3022 does not overlap and not compete with the ACE2 binding site when binds SARS-CoV-2 [24] . Although there is some controversy on the in vivo neutralizing capacity of CR3022 against SARS-CoV-2, the antibody is able to bind to the RBD, conferring in vivo protection [25] [26] [27] . The binding of the CR3022 to the epitope on RBD of SARS-CoV-2 spike was recently studied using structural Molecular Dynamics (MD) simulations [28] . On the other hand, S309 monoclonal antibody was recently identified as a therapeutic agent that potently neutralizes SARS-CoV-2 by recognizing an epitope that contains a fucosylated glycan at position N343 [29] . S309 binds the RBD of SARS-CoV-2 spike in open conformation without competing with ACE2 binding. Currently, an engineered antibody (VIR-7831) based on S309 is under clinical trial assessment in humans (i.e. NCT04545060) [30, 31] . In this work, two IgG-like antibodies have been specifically engineered to detect SARS-CoV-2. After that, in silico studies have been conducted to examine their stability when immobilized on a solid gold surface, which is a usual substrate for virus immunosensors [32] [33] [34] [35] . For this purpose, the crystal structure of the neutralizing antibody CR3022 Fab tightly bounded to SARS-CoV-2 RBD (Code in the Protein Data Bank, PDB: 6YLA) [36] and the crystal structure of the S309 neutralizing antibody Fab fragment in complex with SARS-CoV-2 RBD (Code in the PDB: 6WPT) [37] were combined with the crystal structure of IgG1 B12 antibody (Code in PDB: 1HZH) [38] , which was used to inactivate the human immunodeficiency virus-1 (HIV-1). More specifically, the Fab fragment of CR3022 and the Fv fragment (i.e. the fragment formed by V H and V L , which are held together by non-covalent interactions; see Scheme 1) of S309 were used to replace the corresponding fragments of IgG1 B12, the resulting antibodies being denoted IgG1-CR3022 and IgG1-S309, respectively. It is worth noting that IgG1 is the most abundant IgG subclass in human serum, representing 65% of the total IgG and mediating antibody response against viral pathogens [39] . This molecular engineering process is summarized in Scheme 2. Scheme 2. Process used to engineer IgG1-CR3022 and IgG1-S309 antibodies. [40, 41] was built considering 4536 gold atoms per layer and the crystallographic parameter of gold, a = 2.89 Å. The immobilization of the antibodies on gold was performed according to a well-known strategy that considers the activation with carboxylate groups, functionalizing the surface (Scheme 3a) [42, 43] . Accordingly, the Lys478 residue of IgG1 B12 which is located at the end of the CH 3 heavy chain H of the F c region, was replaced by the covalent linker displayed in Scheme 3b. Both systems were solvated with 162597 explicit water molecules and Na + ions were introduced for charge neutralization, the total amount of particles in the simulated systems being 543221 and 543223 for IgG1-CR3022 and the IgG1-S309, respectively. Details about the construction of the gold surface, the immobilization of the antibodies and the solvation of the tethered systems are provided in the ESI. The constructed models were simulated using AMBER 18 simulation package [44] by applying the force-field parameters described in the ESI. After equilibration and thermalization using the protocols provided in the ESI, 150 ns of NVT conventional MD (cMD) were run for each system. Then, the conformational sampling of the immobilized antibodies was enhanced by performing accelerated MD (aMD) [45] , which artificially reduces energy barriers separating different states of the studied system. More specifically, aMD simulations were conducted along 90 ns using an NVT ensemble and starting from 4 snapshots for each studied system (i.e. the 4 initial structures were taken from cMD simulations at 90, 105, 120 and 140 ns). Accordingly, a total of 360 ns of aMD were produced for each tethered engineered antibody. the two tethered antibodies, which were calculated with respect to the initial models as constructed using the crystal structures of the CR3022 Fab, the S309 Fab and the IgG1 B12 antibody. The two studied systems were equilibrated after 40 ns of production. The averages of RMSD values in the last 110 ns of production take values of 13.5 ± 1.2 Å and 12.8 ± 1.3 Å for IgG1-CR3022 and IgG1-S309, respectively. In both cases, the standard deviation is below 1.5 Å. However, the RMSD values indicate significant distortions with respect to the initial models, which have been attributed to both the fact of combining different crystal structures in the engineered molecules and the restrictions caused by the tethering to the gold surface. Figure 2c , which represents the structure of IgG1-CR3022 and the IgG1-S309 in different colors as function of the RMSF values. As it can be seen, the larger distortions occur in the F c and Fab1 regions for IgG1-CR3022, while they are essentially located in Fab1 for IgG1-S309. The fact that the distortions of the SARS-CoV-2 neutralizing antibodies occur in the Fab1 fragment and such distortions are higher for IgG1-CR3022 than for IgG1-S309 are supported by analyses of the distance between the gold surface and the center of mass of each Fab arm, denoted d Fab1 and d Fab2 in Figure S3a . Inspection of the temporal evolution of d Fab1 along cMD trajectories ( Figure S3b ) indicates that the difference between the average value from the last 110 ns (57 ± 2 and 77 ± 5 Å for IgG1-CR3022 and IgG1-S309, respectively) and the initial value (70 and 72 Å, respectively) is significantly higher for IgG1-CR3022 than for IgG1-S309. Conversely, average d Fab2 values (91 ± 2 and 90 ± 4 Å, respectively) are very similar to the initial values (87 and 92 Å, respectively) for the two antibodies ( Figure S3c ). The effect of the gold surface on the orientation of the engineered antibodies was analyzed considering both cMD and aMD simulations. Firstly, the tilt angle of the engineered antibodies with respect to the constructed surface was examined using the α angle (Figure 3a ), which was defined by a gold atom of the slab at the top, the -carbon atom of the residue used to tether the antibody to the surface (i.e. C  in Scheme 3b), and the -carbon atom of the residue at the top of the F c region (i.e. Pro233 and Pro241 for IgG1-CR3022 and the IgG1-S309, respectively). The evolution of α along the cMD simulation, which is represented in Figure 3b , shows that stabilization is reached after 120 and 95 ns for IgG1-CR3022 and the IgG1-S309, respectively. The tilt angle of the initially constructed models (i.e. 72º and 93º for IgG1-CR3022 and the IgG1-S309, respectively) decreases to α = 54º ± 1º and 50º ± 1º (averaged over the last 30 and 55 ns, respectively). The potential of mean force (PMF) free energy profiles for the tilt angle calculated from aMD samplings are displayed in Figure 3c -d, which only show, for the sake of clarity, the region with G  4.2 kcal/mol. As it can be seen, a significant number of local minima are detected within a G interval of 1.5 kcal/mol for the two engineered antibodies. The global minimum for IgG1-CR3022 appears at α = 54.1 º (green diamond in Figure 3c )), which is very close to the average value obtained by cMD once the trajectory was stabilized (light blue line in Figure 3c ). However, several local minima destabilized by less than 0.1 kcal/mol with respect to the global minimum (i.e. isoenergetic minima) appears at α = 49.5º and 50.7º, suggesting a tilting of around 40º with respect to the ideal orientation (α = 90º). However, due the flatness of the profile in a range of around 10º (i.e. the G is lower than 1.0 kcal/mol for α values comprised between 48º and 57º), this value can be considered as a rough approximation. On the other hand, the PMF profile obtained for IgG1-S309 (Figure 3d ) shows similar trends, with a global minimum at α = 49.4 º (green diamond) and several local minima with G < 0.1 kcal/mol surrounding it. These results indicate that the tilt angle of anchored IgG1-CR3022 and IgG1-S309 are roughly estimated at around 50º-55º, which are quite similar to that obtained for IgG1 protein at room temperature (α= 66º) [46] . The flexibility and distortion of the Fab1 and Fab2 arms with respect to the Fc region were determined by measuring the  and  angles (Figure 4a ), respectively. These were defined by the C  atom of the residue used to tether the antibody to the surface (Scheme respectively, for IgG1-S309, corroborating the cMD results. An important structural feature that deserves discussion is the asymmetry found in the orientation of the two arms, which was defined by the difference between the two angles ( =  -). This asymmetry is significantly higher for IgG1-CR3022 than for IgG1-S309. Thus, cMD simulations predicted  = 53º and 8º for IgG1-CR3022 and IgG1-S309, respectively, similar values being also obtained from aMD when the corresponding global minimum are compared (i.e.  = 49º and 10º for IgG1-CR3022 and IgG1-S309, respectively). However, X-ray crystallography [47] and cryo-electron tomography studies in solution [48] reported that the arrangement of Fab arms relative to the Fc is not symmetrical in IgG1 molecules, which was attributed to their flexibility and dynamics. Indeed, literature shows a great variability in Fab-Fc angles [49] . Detailed analyses of cryo-electron tomography experiments on individual IgG molecules allowed determining the probability distribution of Fab-Fc angles at equilibrium [50] . Results showed that the Fab-Fc angles are approximately uniformly distributed within the two limiting values, 15° and 128°, which reflect the excluded steric constraints imposed by the molecular structure. However, analysis of IgG1 antibodies by individual particle electron tomography suggested that the orientation of the two Fab arms with respect to the Fc tend to be symmetric with  values close to that observed for IgG1-S309 [51] . An important aspect that deserves deeper analysis is the correlation among ,  and  angles, which has been analyzed considering the results from aMD simulations. Figure 5 shows Conversely, the Fab fragments are very similar in all cases, independently of the antibody or the arm orientation (i.e. 24.4-24.5 Å). PMF profiles determined for F c and Fab fragments from aMD simulations, which are displayed in Figure 6d -f, corroborate the structural features obtained from cMD. Additionally, these results provide some interesting findings that deserve consideration. Firstly, the potential well obtained for the Fc of IgG1-S309 is not only shifted with respect to that of IgG1-CR3022 but also is narrower (Figure 6d) . Similarly, the well found in the PMF of the Fab2 for IgG1-S309 is significantly narrower than those obtained for the Fab1 of the same engineered antibody and for the two Fab arms of IgG1-CR3022. These differences, which have been ascribed to the primary structures of the CR3022 and S309 fragments, explain the higher compactness of IgG1-CR3022. On the other hand, Figure S4 displays illustrative 2D PMF maps representing the R g of the whole antibody, and of the F c or the Fab fragments against the ,  and  angles. As it can be seen, the R g values are more restricted for IgG1-CR3022 than for IgG1-S309. Furthermore, the restrictions of the former are practically independently of the ,  and  angles. Other important parameters are those relative to the orientation of the Fab1 and Fab2 branches with respect to the surface, which are described by the φ and θ dihedral angles. Those dihedrals are defined by a gold atom of the superficial slab, the C  atom of the residue used to tether the antibody to the surface (i.e. the residue replacing Lys478 in IgG1 B12), the C  atom of the residue on the top of the Fc region (i.e. Pro233 and Pro241 for IgG1-CR3022 and the IgG1-S309, respectively), and the center of masses of the CDR of Fab1 (φ) or Fab2 (θ), as is schematically sketched in Figure 7a . Analyses of the cMD trajectories indicate that both φ and θ stabilize after 100 ns (Figure 7b-c) . The average values, which were obtained considering the last 50 ns of cMD runs were φ / θ = 80º  6º / -160º  11º and 108º  4º / -143º  2º for IgG1-CR3022 and the IgG1-S309, respectively. This result suggests that the relative disposition between the two arms, which has been defined by  = θ -φ, only exhibits a small dependence on the antibody structure (i.e.  = 120º and 109º, respectively). PMF plots representing the variation of the relative free energy against φ or θ are depicted in Figure 7d -g. Although the position of the global minimum found by aMD is consistent with the average values derived from cMD, the deviation between aMD and cMD is higher than those described for ,  and . However, the regions around the global minimum are very flat and, therefore, such deviations must be considered with caution. In any case, deviations have a relatively small effect on  that increases to 128º and 110º for IgG1-CR3022 and the IgG1-S309, respectively. Finally, structural changes induced by the anchoring of the two designed antibodies to the gold surface were evaluated by examining the secondary structure of the specific binding regions located at the Fc, Fab1 and Fab2. The secondary structure was determined using the DSSP program [54] , which applies a pattern recognition process of hydrogen bond and geometrical features to assign secondary structure to the residues of a protein. The DSSP algorithm was employed considering the snapshots stored during the first 10 ns of the cMD trajectories that, after averaging, allowed to define the initial secondary structure once the conformational tensions associated to the molecular engineering process are eliminated. Besides, the converged secondary structure, which includes changes induced by the immobilization on the gold surface, was obtained by averaging the results derived from the application of the DSSP to the snapshots stored during last 75 ns of the runs. Figure S5 displays the averaged initial and converged secondary structures of the six specific binding regions identified for the Fc fragment of IgG1-CR3022 and IgG1-S309, which have been labeled as Fc_# (where # ranges from 1 to 6). The sequence (primary structure) associated to each of the six Fc_# for each antibody, which is included in Figure S5 , has been labeled using the residue numberings displayed in Figures S1-S2. It is worth mentioning that only structural motifs with average populations higher than 20% have been considered as representative and, therefore, are the only included in Figure S5 . As it can be seen, the secondary structure of the six Fc_# domains was dominated by random, bend, -helix, turn and antiparallel  sheet motifs for the two antibodies. Comparison between the initial and converged secondary structures indicates major conformational rearrangements (i.e. those that imply changes in the populations of the different motifs that add up to more than 30%) for the Fc_2, Fc_3, Fc_5 and Fc_6 domains of IgG1-CR3022 and for the Fc_2, Fc_3 and Fc_6 domains of IgG1-S309. Finally, the secondary structures of the six CDRs of the Fab2 arm, which are labelled as Fab2_# (where # ranges from 1 to 6) in Figure S7 , remained practically unchanged during the MD trajectory. Thus, the only domain that experienced conformational changes involving an accumulated population variation higher than 30% was the Fab2_3 of IgG1-CR3022 ( Figure S7 ). However, detailed analysis of the populations of all identified secondary motifs, which are listed in Figure 8 , indicates that such rearrangements are less important than those observed for Fc and Fab1 domains. More specifically, in the case of Fab2_3, the only relevant conformational change is the transformation of turn and, in less extension, random structures into bend motifs. In this work we presented classical MD simulations of two antibodies, which were designed to detect SARS-CoV-2, immobilized on a gold (111) surface through a covalent linker. The IgG1-CR3022 and IgG1-S309 antibodies were engineered by replacing fragments of the IgG1 B12 antibody by other fragments of CR3022 and S309, respectively. Those fragments were found to bind the RBD of SARS-CoV-2, playing a key role in neutralizing the coronavirus. The choice of the gold surface was based on its frequent utilization as solid support in immunosensors. Previous works showed that the performance of manufactured antibody-based sensors strongly depends on the influence of the solid support and the orientation and structural stability of the immobilized antibody. Accordingly, in this work cMD and aMD simulations were performed to examine and compare the stability of the engineered antibodies when tethered on a solid gold surface. Although cMD simulations showed that, in average, the tilting was 4º higher for IgG1-S309 than for IgG1-CR3022, the rest of the analyzed structural and geometric parameters underwent higher distortions for the latter than for the former. Consistently, the calculated PMF profiles evidenced, not only that the positions of the minima obtained by aMD were very close to those obtained by averaging cMD values, but also that the potential wells were in general narrower for IgG1-S309 than for IgG1-CR3022. These observations are in agreement with the fact that the RMSD and the RMSF were smaller for IgG1-S309 than for IgG1-CR3022 and with the CDR secondary structure analyses at the beginning and at the end of the simulations. Those analyses showed that changes were more pronounced for the latter than for the former. Overall, these results indicate that between the two engineered antibodies, the IgG1-S309 is the most stable from a structural point of view, preserving the conformation of the neutralizing S309 antibody, and the least affected by the tethering to the gold surface. The sensitivity and specificity of immunosensors largely depend on the structure of the antibody fragments used to recognize the target pathogen. The structure of such fragments should not be affected by the molecular engineering process used to design the proposed antibody or by the substrate used for the antibody immobilization. Results derived from this study allow to conclude that, although the IgG1-CR3022 and IgG1-S309 antibodies were engineered using identical molecular recognition principles, the latter is more appropriated for the generation of antibody based sensors for COVID-19 diagnosis. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. IgG subclasses and allotypes: from structure to effector functions AuNPs@Fe 3 O 4 and CMK-8 for IgG anti-Toxocara canis determination Electrochemical immunosensor based on nanoelectrode ensembles for the serological analysis of IgG-type tissue transglutaminase Electrochemical Sensors and Biosensors Based on Nanomaterials and Nanostructures Recent trends in antibody based sensors Current aspects in immunosensors Antibody production, design and use for biosensor-based applications Antibody orientation on biosensor surfaces: a minireview Light assisted antibody immobilization for bio-sensing Tunable control of antibody immobilization using electric field Surface density dependent orientation and immunological recognition of antibody on silicon: TOF-SIMS and surface analysis of two covalent immobilization methods Subunit vaccines against emerging pathogenic human coronaviruses A pneumonia outbreak associated with a new coronavirus of probable bat origin Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor The spike protein of SARS-CoV -a target for vaccine and therapeutic development MERS-CoV spike protein: a key target for antivirals Bispecific IgG neutralizes SARS-CoV-2 variants and prevents escape in mice Convergent antibody responses to SARS-CoV-2 in convalescent individuals SARS-CoV-2 specific antibody and neutralization assays reveal the wide range of the humoral immune response to virus A human monoclonal antibody blocking SARS-CoV-2 infection Humanized single domain antibodies neutralize SARS-CoV-2 by targeting the spike receptor binding domain Neutralization of SARS-CoV-2 by Destruction of the Prefusion Spike A highly conserved cryptic epitope in the receptor binding domains of SARS-CoV-2 and SARS-CoV A human neutralizing antibody targets the receptor-binding site of SARS-CoV-2 Mechanistic Insights to the Binding of Antibody CR3022 Against RBD from SARS-CoV and HCoV-19/SARS-CoV-2: A Computational Study Crossneutralization of SARS-CoV-2 by a human monoclonal SARS-CoV antibody Anti-SARS-CoV-2 neutralizing monoclonal antibodies: clinical pipeline Ultrapotent human antibodies protect against SARS-CoV-2 challenge via multiple mechanisms An electrochemical immunosensor for the corona virus associated with the Middle East respiratory syndrome using an array of gold nanoparticle-modified carbon electrodes Immunosensor-based rapid quantitative detection of Newcastle disease virus antibodies using innovative gold immunochromatographic assay Label-free electrochemical immunosensor for direct detection of Citrus tristeza virus using modified gold electrode Gold nanoparticles and hepatitis B virus Neutralization of SARS-CoV-2 by destruction of the prefusion spike Crossneutralization of SARS-CoV-2 by a human monoclonal SARS-CoV antibody Crystal structure of a neutralizing human IGG against HIV-1: a template for vaccine design IgG subclasses and allotypes: from structure to effector functions Equilibrium morphology of face-centered cubic gold nanoparticles >3 nm and the shape changes induced by temperature Shape control in gold nanoparticles by N-containing ligands: Insights from density functional theory and Wulff constructions Conjugation of antibodies to gold nanorods through Fc portion: Synthesis and molecular specific imaging Design, preparation, and evaluation of a fixed-orientation antibody/gold-nanoparticle conjugate as an immunosensing label Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules Tethering of the IgG1 antibody to amorphous silica for immunosensor development: A molecular dynamics study Crystallographic structure of an intact IgG1 monoclonal antibody Structure and flexibility of individual immunoglobulin G molecules in solution Immunoglobulin structure and function as revealed by electron microscopy Freezing immunoglobulins to see them move 3D structural fluctuation of IgG1 antibody revealed by individual particle electron tomography In-depth analysis of subclass-specific conformational preferences of IgG antibodies Influence of the bispecific antibody IgG subclass on T cell redirection Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features Didac Martí: Investigation, Formal analysis, Methodology, Software, Validation, Visualization. Eduard Martín-Martínez: Investigation, Formal analysis, Methodology, Software, Validation, Visualization. Juan Torras: Conceptualization, Investigation, Supervision, Funding acquisition, Methodology, Validation, Writing -review & editing Funding acquisition, Writing -original draft, Writing -review & editing. Carlos Alemán: Conceptualization, Formal analysis, Funding acquisition, Visualization, Writing -original draft ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work