iNClusive: a database collecting useful information on non-canonical amino acids and their incorporation into proteins for easier genetic code expansion implementation.

The incorporation of non-canonical amino acids (ncAAs) into proteins is a powerful technique used in various research fields. Genetic code expansion (GCE) is the most common way to achieve this: a specific codon is selected to be decoded by a dedicated tRNA orthogonal to the endogenous ones. In the past 30 years, great progress has been made to obtain novel tRNA synthetases (aaRSs) accepting a variety of ncAAs with distinct physicochemical properties, to develop robust in vitro assays or approaches for codon reassignment. This sparked the use of the technique, leading to the accumulation of publications, from which gathering all relevant information can appear daunting. Here we present iNClusive (https://non-canonical-aas.biologie.uni-freiburg.de/), a manually curated, extensive repository using standardized nomenclature that provides organized information on ncAAs successfully incorporated into target proteins as verified by mass spectrometry. Since we focused on tRNA synthetase-based tRNA loading, we provide the sequence of the tRNA and aaRS used for the incorporation. Derived from more than 687 peer-reviewed publications, it currently contains 2432 entries about 466 ncAAs, 569 protein targets, 500 aaRSs and 144 tRNAs. We foresee iNClusive will encourage more researchers to experiment with ncAA incorporation thus contributing to the further development of this exciting technique.

patcHwork: a user-friendly pH sensitivity analysis web server for protein sequences and structures.

pH regulates protein function and interactions by altering the charge of individual residues causing loss or gain of intramolecular noncovalent bonds, which may lead to structural rearrangements. While tools to analyze residue-specific charge distribution of proteins at a given pH exist, currently no tool is available to investigate noncovalent bond changes at two different pH values. To make protein pH sensitivity analysis more accessible, we developed patcHwork, a web server that combines the identification of amino acids undergoing a charge shift with the determination of affected noncovalent bonds at two user-defined pH values. At the sequence-only level, patcHwork applies the Henderson-Hasselbalch equation to determine pH-sensitive residues. When the 3D protein structure is available, patcHwork can be employed to gain mechanistic understanding of the effect of pH. This is achieved using the PDB2PQR and PROPKA tools and noncovalent bond determination algorithms. A user-friendly interface allows visualizing pH-sensitive residues, affected salt bridges, hydrogen bonds and aromatic (pi-pi and cation-pi) interactions. patcHwork can be used to identify patches, a new concept we propose of pH-sensitive residues in close proximity on the protein, which may have a major impact on function. We demonstrate the attractiveness of patcHwork studying experimentally investigated pH-sensitive proteins (https://patchwork.biologie.uni-freiburg.de/).

Engineering AraC to make it responsive to light instead of arabinose.

The L-arabinose-responsive AraC and its cognate PBAD promoter underlie one of the most often used chemically inducible prokaryotic gene expression systems in microbiology and synthetic biology. Here, we change the sensing capability of AraC from L-arabinose to blue light, making its dimerization and the resulting PBAD activation light-inducible. We engineer an entire family of blue light-inducible AraC dimers in Escherichia coli (BLADE) to control gene expression in space and time. We show that BLADE can be used with pre-existing L-arabinose-responsive plasmids and strains, enabling optogenetic experiments without the need to clone. Furthermore, we apply BLADE to control, with light, the catabolism of L-arabinose, thus externally steering bacterial growth with a simple transformation step. Our work establishes BLADE as a highly practical and effective optogenetic tool with plug-and-play functionality-features that we hope will accelerate the broader adoption of optogenetics and the realization of its vast potential in microbiology, synthetic biology and biotechnology.

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study.

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.

Chromatosome Structure and Dynamics from Molecular Simulations.

Chromatosomes are fundamental units of chromatin structure that are formed when a linker histone protein binds to a nucleosome. The positioning of the linker histone on the nucleosome influences the packing of chromatin. Recent simulations and experiments have shown that chromatosomes adopt an ensemble of structures that differ in the geometry of the linker histone-nucleosome interaction. In this article we review the application of Brownian, Monte Carlo, and molecular dynamics simulations to predict the structure of linker histone-nucleosome complexes, to study the binding mechanisms involved, and to predict how this binding affects chromatin fiber structure. These simulations have revealed the sensitivityof the chromatosome structure to variations in DNA and linker histone sequence, as well as to posttranslational modifications, thereby explaining the structural variability observed in experiments. We propose that a concerted application of experimental and computational approaches will reveal the determinants of chromatosome structural variability and how it impacts chromatin packing.

Dependence of Chromatosome Structure on Linker Histone Sequence and Posttranslational Modification.

Linker histone (LH) proteins play a key role in higher-order structuring of chromatin for the packing of DNA in eukaryotic cells and in the regulation of genomic function. The common fruit fly (Drosophila melanogaster) has a single somatic isoform of the LH (H1). It is thus a useful model organism for investigating the effects of the LH on nucleosome compaction and the structure of the chromatosome, the complex formed by binding of an LH to a nucleosome. The structural and mechanistic details of how LH proteins bind to nucleosomes are debated. Here, we apply Brownian dynamics simulations to compare the nucleosome binding of the globular domain of D. melanogaster H1 (gH1) and the corresponding chicken (Gallus gallus) LH isoform, gH5, to identify residues in the LH that critically affect the structure of the chromatosome. Moreover, we investigate the effects of posttranslational modifications on the gH1 binding mode. We find that certain single-point mutations and posttranslational modifications of the LH proteins can significantly affect chromatosome structure. These findings indicate that even subtle differences in LH sequence can significantly shift the chromatosome structural ensemble and thus have implications for chromatin structure and transcriptional regulation.

Conformational selection and dynamic adaptation upon linker histone binding to the nucleosome.

Linker histones are essential for DNA compaction in chromatin. They bind to nucleosomes in a 1:1 ratio forming chromatosomes. Alternative configurations have been proposed in which the globular domain of the linker histone H5 (gH5) is positioned either on- or off-dyad between the nucleosomal and linker DNAs. However, the dynamic pathways of chromatosome assembly remain elusive. Here, we studied the conformational plasticity of gH5 in unbound and off-dyad nucleosome-bound forms with classical and accelerated molecular dynamics simulations. We find that the unbound gH5 converts between open and closed conformations, preferring the closed form. However, the open gH5 contributes to a more rigid chromatosome and restricts the motion of the nearby linker DNA through hydrophobic interactions with thymidines. Moreover, the closed gH5 opens and reorients in accelerated simulations of the chromatosome. Brownian dynamics simulations of chromatosome assembly, accounting for a range of amplitudes of nucleosome opening and different nucleosome DNA sequences, support the existence of both on- and off-dyad binding modes of gH5 and reveal alternative, sequence and conformation-dependent chromatosome configurations. Taken together, these findings suggest that the conformational dynamics of linker histones and nucleosomes facilitate alternative chromatosome configurations through an interplay between induced fit and conformational selection.

SDA 7: A modular and parallel implementation of the simulation of diffusional association software.

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein-protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration-dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object-oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance.

Curve of Spee and its relationship to vertical eruption of teeth among different malocclusion groups.

INTRODUCTION: Our objectives were to assess the depth of the curve of Spee (COS) in different malocclusion groups, to relate this to the eruption of anterior or posterior teeth quantitatively, and to determine whether the depth of the COS is affected by the vertical eruption of anterior or posterior teeth. METHODS: Two hundred conventional lateral cephalograms and 3-dimensional models of untreated patients (70 boys, mean age: 16.4 ± 1.4 years; 130 young women, mean age: 18.1 ± 1.8 years) were included and assigned to 4 malocclusion groups as Class I, Class II Division 1, Class II Division 2, and Class III. The depth of the COS, overjet, and overbite were measured on 3-dimensional models. The perpendicular distance between the incisal tip of the mandibular central incisor (L1-MP), the deepest point of the COS (S-MP), and the distobuccal cusp tip of the mandibular second molar (L7-MP) to the mandibular plane were calculated and proportioned with each other. The Pearson correlation coefficient was calculated, and multiple linear regression analysis was carried out. Also, multivariate analysis of variance was performed at the P <0.05 level. RESULTS: The mesiobuccal cusp of the first molar was the deepest part of the COS in all groups, with a maximum depth of 2.44 ± 0.73 mm in the Class II Division 1 subjects and a minimum depth of 1.76 ± 0.94 in the Class III subjects. The depth of the COS changed as follows: Class II Division 1 > Class II Division 2 > Class I > Class III malocclusion groups. Statistically significant positive correlations were found between the depth of the COS and L1-MP/S-MP (r = 0.541) and L7-MP/S-MP (r = 0.269) in the Class I and Class III subjects, and between the depth of the COS and overjet (r = 0.483) and L7-MP/S-MP (r = 0.289) in the Class II Division 1 subjects. All variables except overjet had positive correlations with the depth of the COS in Class II Division 2 subjects. The multivariate analysis of variance showed statistically significant differences in overjet, overbite, L1-MP/S-MP, L7-MP/S-MP, and the depth of the COS (P <0.001) among the groups. CONCLUSIONS: Although the overjet differed, vertical eruption of the anterior teeth did not differ among the different malocclusion groups and had a significant contribution to the depth of the COS in subjects with Class I and Class III malocclusions.

Development of the curve of Spee in Class II subdivision malocclusion: a longitudinal study.

OBJECTIVE: To compare the depth of curve of Spee (COS) in Class I and II sides of Class II subdivision malocclusion longitudinally and to describe the changes in the COS in relation to dental arch over time. MATERIALS AND METHOD: The study group comprised 17 subjects exhibiting Class II subdivision malocclusion. The depth of COS, intermolar width, distances from mesial anatomic contact points of the first permanent molars to the contact point of the central incisors in Class I and II sides and arch length in mandible were analysed on digital models in three different time points based on cervical vertebral maturation (T1; mean age: 12.4 years, T2; mean age: 15.1 years, and T3: mean age: 19.1 years). Pearson's correlation coefficients were calculated and linear multiple regression analysis with enter method was carried out. RESULTS: No significant differences were found in the depth of COS between Class I and II sides in all time periods. The mesiobuccal cusp of the first molar was the deepest part of COS in both sides and in all time periods, with a maximum depth of 2.37±0.83mm in T1 and a minimum depth of 1.91±0.66mm in T3. Between baseline and final follow-up; the depth of COS, ach length, and distance in Class I side decreased significantly. Arch length had a significant correlation with the depth of COS (r = 0.471) in T3. CONCLUSIONS: As the depth of COS did not differ between either the Class II or I sides, the same mechanics can be used to level the COS in either sides of the mandibular arch.

Int&in: A machine learning-based web server for active split site identification in inteins.

Inteins are proteins that excise themselves out of host proteins and ligate the flanking polypeptides in an auto-catalytic process called protein splicing. In nature, inteins are either contiguous or split. In the case of split inteins, the two fragments must first form a complex for the splicing to occur. Contiguous inteins have previously been artificially split in two fragments because split inteins allow for distinct applications than contiguous ones. Even naturally split inteins have been split at unnatural split sites to obtain fragments with reduced affinity for one another, which are useful to create conditional inteins or to study protein-protein interactions. So far, split sites in inteins have been heuristically identified. We developed Int&in, a web server freely available for academic research (https://intein.biologie.uni-freiburg.de) that runs a machine learning model using logistic regression to predict active and inactive split sites in inteins with high accuracy. The model was trained on a dataset of 126 split sites generated using the gp41-1, Npu DnaE and CL inteins and validated using 97 split sites extracted from the literature. Despite the limited data size, the model, which uses various protein structural features, as well as sequence conservation information, achieves an accuracy of 0.79 and 0.78 for the training and testing sets, respectively. We envision Int&in will facilitate the engineering of novel split inteins for applications in synthetic and cell biology.

ReverseDock: a web server for blind docking of a single ligand to multiple protein targets using AutoDock Vina.

Several platforms exist to perform molecular docking to computationally predict binders to a specific protein target from a library of ligands. The reverse, that is, docking a single ligand to various protein targets, can currently be done by very few web servers, which limits the search to a small set of pre-selected human proteins. However, the possibility to in silico predict which targets a compound identified in a high-throughput drug screen bind would help optimize and reduce the costs of the experimental workflow needed to reveal the molecular mechanism of action of a ligand. Here, we present ReverseDock, a blind docking web server based on AutoDock Vina specifically designed to allow users with no computational expertise to dock a ligand to 100 protein structures of their choice. ReverseDock increases the number and type of proteins a ligand can be docked to, making the task of in silico docking of a ligand to entire families of proteins straightforward. We envision ReverseDock will support researchers by providing the possibility to apply inverse docking computations using web browser. ReverseDock is available at: https://reversedock.biologie.uni-freiburg.de/.

Experimental Characterization of In Silico Red-Shift-Predicted iLOVL470T/Q489K and iLOVV392K/F410V/A426S Mutants.

iLOV is a flavin mononucleotide-binding fluorescent protein used for in vivo cellular imaging similar to the green fluorescent protein. To expand the range of applications of iLOV, spectrally tuned red-shifted variants are desirable to reduce phototoxicity and allow for better tissue penetration. In this report, we experimentally tested two iLOV mutants, iLOVL470T/Q489K and iLOVV392K/F410V/A426S, which were previously computationally proposed by (KhrenovaJ. Phys. Chem. B2017, 121 ( (43), ), pp 10018-10025) to have red-shifted excitation and emission spectra. While iLOVL470T/Q489K is about 20% brighter compared to the WT in vitro, it exhibits a blue shift in contrast to quantum mechanics/molecular mechanics (QM/MM) predictions. Additional optical characterization of an iLOVV392K mutant revealed that V392 is essential for cofactor binding and, accordingly, variants with V392K mutation are unable to bind to FMN. iLOVL470T/Q489K and iLOVV392K/F410V/A426S are expressed at low levels and have no detectable fluorescence in living cells, preventing their utilization in imaging applications.

Hemostatic effects of traditional Inula viscosa and Capsella bursa-pastoris plant mixture extract on rat liver parenchymal bleeding model.

BACKGROUND: Failure to achieve effective bleeding control and problems related to transfusion in liver surgery are the most common causes of post-operative mortality and morbidity. Various methods/drugs including topical hemostatic agents have been em-ployed for bleeding control in liver surgery. This study was aimed to investigate the hemostatic properties of the herb mixture extract of Inula viscosa and Capsella bursa-pastoris (IvCbp) in rat liver laceration model, which have been traditionally used as antiseptic and hemostatic agents public in Hatay/Tukey. METHODS: Thirty rats were divided into three groups equally and blood samples were taken from all rats for preoperative hemoglobin (Hb) measurements. Then, the standard liver resection model was applied to all rats. Sponge for the first rat group, Ankaferd Blood Stopper® Trend-Tech for the second rat group and IvCbp plant extract mixture for the third group were applied to resection areas for 3 minutes. Liver samples of all rats were evaluated in terms of inflammation and necrosis intensity on the 5th post-operative day. RESULTS: Post-operative Hb values were found as 11.0±1.1 g/dL in the sponge group, 11.9±2.0 g/dL in the Ankaferd group, and 14.1±1.2 g/dL in the IvCbp herb mixture group (p<0.001). In the histopathological examination, less necrosis was observed in the herb mixture group compared to the sponge and Ankaferd groups (p=0.001). In addition, no statistically significant necrosis difference was observed between sponge and Ankaferd groups. While less inflammation was observed in the herb mixture group compared to the other groups, Ankaferd group had the highest inflammation score (p<0.001). CONCLUSION: IvCbp herb mixture extract group provide effective hemostatic control, caused less Hb decrease and resulted in less inflammation and necrosis compared to Ankaferd and sponge groups in a rat liver resection model.

Putative second hit rare genetic variants in families with seemingly GBA-associated Parkinson's disease.

Rare variants in the beta-glucocerebrosidase gene (GBA1) are common genetic risk factors for alpha synucleinopathy, which often manifests clinically as GBA-associated Parkinson's disease (GBA-PD). Clinically, GBA-PD closely mimics idiopathic PD, but it may present at a younger age and often aggregates in families. Most carriers of GBA variants are, however, asymptomatic. Moreover, symptomatic PD patients without GBA variant have been reported in families with seemingly GBA-PD. These observations obscure the link between GBA variants and PD pathogenesis and point towards a role for unidentified additional genetic and/or environmental risk factors or second hits in GBA-PD. In this study, we explored whether rare genetic variants may be additional risk factors for PD in two families segregating the PD-associated GBA1 variants c.115+1G>A (ClinVar ID: 93445) and p.L444P (ClinVar ID: 4288). Our analysis identified rare genetic variants of the HSP70 co-chaperone DnaJ homolog subfamily B member 6 (DNAJB6) and lysosomal protein prosaposin (PSAP) as additional factors possibly influencing PD risk in the two families. In comparison to the wild-type proteins, variant DNAJB6 and PSAP proteins show altered functions in the context of cellular alpha-synuclein homeostasis when expressed in reporter cells. Furthermore, the segregation pattern of the rare variants in the genes encoding DNAJB6 and PSAP indicated a possible association with PD in the respective families. The occurrence of second hits or additional PD cosegregating rare variants has important implications for genetic counseling in PD families with GBA1 variant carriers and for the selection of PD patients for GBA targeted treatments.

finDr: A web server for in silico D-peptide ligand identification.

In the rapidly expanding field of peptide therapeutics, the short in vivo half-life of peptides represents a considerable limitation for drug action. D-peptides, consisting entirely of the dextrorotatory enantiomers of naturally occurring levorotatory amino acids (AAs), do not suffer from these shortcomings as they are intrinsically resistant to proteolytic degradation, resulting in a favourable pharmacokinetic profile. To experimentally identify D-peptide binders to interesting therapeutic targets, so-called mirror-image phage display is typically performed, whereby the target is synthesized in D-form and L-peptide binders are screened as in conventional phage display. This technique is extremely powerful, but it requires the synthesis of the target in D-form, which is challenging for large proteins. Here we present finDr, a novel web server for the computational identification and optimization of D-peptide ligands to any protein structure (https://findr.biologie.uni-freiburg.de/). finDr performs molecular docking to virtually screen a library of helical 12-mer peptides extracted from the RCSB Protein Data Bank (PDB) for their ability to bind to the target. In a separate, heuristic approach to search the chemical space of 12-mer peptides, finDr executes a customizable evolutionary algorithm (EA) for the de novo identification or optimization of D-peptide ligands. As a proof of principle, we demonstrate the validity of our approach to predict optimal binders to the pharmacologically relevant target phenol soluble modulin alpha 3 (PSMα3), a toxin of methicillin-resistant Staphylococcus aureus (MRSA). We validate the predictions using in vitro binding assays, supporting the success of this approach. Compared to conventional methods, finDr provides a low cost and easy-to-use alternative for the identification of D-peptide ligands against protein targets of choice without size limitation. We believe finDr will facilitate D-peptide discovery with implications in biotechnology and biomedicine.

Computation of FRAP recovery times for linker histone - chromatin binding on the basis of Brownian dynamics simulations.

Background Fluorescence recovery after photobleaching (FRAP) studies can provide kinetic information about proteins in cells. Single point mutations can significantly affect the binding kinetics of proteins and result in variations in the recovery half time (t50) measured in FRAP experiments. FRAP measurements of linker histone (LH) proteins in the cell nucleus have previously been reported by Brown et al. (2006) and Lele et al. (2006). Methods We performed Brownian dynamics (BD) simulations of the diffusional association of the wild-type and 38 single or double point mutants of the globular domain of mouse linker histone H1.0 (gH1.0) to a nucleosome. From these simulations, we calculated the bimolecular association rate constant (kon), the Gibbs binding free energy (ΔG) and the dissociation rate constant (koff) related to formation of a diffusional encounter complex between the nucleosome and the gH1.0. Results We used these parameters, after application of a correction factor to account for the effects of the crowded environment of the nucleus, to compute FRAP recovery times and curves that are in good agreement with previously published, experimentally measured FRAP recovery time courses. Conclusions Our computational analysis suggests that BD simulations can be used to predict the relative effects of single point mutations on FRAP recovery times related to protein binding. General Significance BD simulations assist in providing a detailed molecular level interpretation of FRAP data.

C-terminal eYFP fusion impairs Escherichia coli MinE function.

The Escherichia coli Min system plays an important role in the proper placement of the septum ring at mid-cell during cell division. MinE forms a pole-to-pole spatial oscillator with the membrane-bound ATPase MinD, resulting in MinD concentration being the lowest at mid-cell. MinC, the direct inhibitor of the septum initiator protein FtsZ, forms a complex with MinD at the membrane, mirroring its polar gradients. Therefore, MinC-mediated FtsZ inhibition occurs away from mid-cell. Min oscillations are often studied in living cells by time-lapse microscopy using fluorescently labelled Min proteins. Here, we show that, despite permitting oscillations to occur in a range of protein concentrations, the enhanced yellow fluorescent protein (eYFP) C-terminally fused to MinE impairs its function. Combining in vivo, in vitro and in silico approaches, we demonstrate that eYFP compromises the ability of MinE to displace MinC from MinD, to stimulate MinD ATPase activity and to directly bind to the membrane. Moreover, we reveal that MinE-eYFP is prone to aggregation. In silico analyses predict that other fluorescent proteins are also likely to compromise several functionalities of MinE, suggesting that the results presented here are not specific to eYFP.

Author Correction: Split intein-mediated selection of cells containing two plasmids using a single antibiotic.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Split intein-mediated selection of cells containing two plasmids using a single antibiotic.

To build or dissect complex pathways in bacteria and mammalian cells, it is often necessary to recur to at least two plasmids, for instance harboring orthogonal inducible promoters. Here we present SiMPl, a method based on rationally designed split enzymes and intein-mediated protein trans-splicing, allowing the selection of cells carrying two plasmids with a single antibiotic. We show that, compared to the traditional method based on two antibiotics, SiMPl increases the production of the antimicrobial non-ribosomal peptide indigoidine and the non-proteinogenic aromatic amino acid para-amino-L-phenylalanine from bacteria. Using a human T cell line, we employ SiMPl to obtain a highly pure population of cells double positive for the two chains of the T cell receptor, TCRα and TCRß, using a single antibiotic. SiMPl has profound implications for metabolic engineering and for constructing complex synthetic circuits in bacteria and mammalian cells.