Rule-based modeling and simulations of the inner kinetochore structure

Background

Combinatorial complexity is a central problem when modeling biochemical reaction networks, since the association of a few components can give rise to a large variation of protein complexes. Available classical modeling approaches are often insufficient for the analysis of very large and complex networks in detail. Recently, we developed a new rule-based modeling approach that facilitates the analysis of spatial and combinatorially complex problems. Here, we explore for the first time how this approach can be applied to a specific biological system, the human kinetochore, which is a multi-protein complex involving over 100 proteins.

Results

Applying our freely available SRSim software to a large data set on kinetochore proteins in human cells, we construct a spatial rule-based simulation model of the human inner kinetochore. The model generates an estimation of the probability distribution of the inner kinetochore 3D architecture and we show how to analyze this distribution using information theory. In our model, the formation of a bridge between CenpA and an H3 containing nucleosome only occurs efficiently for higher protein concentration realized during S-phase but may be not in G1. Above a certain nucleosome distance the protein bridge barely formed pointing towards the importance of chromatin structure for kinetochore complex formation. We define a metric for the distance between structures that allow us to identify structural clusters. Using this modeling technique, we explore different hypothetical chromatin layouts.

Conclusions

Applying a rule-based network analysis to the spatial kinetochore complex geometry allowed us to integrate experimental data on kinetochore proteins, suggesting a 3D model of the human inner kinetochore architecture that is governed by a combinatorial algebraic reaction network. This reaction network can serve as bridge between multiple scales of modeling. Our approach can be applied to other systems beyond kinetochores.     

A decade of high-resolution liquid chromatography of nucleic acids on styrene-divinylbenzene copolymers

The introduction of alkylated, nonporous poly-(styrene-divinylbenzene) microparticles in 1992 enabled the subsequent development of denaturing HPLC that has emerged as the most sensitive screening method for mutations to date. Denaturing HPLC has provided unprecedented insight into human origins and prehistoric migrations, accelerated the cloning of genes involved in mono- and polygenic traits, and facilitated the mutational analysis of more than a hundred candidate genes of human disease. A significant step toward increased sample-throughput and information content was accomplished by the recent introduction of monolithic poly(styrene-divinylbenzene) capillary columns. They have enabled the construction of capillary arrays amenable to multiplex analysis of fluorescent dye-labeled nucleic acids by laser-induced fluorescence detection. Hyphenation of denaturing HPLC with electrospray ionization mass spectrometry, on the other hand, has allowed the direct elucidation of the chemical nature of DNA variation and determination of phase of multiple alleles on a chromosome.     

Mass spectrometric identification of novel posttranslational modification sites in Huntingtin

Huntington's disease (HD) is caused by a CAG triplet repeat expansion in exon 1 of the Huntingtin (Htt) gene, encoding an abnormal expanded polyglutamine (polyQ) tract that confers toxicity to the mutant Htt (mHtt) protein. Recent data suggest that posttranslational modifications of mHtt modulate its cytotoxicity. To further understand the cytotoxic mechanisms of mHtt, we have generated HEK293 cell models stably expressing Strep- and FLAG-tagged Htt containing either 19Q (wild-type Htt), 55Q (mHtt), or 94Q (mHtt) repeats. Following tandem affinity purification, the tagged Htt and associated proteins were subjected to tandem mass spectrometry or 2D nano-LC tandem mass spectrometry and several novel modification sites of mHtt containing 55Q or 94Q were identified. These were phosphorylation sites located at Ser431 and Ser432, and ubiquitination site located at Lys444. The two phosphorylation sites were confirmed by Western blot analysis using phosphorylation site-specific antibodies. In addition, prevention of phosphorylation at the two serine sites altered mHtt toxicity and accumulation. These modifications of mHtt may provide novel therapeutic targets for effective treatment of the disorder.     

Prebiotic processes in planetary atmospheres

Within 40 years of experimental studies in prebiotic chemistry, most of the building blocks of the living systems have been synthesized in plausible conditions of the primitive Earth. The starting ingredients correspond to two complementary classes: volatile organics, and their non volatile oligomers. They may have been formed in the atmosphere on the primitive Earth and/or imported by extra-terrestrial sources. Organic chemistry is involved in meteorites, comets, in the giant planets and several of their satellites. Again this chemistry presents the two complementary aspects. In particular, with a dense reduced atmosphere rich in organic compounds in gas and aerosol phases, Titan appears as a natural laboratory for studying prebiotic chemistry at a planetary scale.     

