Benchmarking the accuracy of structure-based binding affinity predictors on Spike-ACE2 deep mutational interaction set

Since the start of COVID-19 pandemic, a huge effort has been devoted to understanding the Spike (SARS-CoV-2)-ACE2 recognition mechanism. To this end, two deep mutational scanning studies traced the impact of all possible mutations across receptor binding domain (RBD) of Spike and catalytic domain of human ACE2. By concentrating on the interface mutations of these experimental data, we benchmarked six commonly used structure-based binding affinity predictors (FoldX, EvoEF1, MutaBind2, SSIPe, HADDOCK, and UEP). These predictors were selected based on their user-friendliness, accessibility, and speed. As a result of our benchmarking efforts, we observed that none of the methods could generate a meaningful correlation with the experimental binding data. The best correlation is achieved by FoldX (R = -0.51). When we simplified the prediction problem to a binary classification, that is, whether a mutation is enriching or depleting the binding, we showed that the highest accuracy is achieved by FoldX with a 64% success rate. Surprisingly, on this set, simple energetic scoring functions performed significantly better than the ones using extra evolutionary-based terms, as in Mutabind and SSIPe. Furthermore, we demonstrated that recent AI approaches, mmCSM-PPI and TopNetTree, yielded comparable performances to the force field-based techniques. These observations suggest plenty of room to improve the binding affinity predictors in guessing the variant-induced binding profile changes of a host-pathogen system, such as Spike-ACE2. To aid such improvements we provide our benchmarking data at https://github.com/CSB-KaracaLab/RBD-ACE2-MutBench with the option to visualize our mutant models at https://rbd-ace2-mutbench.github.io/.

Dynamic proximity interaction profiling suggests that YPEL2 is involved in cellular stress surveillance

YPEL2 is a member of the evolutionarily conserved YPEL family involved in cellular proliferation, mobility, differentiation, senescence, and death. However, the mechanism by which YPEL2, or YPEL proteins, mediates its effects is largely unknown. Proteins perform their functions in a network of proteins whose identities, amounts, and compositions change spatiotemporally in a lineage-specific manner in response to internal and external stimuli. Here, we explored interaction partners of YPEL2 by using dynamic TurboID-coupled mass spectrometry analyses to infer a function for the protein. Our results using inducible transgene expressions in COS7 cells indicate that proximity interaction partners of YPEL2 are mainly involved in RNA and mRNA metabolic processes, ribonucleoprotein complex biogenesis, regulation of gene silencing by miRNA, and cellular responses to stress. We showed that YPEL2 interacts with the RNA-binding protein ELAVL1 and the selective autophagy receptor SQSTM1. We also found that YPEL2 localizes stress granules in response to sodium arsenite, an oxidative stress inducer, which suggests that YPEL2 participates in stress granule-related processes. Establishing a point of departure in the delineation of structural/functional features of YPEL2, our results suggest that YPEL2 may be involved in stress surveillance mechanisms.

Breaking the conformational ensemble barrier: Ensemble structure modeling challenges in CASP15

For the first time, the 2022 CASP (Critical Assessment of Structure Prediction) community experiment included a section on computing multiple conformations for protein and RNA structures. There was full or partial success in reproducing the ensembles for four of the nine targets, an encouraging result. For protein structures, enhanced sampling with variations of the AlphaFold2 deep learning method was by far the most effective approach. One substantial conformational change caused by a single mutation across a complex interface was accurately reproduced. In two other assembly modeling cases, methods succeeded in sampling conformations near to the experimental ones even though environmental factors were not included in the calculations. An experimentally derived flexibility ensemble allowed a single accurate RNA structure model to be identified. Difficulties included how to handle sparse or low-resolution experimental data and the current lack of effective methods for modeling RNA/protein complexes. However, these and other obstacles appear addressable.

An alternatively spliced PD-L1 isoform PD-L1∆3, and PD-L2 expression in breast cancers: implications for eligibility scoring and immunotherapy response

Targeting PD-1/PD-L1 has shown substantial therapeutic response and unprecedented long-term durable responses in the clinic. However, several challenges persist, encompassing the prediction of treatment effectiveness and patient responses, the emergence of treatment resistance, and the necessity for additional biomarkers. Consequently, we comprehensively explored the often-overlooked isoforms of crucial immunotherapy players, leveraging transcriptomic analysis, structural modeling, and immunohistochemistry (IHC) data. Our investigation has led to the identification of an alternatively spliced isoform of PD-L1 that lacks exon 3 (PD-L1∆3) and the IgV domain required to interact with PD-1. PD-L1∆3 is expressed more than the canonical isoform in a subset of breast cancers and other TCGA tumors. Using the deep learning-based protein modeling tool AlphaFold2, we show the lack of a possible interaction between PD-L1∆3 and PD-1. In addition, we present data on the expression of an additional ligand for PD-1, PD-L2. PD-L2 expression is widespread and positively correlates with PD-L1 levels in breast and other tumors. We report enriched epithelial-mesenchymal transition (EMT) signature in high PD-L2 transcript expressing (PD-L2 > PD-L1) tumors in all breast cancer subtypes, highlighting potential crosstalk between EMT and immune evasion. Notably, the estrogen gene signature is downregulated in ER + breast tumors with high PD-L2. The data on PD-L2 IHC positivity but PD-L1 negativity in breast tumors, together with our results on PD-L1∆3, highlight the need to utilize PD-L2 and PD-L1 isoform-specific antibodies for staining patient tissue sections to offer a more precise prediction of the outcomes of PD-1/PD-L1 immunotherapy.

The impact of AI-based modeling on the accuracy of protein assembly prediction: Insights from CASP15

In CASP15, 87 predictors submitted around 11 000 models on 41 assembly targets. The community demonstrated exceptional performance in overall fold and interface contact predictions, achieving an impressive success rate of 90% (compared to 31% in CASP14). This remarkable accomplishment is largely due to the incorporation of DeepMind's AF2-Multimer approach into custom-built prediction pipelines. To evaluate the added value of participating methods, we compared the community models to the baseline AF2-Multimer predictor. In over 1/3 of cases, the community models were superior to the baseline predictor. The main reasons for this improved performance were the use of custom-built multiple sequence alignments, optimized AF2-Multimer sampling, and the manual assembly of AF2-Multimer-built subcomplexes. The best three groups, in order, are Zheng, Venclovas, and Wallner. Zheng and Venclovas reached a 73.2% success rate over all (41) cases, while Wallner attained 69.4% success rate over 36 cases. Nonetheless, challenges remain in predicting structures with weak evolutionary signals, such as nanobody-antigen, antibody-antigen, and viral complexes. Expectedly, modeling large complexes also remains challenging due to their high memory compute demands. In addition to the assembly category, we assessed the accuracy of modeling interdomain interfaces in the tertiary structure prediction targets. Models on seven targets featuring 17 unique interfaces were analyzed. Best predictors achieved a 76.5% success rate, with the UM-TBM group being the leader. In the interdomain category, we observed that the predictors faced challenges, as in the case of the assembly category, when the evolutionary signal for a given domain pair was weak or the structure was large. Overall, CASP15 witnessed unprecedented improvement in interface modeling, reflecting the AI revolution seen in CASP14.

Genetic complementation screening and molecular docking give new insight on phosphorylation-dependent Mastl kinase activation

Mastl is a mitotic kinase that is essential for error-free chromosome segregation. It is an atypical member of AGC kinase family, possessing a unique non-conserved middle region. The mechanism of Mastl activation has been studied extensively in vitro. Phosphorylation of several residues were identified to be crucial for activation. These sites correspond to T193 and T206 in the activation loop and S861 in the C-terminal tail of mouse Mastl. To date, the significance of these phosphosites was not confirmed in intact mammalian cells. Here, we utilize a genetic complementation approach to determine the essentials of mammalian Mastl kinase activation. We used tamoxifen-inducible conditional knockout mouse embryonic fibroblasts to delete endogenous Mastl and screened various mutants for their ability to complement its loss. S861A mutant was able to complement endogenous Mastl loss. In parallel, we performed computational molecular docking studies to evaluate the significance of this residue for kinase activation. Our in-depth sequence and structure analysis revealed that Mastl pS861 does not belong to a conformational state, where the phosphoresidue contributes to C-tail docking. C-tail of Mastl is relatively short and it lacks a hydrophobic (HF) motif that would otherwise help its anchoring over N-lobe, required for the final steps of kinase activation. Our results show that phosphorylation of Mastl C-tail turn motif (S861) is dispensable for kinase function in cellulo.Communicated by Ramaswamy H. Sarma.

