Aortic disease and cardiomyopathy in patients with a novel DNMT3A gene variant causing Tatton-Brown-Rahman syndrome.

Tatton-Brown-Rahman syndrome (TBRS) is a rare congenital genetic disorder caused by autosomal dominant pathogenic variants in the DNA methyltransferase DNMT3A gene. Typical TBRS clinical features are overgrowth, intellectual disability, and minor facial anomalies. However, since the syndrome was first described in 2014, a widening spectrum of abnormalities is being described. Cardiovascular abnormalities are less commonly reported but can be a major complication of the syndrome. This article describes a family of three individuals diagnosed with TBRS in adulthood and highlights the variable expression of cardiovascular features. A 34-year-old proband presented with progressive aortic dilatation, mitral valve (MV) regurgitation, left ventricular (LV) dilatation, and ventricular arrhythmias. The affected family members (mother and brother) were diagnosed with MV regurgitation, LV dilatation, and arrhythmias. Exome sequencing and computational protein analysis suggested that the novel familial DNMT3A mutation Ser775Tyr is located in the methyltransferase domain, however, distant from the active site or DNA-binding loops. Nevertheless, this bulky substitution may have a significant effect on DNMT3A protein structure, dynamics, and function. Analysis of peripheral blood cfDNA and transcriptome showed shortened mononucleosome fragments and altered gene expression in a number of genes related to cardiovascular health and of yet undescribed function, including several lncRNAs. This highlights the importance of epigenetic regulation by DNMT3A on cardiovascular system development and function. From the clinical perspective, we suggest that new patients diagnosed with congenital DNMT3A variants and TBRS require close examination and follow-up for aortic dilatation and valvular disease because these conditions can progress rapidly. Moreover, personalized treatments, based on the specific DNMT3A variants and the different pathways of their function loss, can be envisioned in the future.

PPI3D: a web server for searching, analyzing and modeling protein-protein, protein-peptide and protein-nucleic acid interactions.

Structure-resolved protein interactions with other proteins, peptides and nucleic acids are key for understanding molecular mechanisms. The PPI3D web server enables researchers to query preprocessed and clustered structural data, analyze the results and make homology-based inferences for protein interactions. PPI3D offers three interaction exploration modes: (i) all interactions for proteins homologous to the query, (ii) interactions between two proteins or their homologs and (iii) interactions within a specific PDB entry. The server allows interactive analysis of the identified interactions in both summarized and detailed manner. This includes protein annotations, structures, the interface residues and the corresponding contact surface areas. In addition, users can make inferences about residues at the interaction interface for the query protein(s) from the sequence alignments and homology models. The weekly updated PPI3D database includes all the interaction interfaces and binding sites from PDB, clustered based on both protein sequence and structural similarity, yielding non-redundant datasets without loss of alternative interaction modes. Consequently, the PPI3D users avoid being flooded with redundant information, a typical situation for intensely studied proteins. Furthermore, PPI3D provides a possibility to download user-defined sets of interaction interfaces and analyze them locally. The PPI3D web server is available at https://bioinformatics.lt/ppi3d.

Prediction of protein assemblies by structure sampling followed by interface-focused scoring.

Proteins often function as part of permanent or transient multimeric complexes, and understanding function of these assemblies requires knowledge of their three-dimensional structures. While the ability of AlphaFold to predict structures of individual proteins with unprecedented accuracy has revolutionized structural biology, modeling structures of protein assemblies remains challenging. To address this challenge, we developed a protocol for predicting structures of protein complexes involving model sampling followed by scoring focused on the subunit-subunit interaction interface. In this protocol, we diversified AlphaFold models by varying construction and pairing of multiple sequence alignments as well as increasing the number of recycles. In cases when AlphaFold failed to assemble a full protein complex or produced unreliable results, additional diverse models were constructed by docking of monomers or subcomplexes. All the models were then scored using a newly developed method, VoroIF-jury, which relies only on structural information. Notably, VoroIF-jury is independent of AlphaFold self-assessment scores and therefore can be used to rank models originating from different structure prediction methods. We tested our protocol in CASP15 and obtained top results, significantly outperforming the standard AlphaFold-Multimer pipeline. Analysis of our results showed that the accuracy of our assembly models was capped mainly by structure sampling rather than model scoring. This observation suggests that better sampling, especially for the antibody-antigen complexes, may lead to further improvement. Our protocol is expected to be useful for modeling and/or scoring protein assemblies.

