PocketOptimizer 2.0: A modular framework for computer-aided ligand-binding design

The ability to design customized proteins to perform specific tasks is of great interest. We are particularly interested in the design of sensitive and specific small molecule ligand-binding proteins for biotechnological or biomedical applications. Computational methods can narrow down the immense combinatorial space to find the best solution and thus provide starting points for experimental procedures. However, success rates strongly depend on accurate modeling and energetic evaluation. Not only intra- but also intermolecular interactions have to be considered. To address this problem, we developed PocketOptimizer, a modular computational protein design pipeline, that predicts mutations in the binding pockets of proteins to increase affinity for a specific ligand. Its modularity enables users to compare different combinations of force fields, rotamer libraries, and scoring functions. Here, we present a much-improved version––PocketOptimizer 2.0. We implemented a cleaner user interface, an extended architecture with more supported tools, such as force fields and scoring functions, a backbone-dependent rotamer library, as well as different improvements in the underlying algorithms. Version 2.0 was tested against a benchmark of design cases and assessed in comparison to the first version. Our results show how newly implemented features such as the new rotamer library can lead to improved prediction accuracy. Therefore, we believe that PocketOptimizer 2.0, with its many new and improved functionalities, provides a robust and versatile environment for the design of small molecule-binding pockets in proteins. It is widely applicable and extendible due to its modular framework. PocketOptimizer 2.0 can be downloaded at https://github.com/Hoecker-Lab/pocketoptimizer .     

Next-Generation Sequencing of Connective Tissue Genes in Patients with Classical Ehlers-Danlos Syndrome

Background: Ehlers-Danlos syndrome (EDS) is a common non-inflammatory, congenital connective tissue disorder. Classical type (cEDS) EDS is one of the more common forms, typically caused by mutations in the COL5A1 and COL5A2 genes, though causative mutations in the COL1A1 gene have also been described. Material and methods: The study group included 59 patients of Polish origin, diagnosed with cEDS. The analysis was performed on genomic DNA (gDNA) with NGS technology, using an Illumina sequencer. Thirty-five genes related to connective tissue were investigated. The pathogenicity of the detected variants was assessed by VarSome. Results: The NGS of 35 genes revealed variants within the COL5A1, COL5A2, COL1A1, and COL1A2 genes for 30 of the 59 patients investigated. Our panel detected no sequence variations for the remaining 29 patients. Discussion: Next-generation sequencing, with an appropriate multigene panel, showed great potential to assist in the diagnosis of EDS and other connective tissue disorders. Our data also show that not all causative genes giving rise to cEDS have been elucidated yet.     

Distinguishing excess mutations and increased cell death based on variant allele frequencies

Tumors often harbor orders of magnitude more mutations than healthy tissues. The increased number of mutations may be due to an elevated mutation rate or frequent cell death and correspondingly rapid cell turnover, or a combination of the two. It is difficult to disentangle these two mechanisms based on widely available bulk sequencing data, where sequences from individual cells are intermixed and, thus, the cell lineage tree of the tumor cannot be resolved. Here we present a method that can simultaneously estimate the cell turnover rate and the rate of mutations from bulk sequencing data. Our method works by simulating tumor growth and finding the parameters with which the observed data can be reproduced with maximum likelihood. Applying this method to a real tumor sample, we find that both the mutation rate and the frequency of death may be high.     

Naive helper T cells from BCG-vaccinated volunteers produce IFN-gamma and IL-5 to mycobacterial antigen-pulsed dendritic cells

Mycobacterium bovis bacillus Calmette-Guérin (BCG) is a live vaccine that has been used in routine vaccination against tuberculosis for nearly 80 years. However, its efficacy is controversial. The failure of BCG vaccination may be at least partially explained by the induction of poor or inappropriate host responses. Dendritic cells (DCs) are likely to play a key role in the induction of immune response to mycobacteria by polarizing the reactivity of T lymphocytes toward a Th1 profile, contributing to the generation of protective cellular immunity against mycobacteria. In this study we aimed to investigate the production of Th1 and Th2 cytokines by naive CD4+ T cells to mycobacterial antigen-pulsed DCs in the group of young, healthy BCG vaccinated volunteers. The response of naive helper T cells was compared with the response of total blood lymphocytes. Our present results clearly showed that circulating naive CD45RA+CD4+ lymphocytes from BCG-vaccinated subjects can become effector helper cells producing IFN-gamma and IL-5 under the stimulation by autologous dendritic cells presenting mycobacterial protein antigen-PPD or infected with live M. bovis BCG bacilli.     

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.     