The Impact of AI-Based Modeling on the Accuracy of Protein Assembly Prediction: Insights from CASP15

In CASP15, 87 predictors submitted around 11,000 models on 41 assembly targets. The community demonstrated exceptional performance in overall fold and interface contact prediction, achieving an impressive success rate of 90% (compared to 31% in CASP14). This remarkable accomplishment is largely due to the incorporation of DeepMind's AF2-Multimer approach into custom-built prediction pipelines. To evaluate the added value of participating methods, we compared the community models to the baseline AF2-Multimer predictor. In over 1/3 of cases the community models were superior to the baseline predictor. The main reasons for this improved performance were the use of custom-built multiple sequence alignments, optimized AF2-Multimer sampling, and the manual assembly of AF2-Multimer-built subcomplexes. The best three groups, in order, are Zheng, Venclovas and Wallner. Zheng and Venclovas reached a 73.2% success rate over all (41) cases, while Wallner attained 69.4% success rate over 36 cases. Nonetheless, challenges remain in predicting structures with weak evolutionary signals, such as nanobody-antigen, antibody-antigen, and viral complexes. Expectedly, modeling large complexes remains also challenging due to their high memory compute demands. In addition to the assembly category, we assessed the accuracy of modeling interdomain interfaces in the tertiary structure prediction targets. Models on seven targets featuring 17 unique interfaces were analyzed. Best predictors achieved the 76.5% success rate, with the UM-TBM group being the leader. In the interdomain category, we observed that the predictors faced challenges, as in the case of the assembly category, when the evolutionary signal for a given domain pair was weak or the structure was large. Overall, CASP15 witnessed unprecedented improvement in interface modeling, reflecting the AI revolution seen in CASP14.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

Rapid and efficient ambient temperature X-ray crystal structure determination at Turkish Light Source

High-resolution biomacromolecular structure determination is essential to better understand protein function and dynamics. Serial crystallography is an emerging structural biology technique which has fundamental limitations due to either sample volume requirements or immediate access to the competitive X-ray beamtime. Obtaining a high volume of well-diffracting, sufficient-size crystals while mitigating radiation damage remains a critical bottleneck of serial crystallography. As an alternative, we introduce the plate-reader module adapted for using a 72-well Terasaki plate for biomacromolecule structure determination at a convenience of a home X-ray source. We also present the first ambient temperature lysozyme structure determined at the Turkish light source (Turkish DeLight). The complete dataset was collected in 18.5 min with resolution extending to 2.39 Å and 100% completeness. Combined with our previous cryogenic structure (PDB ID: 7Y6A), the ambient temperature structure provides invaluable information about the structural dynamics of the lysozyme. Turkish DeLight provides robust and rapid ambient temperature biomacromolecular structure determination with limited radiation damage.

Cr[Formula: see text]AlN and the search for the highest temperature superconductor in the M[Formula: see text]AX family

We have developed a high-throughput computational method to predict the superconducting transition temperature in stable hexagonal M[Formula: see text]AX phases, and applied it to all the known possible choices for M (M: Sc, Ti, V, Cr, Mn, Fe, Y, Zr, Nb, Mo, Lu, Hf and Ta). We combine this with the best candidates for A (A: Al, Cu, Ge and Sn ) and X (X: C and N) from our previous work, and predict T[Formula: see text] for 60 M[Formula: see text]AX-phase materials, 53 of which have never been studied before. From all of these, we identify Cr[Formula: see text]AlN as the best candidate for the highest T[Formula: see text], and confirm its high T[Formula: see text] with more detailed density functional theory electron-phonon coupling calculations. Our detailed calculations predict [Formula: see text] = 14.8 K for Cr[Formula: see text]AlN, which is significantly higher than any [Formula: see text] value known or predicted for any material in the M[Formula: see text]AX family to date.

Differential Occupancy and Regulatory Interactions of KDM6A in Bladder Cell Lines

Epigenetic deregulation is a critical theme which needs further investigation in bladder cancer research. One of the most highly mutated genes in bladder cancer is KDM6A, which functions as an H3K27 demethylase and is one of the MLL3/4 complexes. To decipher the role of KDM6A in normal versus tumor settings, we identified the genomic landscape of KDM6A in normal, immortalized, and cancerous bladder cells. Our results showed differential KDM6A occupancy in the genes involved in cell differentiation, chromatin organization, and Notch signaling depending on the cell type and the mutation status of KDM6A. Transcription factor motif analysis revealed HES1 to be enriched at KDM6A peaks identified in the T24 bladder cancer cell line; moreover, it has a truncating mutation in KDM6A and lacks a demethylase domain. Our co-immunoprecipitation experiments revealed TLE co-repressors and HES1 as potential truncated and wild-type KDM6A interactors. With the aid of structural modeling, we explored how truncated KDM6A could interact with TLE and HES1, as well as RUNX and HHEX transcription factors. These structures provide a solid means of studying the functions of KDM6A independently of its demethylase activity. Collectively, our work provides important contributions to the understanding of KDM6A malfunction in bladder cancer.

PROT-ON: A structure-based detection of designer PROTein interface MutatiONs

The mutation-induced changes across protein-protein interfaces have often been observed to lead to severe diseases. Therefore, several computational tools have been developed to predict the impact of such mutations. Among these tools, FoldX and EvoEF1 stand out as fast and accurate alternatives. Expanding on the capabilities of these tools, we have developed the PROT-ON (PROTein-protein interface mutatiONs) framework, which aims at delivering the most critical protein interface mutations that can be used to design new protein binders. To realize this aim, PROT-ON takes the 3D coordinates of a protein dimer as an input. Then, it probes all possible interface mutations on the selected protein partner with EvoEF1 or FoldX. The calculated mutational energy landscape is statistically analyzed to find the most enriching and depleting mutations. Afterward, these extreme mutations are filtered out according to stability and optionally according to evolutionary criteria. The final remaining mutation list is presented to the user as the designer mutation set. Together with this set, PROT-ON provides several residue- and energy-based plots, portraying the synthetic energy landscape of the probed mutations. The stand-alone version of PROT-ON is deposited at https://github.com/CSB-KaracaLab/prot-on. The users can also use PROT-ON through our user-friendly web service http://proton.tools.ibg.edu.tr:8001/ (runs with EvoEF1 only). Considering its speed and the range of analysis provided, we believe that PROT-ON presents a promising means to estimate designer mutations.

Enzymatic Generation of Double-Modified AdoMet Analogues and Their Application in Cascade Reactions with Different Methyltransferases

Methyltransferases (MTases) have become an important tool for site-specific alkylation and biomolecular labelling. In biocatalytic cascades with methionine adenosyltransferases (MATs), transfer of functional moieties has been realized starting from methionine analogues and ATP. However, the widespread use of S-adenosyl-l-methionine (AdoMet) and the abundance of MTases accepting sulfonium centre modifications limit selective modification in mixtures. AdoMet analogues with additional modifications at the nucleoside moiety bear potential for acceptance by specific MTases. Here, we explored the generation of double-modified AdoMets by an engineered Methanocaldococcus jannaschii MAT (PC-MjMAT), using 19 ATP analogues in combination with two methionine analogues. This substrate screening was extended to cascade reactions and to MTase competition assays. Our results show that MTase targeting selectivity can be improved by using bulky substituents at the N⁶ of adenine. The facile access to >10 new AdoMet analogues provides the groundwork for developing MAT-MTase cascades for orthogonal biomolecular labelling.

Prediction of phonon-mediated superconductivity in new Ti-based M[Formula: see text]AX phases

A high-throughput computational method is used to predict 39 new superconductors in the Ti-based M[Formula: see text]AX phases, and the best candidates are then studied in more detail using density functional theory electron-phonon coupling calculations. The detailed calculations agree with the simple predictions, and Ti[Formula: see text]AlX (X: B, C and N) materials are predicted to have higher values of [Formula: see text] than any currently known hexagonal M[Formula: see text]AX phases. The electronic states at the Fermi level are dominated by the Ti 3d states. The choice of X (X: B, C and N) has a significant impact on the electronic density of states but not on the phonon characteristics. The electron-phonon coupling parameter for Ti[Formula: see text]AlX (X: B, C and N) was determined to be 0.685, 0.743 and 0.775 with a predicted [Formula: see text] of 7.8 K, 10.8 K and 13.0 K, respectively.