The COMER web server for protein analysis by homology.

SUMMARY: Sequence homology is a basic concept in protein evolution, structure and function studies. However, there are not many different tools and services for homology searches being sensitive, accurate and fast at the same time. We present a new web server for protein analysis based on COMER2, a sequence alignment and homology search method that exhibits these characteristics. COMER2 has been upgraded since its last publication to improve its alignment quality and ease of use. We demonstrate how the user can benefit from using it by providing examples of extensive annotation of proteins of unknown function. Among the distinctive features of the web server is the user's ability to submit multiple queries with one click of a button. This and other features allow for transparently running homology searches-in a command-line, programmatic or graphical environment-across multiple databases with multiple queries. They also promote extensive simultaneous protein analysis at the sequence, structure and function levels. AVAILABILITY AND IMPLEMENTATION: The COMER web server is available at https://bioinformatics.lt/comer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PIGN-Related Disease in Two Lithuanian Families: A Report of Two Novel Pathogenic Variants, Molecular and Clinical Characterisation.

Background and Objectives: Pathogenic variants of PIGN are a known cause of multiple congenital anomalies-hypotonia-seizures syndrome 1 (MCAHS1). Many affected individuals have clinical features overlapping with Fryns syndrome and are mainly characterised by developmental delay, congenital anomalies, hypotonia, seizures, and specific minor facial anomalies. This study investigates the clinical and molecular data of three individuals from two unrelated families, the clinical features of which were consistent with a diagnosis of MCAHS1. Materials and Methods: Next-generation sequencing (NGS) technology was used to identify the changes in the DNA sequence. Sanger sequencing of gDNA of probands and their parents was used for validation and segregation analysis. Bioinformatics tools were used to investigate the consequences of pathogenic or likely pathogenic PIGN variants at the protein sequence and structure level. Results: The analysis of NGS data and segregation analysis revealed a compound heterozygous NM_176787.5:c.[1942G>T];[1247_1251del] PIGN genotype in family 1 and NG_033144.1(NM_176787.5):c.[932T>G];[1674+1G>C] PIGN genotype in family 2. In silico, c.1942G>T (p.(Glu648Ter)), c.1247_1251del (p.(Glu416GlyfsTer22)), and c.1674+1G>C (p.(Glu525AspfsTer68)) variants are predicted to result in a premature termination codon that leads to truncated and functionally disrupted protein causing the phenotype of MCAHS1 in the affected individuals. Conclusions: PIGN-related disease represents a wide spectrum of phenotypic features, making clinical diagnosis inaccurate and complicated. The genetic testing of every individual with this phenotype provides new insights into the origin and development of the disease.

Donor Splice Site Variant in SLC9A6 Causes Christianson Syndrome in a Lithuanian Family: A Case Report.

Background and Objectives: The pathogenic variants of SLC9A6 are a known cause of a rare, X-linked neurological disorder called Christianson syndrome (CS). The main characteristics of CS are developmental delay, intellectual disability, and neurological findings. This study investigated the genetic basis and explored the molecular changes that led to CS in two male siblings presenting with intellectual disability, epilepsy, behavioural problems, gastrointestinal dysfunction, poor height, and weight gain. Materials and Methods: Next-generation sequencing of a tetrad was applied to identify the DNA changes and Sanger sequencing of proband's cDNA was used to evaluate the impact of a splice site variant on mRNA structure. Bioinformatical tools were used to investigate SLC9A6 protein structure changes. Results: Sequencing and bioinformatical analysis revealed a novel donor splice site variant (NC_000023.11(NM_001042537.1):c.899 + 1G > A) that leads to a frameshift and a premature stop codon. Protein structure modelling showed that the truncated protein is unlikely to form any functionally relevant SLC9A6 dimers. Conclusions: Molecular and bioinformatical analysis revealed the impact of a novel donor splice site variant in the SLC9A6 gene that leads to truncated and functionally disrupted protein causing the phenotype of CS in the affected individuals.

Modeling of protein complexes in CASP14 with emphasis on the interaction interface prediction.

