Computationally guided AAV engineering for enhanced gene delivery.

Gene delivery vehicles based on adeno-associated viruses (AAVs) are enabling increasing success in human clinical trials, and they offer the promise of treating a broad spectrum of both genetic and non-genetic disorders. However, delivery efficiency and targeting must be improved to enable safe and effective therapies. In recent years, considerable effort has been invested in creating AAV variants with improved delivery, and computational approaches have been increasingly harnessed for AAV engineering. In this review, we discuss how computationally designed AAV libraries are enabling directed evolution. Specifically, we highlight approaches that harness sequences outputted by next-generation sequencing (NGS) coupled with machine learning (ML) to generate new functional AAV capsids and related regulatory elements, pushing the frontier of what vector engineering and gene therapy may achieve.

Ancestral ribonucleases back in motion for evolutionary-dynamics guided protein design.

The dynamics behavior of a protein is essential for its functionality. Here, Doucet et al. demonstrate how the evolutionary analysis of conformational pathways within a protein family serves to identify common core scaffolds that accommodate branch-specific functional regions controlled by flexibility switches, offering a model for evolutionary-dynamics based protein design.

Metamorphic protein folding as evolutionary adaptation.

Metamorphic proteins switch reversibly between multiple distinct, stable structures, often with different functions. It was previously hypothesized that metamorphic proteins arose as intermediates in the evolution of a new fold - rare and transient exceptions to the 'one sequence, one fold' paradigm. However, as described herein, mounting evidence suggests that metamorphic folding is an adaptive feature, preserved and optimized over evolutionary time as exemplified by the NusG family and the chemokine XCL1. Analysis of extant protein families and resurrected protein ancestors demonstrates that large regions of sequence space are compatible with metamorphic folding. As a category that enhances biological fitness, metamorphic proteins are likely to employ fold switching to perform important biological functions and may be more common than previously thought.

Evolution and diversification of carboxylesterase-like [4+2] cyclases in aspidosperma and iboga alkaloid biosynthesis.

Monoterpene indole alkaloids (MIAs) are a large and diverse class of plant natural products, and their biosynthetic construction has been a subject of intensive study for many years. The enzymatic basis for the production of aspidosperma and iboga alkaloids, which are produced exclusively by members of the Apocynaceae plant family, has recently been discovered. Three carboxylesterase (CXE)-like enzymes from Catharanthus roseus and Tabernanthe iboga catalyze regio- and enantiodivergent [4+2] cycloaddition reactions to generate the aspidosperma (tabersonine synthase, TS) and iboga (coronaridine synthase, CorS; catharanthine synthase, CS) scaffolds from a common biosynthetic intermediate. Here, we use a combined phylogenetic and biochemical approach to investigate the evolution and functional diversification of these cyclase enzymes. Through ancestral sequence reconstruction, we provide evidence for initial evolution of TS from an ancestral CXE followed by emergence of CorS in two separate lineages, leading in turn to CS exclusively in the Catharanthus genus. This progression from aspidosperma to iboga alkaloid biosynthesis is consistent with the chemotaxonomic distribution of these MIAs. We subsequently generate and test a panel of chimeras based on the ancestral cyclases to probe the molecular basis for differential cyclization activity. Finally, we show through partial heterologous reconstitution of tabersonine biosynthesis using non-pathway enzymes how aspidosperma alkaloids could have first appeared as "underground metabolites" via recruitment of promiscuous enzymes from common protein families. Our results provide insight into the evolution of biosynthetic enzymes and how new secondary metabolic pathways can emerge through small but important sequence changes following co-option of preexisting enzymatic functions.

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction.

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes' innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Ancestral sequence reconstruction dissects structural and functional differences among eosinophil ribonucleases.

Evolutionarily conserved structural folds can give rise to diverse biological functions, yet predicting atomic-scale interactions that contribute to the emergence of novel activities within such folds remains challenging. Pancreatic-type ribonucleases illustrate this complexity, sharing a core structure that has evolved to accommodate varied functions. In this study, we used ancestral sequence reconstruction to probe evolutionary and molecular determinants that distinguish biological activities within eosinophil members of the RNase 2/3 subfamily. Our investigation unveils functional, structural, and dynamical behaviors that differentiate the evolved ancestral ribonuclease (AncRNase) from its contemporary eosinophil RNase orthologs. Leveraging the potential of ancestral reconstruction for protein engineering, we used AncRNase predictions to design a minimal 4-residue variant that transforms human RNase 2 into a chimeric enzyme endowed with the antimicrobial and cytotoxic activities of RNase 3 members. This work provides unique insights into mutational and evolutionary pathways governing structure, function, and conformational states within the eosinophil RNase subfamily, offering potential for targeted modulation of RNase-associated functions.