Novel biallelic variants affecting the OTU domain of the gene OTUD6B associate with severe intellectual disability syndrome and molecular dynamics simulations

Intellectual developmental disorder with dysmorphic facies, seizures, and distal limb anomalies (IDDFSDA) is an autosomal recessive multisystem disorder caused by compound heterozygous or homozygous variants in the gene OTUD6B. Herein, we describe novel pathogenic compound heterozygous variants in OTUD6B identified via whole-exome sequencing in an index case exhibited the severe IDDFSDA phenotype. The potential pathogenicity of the novel frameshift and missense variants in the index case was investigated using in silico tools. The truncating frameshift variant in one allele was predicted to undergo degradation via nonsense-mediated decay of the mRNA molecule. To predict the severity of the damage to the protein caused by the missense variant in the other allele and its effects on phenotypic severity was further investigated together with a previously reported first homozygous missense variant in the same domain in another patient with a less severe IDDFSDA phenotype using structural modeling and molecular dynamics (MD) simulations for the first time. Based on these analyzes, it is anticipated that Tyr216Cys in the earlier reported case with less severe IDDFSDA will lead to localized destabilization, whereas Ile274Arg in the presented index case with the severe IDDFSDA phenotype will lead to significant distortion in the overall fold of OTUD6B. Our findings suggest that compound LOF and ultrarare missense variants may be contribute to the underlying variability expressivity associated with this disorder. In conclusion, our findings support that the clinical severity could be related with the predicted functional severity of the variations in OTUD6B. However, additional functional studies are required.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

A C-term truncated EIF2Bγ protein encoded by an intronically polyadenylated isoform introduces unfavorable EIF2Bγ-EIF2γ interactions

Eukaryotic translation initiates upon recruitment of the EIF2-GTP·Met-tRNA_i ternary complex (TC) to the ribosomes. EIF2 (α, β, γ subunits) is a GTPase. The GDP to GTP exchange within EIF2 is facilitated by the guanine nucleotide exchange factor EIF2B (α-ε subunits). During stress-induced conditions, phosphorylation of the α-subunit of EIF2 turns EIF2 into an inhibitor of EIF2B. In turn, inhibition of EIF2B decreases TC formation and triggers the internal stress response (ISR), which determines the cell fate. Deregulated ISR has been linked to neurodegenerative disorders and cancer, positioning EIF2B as a promising therapeutic target. Hence, a better understanding of the mechanisms/factors that regulate EIF2B activity is required. Here, combining transcript and protein level analyses, we describe an intronically polyadenylated (IPA) transcript of EIF2B's γ-subunit. We show that the IPA mRNA isoform is translated into a C-terminus truncated protein. Using structural modeling, we predict that the truncated EIF2Bγ protein has unfavorable interactions with EIF2γ, leading to a potential decrease in the stability of the nonproductive EIF2:EIF2B complex. While we discovered and confirmed the IPA mRNA isoform in breast cancer cells, the expression of this isoform is not cancer-specific and is widely present in normal tissues. Overall, our data show that a truncated EIF2Bγ protein co-exists with the canonical protein and is an additional player to regulate the equilibrium between productive and nonproductive states of the EIF2:EIF2B complex. These results may have implications in stress-induced translation control in normal and disease states. Our combinatorial approach demonstrates the need to study noncanonical mRNA and protein isoforms to understand protein interactions and intricate molecular mechanisms.

Assessment of the CASP14 assembly predictions

In CASP14, 39 research groups submitted more than 2500 3D models on 22 protein complexes. In general, the community performed well in predicting the fold of the assemblies (for 80% of the targets), although it faced significant challenges in reproducing the native contacts. This is especially the case for the complexes without whole-assembly templates. The leading predictor, BAKER-experimental, used a methodology combining classical techniques (template-based modeling, protein docking) with deep learning-based contact predictions and a fold-and-dock approach. The Venclovas team achieved the runner-up position with template-based modeling and docking. By analyzing the target interfaces, we showed that the complexes with depleted charged contacts or dominating hydrophobic interactions were the most challenging ones to predict. We also demonstrated that if AlphaFold2 predictions were at hand, the interface prediction challenge could be alleviated for most of the targets. All in all, it is evident that new approaches are needed for the accurate prediction of assemblies, which undoubtedly will expand on the significant improvements in the tertiary structure prediction field.

Investigation of diagnostic use of serum anti-Müllerian hormone concentration in dioestrus and anoestrus bitches before and after ovariohysterectomy and the relationship with ovarian follicle numbers

The aims of the study were to (1) compare the serum concentration of anti-Müllerian hormone (AMH) with the number of follicles in ovaries and (2) determine the serum AMH con- centration before and after ovariohysterectomy in dioestrus and anoestrus bitches. Sixteen bitches were divided into two groups: Group I (n=8) consisted of dioestrus and group II (n=8) anoestrus bitches. The blood samples for AMH assesment were taken before ovariohysterectomy (day 0) and on day 1, 5 and 10. Both in group I and II, serum AMH concentrations on day 1 and 5 were significantly different compared to day 0 (p⟨0.05). However, the concentrations at day 10 were under the minimum detectable concentration (1.0 ng/mL) and this finding revealed that ovaries are the only source of AMH synthesis. Follicle counts were not statistically different between the groups (p>0.05). Significantly positive correlation in serum AMH with secondary follicle num- bers (r=.942, p⟨0.01), as well as negative correlation with antral follicle numbers (r=-.765, p⟨0.05) were determined in the group I. In the group II, positive correlations between serum AMH concentration and secondary follicle numbers (r=.960, p⟨0.01) and early antral follicles (r=.726, p⟨0.05) were noted. Assesment of AMH concentration seems to not only provide the diagnosis of the presence of ovaries but also correlate with the number of secondary follicles in young dioestrus and anoestrus bitches.

Predicting the Specificity- Determining Positions of Receptor Tyrosine Kinase Axl

Owing to its clinical significance, modulation of functionally relevant amino acids in protein-protein complexes has attracted a great deal of attention. To this end, many approaches have been proposed to predict the partner-selecting amino acid positions in evolutionarily close complexes. These approaches can be grouped into sequence-based machine learning and structure-based energy-driven methods. In this work, we assessed these methods’ ability to map the specificity-determining positions of Axl, a receptor tyrosine kinase involved in cancer progression and immune system diseases. For sequence-based predictions, we used SDPpred, Multi-RELIEF, and Sequence Harmony. For structure-based predictions, we utilized HADDOCK refinement and molecular dynamics simulations. As a result, we observed that (i) sequence-based methods overpredict partner-selecting residues of Axl and that (ii) combining Multi-RELIEF with HADDOCK-based predictions provides the key Axl residues, covered by the extensive molecular dynamics simulations. Expanding on these results, we propose that a sequence-structure-based approach is necessary to determine specificity-determining positions of Axl, which can guide the development of therapeutic molecules to combat Axl misregulation.

Autosomal recessive variants in TUBGCP2 alter the γ-tubulin ring complex leading to neurodevelopmental disease

Microtubules help building the cytoskeleton of neurons and other cells. Several components of the gamma-tubulin (γ-tubulin) complex have been previously reported in human neurodevelopmental diseases. We describe two siblings from a consanguineous Turkish family with dysmorphic features, developmental delay, brain malformation, and epilepsy carrying a homozygous mutation (p.Glu311Lys) in TUBGCP2 encoding the γ-tubulin complex 2 (GCP2) protein. This variant is predicted to disrupt the electrostatic interaction of GCP2 with GCP3. In primary fibroblasts carrying the variant, we observed a faint delocalization of γ-tubulin during the cell cycle but normal GCP2 protein levels. Through mass spectrometry, we observed dysregulation of multiple proteins involved in the assembly and organization of the cytoskeleton and the extracellular matrix, controlling cellular adhesion and of proteins crucial for neuronal homeostasis including axon guidance. In summary, our functional and proteomic studies link TUBGCP2 and the γ-tubulin complex to the development of the central nervous system in humans.

How Far Are We from the Rapid Prediction of Drug Resistance Arising Due to Kinase Mutations?

