Structure-based Profiling Research Articles

Structure-Based Profiling of Potential Phytomolecules with AKT1 a Key Cancer Drug Target

Identifying cancer biomarkers is imperative, as upregulated genes offer a better microenvironment for the tumor; hence, targeted inhibition is preferred. The theme of our study is to predict molecular interactions between cancer biomarker proteins and selected natural compounds. We identified an overexpressed potential molecular target (AKT1) and computationally evaluated its inhibition by four dietary ligands (isoliquiritigenin, shogaol, tehranolide, and theophylline). The three-dimensional structures of protein and phytochemicals were retrieved from the RCSB PDB database (4EKL) and NCBI's PubChem, respectively. Rational structure-based docking studies were performed using AutoDock. Results were analyzed based primarily on the estimated free binding energy (kcal/mol), hydrogen bonds, and inhibition constant, Ki, to identify the most effective anti-cancer phytomolecule. Toxicity and drug-likeliness prediction were performed using OSIRIS and SwissADME. Amongst the four phytocompounds, tehranolide has better potential to suppress the expression of AKT1 and could be used for anti-cancer drug development, as inhibition of AKT1 is directly associated with the inhibition of growth, progression, and metastasis of the tumor. Docking analyses reveal that tehranolide has the most efficiency in inhibiting AKT1 and has the potential to be used for the therapeutic management of cancer. Natural compounds targeting cancer biomarkers offer less rejection, minimal toxicity, and fewer side effects.

BindProfX: Assessing Mutation-Induced Binding Affinity Change by Protein Interface Profiles with Pseudo-Counts

Understanding how gene-level mutations affect the binding affinity of protein–protein interactions is a key issue of protein engineering. Due to the complexity of the problem, using physical force field to predict the mutation-induced binding free-energy change remains challenging. In this work, we present a renewed approach to calculate the impact of gene mutations on the binding affinity through the structure-based profiling of protein–protein interfaces, where the binding free-energy change (ΔΔG) is counted as the logarithm of relative probability of mutant amino acids over wild-type ones in the interface alignment matrix; three pseudo-counts are introduced to alleviate the limit of the current interface library. Compared with a previous profile score that was based on the log-odds likelihood calculation, the correlation between predicted and experimental ΔΔG of single-site mutations is increased in this approach from 0.33 to 0.68. The structure-based profile score is found complementary to the physical potentials, where a linear combination of the profile score with the FoldX potential could increase the ΔΔG correlation from 0.46 to 0.74. It is also shown that the profile score is robust for counting the coupling effect of multiple individual mutations. For the mutations involving more than two mutation sites where the correlation between FoldX and experimental data vanishes, the profile-based calculation retains a strong correlation with the experimental measurements.

ResQ: An Approach to Unified Estimation of B-Factor and Residue-Specific Error in Protein Structure Prediction

Computer-based structure prediction becomes a major tool to provide large-scale structure models for annotating biological function of proteins. Information of residue-level accuracy and thermal mobility (or B-factor), which is critical to decide how biologists utilize the predicted models, is however missed in most structure prediction pipelines. We developed ResQ for unified residue-level model quality and B-factor estimations by combining local structure assembly variations with sequence-based and structure-based profiling. ResQ was tested on 635 non-redundant proteins with structure models generated by I-TASSER, where the average difference between estimated and observed distance errors is 1.4Å for the confidently modeled proteins. ResQ was further tested on structure decoys from CASP9-11 experiments, where the error of local structure quality prediction is consistently lower than or comparable to other state-of-the-art predictors. Finally, ResQ B-factor profile was used to assist molecular replacement, which resulted in successful solutions on several proteins that could not be solved from constant B-factor settings.

Marker-trait association study for sucrose and yield contributing traits in sugarcane (Saccharum spp. hybrid)