Polychromatic (eight-color) slide-based cytometry for the phenotyping of leukocyte, NK, and NKT subsets.

BACKGROUND: Natural killer (NK) and NK T (NKT) cells are important in innate immune defense. Their unequivocal identification requires at least four antigens. Based on the expression of additional antigens, they can be further divided into functional subsets. For more accurate immunophenotyping and to describe multiple expression patterns of leukocyte subsets, an increased number of measurable colors is necessary. To take advantage of the technologic features offered by slide-based cytometry, repeated analysis was combined with sequential optical-filter changing. METHODS: Human peripheral blood leukocytes from healthy adult volunteers were labeled with antibodies by direct or indirect staining. Tandem dyes of Cy7 (phycoerythrin [PE]-/allophycocyanin [APC]-Cy7), Cy5.5 (PE-/APC-Cy5.5), and PE-Cy5 and the fluorochromes fluorescein isothiocyanate (FITC), PE, and APC were tested alone and in combinations. Optical filters of the laser scanning cytometer were 555 DRLP/BP 530/30 nm for photomultiplier tube (PMT) 1/FITC, 605 DRLP/BP 580/30 nm for PMT 2/PE, 740 DCXR/BP 670/20 nm for PMT 3/Cy5/APC, and BP 810/90 nm for PMT 4/Cy7. Filter PMT 3 was replaced for detection of PE/Cy5.5 and APC/Cy5.5 by 740 LP/BP 710/20 nm and the sample was remeasured. Both data files were merged into one to combine the different information on a single-cell basis. The combination of eight antibodies against CD3, CD4, CD8, CD14, CD16, CD19, CD45, and CD56 was used to characterize NK and NKT cells and their subsets. RESULTS: In this way Cy5.5 is measurable at 488-nm and 633-nm excitation. Further, with the two different filters it is possible to distinguish Cy5 from Cy5.5 in the same detection channel (PMT 3). With this method we identified NK and NKT cells, subsets of NK (CD3-16+56+, CD3-16+56-, CD3-16-56+) and NKT (CD3+16+56+, CD3+16-56+) and their CD4+8-, CD4-8+, CD4-8- and CD4+8+ subsets. CONCLUSION: With our adaptations it is possible to discriminate tandem conjugates of Cy5, Cy5.5, and Cy7 for eight-color immunophenotyping. Using this method, novel rare subsets of NK and NKT cells that are CD4/CD8 double positive are reported for the first time.     

Slow Oscillation Amplitudes and Up-State Lengths Relate to Memory Improvement

There is growing evidence of the active involvement of sleep in memory consolidation. Besides hippocampal sharp wave-ripple complexes and sleep spindles, slow oscillations appear to play a key role in the process of sleep-associated memory consolidation. Furthermore, slow oscillation amplitude and spectral power increase during the night after learning declarative and procedural memory tasks. However, it is unresolved whether learning-induced changes specifically alter characteristics of individual slow oscillations, such as the slow oscillation up-state length and amplitude, which are believed to be important for neuronal replay. 24 subjects (12 men) aged between 20 and 30 years participated in a randomized, within-subject, multicenter study. Subjects slept on three occasions for a whole night in the sleep laboratory with full polysomnography. Whereas the first night only served for adaptation purposes, the two remaining nights were preceded by a declarative word-pair task or by a non-learning control task. Slow oscillations were detected in non-rapid eye movement sleep over electrode Fz. Results indicate positive correlations between the length of the up-state as well as the amplitude of both slow oscillation phases and changes in memory performance from pre to post sleep. We speculate that the prolonged slow oscillation up-state length might extend the timeframe for the transfer of initial hippocampal to long-term cortical memory representations, whereas the increase in slow oscillation amplitudes possibly reflects changes in the net synaptic strength of cortical networks.     

Laser scanning cytometry in human brain slices.

BACKGROUND: The Laser Scanning Cytometry (LSC) offers quantitative fluorescence analysis of cell suspensions and tissue sections. METHODS: We adapted this technique to immunohistochemical labelled human brain slices. RESULTS: We were able to identify neurons according to their labelling and to display morphological structures such as the lamination of the entorhinal cortex. Further, we were able to distinguish between neurons with and without cyclin B1 expression and we could assign the expression of cyclin B1 to the cell islands of layer II and the pyramidal neurons of layer V of the entorhinal cortex in Alzheimer's disease effected brain. In addition, we developed a method depicting the three-dimensional distribution of the cells in intact tissue sections. CONCLUSIONS: In this pilot experiments we could demonstrate the power of the LSC for the analysis of human brain sections.