In all living organisms, protein kinases regulate various cell signaling events through phosphorylation. The phosphorylation occurs upon transferring an ATP's terminal phosphate to a target residue. Because of the central role of protein kinases in several proliferative pathways, point mutations occurring within the kinase's ATP-binding site can lead to a constitutively active enzyme, and ultimately, to cancer. A select set of these point mutations can also make the enzyme drug resistant toward the available kinase inhibitors. Because of technical and economical limitations, rapid experimental exploration of the impact of these mutations remains to be a challenge. This underscores the importance of kinase-ligand binding affinity prediction tools that are poised to measure the efficacy of inhibitors in the presence of kinase mutations. To this end, here, we compare the performances of six web-based scoring tools (DSX-ONLINE, KDEEP, HADDOCK2.2, PDBePISA, Pose&Rank, and PRODIGY-LIG) in assessing the impact of kinase mutations on their interactions with their inhibitors. This assessment is carried out on a new structure-based BINDKIN benchmark we compiled. BINDKIN contains wild-type and mutant structure pairs of kinase-inhibitor complexes, together with their corresponding experimental binding affinities (in the form of IC₅₀, K _d, and K _i). The performance of various web servers over BINDKIN shows that they cannot predict the binding affinities (ΔGs) of wild-type and mutant cases directly. Still, they could catch whether a mutation improves or worsens the ligand binding (ΔΔGs) where the highest Pearson's R correlation coefficient is reached by DSX-ONLINE over the K _i dataset. When homology models are used instead of K _i-associated crystal structures, DSX-ONLINE loses its predictive capacity. These results highlight that there is room to improve the available scoring functions to estimate the impact of protein kinase point mutations on inhibitor binding. The BINDKIN benchmark with all related results is freely accessible online (https://github.com/CSB-KaracaLab/BINDKIN).

Primary Unilateral Small Cell Neuroendocrine Carcinoma of the Kidney in a Dog

Primary small cell carcinomas are rare in domestic animals. A mass measuring 15 × 20 × 9 cm was detected in the left abdominal cavity of a 7.5-year-old female golden retriever. The cut surface of the excised mass showed a tumour replacing the left kidney. Microscopically, the mass was composed of polymorphic, small basophilic cells with a high nuclear to cytoplasmic ratio and round, oval or short slender fusiform nuclei with condensed or finely granular chromatin, absent or inconspicuous nucleoli, and scant, faintly eosinophilic cytoplasm with poorly defined cytoplasmic borders. Immunohistochemically, most of the neoplastic cells were immunoreactive for thyroid transcription factor 1 and CD56, moderately positive for vimentin and weakly or sparsely labelled for chromogranin A, synaptophysin, Wilms' tumour 1 protein, neuron-specific enolase, pan-cytokeratin (CK) AE1/AE3 and epithelial membrane antigen. The tumour cells were negative for glial fibrillary acidic protein (GFAP), ionized calcium-binding adapter molecule 1, CK7, CK20, CD3, CD45 and CD99. These findings indicated a neuroendocrine origin of the tumour. To the best of author's knowledge, this is the first report of a small cell neuroendocrine carcinoma originating as a primary tumour in the kidney of a dog.

More Favorable Palmitic Acid Over Palmitoleic Acid Modification of Wnt3 Ensures Its Localization and Activity in Plasma Membrane Domains

While the lateral organization of plasma membrane components has been shown to control binding of Wnt ligands to their receptors preferentially in the ordered membrane domains, the role of posttranslational lipid modification of Wnt on this selective binding is unknown. Here, we identify that the canonical Wnt is presumably acylated by palmitic acid, a saturated 16-carbon fatty acid, at a conserved serine residue. Acylation of Wnt3 is dispensable for its secretion and binding to Fz8 while it is essential for Wnt3's proper binding and domain-like diffusion in the ordered membrane domains. We further unravel that non-palmitoylated Wnt3 is unable to activate Wnt/β-catenin signaling either in zebrafish embryos or in mammalian cells. Based on these results, we propose that the lipidation of canonical Wnt, presumably by a saturated fatty acid, determines its competence in interacting with the receptors in the appropriate domains of the plasma membrane, ultimately keeping the signaling activity under control.

In Silico Phylogenetic and Structural Analyses of Plant Endogenous Danger Signaling Molecules upon Stress

The plant innate immune system has two major branches, the pathogen-triggered immunity and the effector-triggered immunity (ETI). The effectors are molecules released by plant attackers to evade host immunity. In addition to the foreign intruders, plants possess endogenous instigators produced in response to general cellular injury termed as damage-associated molecular patterns (DAMPs). In plants, DAMPs or alarmins are released by damaged, stressed, or dying cells following abiotic stress such as radiation, oxidative and drought stresses. In turn, a cascade of downstream signaling events is initiated leading to the upregulation of defense or response-related genes. In the present study, we have investigated more thoroughly the conservation status of the molecular mechanisms implicated in the danger signaling primarily in plants. Towards this direction, we have performed in silico phylogenetic and structural analyses of the associated biomolecules in taxonomically diverse plant species. On the basis of our results, the defense mechanisms appear to be largely conserved within the plant kingdom. Of note, the sequence and/or function of several components of these mechanisms was found to be conserved in animals, as well. At the same time, the molecules involved in plant defense were found to form a dense protein-protein interaction (PPi) network, suggesting a crosstalk between the various defense mechanisms to a variety of stresses, like oxidative stress.

Novel mutations in KMT2B offer pathophysiological insights into childhood-onset progressive dystonia

Rapid progress has recently been made in the elucidation of the genetic basis of childhood-onset inherited generalized dystonia (IGD) due to the implementation of genomic sequencing methodologies. We identified four patients with childhood-onset IGD harboring novel disease-causing mutations in lysine-specific histone methyltransferase 2B gene (KMT2B) by whole-exome sequencing. The main focus of this paper is to gain novel pathophysiological insights through understanding the molecular consequences of these mutations. The disease course is mostly progressive, evolving from lower limbs into generalized dystonia, which could be associated with dysarthria, dysphonia, intellectual disability, orofacial dyskinesia, and sometimes distinct dysmorphic facial features. In two patients, motor performances improved after bilateral implantation of deep brain stimulation in the globus pallidus internus (GPi-DBS). Pharmacotherapy with trihexyphenidyl reduced dystonia in two patients. We discovered three novel KMT2B mutations. Our analyses revealed that the mutation in patient 1 (c.7463A > G, p.Y2488C) is localized in the highly conserved FYRC domain of KMT2B. This mutation holds the potential to alter the inter-domain FYR interactions, which could lead to KMT2B instability. The mutations in patients 2 and 3 (c.3596_3697insC, p.M1202Dfs*22; c.4229delA, p.Q1410Rfs*12) lead to predicted unstable transcripts, likely to be subject to degradation by non-sense-mediated decay. Childhood-onset progressive dystonia with orofacial involvement is one of the main clinical manifestations of KMT2B mutations. In all, 26% (18/69) of the reported cases have T2 signal alterations of the globus pallidus internus, mostly at a younger age. Anticholinergic medication and GPi-DBS are promising treatment options and shall be considered early.

Transposase-DNA Complex Structures Reveal Mechanisms for Conjugative Transposition of Antibiotic Resistance

Conjugative transposition drives the emergence of multidrug resistance in diverse bacterial pathogens, yet the mechanisms are poorly characterized. The Tn1549 conjugative transposon propagates resistance to the antibiotic vancomycin used for severe drug-resistant infections. Here, we present four high-resolution structures of the conserved Y-transposase of Tn1549 complexed with circular transposon DNA intermediates. The structures reveal individual transposition steps and explain how specific DNA distortion and cleavage mechanisms enable DNA strand exchange with an absolute minimum homology requirement. This appears to uniquely allow Tn916-like conjugative transposons to bypass DNA homology and insert into diverse genomic sites, expanding gene transfer. We further uncover a structural regulatory mechanism that prevents premature cleavage of the transposon DNA before a suitable target DNA is found and generate a peptide antagonist that interferes with the transposase-DNA structure to block transposition. Our results reveal mechanistic principles of conjugative transposition that could help control the spread of antibiotic resistance genes.

•

Antibiotic resistance-carrying conjugative transposon integrase structure revealed

•

DNA distortion and special cleavage site allow insertion into diverse genomic sites

•

Key structural features are conserved among numerous conjugative transposons

•

Structures uncover auto-inhibition, allowing transposition antagonist design

Structural snapshots of Xer recombination reveal activation by synaptic complex remodeling and DNA bending

