Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus.

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho.

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ∼5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 µM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1-KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

First Sprayable Double-Stranded RNA-Based Biopesticide Product Targets Proteasome Subunit Beta Type-5 in Colorado Potato Beetle (Leptinotarsa decemlineata).

Colorado potato beetle (CPB, Leptinotarsa decemlineata) is a major pest of potato and other solanaceous vegetables in the Northern Hemisphere. The insect feeds on leaves and can completely defoliate crops. Because of the repeated use of single insecticide classes without rotating active ingredients, many chemicals are no longer effective in controlling CPB. Ledprona is a sprayable double-stranded RNA biopesticide with a new mode of action that triggers the RNA interference pathway. Laboratory assays with second instar larvae fed Ledprona showed a dose-response where 25×10-6g/L of dsPSMB5 caused 90% mortality after 6days of initial exposure. We also showed that exposure to Ledprona for 6h caused larval mortality and decreased target messenger RNA (mRNA) expression. Decrease in PSMB5 protein levels was observed after 48h of larval exposure to Ledprona. Both PSMB5 mRNA and protein levels did not recover over time. Ledprona efficacy was demonstrated in a whole plant greenhouse trial and performed similarly to spinosad. Ledprona, currently pending registration at EPA, represents a new biopesticide class integrated pest management and insecticide resistance management programs directed against CPB.

Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes.

Protein-protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. In this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein-protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.

Computational Methods for Predicting Functions at the mRNA Isoform Level.

Multiple mRNA isoforms of the same gene are produced via alternative splicing, a biological mechanism that regulates protein diversity while maintaining genome size. Alternatively spliced mRNA isoforms of the same gene may sometimes have very similar sequence, but they can have significantly diverse effects on cellular function and regulation. The products of alternative splicing have important and diverse functional roles, such as response to environmental stress, regulation of gene expression, human heritable, and plant diseases. The mRNA isoforms of the same gene can have dramatically different functions. Despite the functional importance of mRNA isoforms, very little has been done to annotate their functions. The recent years have however seen the development of several computational methods aimed at predicting mRNA isoform level biological functions. These methods use a wide array of proteo-genomic data to develop machine learning-based mRNA isoform function prediction tools. In this review, we discuss the computational methods developed for predicting the biological function at the individual mRNA isoform level.

Combinatorial biosynthesis and the basis for substrate promiscuity in class I diterpene synthases.

Terpene synthases are capable of mediating complex reactions, but fundamentally simply catalyze lysis of allylic diphosphate esters with subsequent deprotonation. Even with the initially generated tertiary carbocation this offers a variety of product outcomes, and deprotonation further can be preceded by the addition of water. This is particularly evident with labdane-related diterpenes (LRDs) where such lysis follows bicyclization catalyzed by class II diterpene cyclases (DTCs) that generates preceding structural variation. Previous investigation revealed that two diterpene synthases (DTSs), one bacterial and the other plant-derived, exhibit extreme substrate promiscuity, but yet still typically produce exo-ene or tertiary alcohol LRD derivatives, respectively (i.e., demonstrating high catalytic specificity), enabling rational combinatorial biosynthesis. Here two DTSs that produce either cis or trans endo-ene LRD derivatives, also plant and bacterial (respectively), were examined for their potential analogous utility. Only the bacterial trans-endo-ene forming DTS was found to exhibit significant substrate promiscuity (with moderate catalytic specificity). This further led to investigation of the basis for substrate promiscuity, which was found to be more closely correlated with phylogenetic origin than reaction complexity. Specifically, bacterial DTSs exhibited significantly more substrate promiscuity than those from plants, presumably reflecting their distinct evolutionary context. In particular, plants typically have heavily elaborated LRD metabolism, in contrast to the rarity of such natural products in bacteria, and the lack of potential substrates presumably alleviates selective pressure against such promiscuity. Regardless of such speculation, this work provides novel biosynthetic access to almost 19 LRDs, demonstrating the power of the combinatorial approach taken here.

Coupling dynamics and evolutionary information with structure to identify protein regulatory and functional binding sites.

Binding sites in proteins can be either specifically functional binding sites (active sites) that bind specific substrates with high affinity or regulatory binding sites (allosteric sites), that modulate the activity of functional binding sites through effector molecules. Owing to their significance in determining protein function, the identification of protein functional and regulatory binding sites is widely acknowledged as an important biological problem. In this work, we present a novel binding site prediction method, Active and Regulatory site Prediction (AR-Pred), which supplements protein geometry, evolutionary, and physicochemical features with information about protein dynamics to predict putative active and allosteric site residues. As the intrinsic dynamics of globular proteins plays an essential role in controlling binding events, we find it to be an important feature for the identification of protein binding sites. We train and validate our predictive models on multiple balanced training and validation sets with random forest machine learning and obtain an ensemble of discrete models for each prediction type. Our models for active site prediction yield a median area under the curve (AUC) of 91% and Matthews correlation coefficient (MCC) of 0.68, whereas the less well-defined allosteric sites are predicted at a lower level with a median AUC of 80% and MCC of 0.48. When tested on an independent set of proteins, our models for active site prediction show comparable performance to two existing methods and gains compared to two others, while the allosteric site models show gains when tested against three existing prediction methods. AR-Pred is available as a free downloadable package at https://github.com/sambitmishra0628/AR-PRED_source.

Comparisons of Protein Dynamics from Experimental Structure Ensembles, Molecular Dynamics Ensembles, and Coarse-Grained Elastic Network Models.

Predicting protein motions is important for bridging the gap between protein structure and function. With growing numbers of structures of the same or closely related proteins becoming available, it is now possible to understand more about the intrinsic dynamics of a protein with principal component analysis (PCA) of the motions apparent within ensembles of experimental structures. In this paper, we compare the motions extracted from experimental ensembles of 50 different proteins with the modes of motion predicted by several types of coarse-grained elastic network models (ENMs) which additionally take into account more details of either the protein geometry or the amino acid specificity. We further compare the structural variations in the experimental ensembles with the motions sampled in molecular dynamics (MD) simulations for a smaller subset of 17 proteins with available trajectories. We find that the correlations between the motions extracted from MD trajectories and experimental structure ensembles are slightly different than those for the ENMs, possibly reflecting potential sampling biases. We find that there are small gains in the predictive power of the ENMs in reproducing motions present in either experimental or MD ensembles by accounting for the protein geometry rather than the amino acid specificity of the interactions.