Characterization of a novel Obg-like ATPase in the protozoan Trypanosoma cruzi

We characterized a gene encoding an YchF-related protein, TcYchF, potentially associated with the protein translation machinery of Trypanosoma cruzi. YchF belongs to the translation factor-related (TRAFAC) class of P-loop NTPases. The coding region of the gene is 1185 bp long and encodes a 44.3 kDa protein. BlastX searches showed TcYchF to be very similar (45-86%) to putative GTP-binding proteins from eukaryotes, including some species of trypanosomatids (Leishmania major and Trypanosoma brucei). A lower but significant level of similarity (38-43%) was also found between the predicted sequences of TcYchF and bacterial YyaF/YchF GTPases of the Spo0B-associated GTP-binding protein (Obg) family. Some of the most important features of the G domain of this family of GTPases are conserved in TcYchF. However, we found that TcYchF preferentially hydrolyzed ATP rather than GTP. The function of YyaF/YchF is unknown, but other members of the Obg family are known to be associated with ribosomal subunits. Immunoblots of the polysome fraction from sucrose gradients showed that TcYchF was associated with ribosomal subunits and polysomes. Immunoprecipitation assays showed that TcYchF was also associated with the proteasome of T. cruzi. Furthermore, inactivation of the T. brucei homolog of TcYchF by RNA interference inhibited the growth of procyclic forms of the parasite. These data suggest that this protein plays an important role in the translation machinery of trypanosomes.     

The mutational dynamics of the SARS-CoV-2 virus in serial passages in vitro

Since its outbreak in 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) keeps surprising the medical community by evolving diverse immune escape mutations in a rapid and effective manner. To gain deeper insight into mutation frequency and dynamics, we isolated ten ancestral strains of SARS-CoV-2 and performed consecutive serial incubation in ten replications in a suitable and common cell line and subsequently analysed them using RT-qPCR and whole genome sequencing. Along those lines we hoped to gain fundamental insights into the evolutionary capacity of SARS-CoV-2 in vitro. Our results identified a series of adaptive genetic changes, ranging from unique convergent substitutional mutations and hitherto undescribed insertions. The region coding for spike proved to be a mutational hotspot, evolving a number of mutational changes including the already known substitutions at positions S:484 and S:501. We discussed the evolution of all specific adaptations as well as possible reasons for the seemingly inhomogeneous potential of SARS-CoV-2 in the adaptation to cell culture. The combination of serial passage in vitro with whole genome sequencing uncovers the immense mutational potential of some SARS-CoV-2 strains. The observed genetic changes of SARS-CoV-2 in vitro could not be explained solely by selectively neutral mutations but possibly resulted from the action of directional selection accumulating favourable genetic changes in the evolving variants, along the path of increasing potency of the strain. Competition among a high number of quasi-species in the SARS-CoV-2 in vitro population gene pool may reinforce directional selection and boost the speed of evolutionary change.     

Kinetic analysis of the zinc-dependent deacetylase in the lipid A biosynthetic pathway

The first committed step of lipid A biosynthesis in Gram-negative bacteria is catalyzed by the zinc-dependent hydrolase LpxC that removes an acetate from the nitrogen at the 2' '-position of UDP-3-O-acyl-N-acetylglucosamine. Recent structural characterization by both NMR and X-ray crystallography provides many important details about the active site environment of LpxC from Aquifex aeolicus, a heat-stable orthologue that displays 32% sequence identity to LpxC from Escherichia coli. The detailed reaction mechanism and specific roles of active site residues for LpxC from A. aeolicus are further analyzed here. The pH dependencies of k(cat)/K(M) and k(cat) for the deacetylation of the substrate UDP-3-O-[(R)-3-hydroxymyristoyl]-GlcNAc are both bell-shaped. The ascending acidic limb (pK(1)) was fitted to 6.1 +/- 0.2 for k(cat) and 5.7 +/- 0.2 for k(cat)/K(M). The descending basic limb (pK(2)) was fitted to 8.0 +/- 0.2 for k(cat) and 8.4 +/- 0.2 for k(cat)/K(M). The pH dependence of the E73A mutant exhibits loss of the acidic limb, and the mutant retains only 0.15% activity versus the wild type. The pH dependencies of the other active site mutants H253A, K227A, H253A/K227A, and D234N remain bell-shaped, although their significantly lower activities (0.25%, 0.05%, 0.007%, and 0.57%, respectively) suggest that they contribute significantly to catalysis. Our cumulative data support a mechanism for LpxC wherein Glu73 serves as the general base for deprotonation and activation of the zinc-bound water.     