Bacterial Xer site-specific recombinases play an essential genome maintenance role by unlinking chromosome multimers, but their mechanism of action has remained structurally uncharacterized. Here, we present two high-resolution structures of Helicobacter pylori XerH with its recombination site DNA dif_H, representing pre-cleavage and post-cleavage synaptic intermediates in the recombination pathway. The structures reveal that activation of DNA strand cleavage and rejoining involves large conformational changes and DNA bending, suggesting how interaction with the cell division protein FtsK may license recombination at the septum. Together with biochemical and in vivo analysis, our structures also reveal how a small sequence asymmetry in dif_H defines protein conformation in the synaptic complex and orchestrates the order of DNA strand exchanges. Our results provide insights into the catalytic mechanism of Xer recombination and a model for regulation of recombination activity during cell division.

DOI:http://dx.doi.org/10.7554/eLife.19706.001

Similar to humans, bacteria store their genetic material in the form of DNA and arrange it into structures called chromosomes. In fact, most bacteria have a single circular chromosome. Bacteria multiply by simply dividing in two, and before that happens they must replicate their DNA so that each of the newly formed cells receives one copy of the chromosome.

Occasionally, mistakes during the DNA replication process can cause the two chromosomes to become tangled with each other; this prevents them from separating into the newly formed cells. For instance, the chromosomes can become physically connected like links in a chain, or merge into one long string. This kind of tangling can result in cell death, so bacteria encode enzymes called Xer recombinases that can untangle chromosomes. These enzymes separate the chromosomes by cutting and rejoining the DNA strands in a process known as Xer recombination.

Although Xer recombinases have been studied in quite some detail, many questions remain unanswered about how they work. How do Xer recombinases interact with DNA? How do they ensure they only work on tangled chromosomes? And how does a protein called FtsK ensure that Xer recombination takes place at the correct time and place?

Bebel et al. have now studied the Xer recombinase from a bacterium called Helicobacter pylori, which causes stomach ulcers, using a technique called X-ray crystallography. This enabled the three-dimensional structure of the Xer recombinase to be visualized as it interacted with DNA to form a Xer-DNA complex. Structures of the enzyme before and after it cut the DNA show that Xer-DNA complexes first assemble in an inactive state and are then activated by large conformational changes that make the DNA bend. Bebel et al. propose that the FtsK protein might trigger these changes and help to bend the DNA as it activates Xer recombination.

Further work showed that the structures could be used to model and understand Xer recombinases from other species of bacteria. The next step is to analyze how FtsK activates Xer recombinases and to see if this process is universal amongst bacteria. Understanding how this process can be interrupted could help to develop new drugs that can kill harmful bacteria.

DOI:http://dx.doi.org/10.7554/eLife.19706.002

Sense and simplicity in HADDOCK scoring: Lessons from CASP-CAPRI round 1

Our information-driven docking approach HADDOCK is a consistent top predictor and scorer since the start of its participation in the CAPRI community-wide experiment. This sustained performance is due, in part, to its ability to integrate experimental data and/or bioinformatics information into the modelling process, and also to the overall robustness of the scoring function used to assess and rank the predictions. In the CASP-CAPRI Round 1 scoring experiment we successfully selected acceptable/medium quality models for 18/14 of the 25 targets - a top-ranking performance among all scorers. Considering that for only 20 targets acceptable models were generated by the community, our effective success rate reaches as high as 90% (18/20). This was achieved using the standard HADDOCK scoring function, which, thirteen years after its original publication, still consists of a simple linear combination of intermolecular van der Waals and Coulomb electrostatics energies and an empirically derived desolvation energy term. Despite its simplicity, this scoring function makes sense from a physico-chemical perspective, encoding key aspects of biomolecular recognition. In addition to its success in the scoring experiment, the HADDOCK server takes the first place in the server prediction category, with 16 successful predictions. Much like our scoring protocol, because of the limited time per target, the predictions relied mainly on either an ab initio center-of-mass and symmetry restrained protocol, or on a template-based approach whenever applicable. These results underline the success of our simple but sensible prediction and scoring scheme. Proteins 2017; 85:417-423. © 2016 Wiley Periodicals, Inc.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

Human genome meeting 2016 : Houston, TX, USA. 28 February - 2 March 2016

