Static Homology Research Articles

Embedding Dynamics in Intrinsic Physicochemical Profiles of Market-Stage Antibody-Based Biotherapeutics.

Adequate stability, manufacturability, and safety are crucial to bringing an antibody-based biotherapeutic to the market. Following the concept of holistic in silico developability, we introduce a physicochemical description of 91 market-stage antibody-based biotherapeutics based on orthogonal molecular properties of variable regions (Fvs) embedded in different simulation environments, mimicking conditions experienced by antibodies during manufacturing, formulation, and in vivo. In this work, the evaluation of molecular properties includes conformational flexibility of the Fvs using molecular dynamics (MD) simulations. The comparison between static homology models and simulations shows that MD significantly affects certain molecular descriptors like surface molecular patches. Moreover, the structural stability of a subset of Fv regions is linked to changes in their specific molecular interactions with ions in different experimental conditions. This is supported by the observation of differences in protein melting temperatures upon addition of NaCl. A DEvelopability Navigator In Silico (DENIS) is proposed to compare mAb candidates for their similarity with market-stage biotherapeutics in terms of physicochemical properties and conformational stability. Expanding on our previous developability guidelines (Ahmed et al. Proc. Natl. Acad. Sci. 2021, 118 (37), e2020577118), the hydrodynamic radius and the protein strand ratio are introduced as two additional descriptors that enable a more comprehensive in silico characterization of biotherapeutic drug candidates. Test cases show how this approach can facilitate identification and optimization of intrinsically developable lead candidates. DENIS represents an advanced computational tool to progress biotherapeutic drug candidates from discovery into early development by predicting drug properties in different aqueous environments.

A revised phylogeny of the New Caledonian endemic genus Troglosiro (Opiliones : Cyphophthalmi : Troglosironidae) with the description of four new species

The Cyphophthalmi genus Troglosiro (the only genus of the family Troglosironidae) is endemic to New Caledonia, representing one of the oldest lineages of this emerged part of Zealandia. Its species are short-range endemics, many known from single localities. Here we examined the phylogenetic relationships of Troglosironidae using standard Sanger-sequenced markers (nuclear 18S rRNA, 28S rRNA, and mitochondrial 16S rRNA and cytochrome c oxidase subunit I) and a combination of phylogenetic methods, including parsimony under Direct Optimization and maximum likelihood with static homology. We also applied a diversity of species delimitation methods, including distance-based, topology-based and unsupervised machine learning to evaluate previous species designations. Finally, we used a combination of genetic and morphological information to describe four new species – T. dogny sp. nov., T. pin sp. nov., T. pseudojuberthiei sp. nov. and T. sharmai sp. nov. – and discuss them in the broader context of the phylogeny and biogeographic history of the family. A key to the species of Troglosiro is also provided. urn:lsid:zoobank.org:pub:93541314-8309-468C-BB77-B34C3A81137E

Conformational Ensembles of Antibodies Determine Their Hydrophobicity

A major challenge in the development of antibody biotherapeutics is their tendency to aggregate. One root cause for aggregation is exposure of hydrophobic surface regions to the solvent. Many current techniques predict the relative aggregation propensity of antibodies via precalculated scales for the hydrophobicity or aggregation propensity of single amino acids. However, those scales cannot describe the nonadditive effects of a residue’s surrounding on its hydrophobicity. Therefore, they are inherently limited in their ability to describe the impact of subtle differences in molecular structure on the overall hydrophobicity. Here, we introduce a physics-based approach to describe hydrophobicity in terms of the hydration free energy using grid inhomogeneous solvation theory (GIST). We apply this method to assess the effects of starting structures, conformational sampling, and protonation states on the hydrophobicity of antibodies. Our results reveal that high-quality starting structures, i.e., crystal structures, are crucial for the prediction of hydrophobicity and that conformational sampling can compensate errors introduced by the starting structure. On the other hand, sampling of protonation states only leads to good results when combined with high-quality structures, whereas it can even be detrimental otherwise. We conclude by pointing out that a single static homology model may not be adequate for predicting hydrophobicity.

Phenotypic characters of static homology increase phylogenetic stability under direct optimization of otherwise dynamic homology characters.