Linkage disequilibrium (LD)-based marker-trait association (MTA) was used to identify markers for sucrose and yield contributing traits in a panel of 108 sugarcane genotypes from sub-tropical India. Population structure (Q), kinship (K), and MTA study exploited a set of 989 marker loci generated from 123 genomic- and expressed sequence tag-SSR primers. The mixed linear model (MLM) coupled with a modified algorithm for population structure (Q) analysis was able to control both type I and type II errors and provided a deeper understanding of the genetics, population stratification and its manifestations on LD in the sugarcane genome. Significant associations were identified for four markers with cane diameter, seven markers each with cane length and number of millable canes (NMCs), eleven markers with number of nodes, six with sucrose per cent, and five markers with average cane weight. A total of 15 markers stable for all the 3 years of study explained 57 % trait variation for NMCs, 34 % for cane width, 27 % for cane length, 20 % for sucrose content, and 19 % for number of nodes. The frequent deviation of structure-based profiles from pedigree-based grouping in this complex heterozygous system reinforced the importance of using genotypic data for selection and breeding. The results contribute to a deeper insight of the complex genome and the identified MTAs could be exploited to fine-tune marker-assisted breeding programmes in genetically complex sugarcane crop.

Convergent evolution in structural elements of proteins investigated using cross profile analysis

BackgroundEvolutionary relations of similar segments shared by different protein folds remain controversial, even though many examples of such segments have been found. To date, several methods such as those based on the results of structure comparisons, sequence-based classifications, and sequence-based profile-profile comparisons have been applied to identify such protein segments that possess local similarities in both sequence and structure across protein folds. However, to capture more precise sequence-structure relations, no method reported to date combines structure-based profiles, and sequence-based profiles based on evolutionary information. The former are generally regarded as representing the amino acid preferences at each position of a specific conformation of protein segment. They might reflect the nature of ancient short peptide ancestors, using the results of structural classifications of protein segments.ResultsThis report describes the development and use of "Cross Profile Analysis" to compare sequence-based profiles and structure-based profiles based on amino acid occurrences at each position within a protein segment cluster. Using systematic cross profile analysis, we found structural clusters of 9-residue and 15-residue segments showing remarkably strong correlation with particular sequence profiles. These correlations reflect structural similarities among constituent segments of both sequence-based and structure-based profiles. We also report previously undetectable sequence-structure patterns that transcend protein family and fold boundaries, and present results of the conformational analysis of the deduced peptide of a segment cluster. These results suggest the existence of ancient short-peptide ancestors.ConclusionsCross profile analysis reveals the polyphyletic and convergent evolution of β-hairpin-like structures, which were verified both experimentally and computationally. The results presented here give us new insights into the evolution of short protein segments.

SP5: improving protein fold recognition by using torsion angle profiles and profile-based gap penalty model.

How to recognize the structural fold of a protein is one of the challenges in protein structure prediction. We have developed a series of single (non-consensus) methods (SPARKS, SP2, SP3, SP4) that are based on weighted matching of two to four sequence and structure-based profiles. There is a robust improvement of the accuracy and sensitivity of fold recognition as the number of matching profiles increases. Here, we introduce a new profile-profile comparison term based on real-value dihedral torsion angles. Together with updated real-value solvent accessibility profile and a new variable gap-penalty model based on fractional power of insertion/deletion profiles, the new method (SP5) leads to a robust improvement over previous SP method. There is a 2% absolute increase (5% relative improvement) in alignment accuracy over SP4 based on two independent benchmarks. Moreover, SP5 makes 7% absolute increase (22% relative improvement) in success rate of recognizing correct structural folds, and 32% relative improvement in model accuracy of models within the same fold in Lindahl benchmark. In addition, modeling accuracy of top-1 ranked models is improved by 12% over SP4 for the difficult targets in CASP 7 test set. These results highlight the importance of harnessing predicted structural properties in challenging remote-homolog recognition. The SP5 server is available at http://sparks.informatics.iupui.edu.

Relational Database Driven Two‐Dimensional Chemical Graph Analysis