The goal of CASP experiments is to monitor the progress in the protein structure prediction field. During the 14th CASP edition we aimed to test our capabilities of predicting structures of protein complexes. Our protocol for modeling protein assemblies included both template-based modeling and free docking. Structural templates were identified using sensitive sequence-based searches. If sequence-based searches failed, we performed structure-based template searches using selected CASP server models. In the absence of reliable templates we applied free docking starting from monomers generated by CASP servers. We evaluated and ranked models of protein complexes using an improved version of our protein structure quality assessment method, VoroMQA, taking into account both interaction interface and global structure scores. If reliable templates could be identified, generally accurate models of protein assemblies were generated with the exception of an antibody-antigen interaction. The success of free docking mainly depended on the accuracy of initial subunit models and on the scoring of docking solutions. To put our overall results in perspective, we analyzed our performance in the context of other CASP groups. Although the subunits in our assembly models often were not of the top quality, these models had, overall, the best-predicted intersubunit interfaces according to several accuracy measures. We attribute our relative success primarily to the emphasis on the interaction interface when modeling and scoring.

Mapping of Recognition Sites of Monoclonal Antibodies Responsible for the Inhibition of Pneumolysin Functional Activity.

The pathogenicity of many bacteria, including Streptococcus pneumoniae, depends on pore-forming toxins (PFTs) that cause host cell lysis by forming large pores in cholesterol-containing cell membranes. Therefore, PFTs-neutralising antibodies may provide useful tools for reducing S. pneumoniae pathogenic effects. This study aimed at the development and characterisation of monoclonal antibodies (MAbs) with neutralising activity to S. pneumoniae PFT pneumolysin (PLY). Five out of 10 produced MAbs were able to neutralise the cytolytic activity of PLY on a lung epithelial cell line. Epitope mapping with a series of recombinant overlapping PLY fragments revealed that neutralising MAbs are directed against PLY loops L1 and L3 within domain 4. The epitopes of MAbs 3A9, 6E5 and 12F11 located at L1 loop (aa 454-471) were crucial for PLY binding to the immobilised cholesterol. In contrast, the MAb 12D10 recognising L3 (aa 403-423) and the MAb 3F3 against the conformational epitope did not interfere with PLY-cholesterol interaction. Due to conformation-dependent binding, the approach to use overlapping peptides for fine epitope mapping of the neutralising MAbs was unsuccessful. Therefore, the epitopes recognised by the MAbs were analysed using computational methods. This study provides new data on PLY sites involved in functional activity.

Template-Based Modeling of Protein Complexes Using the PPI3D Web Server.

There is a large gap between the numbers of known protein-protein interactions and the corresponding experimentally solved structures of protein complexes. Fortunately, this gap can be in part bridged by computational structure modeling methods. Currently, template-based modeling is the most accurate means to predict both individual protein structures and protein complexes. One of the major issues in template-based modeling is to identify homologous structures that could be utilized as templates. To simplify this task, we have developed the PPI3D web server. The server is not only able to search for homologous protein complexes, but also provides means to analyze identified interactions and to model protein complexes. In recent CASP and CAPRI experiments, PPI3D proved to be a useful tool for homology modeling of multimeric proteins. In this chapter, we provide a brief description of the PPI3D web server capabilities and how to use the server for modeling of protein complexes.

Template-based modeling of diverse protein interactions in CAPRI rounds 38-45.

Structures of proteins complexed with other proteins, peptides, or ligands are essential for investigation of molecular mechanisms. However, the experimental structures of protein complexes of interest are often not available. Therefore, computational methods are widely used to predict these structures, and, of those methods, template-based modeling is the most successful. In the rounds 38-45 of the Critical Assessment of PRediction of Interactions (CAPRI), we applied template-based modeling for 9 of 11 protein-protein and protein-peptide interaction targets, resulting in medium and high-quality models for six targets. For the protein-oligosaccharide docking targets, we used constraints derived from template structures, and generated models of at least acceptable quality for most of the targets. Apparently, high flexibility of oligosaccharide molecules was the main cause preventing us from obtaining models of higher quality. We also participated in the CAPRI scoring challenge, the goal of which was to identify the highest quality models from a large pool of decoys. In this experiment, we tested VoroMQA, a scoring method based on interatomic contact areas. The results showed VoroMQA to be quite effective in scoring strongly binding and obligatory protein complexes, but less successful in the case of transient interactions. We extensively used manual intervention in both CAPRI modeling and scoring experiments. This oftentimes allowed us to select the correct templates from available alternatives and to limit the search space during the model scoring.