Direct optimization of unaligned sequence characters provides a natural framework to explore the sensitivity of phylogenetic hypotheses to variation in analytical parameters. Phenotypic data, when combined into such analyses, are typically analyzed with static homology correspondences unlike the dynamic homology sequence data. Static homology characters may be expected to constrain the direct optimization and thus, potentially increase the similarity of phylogenetic hypotheses under different cost sets. However, whether a total-evidence approach increases the phylogenetic stability or not remains empirically largely unexplored. Here, I studied the impact of static homology data on sensitivity using six empirical data sets composed of several molecular markers and phenotypic data. The inclusion of static homology phenotypic data increased the average stability of phylogenetic hypothesis in five out of the six data sets. To investigate if any static homology characters would have similar effect, the analyses were repeated with randomized phenotypic data, and with one of the molecular markers fixed as static homology characters. These analyses had, on average, almost no effect on the phylogenetic stability, although the randomized phenotypic data sometimes resulted in even higher stability than empirical phenotypic data. The impact was related to the strength of the phylogenetic signal in the phenotypic data: higher average jackknife support of the phenotypic tree correlated with stronger stabilizing effect in the total-evidence analysis. Phenotypic data with a strong signal made the total-evidence trees topologically more similar to the phenotypic trees, thus, they constrained the dynamic homology correspondences of the sequence data. Characters that increase phylogenetic stability are particularly valuable for phylogenetic inference. These results indicate an important role and additive value of phenotypic data in increasing the stability of phylogenetic hypotheses in total-evidence analyses.

Molecular phylogeny of the genus Bolusiella (Orchidaceae, Angraecinae)

Recent molecular studies have suggested the monophyly of Bolusiella, a small orchid genus comprising five species and one subspecies from Continental Africa, but sampling has been limited. Using the species delimitation presented in the recent taxonomic revision of the genus, this study aimed to confirm the monophyly of Bolusiella and assess the interspecific relationships using a comprehensive sampling and various analytical methods. DNA sequences of one nuclear spacer region (ITS-1) and five plastid regions (matK, rps16, trnL–trnF, trnC–petN, and ycf1) from 20 specimens representing all five species of the genus were analyzed using static homology, dynamic homology, and Bayesian methods. The monophyly of both the genus Bolusiella and each of its five species was confirmed, corroborating the previously published taxonomic revision. The use of dynamic homology methods was not conclusive for this particular group. The results of the total evidence analysis (combining all six sequence regions) using the dynamic homology approach yielded a slightly different hypothesis regarding interspecific relationships (namely the exchange of B. talbotii and Bolusiella iridifolia as the earliest diverging lineage), probably because the nodes in question are supported by a small subset of conflicting characters, compared to the hypotheses resulting from the static homology and Bayesian methods, which are congruent with the results of previous studies.

Hazard assessment of benchmark metallic nanomaterials in alveolar epithelial cells

A major challenge in the development of antibody biotherapeutics is their tendency to aggregate. One root cause for aggregation is exposure of hydrophobic surface regions to the solvent. Many current techniques predict the relative aggregation propensity of antibodies via precalculated scales for the hydrophobicity or aggregation propensity of single amino acids. However, those scales cannot describe the nonadditive effects of a residue’s surrounding on its hydrophobicity. Therefore, they are inherently limited in their ability to describe the impact of subtle differences in molecular structure on the overall hydrophobicity. Here, we introduce a physics-based approach to describe hydrophobicity in terms of the hydration free energy using grid inhomogeneous solvation theory (GIST). We apply this method to assess the effects of starting structures, conformational sampling, and protonation states on the hydrophobicity of antibodies. Our results reveal that high-quality starting structures, i.e., crystal structures, are crucial for the prediction of hydrophobicity and that conformational sampling can compensate errors introduced by the starting structure. On the other hand, sampling of protonation states only leads to good results when combined with high-quality structures, whereas it can even be detrimental otherwise. We conclude by pointing out that a single static homology model may not be adequate for predicting hydrophobicity.

Homology Modeling of the CheW Coupling Protein of the Chemotaxis Signaling Complex

