DockNet: high-throughput protein-protein interface contact prediction.

MOTIVATION: Over 300 000 protein-protein interaction (PPI) pairs have been identified in the human proteome and targeting these is fast becoming the next frontier in drug design. Predicting PPI sites, however, is a challenging task that traditionally requires computationally expensive and time-consuming docking simulations. A major weakness of modern protein docking algorithms is the inability to account for protein flexibility, which ultimately leads to relatively poor results. RESULTS: Here, we propose DockNet, an efficient Siamese graph-based neural network method which predicts contact residues between two interacting proteins. Unlike other methods that only utilize a protein's surface or treat the protein structure as a rigid body, DockNet incorporates the entire protein structure and places no limits on protein flexibility during an interaction. Predictions are modeled at the residue level, based on a diverse set of input node features including residue type, surface accessibility, residue depth, secondary structure, pharmacophore and torsional angles. DockNet is comparable to current state-of-the-art methods, achieving an area under the curve (AUC) value of up to 0.84 on an independent test set (DB5), can be applied to a variety of different protein structures and can be utilized in situations where accurate unbound protein structures cannot be obtained. AVAILABILITY AND IMPLEMENTATION: DockNet is available at https://github.com/npwilliams09/docknet and an easy-to-use webserver at https://biosig.lab.uq.edu.au/docknet. All other data underlying this article are available in the article and in its online supplementary material. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Functional genomics for curation of variants in telomere biology disorder associated genes: A systematic review.

PURPOSE: Patients with an underlying telomere biology disorder (TBD) have variable clinical presentations, and they can be challenging to diagnose clinically. A genomic diagnosis for patients presenting with TBD is vital for optimal treatment. Unfortunately, many variants identified during diagnostic testing are variants of uncertain significance. This complicates management decisions, delays treatment, and risks nonuptake of potentially curative therapies. Improved application of functional genomic evidence may reduce variants of uncertain significance classifications. METHODS: We systematically searched the literature for published functional assays interrogating TBD gene variants. When possible, established likely benign/benign and likely pathogenic/pathogenic variants were used to estimate the assay sensitivity, specificity, positive predictive value, negative predictive value, and odds of pathogenicity. RESULTS: In total, 3131 articles were screened and 151 met inclusion criteria. Sufficient data to enable a PS3/BS3 recommendation were available for TERT variants only. We recommend that PS3 and BS3 can be applied at a moderate and supportive level, respectively. PS3/BS3 application was limited by a lack of assay standardization and limited inclusion of benign variants. CONCLUSION: Further assay standardization and assessment of benign variants are required for optimal use of the PS3/BS3 criterion for TBD gene variant classification.

A novel small molecule inhibitor of human Drp1.

Mitochondrial dynamin-related protein 1 (Drp1) is a large GTPase regulator of mitochondrial dynamics and is known to play an important role in numerous pathophysiological processes. Despite being the most widely used Drp1 inhibitor, the specificity of Mdivi-1 towards human Drp1 has not been definitively proven and there have been numerous issues reported with its use including off-target effects. In our hands Mdivi-1 showed varying binding affinities toward human Drp1, potentially impacted by compound aggregation. Herein, we sought to identify a novel small molecule inhibitor of Drp1. From an initial virtual screening, we identified DRP1i27 as a compound which directly bound to the human isoform 3 of Drp1 via surface plasmon resonance and microscale thermophoresis. Importantly, DRP1i27 was found to have a dose-dependent increase in the cellular networks of fused mitochondria but had no effect in Drp1 knock-out cells. Further analogues of this compound were identified and screened, though none displayed greater affinity to human Drp1 isoform 3 than DRP1i27. To date, this is the first small molecule inhibitor shown to directly bind to human Drp1.

Structural Analysis of the Complex of Human Transthyretin with 3',5'-Dichlorophenylanthranilic Acid at 1.5 Å Resolution.