O1 The metabolomics approach to autism: identification of biomarkers for early detection of autism spectrum disorder A. K. Srivastava, Y. Wang, R. Huang, C. Skinner, T. Thompson, L. Pollard, T. Wood, F. Luo, R. Stevenson O2 Phenome-wide association study for smoking- and drinking-associated genes in 26,394 American women with African, Asian, European, and Hispanic descents R. Polimanti, J. Gelernter O3 Effects of prenatal environment, genotype and DNA methylation on birth weight and subsequent postnatal outcomes: findings from GUSTO, an Asian birth cohort X. Lin, I. Y. Lim, Y. Wu, A. L. Teh, L. Chen, I. M. Aris, S. E. Soh, M. T. Tint, J. L. MacIsaac, F. Yap, K. Kwek, S. M. Saw, M. S. Kobor, M. J. Meaney, K. M. Godfrey, Y. S. Chong, J. D. Holbrook, Y. S. Lee, P. D. Gluckman, N. Karnani, GUSTO study group O4 High-throughput identification of specific qt interval modulating enhancers at the SCN5A locus A. Kapoor, D. Lee, A. Chakravarti O5 Identification of extracellular matrix components inducing cancer cell migration in the supernatant of cultivated mesenchymal stem cells C. Maercker, F. Graf, M. Boutros O6 Single cell allele specific expression (ASE) IN T21 and common trisomies: a novel approach to understand DOWN syndrome and other aneuploidies G. Stamoulis, F. Santoni, P. Makrythanasis, A. Letourneau, M. Guipponi, N. Panousis, M. Garieri, P. Ribaux, E. Falconnet, C. Borel, S. E. Antonarakis O7 Role of microRNA in LCL to IPSC reprogramming S. Kumar, J. Curran, J. Blangero O8 Multiple enhancer variants disrupt gene regulatory network in Hirschsprung disease S. Chatterjee, A. Kapoor, J. Akiyama, D. Auer, C. Berrios, L. Pennacchio, A. Chakravarti O9 Metabolomic profiling for the diagnosis of neurometabolic disorders T. R. Donti, G. Cappuccio, M. Miller, P. Atwal, A. Kennedy, A. Cardon, C. Bacino, L. Emrick, J. Hertecant, F. Baumer, B. Porter, M. Bainbridge, P. Bonnen, B. Graham, R. Sutton, Q. Sun, S. Elsea O10 A novel causal methylation network approach to Alzheimer’s disease Z. Hu, P. Wang, Y. Zhu, J. Zhao, M. Xiong, David A Bennett O11 A microRNA signature identifies subtypes of triple-negative breast cancer and reveals MIR-342-3P as regulator of a lactate metabolic pathway A. Hidalgo-Miranda, S. Romero-Cordoba, S. Rodriguez-Cuevas, R. Rebollar-Vega, E. Tagliabue, M. Iorio, E. D’Ippolito, S. Baroni O12 Transcriptome analysis identifies genes, enhancer RNAs and repetitive elements that are recurrently deregulated across multiple cancer types B. Kaczkowski, Y. Tanaka, H. Kawaji, A. Sandelin, R. Andersson, M. Itoh, T. Lassmann, the FANTOM5 consortium, Y. Hayashizaki, P. Carninci, A. R. R. Forrest O13 Elevated mutation and widespread loss of constraint at regulatory and architectural binding sites across 11 tumour types C. A. Semple O14 Exome sequencing provides evidence of pathogenicity for genes implicated in colorectal cancer E. A. Rosenthal, B. Shirts, L. Amendola, C. Gallego, M. Horike-Pyne, A. Burt, P. Robertson, P. Beyers, C. Nefcy, D. Veenstra, F. Hisama, R. Bennett, M. Dorschner, D. Nickerson, J. Smith, K. Patterson, D. Crosslin, R. Nassir, N. Zubair, T. Harrison, U. Peters, G. Jarvik, NHLBI GO Exome Sequencing Project O15 The tandem duplicator phenotype as a distinct genomic configuration in cancer F. Menghi, K. Inaki, X. Woo, P. Kumar, K. Grzeda, A. Malhotra, H. Kim, D. Ucar, P. Shreckengast, K. Karuturi, J. Keck, J. Chuang, E. T. Liu O16 Modeling genetic interactions associated with molecular subtypes of breast cancer B. Ji, A. Tyler, G. Ananda, G. Carter O17 Recurrent somatic mutation in the MYC associated factor X in brain tumors H. Nikbakht, M. Montagne, M. Zeinieh, A. Harutyunyan, M. Mcconechy, N. Jabado, P. Lavigne, J. Majewski O18 Predictive biomarkers to metastatic pancreatic cancer treatment J. B. Goldstein, M. Overman, G. Varadhachary, R. Shroff, R. Wolff, M. Javle, A. Futreal, D. Fogelman O19 DDIT4 gene expression as a prognostic marker in several malignant tumors L. Bravo, W. Fajardo, H. Gomez, C. Castaneda, C. Rolfo, J. A. Pinto O20 Spatial organization of the genome and genomic alterations in human cancers K. C. Akdemir, L. Chin, A. Futreal, ICGC PCAWG Structural Alterations Group O21 Landscape of targeted therapies in solid tumors S. Patterson, C. Statz, S. Mockus O22 Genomic analysis reveals novel drivers and progression pathways in skin basal cell carcinoma S. N. Nikolaev, X. I. Bonilla, L. Parmentier, B. King, F. Bezrukov, G. Kaya, V. Zoete, V. Seplyarskiy, H. Sharpe, T. McKee, A. Letourneau, P. Ribaux, K. Popadin, N. Basset-Seguin, R. Ben Chaabene, F. Santoni, M. Andrianova, M. Guipponi, M. Garieri, C. Verdan, K. Grosdemange, O. Sumara, M. Eilers, I. Aifantis, O. Michielin, F. de Sauvage, S. Antonarakis O23 Identification of differential biomarkers of hepatocellular carcinoma and cholangiocarcinoma via transcriptome microarray meta-analysis S. Likhitrattanapisal O24 Clinical validity and actionability of multigene tests for hereditary cancers in a large multi-center study S. Lincoln, A. Kurian, A. Desmond, S. Yang, Y. Kobayashi, J. Ford, L. Ellisen O25 Correlation with tumor ploidy status is essential for correct determination of genome-wide copy number changes by SNP array T. L. Peters, K. R. Alvarez, E. F. Hollingsworth, D. H. Lopez-Terrada O26 Nanochannel based next-generation mapping for interrogation of clinically relevant structural variation A. Hastie, Z. Dzakula, A. W. Pang, E. T. Lam, T. Anantharaman, M. Saghbini, H. Cao, BioNano Genomics O27 Mutation spectrum in a pulmonary arterial hypertension (PAH) cohort and identification of associated truncating mutations in TBX4 C. Gonzaga-Jauregui, L. Ma, A. King, E. Berman Rosenzweig, U. Krishnan, J. G. Reid, J. D. Overton, F. Dewey, W. K. Chung O28 NORTH CAROLINA macular dystrophy (MCDR1): mutations found affecting PRDM13 K. Small, A. DeLuca, F. Cremers, R. A. Lewis, V. Puech, B. Bakall, R. Silva-Garcia, K. Rohrschneider, M. Leys, F. S. Shaya, E. Stone O29 PhenoDB and genematcher, solving unsolved whole exome sequencing data N. L. Sobreira, F. Schiettecatte, H. Ling, E. Pugh, D. Witmer, K. Hetrick, P. Zhang, K. Doheny, D. Valle, A. Hamosh O30 Baylor-Johns Hopkins Center for Mendelian genomics: a four year review S. N. Jhangiani, Z. Coban Akdemir, M. N. Bainbridge, W. Charng, W. Wiszniewski, T. Gambin, E. Karaca, Y. Bayram, M. K. Eldomery, J. Posey, H. Doddapaneni, J. Hu, V. R. Sutton, D. M. Muzny, E. A. Boerwinkle, D. Valle, J. R. Lupski, R. A. Gibbs O31 Using read overlap assembly to accurately identify structural genetic differences in an ashkenazi jewish trio S. Shekar, W. Salerno, A. English, A. Mangubat, J. Bruestle O32 Legal interoperability: a sine qua non for international data sharing A. Thorogood, B. M. Knoppers, Global Alliance for Genomics and Health - Regulatory and Ethics Working Group O33 High throughput screening platform of competent sineups: that can enhance translation activities of therapeutic target H. Takahashi, K. R. Nitta, A. Kozhuharova, A. M. Suzuki, H. Sharma, D. Cotella, C. Santoro, S. Zucchelli, S. Gustincich, P. Carninci O34 The undiagnosed diseases network international (UDNI): clinical and laboratory research to meet patient needs J. J. Mulvihill, G. Baynam, W. Gahl, S. C. Groft, K. Kosaki, P. Lasko, B. Melegh, D. Taruscio O36 Performance of computational algorithms in pathogenicity predictions for activating variants in oncogenes versus loss of function mutations in tumor suppressor genes R. Ghosh, S. Plon O37 Identification and electronic health record incorporation of clinically actionable pharmacogenomic variants using prospective targeted sequencing S. Scherer, X. Qin, R. Sanghvi, K. Walker, T. Chiang, D. Muzny, L. Wang, J. Black, E. Boerwinkle, R. Weinshilboum, R. Gibbs O38 Melanoma reprogramming state correlates with response to CTLA-4 blockade in metastatic melanoma T. Karpinets, T. Calderone, K. Wani, X. Yu, C. Creasy, C. Haymaker, M. Forget, V. Nanda, J. Roszik, J. Wargo, L. Haydu, X. Song, A. Lazar, J. Gershenwald, M. Davies, C. Bernatchez, J. Zhang, A. Futreal, S. Woodman O39 Data-driven refinement of complex disease classification from integration of heterogeneous functional genomics data in GeneWeaver E. J. Chesler, T. Reynolds, J. A. Bubier, C. Phillips, M. A. Langston, E. J. Baker O40 A general statistic framework for genome-based disease risk prediction M. Xiong, L. Ma, N. Lin, C. Amos O41 Integrative large-scale causal network analysis of imaging and genomic data and its application in schizophrenia studies N. Lin, P. Wang, Y. Zhu, J. Zhao, V. Calhoun, M. Xiong O42 Big data and NGS data analysis: the cloud to the rescue O. Dobretsberger, M. Egger, F. Leimgruber O43 Cpipe: a convergent clinical exome pipeline specialised for targeted sequencing S. Sadedin, A. Oshlack, Melbourne Genomics Health Alliance O44 A Bayesian classification of biomedical images using feature extraction from deep neural networks implemented on lung cancer data V. A. A. Antonio, N. Ono, Clark Kendrick C. Go O45 MAV-SEQ: an interactive platform for the Management, Analysis, and Visualization of sequence data Z. Ahmed, M. Bolisetty, S. Zeeshan, E. Anguiano, D. Ucar O47 Allele specific enhancer in EPAS1 intronic regions may contribute to high altitude adaptation of Tibetans C. Zeng, J. Shao O48 Nanochannel based next-generation mapping for structural variation detection and comparison in trios and populations H. Cao, A. Hastie, A. W. Pang, E. T. Lam, T. Liang, K. Pham, M. Saghbini, Z. Dzakula O49 Archaic introgression in indigenous populations of Malaysia revealed by whole genome sequencing Y. Chee-Wei, L. Dongsheng, W. Lai-Ping, D. Lian, R. O. Twee Hee, Y. Yunus, F. Aghakhanian, S. S. Mokhtar, C. V. Lok-Yung, J. Bhak, M. Phipps, X. Shuhua, T. Yik-Ying, V. Kumar, H. Boon-Peng O50 Breast and ovarian cancer prevention: is it time for population-based mutation screening of high risk genes? I. Campbell, M.-A. Young, P. James, Lifepool O53 Comprehensive coverage from low DNA input using novel NGS library preparation methods for WGS and WGBS C. Schumacher, S. Sandhu, T. Harkins, V. Makarov O54 Methods for large scale construction of robust PCR-free libraries for sequencing on Illumina HiSeqX platform H. DoddapaneniR. Glenn, Z. Momin, B. Dilrukshi, H. Chao, Q. Meng, B. Gudenkauf, R. Kshitij, J. Jayaseelan, C. Nessner, S. Lee, K. Blankenberg, L. Lewis, J. Hu, Y. Han, H. Dinh, S. Jireh, K. Walker, E. Boerwinkle, D. Muzny, R. Gibbs O55 Rapid capture methods for clinical sequencing J. Hu, K. Walker, C. Buhay, X. Liu, Q. Wang, R. Sanghvi, H. Doddapaneni, Y. Ding, N. Veeraraghavan, Y. Yang, E. Boerwinkle, A. L. Beaudet, C. M. Eng, D. M. Muzny, R. A. Gibbs O56 A diploid personal human genome model for better genomes from diverse sequence data K. C. C. Worley, Y. Liu, D. S. T. Hughes, S. C. Murali, R. A. Harris, A. C. English, X. Qin, O. A. Hampton, P. Larsen, C. Beck, Y. Han, M. Wang, H. Doddapaneni, C. L. Kovar, W. J. Salerno, A. Yoder, S. Richards, J. Rogers, J. R. Lupski, D. M. Muzny, R. A. Gibbs O57 Development of PacBio long range capture for detection of pathogenic structural variants Q. Meng, M. Bainbridge, M. Wang, H. Doddapaneni, Y. Han, D. Muzny, R. Gibbs O58 Rhesus macaques exhibit more non-synonymous variation but greater impact of purifying selection than humans R. A. Harris, M. Raveenedran, C. Xue, M. Dahdouli, L. Cox, G. Fan, B. Ferguson, J. Hovarth, Z. Johnson, S. Kanthaswamy, M. Kubisch, M. Platt, D. Smith, E. Vallender, R. Wiseman, X. Liu, J. Below, D. Muzny, R. Gibbs, F. Yu, J. Rogers O59 Assessing RNA structure disruption induced by single-nucleotide variation J. Lin, Y. Zhang, Z. Ouyang P1 A meta-analysis of genome-wide association studies of mitochondrial dna copy number A. Moore, Z. Wang, J. Hofmann, M. Purdue, R. Stolzenberg-Solomon, S. Weinstein, D. Albanes, C.-S. Liu, W.-L. Cheng, T.-T. Lin, Q. Lan, N. Rothman, S. Berndt P2 Missense polymorphic genetic combinations underlying down syndrome susceptibility E. S. Chen P4 The evaluation of alteration of ELAM-1 expression in the endometriosis patients H. Bahrami, A. Khoshzaban, S. Heidari Keshal P5 Obesity and the incidence of apolipoprotein E polymorphisms in an assorted population from Saudi Arabia population K. K. R. Alharbi P6 Genome-associated personalized antithrombotical therapy for patients with high risk of thrombosis and bleeding M. Zhalbinova, A. Akilzhanova, S. Rakhimova, M. Bekbosynova, S. Myrzakhmetova P7 Frequency of Xmn1 polymorphism among sickle cell carrier cases in UAE population M. Matar P8 Differentiating inflammatory bowel diseases by using genomic data: dimension of the problem and network organization N. Mili, R. Molinari, Y. Ma, S. Guerrier P9 Vulnerability of genetic variants to the risk of autism among Saudi children N. Elhawary, M. Tayeb, N. Bogari, N. Qotb P10 Chromatin profiles from ex vivo purified dopaminergic neurons establish a promising model to support studies of neurological function and dysfunction S. A. McClymont, P. W. Hook, L. A. Goff, A. McCallion P11 Utilization of a sensitized chemical mutagenesis screen to identify genetic modifiers of retinal dysplasia in homozygous Nr2e3^rd7 mice Y. Kong, J. R. Charette, W. L. Hicks, J. K. Naggert, L. Zhao, P. M. Nishina P12 Ion torrent next generation sequencing of recessive polycystic kidney disease in Saudi patients B. M. Edrees, M. Athar, F. A. Al-Allaf, M. M. Taher, W. Khan, A. Bouazzaoui, N. A. Harbi, R. Safar, H. Al-Edressi, A. Anazi, N. Altayeb, M. A. Ahmed, K. Alansary, Z. Abduljaleel P13 Digital expression profiling of Purkinje neurons and dendrites in different subcellular compartments A. Kratz, P. Beguin, S. Poulain, M. Kaneko, C. Takahiko, A. Matsunaga, S. Kato, A. M. Suzuki, N. Bertin, T. Lassmann, R. Vigot, P. Carninci, C. Plessy, T. Launey P14 The evolution of imperfection and imperfection of evolution: the functional and functionless fractions of the human genome D. Graur P16 Species-independent identification of known and novel recurrent genomic entities in multiple cancer patients J. Friis-Nielsen, J. M. Izarzugaza, S. Brunak P18 Discovery of active gene modules which are densely conserved across multiple cancer types reveal their prognostic power and mutually exclusive mutation patterns B. S. Soibam P19 Whole exome sequencing of dysplastic leukoplakia tissue indicates sequential accumulation of somatic mutations from oral precancer to cancer D. Das, N. Biswas, S. Das, S. Sarkar, A. Maitra, C. Panda, P. Majumder P21 Epigenetic mechanisms of carcinogensis by hereditary breast cancer genes J. J. Gruber, N. Jaeger, M. Snyder P22 RNA direct: a novel RNA enrichment strategy applied to transcripts associated with solid tumors K. Patel, S. Bowman, T. Davis, D. Kraushaar, A. Emerman, S. Russello, N. Henig, C. Hendrickson P23 RNA sequencing identifies gene mutations for neuroblastoma K. Zhang P24 Participation of SFRP1 in the modulation of TMPRSS2-ERG fusion gene in prostate cancer cell lines M. Rodriguez-Dorantes, C. D. Cruz-Hernandez, C. D. P. Garcia-Tobilla, S. Solorzano-Rosales P25 Targeted Methylation Sequencing of Prostate Cancer N. Jäger, J. Chen, R. Haile, M. Hitchins, J. D. Brooks, M. Snyder P26 Mutant TPMT alleles in children with acute lymphoblastic leukemia from México City and Yucatán, Mexico S. Jiménez-Morales, M. Ramírez, J. Nuñez, V. Bekker, Y. Leal, E. Jiménez, A. Medina, A. Hidalgo, J. Mejía P28 Genetic modifiers of Alström syndrome J. Naggert, G. B. Collin, K. DeMauro, R. Hanusek, P. M. Nishina P31 Association of genomic variants with the occurrence of angiotensin-converting-enzyme inhibitor (ACEI)-induced coughing among Filipinos E. M. Cutiongco De La Paz, R. Sy, J. Nevado, P. Reganit, L. Santos, J. D. Magno, F. E. Punzalan , D. Ona , E. Llanes, R. L. Santos-Cortes , R. Tiongco, J. Aherrera, L. Abrahan, P. Pagauitan-Alan; Philippine Cardiogenomics Study Group P32 The use of “humanized” mouse models to validate disease association of a de novo GARS variant and to test a novel gene therapy strategy for Charcot-Marie-Tooth disease type 2D K. H. Morelli, J. S. Domire, N. Pyne, S. Harper, R. Burgess P34 Molecular regulation of chondrogenic human induced pluripotent stem cells M. A. Gari, A. Dallol, H. Alsehli, A. Gari, M. Gari, A. Abuzenadah P35 Molecular profiling of hematologic malignancies: implementation of a variant assessment algorithm for next generation sequencing data analysis and clinical reporting M. Thomas, M. Sukhai, S. Garg, M. Misyura, T. Zhang, A. Schuh, T. Stockley, S. Kamel-Reid P36 Accessing genomic evidence for clinical variants at NCBI S. Sherry, C. Xiao, D. Slotta, K. Rodarmer, M. Feolo, M. Kimelman, G. Godynskiy, C. O’Sullivan, E. Yaschenko P37 NGS-SWIFT: a cloud-based variant analysis framework using control-accessed sequencing data from DBGAP/SRA C. Xiao, E. Yaschenko, S. Sherry P38 Computational assessment of drug induced hepatotoxicity through gene expression profiling C. Rangel-Escareño, H. Rueda-Zarate P40 Flowr: robust and efficient pipelines using a simple language-agnostic approach;ultraseq; fast modular pipeline for somatic variation calling using flowr S. Seth, S. Amin, X. Song, X. Mao, H. Sun, R. G. Verhaak, A. Futreal, J. Zhang P41 Applying “Big data” technologies to the rapid analysis of heterogenous large cohort data S. J. Whiite, T. Chiang, A. English, J. Farek, Z. Kahn, W. Salerno, N. Veeraraghavan, E. Boerwinkle, R. Gibbs P42 FANTOM5 web resource for the large-scale genome-wide transcription start site activity profiles of wide-range of mammalian cells T. Kasukawa, M. Lizio, J. Harshbarger, S. Hisashi, J. Severin, A. Imad, S. Sahin, T. C. Freeman, K. Baillie, A. Sandelin, P. Carninci, A. R. R. Forrest, H. Kawaji, The FANTOM Consortium P43 Rapid and scalable typing of structural variants for disease cohorts W. Salerno, A. English, S. N. Shekar, A. Mangubat, J. Bruestle, E. Boerwinkle, R. A. Gibbs P44 Polymorphism of glutathione S-transferases and sulphotransferases genes in an Arab population A. H. Salem, M. Ali, A. Ibrahim, M. Ibrahim P46 Genetic divergence of CYP3A5*3 pharmacogenomic marker for native and admixed Mexican populations J. C. Fernandez-Lopez, V. Bonifaz-Peña, C. Rangel-Escareño, A. Hidalgo-Miranda, A. V. Contreras P47 Whole exome sequence meta-analysis of 13 white blood cell, red blood cell, and platelet traits L. Polfus, CHARGE and NHLBI Exome Sequence Project Working Groups P48 Association of adipoq gene with type 2 diabetes and related phenotypes in african american men and women: The jackson heart study S. Davis, R. Xu, S. Gebeab, P Riestra, A Gaye, R. Khan, J. Wilson, A. Bidulescu P49 Common variants in casr gene are associated with serum calcium levels in koreans S. H. Jung, N. Vinayagamoorthy, S. H. Yim, Y. J. Chung P50 Inference of multiple-wave population admixture by modeling decay of linkage disequilibrium with multiple exponential functions Y. Zhou, S. Xu P51 A Bayesian framework for generalized linear mixed models in genome-wide association studies X. Wang, V. Philip, G. Carter P52 Targeted sequencing approach for the identification of the genetic causes of hereditary hearing impairment A. A. Abuzenadah, M. Gari, R. Turki, A. Dallol P53 Identification of enhancer sequences by ATAC-seq open chromatin profiling A. Uyar, A. Kaygun, S. Zaman, E. Marquez, J. George, D. Ucar P54 Direct enrichment for the rapid preparation of targeted NGS libraries C. L. Hendrickson, A. Emerman, D. Kraushaar, S. Bowman, N. Henig, T. Davis, S. Russello, K. Patel P56 Performance of the Agilent D5000 and High Sensitivity D5000 ScreenTape assays for the Agilent 4200 Tapestation System R. Nitsche, L. Prieto-Lafuente P57 ClinVar: a multi-source archive for variant interpretation M. Landrum, J. Lee, W. Rubinstein, D. Maglott P59 Association of functional variants and protein physical interactions of human MUTY homolog linked with familial adenomatous polyposis and colorectal cancer syndrome Z. Abduljaleel, W. Khan, F. A. Al-Allaf, M. Athar , M. M. Taher, N. Shahzad P60 Modification of the microbiom constitution in the gut using chicken IgY antibodies resulted in a reduction of acute graft-versus-host disease after experimental bone marrow transplantation A. Bouazzaoui, E. Huber, A. Dan, F. A. Al-Allaf, W. Herr, G. Sprotte, J. Köstler, A. Hiergeist, A. Gessner, R. Andreesen, E. Holler P61 Compound heterozygous mutation in the LDLR gene in Saudi patients suffering severe hypercholesterolemia F. Al-Allaf, A. Alashwal, Z. Abduljaleel, M. Taher, A. Bouazzaoui, H. Abalkhail, A. Al-Allaf, R. Bamardadh, M. Athar

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.