Homology models of the E. coli and T. maritima chemotaxis protein CheW were constructed to assess the quality of structural predictions and their applicability in chemotaxis research: i) a model of E. coli CheW was constructed using the T. maritima CheW NMR structure as a template, and ii) a model of T. maritima CheW was constructed using the E. coli CheW NMR structure as a template. The conformational space accessible to the homology models and to the NMR structures was investigated using molecular dynamics and Monte Carlo simulations. The results show that even though static homology models of CheW may be partially structurally different from their corresponding experimentally determined structures, the conformational space they can access through their dynamic variations can be similar, for specific regions of the protein, to that of the experimental NMR structures. When CheW homology models are allowed to explore their local accessible conformational space, modeling can provide a rational path to predicting CheW interactions with the MCP and CheA proteins of the chemotaxis complex. Homology models of CheW (and potentially, of other chemotaxis proteins) should be seen as snapshots of an otherwise larger ensemble of accessible conformational space.

Phylogenetics of scolopendromorph centipedes: can denser taxon sampling improve an artificial classification?

Previous phylogenetic analyses of the centipede order Scolopendromorpha indicated a fundamental division into blind and ocellate clades. These analyses corroborated the monophyly of most families and tribes but suggested that several species-rich, cosmopolitan genera in traditional and current classifications are polyphyletic. Denser taxon sampling is applied to a dataset of 122 morphological characters and sequences for four nuclear and mitochondrial loci. Phylogenetic analyses including 98 species and subspecies of Scolopendromorpha employ parsimony under dynamic and static homology schemes as well as maximum likelihood and Bayesian inference of multiple sequence alignments. The monotypic Australian genera Notiasemus and Kanparka nest within Cormocephalus and Scolopendra, respectively, and the New Caledonian Campylostigmus is likewise a clade within Cormocephalus. New World Scolopendra are more closely related to Hemiscolopendra and Arthrorhabdus than to Scolopendra s.s., which is instead closely allied to Asanada; the tribe Asanadini nests within Scolopendrini for molecular and combined datasets. The generic classification of Otostigmini has a poor fit to phylogenetic relationships, although nodal support within this tribe is weak. New synonymies are proposed for Ectonocryptopinae Shelley & Mercurio, 2005 (= Newportiinae Pocock, 1896), Asanadini Verhoeff, 1907 (= Scolopendrini Leach, 1814), and Kanparka Waldock & Edgecombe, 2012 (= Scolopendra Linnaeus, 1758). Scolopendrid systematics largely depicts incongruence between phylogeny and classification rather than between morphology and molecules.

Molecular phylogeny of the New World Dipsadidae (Serpentes: Colubroidea): a reappraisal.

We present a phylogenetic analysis of the New World dipsadids based on an expanded data matrix that includes 246 terminal taxa including 196 dipsadids. The species are sampled for eight genes (12S, 16S, cytb, nd2, nd4, bdnf, c-mos, rag2). The data are explored using two distinct optimality procedures-maximum parsimony and maximum likelihood-and two alignment strategies-dynamic homology and static homology. Two previously unsampled dipsadid genera, Sordellina and Rhachidelus, are now included in the analysis. The definitions of the genera, Erythrolamprus, Clelia, Hypsirhynchus, Philodryas and Phimophis, and the tribes Alsophiini, Echinantherini and Conophiini, are revised. In order to maintain monophyly, the genus Umbrivaga is synonymized with Erythrolamprus, and two new genera are erected to accommodate Phimophis iglesiasi and Clelia rustica, as well as their closely related species. The West Indian genera Schwartzophis, Darlingtonia, Antillophis and Ocyophis are resurrected. © The Willi Hennig Society 2012.

Study of Human Orexin-1 and -2 G-Protein-Coupled Receptors with Novel and Published Antagonists by Modeling, Molecular Dynamics Simulations, and Site-Directed Mutagenesis