Structural modeling of protein complexes: Current capabilities and challenges.

Proteins frequently interact with each other, and the knowledge of structures of the corresponding protein complexes is necessary to understand how they function. Computational methods are increasingly used to provide structural models of protein complexes. Not surprisingly, community-wide Critical Assessment of protein Structure Prediction (CASP) experiments have recently started monitoring the progress in this research area. We participated in CASP13 with the aim to evaluate our current capabilities in modeling of protein complexes and to gain a better understanding of factors that exert the largest impact on these capabilities. To model protein complexes in CASP13, we applied template-based modeling, free docking and hybrid techniques that enabled us to generate models of the topmost quality for 27 of 42 multimers. If templates for protein complexes could be identified, we modeled the structures with reasonable accuracy by straightforward homology modeling. If only partial templates were available, it was nevertheless possible to predict the interaction interfaces correctly or to generate acceptable models for protein complexes by combining template-based modeling with docking. If no templates were available, we used rigid-body docking with limited success. However, in some free docking models, despite the incorrect subunit orientation and missed interface contacts, the approximate location of protein binding sites was identified correctly. Apparently, our overall performance in docking was limited by the quality of monomer models and by the imperfection of scoring methods. The impact of human intervention on our results in modeling of protein complexes was significant indicating the need for improvements of automatic methods.

SKEMPI 2.0: an updated benchmark of changes in protein-protein binding energy, kinetics and thermodynamics upon mutation.

Motivation: Understanding the relationship between the sequence, structure, binding energy, binding kinetics and binding thermodynamics of protein-protein interactions is crucial to understanding cellular signaling, the assembly and regulation of molecular complexes, the mechanisms through which mutations lead to disease, and protein engineering. Results: We present SKEMPI 2.0, a major update to our database of binding free energy changes upon mutation for structurally resolved protein-protein interactions. This version now contains manually curated binding data for 7085 mutations, an increase of 133%, including changes in kinetics for 1844 mutations, enthalpy and entropy changes for 443 mutations, and 440 mutations, which abolish detectable binding. Availability and implementation: The database is available as supplementary data and at https://life.bsc.es/pid/skempi2/. Supplementary information: Supplementary data are available at Bioinformatics online.

Exogenous interleukin-1α signaling negatively impacts acquired chemoresistance and alters cell adhesion molecule expression pattern in colorectal carcinoma cells HCT116.

Proinflammatory cytokine and chemokine signaling from the tumor microenvironment is thought to be crucial for developing and sustaining colorectal cancer by regulating a multitude of pathways associated with a variety of cellular mechanisms. Among these pathways there is acquired chemoresistance, which is usually a major obstacle in the way towards successful chemotherapeutic treatment of advanced colorectal cancer cases. Despite of an emerging body of data published on the role of cytokine signaling network in cancer, little is known about the effects of the upstream cytokine interleukin-1α (IL-1α) signaling to the cancer cells. In this study we have shown that the increase in exogenous IL-1α signaling increases chemosensitivity of both chemosensitive and chemoresistant colorectal cancer cell lines, treated with a widely used cytotoxic antimetabolite 5-fluorouracil (5-FU). This was a result of increased cell death but not of the changes in 5-FU-induced cell cycle arrest. Noticeably, combined exogenous IL-1α and 5-FU treatment had significant effects on the expression of cell adhesion molecules, suggesting a decrease in adhesion-dependent chemoresistance and, on the other hand, an increase in metastatic potential of the cells. These results lead to a conclusion that modulation of IL-1 receptor activity could have applications as a part of combination therapy for advanced and highly metastatic colorectal cancers.

Modeling of protein complexes in CAPRI Round 37 using template-based approach combined with model selection.

We participated in Round 37 of the Critical Assessment of PRediction of Interactions (CAPRI), held jointly with the 12th edition of the Critical Assessment of protein Structure Prediction (CASP12), having two major objectives. First, we intended to test the utility of our PPI3D web server in finding and selecting templates for comparative modeling of structures of protein complexes. Our second aim was to evaluate the ability of our model accuracy estimation method VoroMQA to score and rank structural models for protein-protein interactions. Using sequence search in PPI3D and HHpred servers we identified multimeric templates for 7 of 11 CAPRI targets, and models of at least acceptable quality were constructed for 6 of them. The clustering and visual analysis features implemented in the PPI3D software were instrumental in detecting alternative protein-protein interaction interfaces among the identified templates. When a single binding mode was observed for homologous proteins, the structural modeling of the protein complex was fairly straightforward, whereas choosing the correct interaction template from several alternatives turned out to be a difficult task requiring manual intervention. The combination of full structure and interaction interface VoroMQA scores effectively ranked structural models of protein complexes and selected models of better quality from the CAPRI Scoring sets. The overall results show possible uses of PPI3D and VoroMQA in structural modeling of protein-protein interactions and suggest ways for further improvements of both methods.