EVALUATION OF PREIMPLANTATION GENETIC ANEUPLOIDY SCREENING CASES AT A REFERENCE GENETICS CENTER: 10 YEARS' EXPERIENCE

The aim of this study is to review and evaluate our preimplantation genetic screening (PGS) records in terms of their demographic data, indications, cytogenetic results, pregnancy outcomes and discuss these findings in different aspects. PGS was performed in a total of 84 couples (87 cycles) between the period 2005 to 2015. Biopsied blastomeres from embryos on day 3 were fixed and fluorescence in situ hybridization was carried out for chromosomes 13, 16, 18, 21, 22, X and Y depending on the indication. The diagnostic and clinical data were retrospectively evaluated. A total of 450 blastomeres were biopsied. Ninety-eight of them were found to be suitable for transfer. They were transferred to 72 patients in 75 cycles resulting in 23 pregnancies and 20 healthy births. The most common indication was unexplained infertility. The implantation rate was calculated as 23.4% whereas the take-home baby rate was 26.6% per transfer. The highest rate of healthy living births is achieved in patients having low grade maternal mosaic sex chromosomal aneuploidy. All living births achieved by PGS had normal chromosomal structure which we can propose it as an alternative test for couples at risk to select normal embryos to improve the outcomes of assisted reproductive procedures and to avoid the transfer of chromosomally unbalanced and multiple embryos.

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.