Emergence of an Orphan Nitrogenase Protein Following Atmospheric Oxygenation.

Molecular innovations within key metabolisms can have profound impacts on element cycling and ecological distribution. Yet, much of the molecular foundations of early evolved enzymes and metabolisms are unknown. Here, we bring one such mystery to relief by probing the birth and evolution of the G-subunit protein, an integral component of certain members of the nitrogenase family, the only enzymes capable of biological nitrogen fixation. The G-subunit is a Paleoproterozoic-age orphan protein that appears more than 1 billion years after the origin of nitrogenases. We show that the G-subunit arose with novel nitrogenase metal dependence and the ecological expansion of nitrogen-fixing microbes following the transition in environmental metal availabilities and atmospheric oxygenation that began ∼2.5 billion years ago. We identify molecular features that suggest early G-subunit proteins mediated cofactor or protein interactions required for novel metal dependency, priming ancient nitrogenases and their hosts to exploit these newly diversified geochemical environments. We further examined the degree of functional specialization in G-subunit evolution with extant and ancestral homologs using laboratory reconstruction experiments. Our results indicate that permanent recruitment of the orphan protein depended on the prior establishment of conserved molecular features and showcase how contingent evolutionary novelties might shape ecologically important microbial innovations.

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction.

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels.

Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.

Ancient whale rhodopsin reconstructs dim-light vision over a major evolutionary transition: Implications for ancestral diving behavior.

Cetaceans are fully aquatic mammals that descended from terrestrial ancestors, an iconic evolutionary transition characterized by adaptations for underwater foraging via breath-hold diving. Although the evolutionary history of this specialized behavior is challenging to reconstruct, coevolving sensory systems may offer valuable clues. The dim-light visual pigment, rhodopsin, which initiates phototransduction in the rod photoreceptors of the eye, has provided insight into the visual ecology of depth in several aquatic vertebrate lineages. Here, we use ancestral sequence reconstruction and protein resurrection experiments to quantify light-activation metrics in rhodopsin pigments from ancestors bracketing the cetacean terrestrial-to-aquatic transition. By comparing multiple reconstruction methods on a broadly sampled cetartiodactyl species tree, we generated highly robust ancestral sequence estimates. Our experimental results provide direct support for a blue-shift in spectral sensitivity along the branch separating cetaceans from terrestrial relatives. This blue-shift was 14 nm, resulting in a deep-sea signature (λmax = 486 nm) similar to many mesopelagic-dwelling fish. We also discovered that the decay rates of light-activated rhodopsin increased in ancestral cetaceans, which may indicate an accelerated dark adaptation response typical of deeper-diving mammals. Because slow decay rates are thought to help sequester cytotoxic photoproducts, this surprising result could reflect an ecological trade-off between rod photoprotection and dark adaptation. Taken together, these ancestral shifts in rhodopsin function suggest that some of the first fully aquatic cetaceans could dive into the mesopelagic zone (>200 m). Moreover, our reconstructions indicate that this behavior arose before the divergence of toothed and baleen whales.

Template switching in DNA replication can create and maintain RNA hairpins.

The evolutionary origin of RNA stem structures and the preservation of their base pairing under a spontaneous and random mutation process have puzzled theoretical evolutionary biologists. DNA replication-related template switching is a mutation mechanism that creates reverse-complement copies of sequence regions within a genome by replicating briefly along either the complementary or nascent DNA strand. Depending on the relative positions and context of the four switch points, this process may produce a reverse-complement repeat capable of forming the stem of a perfect DNA hairpin or fix the base pairing of an existing stem. Template switching is typically thought to trigger large structural changes, and its possible role in the origin and evolution of RNA genes has not been studied. Here, we show that the reconstructed ancestral histories of RNA genes contain mutation patterns consistent with the DNA replication-related template switching. In addition to multibase compensatory mutations, the mechanism can explain complex sequence changes, although mutations breaking the structure rarely get fixed in evolution. Our results suggest a solution for the long-standing dilemma of RNA gene evolution and demonstrate how template switching can both create perfect stems with a single mutation event and help maintaining the stem structure over time. Interestingly, template switching also provides an elegant explanation for the asymmetric base pair frequencies within RNA stems.

Vanillyl alcohol oxidase from Diplodia corticola: Residues Ala420 and Glu466 allow for efficient catalysis of syringyl derivatives.