Direct peptide profiling of lateral cell groups of the antennal lobes of Manduca sexta reveals specific composition and changes in neuropeptide expression during development

The paired antennal lobes are the first integration centers for odor information in the insect brain. In the sphinx moth Manduca sexta, like in other holometabolous insects, they are formed during metamorphosis. To further understand mechanisms involved in the formation of this particularly well investigated brain area, we performed a direct peptide profiling of a well defined cell group (the lateral cell group) of the antennal lobe throughout development by MALDI-TOF mass spectrometry. Although the majority of the about 100 obtained ion signals represent still unknown substances, this first peptidomic characterization of this cell group indicated the occurrence of 12 structurally known neuropeptides. Among these peptides are helicostatin 1, cydiastatins 2, 3, and 4, M. sexta-allatotropin (Mas-AT), M. sexta-FLRFamide (Mas-FLRFamide) I, II, and III, nonblocked Mas-FLRFamide I, and M. sexta-myoinhibitory peptides (Mas-MIPs) III, V, and VI. The identity of two of the allatostatins (cydiastatins 3 and 4) and Mas-AT were confirmed by tandem mass spectrometry (MALDI-TOF/TOF). During development of the antennal lobe, number and frequency of ion signals including those representing known peptides generally increased at the onset of glomeruli formation at pupal Stage P7/8, with cydiastatin 2, helicostatin 1, and Mas-MIP V being the exceptions. Cydiastatin 2 showed transient occurrence mainly during the period of glomerulus formation, helicostatin 1 was restricted to late pupae and adults, while Mas-MIP V occurred exclusively in adult antennal lobes. The power of the applied direct mass spectrometric profiling lies in the possibility of chemically identifying neuropeptides of a given cell population in a fast and reliable manner, at any developmental stage in single specimens. The identification of neuropeptides in the antennal lobes now allows to specifically address the function of these signaling molecules during the formation of the antennal lobe network.     

Identification of residues in the human guanylate-binding protein 1 critical for nucleotide binding and cooperative GTP hydrolysis

The guanylate-binding proteins (GBPs) form a group of interferon-gamma inducible GTP-binding proteins which belong to the family of dynamin-related proteins. Like other members of this family, human guanylate-binding protein 1 (hGBP1) shows nucleotide-dependent oligomerisation that stimulates the GTPase activity of the protein. A unique feature of the GBPs is their ability to hydrolyse GTP to GDP and GMP. In order to elucidate the relationship between these findings, we designed point mutants in the phosphate-binding loop (P-loop) as well as in the switch I and switch II regions of the protein based on the crystal structure of hGBP1. These mutant proteins were analysed for their interaction with guanine nucleotides labeled with a fluorescence dye and for their ability to hydrolyse GTP in a cooperative manner. We identified mutations of amino acid residues that decrease GTPase activity by orders of magnitude a part of which are conserved in GTP-binding proteins. In addition, mutants in the P-loop were characterized that strongly impair binding of nucleotide. In consequence, together with altered GTPase activity and given cellular nucleotide concentrations this results in hGBP1 mutants prevailingly resting in the nucleotide-free (K51A and S52N) or the GTP bound form (R48A), respectively. Using size-exclusion chromatography and analytical ultracentrifugation we addressed the impact on protein oligomerisation. In summary, mutants of hGBP1 were identified and biochemically characterized providing hGBP1 locked in defined states in order to investigate their functional role in future cell biology studies.     

SMN-mediated assembly of RNPs: a complex story

Although many RNA-protein complexes or ribonucleoproteins (RNPs) assemble spontaneously in vitro, little is known about how they form in the environment of a living cell. Insight into RNP assembly has come unexpectedly from functional analyses of the survival motor neuron (SMN) protein, a gene product that is affected in the neuromuscular disease spinal muscular atrophy. These studies show that the assembly of spliceosomal U-rich small nuclear RNPs in vivo depends on the activity of two large protein complexes, one of which contains the SMN protein. These complexes might also facilitate the assembly of other cellular RNPs.