TWINS WITH KLEEFSTRA SYNDROME DUE TO CHROMOSOME 9q34.3 MICRODELETION

Kleefstra or 9q subtelomeric deletion syndrome (9qSTDS) is a rare microdeletion syndrome. The most prominent phenotypic features include hypotonia, developmental retardation, as well as typical dysmorphic face. It has been shown that terminal deletions of the chromosome 9q34.3 region, or EHMT1 gene mutations, lead to Kleefstra syndrome. We present 16-month-old twin sisters, one of whom had originally been referred for Down syndrome screening due to hypotonia, growth and development retardation, dysmorphic facial signs, and accompanying congenital heart disease. They were subsequently diagnosed as Kleefstra syndrome based on subtelomeric FISH analysis. In conclusion, Kleefstra syndrome should be considered in the differential diagnosis of Down syndrome because it presents with very similar phenotypic features.

Information-driven structural modelling of protein-protein interactions

Protein-protein docking aims at predicting the three-dimensional structure of a protein complex starting from the free forms of the individual partners. As assessed in the CAPRI community-wide experiment, the most successful docking algorithms combine pure laws of physics with information derived from various experimental or bioinformatics sources. Of these so-called "information-driven" approaches, HADDOCK stands out as one of the most successful representatives. In this chapter, we briefly summarize which experimental information can be used to drive the docking prediction in HADDOCK, and then focus on the docking protocol itself. We discuss and illustrate with a tutorial example a "classical" protein-protein docking prediction, as well as more recent developments for modelling multi-body systems and large conformational changes.

Blind prediction of interfacial water positions in CAPRI

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

A twin sibling with Prader-Willi syndrome caused by uniparental disomy conceived after in vitro fertilization

The use of assisted reproductive technologies (ART) has increased gradually in the treatment of infertility worldwide. On the other hand ART has been found to be associated with an increased risk of congenital malformations including imprinting defects as well. Although a number of imprinting syndromes have been reported to be related with ART, no case with uniparental disomy (UPD) caused Prader-Willi syndrome (PWS) [OMIM ID: 176270] has been reported in the literature. Here we present a dizygotic twin in which one of them was born with maternal UPD15 following ART. The proband was a 2-year-old boy who had feeding difficulties, generalized hypotonia, frontal bossing, broad forehead, small hands and feet. Laboratory investigations revealed minimal dilatation in 3rd and 4th ventricles and corpus callosum hypoplasia in magnetic resonance imaging, supravalvular pulmonary stenosis in echocardiography and pelvicaliectasia in the USG examinations. Methylation and microsatellite markers analyses showed maternal UPD for chromosome 15. Here we report, for the first time UPD caused PWS patient born after ART.

Defining the limits of homology modeling in information-driven protein docking

Information-driven docking is currently one of the most successful approaches to obtain structural models of protein interactions as demonstrated in the latest round of CAPRI. While various experimental and computational techniques can be used to retrieve information about the binding mode, the availability of three-dimensional structures of the interacting partners remains a limiting factor. Fortunately, the wealth of structural information gathered by large-scale initiatives allows for homology-based modeling of a significant fraction of the protein universe. Defining the limits of information-driven docking based on such homology models is therefore highly relevant. Here we show, using previous CAPRI targets, that out of a variety of measures, the global sequence identity between template and target is a simple but reliable predictor of the achievable quality of the docking models. This indicates that a well-defined overall fold is critical for the interaction. Furthermore, the quality of the data at our disposal to characterize the interaction plays a determinant role in the success of the docking. Given reliable interface information we can obtain acceptable predictions even at low global sequence identity. These results, which define the boundaries between trustworthy and unreliable predictions, should guide both experts and nonexperts in defining the limits of what is achievable by docking. This is highly relevant considering that the fraction of the interactome amenable for docking is only bound to grow as the number of experimentally solved structures increases.