The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. Site-directed mutagenesis (SDM) and domain exchange (chimera) studies have provided important insight into key features of the OX1 and OX2 binding sites. However, the precise determinants of antagonist binding and selectivity are still not fully known. In this work, we used homology modeling of OX receptors to direct further SDM studies. These SDM studies were followed by molecular dynamics (MD) simulations to rationalize the full scope of the SDM data and to explain the role of each mutated residue in the binding and selectivity of a set of OX antagonists: Almorexant (dual OX1 and OX2 antagonist), SB-674042 (OX1 selective antagonist), EMPA (OX2 selective antagonist), and others. Our primary interest was focused on transmembrane helix 3 (TM3), which is identified as being of great importance for the selectivity of OX antagonists. These studies revealed conformational differences between the TM3 helices of OX1 and OX2, resulting from differences in amino acid sequences of the OX receptors that affect key interhelical interactions formed between TM3 and neighboring TM domains. The MD simulation protocol used here, which was followed by flexible docking studies, went beyond the use of static models and allowed for a more detailed exploration of the OX structures. In this work, we have demonstrated how even small differences in the amino acid sequences of GPCRs can lead to significant differences in structure, antagonist binding affinity, and selectivity of these receptors. The MD simulations allowed refinement of the OX receptor models to a degree that was not possible with static homology modeling alone and provided a deeper rationalization of the SDM data obtained. To validate these findings and to demonstrate that they can be usefully applied to the design of novel, very selective OX antagonists, we show here two examples of antagonists designed in house: EP-109-0092 (OX1 selective) and EP-009-0513 (OX2 selective).

An extraordinary new genus of spiders from Western Australia with an expanded hypothesis on the phylogeny of Tetragnathidae (Araneae)

We describe Pinkfloydia Hormiga & Dimitrov gen. nov., a new genus of tetragnathid spiders from Western Australia and study its phylogenetic placement. The taxon sampling from our previous cladistic studies was expanded, with the inclusion of representatives of additional tetragnathid genera and outgroup taxa. Sequences from six genetic markers, 12S, 16S, 18S, 28S, cytochrome c oxidase subunit 1, and histone 3, along with morphological and behavioural data were used to infer tetragnathid relationships. These data were analysed using parsimony (under both static homology and dynamic optimization) and Bayesian methods. Our results indicate that Pinkfloydia belongs to the ‘Nanometa’ clade. We also propose a revised set of synapomorphies to define this lineage. Based on the new evidence presented here we propose a revised hypothesis for the intrafamilial relationships of Tetragnathidae and show that Mimetidae is most likely the sister group of Tetragnathidae. The single species in this genus so far, Pinkfloydia harveii Dimitrov& Hormiga sp. nov., is described in detail and its web architecture documented and illustrated. © 2011 The Linnean Society of London, Zoological Journal of the Linnean Society, 2011, 161, 735–768.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence.

BackgroundMites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei.ResultsThe phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infra-ordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes.ConclusionsPrevious studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Assessing the molluscan hypothesis Serialia (Monoplacophora+Polyplacophora) using novel molecular data

A consensus on molluscan relationships has yet to be achieved, largely because of conflicting morphological and molecular hypotheses. Monoplacophora show marked seriality of ctenidia, atria, muscles and nephridia and this has been interpreted as plesiomorphic for Mollusca, reflecting a segmented ancestry. More recently this seriality, also partly seen in Polyplacophora, has been seen as a derived condition. Analysis of the first published monoplacophoran DNA sequence from Laevilipilina antarctica Warén & Hain, 1992 [Giribet, G., Okusu, A., Lindgren, A.R., Huff, S., Schrödl, M., Nishiguchi, M.K., 2006. Evidence for a clade composed of molluscs with serially repeated structures: Monoplacophorans are related to chitons. Proc. Natl. Acad. Sci. USA 103, 7723-7728. 10.1073/pnas.0602578103], showed Monoplacophora inside Polyplacophora. These taxa were then grouped under the name Serialia, reflecting the hypothesis that their seriality is a synapomorphy. Subsequent examination revealed that part of the L. antarctica published sequence was the result of contamination with Polyplacophora (Giribet, Supplementary Material S1). We collected and sequenced another monoplacophoran, Laevipilina hyalina McLean, 1979, resulting in the first multi-gene dataset representing all molluscan classes. Our analyses did not show unambiguous support for Serialia. Model-based approaches strongly supported Serialia as a clade, however, parsimony analyses under dynamic and static homology did not resolve the position of Monoplacophora. Although our study provides support for Serialia and none for Conchifera, it appears that further resolution of molluscan relationships will require large increases of data.