Human transthyretin (hTTR) can form amyloid deposits that accumulate in nerves and organs, disrupting cellular function. Molecules such as tafamidis that bind to and stabilize the TTR tetramer can reduce such amyloid formation. Here, we studied the interaction of VCP-6 (2-((3,5-dichlorophenyl)amino)benzoic acid) with hTTR. VCP-6 binds to hTTR with 5 times the affinity of the cognate ligand, thyroxine (T4). The structure of the hTTR:VCP-6 complex was determined by X-ray crystallography at 1.52 Šresolution. VCP-6 binds deeper in the binding channel than T4 with the 3',5'-dichlorophenyl ring binding in the 'forward' mode towards the channel centre. The dichlorophenyl ring lies along the 2-fold axis coincident with the channel centre, while the 2-carboxylatephenylamine ring of VCP-6 is symmetrically displaced from the 2-fold axis, allowing the 2-carboxylate group to form a tight intermolecular hydrogen bond with Nζ of Lys15 and an intramolecular hydrogen bond with the amine of VCP-6, stabilizing its conformation and explaining the greater affinity of VCP-6 compared to T4. This arrangement maintains optimal halogen bonding interactions in the binding sites, via chlorine atoms rather than iodine of the thyroid hormone, thereby explaining why the dichloro substitution pattern is a stronger binder than either the diiodo or dibromo analogues.

Functional interaction between compound heterozygous TERT mutations causes severe telomere biology disorder.

Telomere biology disorders (TBDs) are a spectrum of multisystem inherited disorders characterized by bone marrow failure, resulting from mutations in the genes encoding telomerase or other proteins involved in maintaining telomere length and integrity. Pathogenicity of variants in these genes can be hard to evaluate, because TBD mutations show highly variable penetrance and genetic anticipation related to inheritance of shorter telomeres with each generation. Thus, detailed functional analysis of newly identified variants is often essential. Herein, we describe a patient with compound heterozygous variants in the TERT gene, which encodes the catalytic subunit of telomerase, hTERT. This patient had the extremely severe Hoyeraal-Hreidarsson form of TBD, although his heterozygous parents were clinically unaffected. Molecular dynamic modeling and detailed biochemical analyses demonstrate that one allele (L557P) affects association of hTERT with its cognate RNA component hTR, whereas the other (K1050E) affects the binding of telomerase to its DNA substrate and enzyme processivity. Unexpectedly, the data demonstrate a functional interaction between the proteins encoded by the two alleles, with wild-type hTERT rescuing the effect of K1050E on processivity, whereas L557P hTERT does not. These data contribute to the mechanistic understanding of telomerase, indicating that RNA binding in one hTERT molecule affects the processivity of telomere addition by the other molecule. This work emphasizes the importance of functional characterization of TERT variants to reach a definitive molecular diagnosis for patients with TBD, and, in particular, it illustrates the importance of analyzing the effects of compound heterozygous variants in combination, to reveal interallelic effects.

A structural model of the human plasminogen and Aspergillus fumigatus enolase complex.

The metabolic enzyme, enolase, plays a crucial role in the cytoplasm where it maintains cellular energy production within the process of glycolysis. The main role of enolase in glycolysis is to convert 2-phosphoglycerate to phosphoenolpyruvate; however, enolase can fulfill roles that deviate from this function. In pathogenic bacteria and fungi, enolase is also located on the cell surface where it functions as a virulence factor. Surface-expressed enolase is a receptor for human plasma proteins, including plasminogen, and this interaction facilitates nutrient acquisition and tissue invasion. A novel approach to developing antifungal drugs is to inhibit the formation of this complex. To better understand the structure of enolase and the interactions that may govern complex formation, we have solved the first X-ray crystal structure of enolase from Aspergillus fumigatus (2.0 Å) and have shown that it preferentially adopts a dimeric quaternary structure using native mass spectrometry. Two additional X-ray crystal structures of A. fumigatus enolase bound to the endogenous substrate 2-phosphoglycerate and product phosphoenolpyruvate were determined and kinetic characterization was carried out to better understand the details of its canonical function. From these data, we have produced a model of the A. fumigatus enolase and human plasminogen complex to provide structural insights into the mechanisms of virulence and aid future development of small molecules or peptidomimetics for antifungal drug design.

Analysis of physicochemical properties of protein-protein interaction modulators suggests stronger alignment with the "rule of five".