Vanillyl alcohol oxidases (VAOs) belong to the 4-phenol oxidases family and are found predominantly in lignin-degrading ascomycetes. Systematical investigation of the enzyme family at the sequence level resulted in discovery and characterization of the second recombinantly produced VAO member, DcVAO, from Diplodia corticola. Remarkably high activities for 2,6-substituted substrates like 4-allyl-2,6-dimethoxy-phenol (3.5 ± 0.02 U mg-1) or 4-(hydroxymethyl)-2,6-dimethoxyphenol (6.3 ± 0.5 U mg-1) were observed, which could be attributed to a Phe to Ala exchange in the catalytic center. In order to rationalize this rare substrate preference among VAOs, we resurrected and characterized three ancestral enzymes and performed mutagenesis analyses. The results indicate that a Cys/Glu exchange was required to retain activity for É£-hydroxylations and shifted the acceptance towards benzyl ethers (up to 4.0 ± 0.1 U mg-1). Our findings contribute to the understanding of the functionality of VAO enzyme group, and with DcVAO, we add a new enzyme to the repertoire of ether cleaving biocatalysts.

Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation.

Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.

Enhancing PET Degrading Enzymes: A Combinatory Approach.

Plastic waste has become a substantial environmental issue. A potential strategy to mitigate this problem is to use enzymatic hydrolysis of plastics to depolymerize post-consumer waste and allow it to be reused. Over the last few decades, the use of enzymatic PET-degrading enzymes has shown promise as a great solution for creating a circular plastic waste economy. PsPETase from Piscinibacter sakaiensis has been identified as an enzyme with tremendous potential for such applications. But to improve its efficiency, enzyme engineering has been applied aiming at enhancing its thermal stability, enzymatic activity, and ease of production. Here, we combine different strategies such as structure-based rational design, ancestral sequence reconstruction and machine learning to engineer a more highly active Combi-PETase variant with a melting temperature of 70 °C and optimal performance at 60 °C. Furthermore, this study demonstrates that these approaches, commonly used in other works of enzyme engineering, are most effective when utilized in combination, enabling the improvement of enzymes for industrial applications.

Engineering the thermal stability of a polyphosphate kinase by ancestral sequence reconstruction to expand the temperature boundary for an industrially applicable ATP regeneration system.

ATP-dependent energy-consuming enzymatic reactions are widely used in cell-free biocatalysis. However, the direct addition of large amounts of expensive ATP can greatly increase cost, and enzymatic production is often difficult to achieve as a result. Although a polyphosphate kinase (PPK)-polyphosphate-based ATP regeneration system has the potential to solve this challenge, the generally poor thermal stability of PPKs limits the widespread use of this method. In this paper, we evaluated the thermal stability of a PPK from Sulfurovum lithotrophicum (SlPPK2). After directed evolution and computation-supported design, we found that SlPPK2 is very recalcitrant and cannot acquire beneficial mutations. Inspired by the usually outstanding stability of ancestral enzymes, we reconstructed the ancestral sequence of the PPK family and used it as a guide to construct three heat-stable variants of SlPPK2, of which the L35F/T144S variant has a half-life of more than 14 h at 60°C. Molecular dynamics simulations were performed on all enzymes to analyze the reasons for the increased thermal stability. The results showed that mutations at these two positions act synergistically from the interior and surface of the protein, leading to a more compact structure. Finally, the robustness of the L35F/T144S variant was verified in the synthesis of nucleotides at high temperature. In practice, the use of this high-temperature ATP regeneration system can effectively avoid byproduct accumulation. Our work extends the temperature boundary of ATP regeneration and has great potential for industrial applications.IMPORTANCEATP regeneration is an important basic applied study in the field of cell-free biocatalysis. Polyphosphate kinase (PPK) is an enzyme tool widely used for energy regeneration during enzymatic reactions. However, the thermal stability of the PPKs reported to date that can efficiently regenerate ATP is usually poor, which greatly limits their application. In this study, the thermal stability of a difficult-to-engineer PPK from Sulfurovum lithotrophicum was improved, guided by an ancestral sequence reconstruction strategy. The optimal variant has a 4.5-fold longer half-life at 60°C than the wild-type enzyme, thus enabling the extension of the temperature boundary for ATP regeneration. The ability of this variant to regenerate ATP was well demonstrated during high-temperature enzymatic production of nucleotides.

Improving the Quality of Wheat Flour Bread by a Thermophilic Xylanase with Ultra Activity and Stability Reconstructed by Ancestral Sequence and Computational-Aided Analysis.