Welcome back New Zealand: regional biogeography and Gondwanan origin of three endemic genera of mite harvestmen (Arachnida, Opiliones, Cyphophthalmi)

AbstractAim Biogeographers have long been intrigued by New Zealand’s biota due to its unique combination of typical ‘continental’ and ‘island’ characteristics. The New Zealand plateau rifted from the former supercontinent Gondwana c. 80 Ma, and has been isolated from other land masses ever since. Therefore, the flora and fauna of New Zealand include lineages that are Gondwanan in origin, but also include a very large number of endemics. In this study, we analyse the evolutionary relationships of three genera of mite harvestmen (Arachnida, Opiliones, Cyphophthalmi) endemic to New Zealand, both to each other and to their temperate Gondwanan relatives found in Australia, Chile, Sri Lanka and South Africa.Location New Zealand (North Island, South Island and Stewart Island).Methods A total of 94 specimens of the family Pettalidae in the suborder Cyphophthalmi were studied, representing 31 species and subspecies belonging to three endemic genera from New Zealand (Aoraki, Neopurcellia and Rakaia) plus six other members of the family from Chile, South Africa, Sri Lanka and Australia. The phylogeny of these taxa was constructed using morphological and molecular data from five nuclear and mitochondrial genes (18S rRNA, 28S rRNA, 16S rRNA, cytochrome c oxidase subunit I and histone H3, totalling c. 5 kb), which were analysed using dynamic as well as static homology under a variety of optimality criteria.Results The results showed that each of the three New Zealand cyphophthalmid genera is monophyletic, and occupies a distinct geographical region within the archipelago, grossly corresponding to palaeogeographical regions. All three genera of New Zealand mite harvestmen fall within the family Pettalidae with a classic temperate Gondwanan distribution, but they do not render any other genera paraphyletic.Main conclusions Our study shows that New Zealand’s three genera of mite harvestmen are unequivocally related to other members of the temperate Gondwanan family Pettalidae. Monophyly of each genus contradicts the idea of recent dispersal to New Zealand. Within New Zealand, striking biogeographical patterns are apparent in this group of short‐range endemics, particularly in the South Island. These patterns are interpreted in the light of New Zealand’s turbulent geological history and present‐day patterns of forest cover.

Phylogenetic relationships of the spider family Tetragnathidae (Araneae, Araneoidea) based on morphological and DNA sequence data.

The monophyly of Tetragnathidae including the species composition of the family (e.g., Are Nephila and their relatives part of this lineage?), the phylogenetic relationships of its various lineages, and the exact placement of Tetragnathidae within Araneoidea have been three recalcitrant problems in spider systematics. Most studies on tetragnathid phylogeny have focused on morphological and behavioral data, but little molecular work has been published to date. To address these issues we combine previous morphological and behavioral data with novel molecular data including nuclear ribosomal RNA genes 18S and 28S, mitochondrial ribosomal RNA genes 12S and 16S and protein-coding genes from the mitochondrion [cytochrome c oxidase subunit I (COI)] and from the nucleus (histone H3), totaling ca. 6.3 kb of sequence data per taxon. These data were analyzed using direct optimization and static homology using both parsimony and Bayesian methods. Our results indicate monophyly of Tetragnathidae, Tetragnathinae, Leucauginae, the "Nanometa clade" and the subfamily Metainae, which, with the exception of the later subfamily, received high nodal support. Morphological synapomorphies that support these clades are also discussed. The position of tetragnathids with respect to the rest of the araneoid spiders remains largely unresolved but tetragnathids and nephilids were never recovered as sister taxa. The combined dataset suggests that Nephilidae is sister to Araneidae; furthermore, the sister group of Nephila is the clade composed by Herennia plus Nephilengys and this pattern has clear implications for understanding the comparative biology of the group. Tetragnathidae is most likely sister to some members of the "reduced piriform clade" and nephilids constitute the most-basal lineage of araneids.