This manuscript presents a method for pre-computing and storing molecular features or ''scaffolds'' that can be used for rapid clustering of diverse compound sets within the context of a relational database based on hierarchies of scaffold structures. In addition, a method for rapid structure-based profiling of a large compound collection is demonstrated. Pre-organizing compounds by shared structural features in this way facilitates the merger of chemical features and biological data within a relational database.

Evolutionary profiles from the QR factorization of multiple sequence alignments

We present an algorithm to generate complete evolutionary profiles that represent the topology of the molecular phylogenetic tree of the homologous group. The method, based on the multidimensional QR factorization of numerically encoded multiple sequence alignments, removes redundancy from the alignments and orders the protein sequences by increasing linear dependence, resulting in the identification of a minimal basis set of sequences that spans the evolutionary space of the homologous group of proteins. We observe a general trend that these smaller, more evolutionarily balanced profiles have comparable and, in many cases, better performance in database searches than conventional profiles containing hundreds of sequences, constructed in an iterative and computationally intensive procedure. For more diverse families or superfamilies, with sequence identity <30%, structural alignments, based purely on the geometry of the protein structures, provide better alignments than pure sequence-based methods. Merging the structure and sequence information allows the construction of accurate profiles for distantly related groups. These structure-based profiles outperformed other sequence-based methods for finding distant homologs and were used to identify a putative class II cysteinyl-tRNA synthetase (CysRS) in several archaea that eluded previous annotation studies. Phylogenetic analysis showed the putative class II CysRSs to be a monophyletic group and homology modeling revealed a constellation of active site residues similar to that in the known class I CysRS.

On the Role of Structural Information in Remote Homology Detection and Sequence Alignment: New Methods Using Hybrid Sequence Profiles

Structural alignments often reveal relationships between proteins that cannot be detected using sequence alignment alone. However, profile search methods based entirely on structural alignments alone have not been found to be effective in finding remote homologs. Here, we explore the role of structural information in remote homolog detection and sequence alignment. To this end, we develop a series of hybrid multidimensional alignment profiles that combine sequence, secondary and tertiary structure information into hybrid profiles. Sequence-based profiles are profiles whose position-specific scoring matrix is derived from sequence alignment alone; structure-based profiles are those derived from multiple structure alignments. We compare pure sequence-based profiles to pure structure-based profiles, as well as to hybrid profiles that use combined sequence-and-structure-based profiles, where sequence-based profiles are used in loop/motif regions and structural information is used in core structural regions. All of the hybrid methods offer significant improvement over simple profile-to-profile alignment. We demonstrate that both sequence-based and structure-based profiles contribute to remote homology detection and alignment accuracy, and that each contains some unique information. We discuss the implications of these results for further improvements in amino acid sequence and structural analysis.

Secondary structure‐based profiles: Use of structure‐conserving scoring tables in searching protein sequence databases for structural similarities

The profile method, for detecting distantly related proteins by sequence comparison, has been extended to incorporate secondary structure information from known X-ray structures. The sequence of a known structure is aligned to sequences of other members of a given folding class. From the known structure, the secondary structure (alpha-helix, beta-strand or "other") is assigned to each position of the aligned sequences. As in the standard profile method, a position-dependent scoring table, termed a profile, is calculated from the aligned sequences. However, rather than using the standard Dayhoff mutation table in calculating the profile, we use distinct amino acid mutation tables for residues in alpha-helices, beta-strands or other secondary structures to calculate the profile. In addition, we also distinguish between internal and external residues. With this new secondary structure-based profile method, we created a profile for eight-stranded, antiparallel beta barrels of the insecticyanin folding class. It is based on the sequences of retinol-binding protein, insecticyanin and beta-lactoglobulin. Scanning the sequence database with this profile, it was possible to detect the sequence of avidin. The structure of streptavidin is known, and it appears to be distantly related to the antiparallel beta barrels. Also detected is the sequence of complement component C8, which we therefore predict to be a member of this folding class.