The PPI3D web server for searching, analyzing and modeling protein-protein interactions in the context of 3D structures.

Summary: The PPI3D web server is focused on searching and analyzing the structural data on protein-protein interactions. Reducing the data redundancy by clustering and analyzing the properties of interaction interfaces using Voronoi tessellation makes this software a highly effective tool for addressing different questions related to protein interactions. Availability and Implementation: The server is freely accessible at http://bioinformatics.lt/software/ppi3d/ . Contact: ceslovas.venclovas@bti.vu.lt. Supplementary information: Supplementary data are available at Bioinformatics online.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data.

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

QSAR analysis of blood-brain distribution: the influence of plasma and brain tissue binding.

The extent of brain delivery expressed as steady-state brain/blood distribution ratio (log BB) is the most frequently used parameter for characterizing central nervous system exposure of drugs and drug candidates. The aim of the current study was to propose a physicochemical QSAR model for log BB prediction. Model development involved the following steps: (i) A data set consisting of 470 experimental log BB values determined in rodents was compiled and verified to ensure that selected data represented drug disposition governed by passive diffusion across blood-brain barrier. (ii) Available log BB values were corrected for unbound fraction in plasma to separate the influence of drug binding to brain and plasma constituents. (iii) The resulting ratios of total brain to unbound plasma concentrations reflecting brain tissue binding were described by a nonlinear ionization-specific model in terms of octanol/water log P and pK(a). The results of internal and external validation demonstrated good predictive power of the obtained model as both log BB and brain tissue binding strength were predicted with residual mean square error of 0.4 log units. The statistical parameters were similar among training and validation sets, indicating that the model is not likely to be overfitted.

Trainable structure-activity relationship model for virtual screening of CYP3A4 inhibition.

A new structure-activity relationship model predicting the probability for a compound to inhibit human cytochrome P450 3A4 has been developed using data for >800 compounds from various literature sources and tested on PubChem screening data. Novel GALAS (Global, Adjusted Locally According to Similarity) modeling methodology has been used, which is a combination of baseline global QSAR model and local similarity based corrections. GALAS modeling method allows forecasting the reliability of prediction thus defining the model applicability domain. For compounds within this domain the statistical results of the final model approach the data consistency between experimental data from literature and PubChem datasets with the overall accuracy of 89%. However, the original model is applicable only for less than a half of PubChem database. Since the similarity correction procedure of GALAS modeling method allows straightforward model training, the possibility to expand the applicability domain has been investigated. Experimental data from PubChem dataset served as an example of in-house high-throughput screening data. The model successfully adapted itself to both data classified using the same and different IC50 threshold compared with the training set. In addition, adjustment of the CYP3A4 inhibition model to compounds with a novel chemical scaffold has been demonstrated. The reported GALAS model is proposed as a useful tool for virtual screening of compounds for possible drug-drug interactions even prior to the actual synthesis.

Probabilistic prediction of the human CYP3A4 and CYP2D6 metabolism sites.

This article briefly introduces the results of in silico prediction of the most probable metabolism sites for the human cytochrome P450 3A4 and 2D6 isoforms. Ligand-based QSAR models have been developed using a novel GALAS modeling approach, and provide probabilities of being a target of CYP3A4 or CYP2D6 for any atom in a molecule. The GALAS-model development methodology allows evaluation of the reliability of predictions in the form of estimated prediction Reliability Indices (RIs). For all the models considered in this study, the number of misclassifications and inconclusive results was reduced significantly when only predictions of high quality (RI>0.5) were taken into account, demonstrating that RI reflects accuracy of prediction. The applicability domain of regioselectivity models is shown to be easily expandable to cover compound classes of interest to the user. The results obtained so far show promising perspectives for the utilization of the GALAS modeling in the analysis of regioselectivity for other important biotransformation enzymes--a work currently in progress.