Despite the important roles played by protein-protein interactions (PPIs) in disease, they have been long considered as 'undruggable'. However, recent advances have suggested that PPIs are druggable but may not follow conventional rules of 'drug ability'. Here we explore which physicochemical parameters are essential for a PPI modulator to be a clinical drug by analysing the physicochemical properties of small-molecule PPI modulators in the market, in clinical trials, and published. Our analysis reveals that those compounds currently on the market have a larger range of values for most of the physicochemical parameters, whereas those in clinical trials fit much more stringently to standard drug-like parameters. This observation was particularly true for molecular weight, clog P and topological polar surface area, where aside from a few outliers, most of the compounds in clinical trials fit within standard drug-like parameters. This implies that the newer PPI modulators are more drug-like than those currently on the market, suggesting that designing new PPI-specific screening libraries should remain within standard drug-like parameters in order to obtain a clinical candidate. Taken together, our analysis has important implications for designing future drug discovery campaigns aimed at targeting PPIs.

Simplified heavy-atom derivatization of protein structures via co-crystallization with the MAD tetragon tetrabromoterephthalic acid.

The phase problem is a persistent bottleneck that impedes the structure-determination pipeline and must be solved to obtain atomic resolution crystal structures of macromolecules. Although molecular replacement has become the predominant method of solving the phase problem, many scenarios still exist in which experimental phasing is needed. Here, a proof-of-concept study is presented that shows the efficacy of using tetrabromoterephthalic acid (B4C) as an experimental phasing compound. Incorporating B4C into the crystal lattice using co-crystallization, the crystal structure of hen egg-white lysozyme was solved using MAD phasing. The strong anomalous signal generated by its four Br atoms coupled with its compatibility with commonly used crystallization reagents render B4C an effective experimental phasing compound that can be used to overcome the phase problem.

Derivatization of Protein Crystals with I3C using Random Microseed Matrix Screening.

Protein structure elucidation using X-ray crystallography requires both high quality diffracting crystals and computational solution of the diffraction phase problem. Novel structures that lack a suitable homology model are often derivatized with heavy atoms to provide experimental phase information. The presented protocol efficiently generates derivatized protein crystals by combining random microseeding matrix screening with derivatization with a heavy atom molecule I3C (5-amino-2,4,6-triiodoisophthalic acid). By incorporating I3C into the crystal lattice, the diffraction phase problem can be efficiently solved using single wavelength anomalous dispersion (SAD) phasing. The equilateral triangle arrangement of iodine atoms in I3C allows for rapid validation of a correct anomalous substructure. This protocol will be useful to structural biologists who solve macromolecular structures using crystallography-based techniques with interest in experimental phasing.

Targeting Unconventional Pathways in Pursuit of Novel Antifungals.

The impact of invasive fungal infections on human health is a serious, but largely overlooked, public health issue. Commonly affecting the immunocompromised community, fungal infections are predominantly caused by species of Candida, Cryptococcus, and Aspergillus. Treatments are reliant on the aggressive use of pre-existing antifungal drug classes that target the fungal cell wall and membrane. Despite their frequent use, these drugs are subject to unfavorable drug-drug interactions, can cause undesirable side-effects and have compromised efficacy due to the emergence of antifungal resistance. Hence, there is a clear need to develop novel classes of antifungal drugs. A promising approach involves exploiting the metabolic needs of fungi by targeted interruption of essential metabolic pathways. This review highlights potential antifungal targets including enolase, a component of the enolase-plasminogen complex, and enzymes from the mannitol biosynthesis and purine nucleotide biosynthesis pathways. There has been increased interest in the enzymes that comprise these particular pathways and further investigation into their merits as antifungal targets and roles in fungal survival and virulence are warranted. Disruption of these vital processes by targeting unconventional pathways with small molecules or antibodies may serve as a promising approach to discovering novel classes of antifungals.

Combining random microseed matrix screening and the magic triangle for the efficient structure solution of a potential lysin from bacteriophage P68.

Two commonly encountered bottlenecks in the structure determination of a protein by X-ray crystallography are screening for conditions that give high-quality crystals and, in the case of novel structures, finding derivatization conditions for experimental phasing. In this study, the phasing molecule 5-amino-2,4,6-triiodoisophthalic acid (I3C) was added to a random microseed matrix screen to generate high-quality crystals derivatized with I3C in a single optimization experiment. I3C, often referred to as the magic triangle, contains an aromatic ring scaffold with three bound I atoms. This approach was applied to efficiently phase the structures of hen egg-white lysozyme and the N-terminal domain of the Orf11 protein from Staphylococcus phage P68 (Orf11 NTD) using SAD phasing. The structure of Orf11 NTD suggests that it may play a role as a virion-associated lysin or endolysin.