Xylanase is an essential component used to hydrolyze the xylan in wheat flour to enhance the quality of bread. Presently, cold-activated xylanase is popularly utilized to aid in the development of dough. In this study, ancestral sequence reconstruction and molecular docking of xylanase and wheat xylan were used to enhance the activity and stability of a thermophilic xylanase. The results indicated that the ancestral enzyme TmxN3 exhibited significantly improved activity and thermal stability. The Vmax increased by 2.7 times, and the catalytic efficiency (Kcat/Km) increased by 1.7 times in comparison to TmxB. After being incubated at 100 °C for 120 min, it still retained 87.3% of its activity, and the half-life in 100 °C was 330 min, while the wild type xylanase was only 55 min. This resulted in an improved shelf life of bread, while adding TmxN3 considerably enhanced its quality with excellent volume and reduced hardness, chewiness, and gumminess. The results showed that the hardness was reduced by 55.2%, the chewiness was reduced by 40.11%, and the gumminess was reduced by 53.52%. To facilitate its industrial application, we further optimized the production conditions in a 5L bioreactor, and the xylanase activity reached 1.52 × 106 U/mL culture.

Engineering Biocatalysts for the C-H Activation of Fatty Acids by Ancestral Sequence Reconstruction.

Selective, one-step C-H activation of fatty acids from biomass is an attractive concept in sustainable chemistry. Biocatalysis has shown promise for generating high-value hydroxy acids, but to date enzyme discovery has relied on laborious screening and produced limited hits, which predominantly oxidise the subterminal positions of fatty acids. Herein we show that ancestral sequence reconstruction (ASR) is an effective tool to explore the sequence-activity landscape of a family of multidomain, self-sufficient P450 monooxygenases. We resurrected 11 catalytically active CYP116B ancestors, each with a unique regioselectivity fingerprint that varied from subterminal in the older ancestors to mid-chain in the lineage leading to the extant, P450-TT. In lineages leading to extant enzymes in thermophiles, thermostability increased from ancestral to extant forms, as expected if thermophily had arisen de novo. Our studies show that ASR can be applied to multidomain enzymes to develop active, self-sufficient monooxygenases as regioselective biocatalysts for fatty acid hydroxylation.

Engineering functional thermostable proteins using ancestral sequence reconstruction.

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Single mutations toggle the substrate selectivity of multifunctional Camptotheca secologanic acid synthases.

Terpene indole alkaloids (TIAs) are plant-derived specialized metabolites with widespread use in medicine. Species-specific pathways derive various TIAs from common intermediates, strictosidine or strictosidinic acid, produced by coupling tryptamine with secologanin or secologanic acid. The penultimate reaction in this pathway is catalyzed by either secologanin synthase (SLS) or secologanic acid synthase (SLAS) according to whether plants produce secologanin from loganin or secologanic acid from loganic acid. Previous work has identified SLSs and SLASs from different species, but the determinants of selectivity remain unclear. Here, combining molecular modeling, ancestral sequence reconstruction, and biochemical methodologies, we identified key residues that toggle SLS and SLAS selectivity in two CYP72A (cytochrome P450) subfamily enzymes from Camptotheca acuminata. We found that the positions of foremost importance are in substrate recognition sequence 1 (SRS1), where mutations to either of two adjacent histidine residues switched selectivity; His131Phe selects for and increases secologanin production whereas His132Asp selects for secologanic acid production. Furthermore, a change in SRS3 in the predicted substrate entry channel (Arg/Lys270Thr) and another in SRS4 at the start of the I-helix (Ser324Glu) decreased enzyme activity toward either substrate. We propose that the Camptotheca SLASs have maintained the broadened activities found in a common asterid ancestor, even as the Camptotheca lineage lost its ability to produce loganin while the campanulid and lamiid lineages specialized to produce secologanin by acquiring mutations in SRS1. The identification here of the residues essential for the broad substrate scope of SLASs presents opportunities for more tailored heterologous production of TIAs.

Early Nitrogenase Ancestors Encompassed Novel Active Site Diversity.

Ancestral sequence reconstruction (ASR) infers predicted ancestral states for sites within sequences and can constrain the functions and properties of ancestors of extant protein families. Here, we compare the likely sequences of inferred nitrogenase ancestors to extant nitrogenase sequence diversity. We show that the most-likely combinations of ancestral states for key substrate channel residues are not represented in extant sequence space, and rarely found within a more broadly defined physiochemical space-supporting that the earliest ancestors of extant nitrogenases likely had alternative substrate channel composition. These differences may indicate differing environmental selection pressures acting on nitrogenase substrate specificity in ancient environments. These results highlight ASR's potential as an in silico tool for developing hypotheses about ancestral enzyme functions, as well as improving hypothesis testing through more targeted in vitro and in vivo experiments.