Integrating Bonded and Nonbonded Potentials in the Knowledge-Based Scoring Function for Protein Structure Prediction

Illuminating the "Twilight Zone": Advances in Difficult Protein Modeling

Homology modeling was long considered a method of choice in tertiary protein structure prediction. However, it used to provide models of acceptable quality only when templates with appreciable sequence identity with a target could be found. The threshold value was long assumed to be around 20-30%. Below this level, obtained sequence identity was getting dangerously close to values that can be obtained by chance, after aligning any random, unrelated sequences. In these cases, other approaches, including ab initio folding simulations or fragment assembly, were usually employed. The most recent editions of the CASP and CAMEO community-wide modeling methods assessment have brought some surprising outcomes, proving that much more clues can be inferred from protein sequence analyses than previously thought. In this chapter, we focus on recent advances in the field of difficult protein modeling, pushing the threshold deep into the "twilight zone", with particular attention devoted to improvements in applications of machine learning and model evaluation.

Structure Prediction and Potential Inhibitors Docking of Enterovirus 2C Proteins

Human enterovirus infections are mostly asymptomatic and occasionally could be severe and life-threatening. The conserved non-structural 2C from enteroviruses protein is a promising target in antiviral therapies against human enteroviruses. Understanding of 2C-drug interactions is crucial for developing the potential antiviral agents. While functions of enterovirus 2C proteins have been widely studied, three-dimensional structure information of 2C is limited. In this study, the structures of 2C proteins from 20 enteroviruses were simulated and reconstructed using I-TASSER programs. Subsequent docking studies of the known 22 antiviral inhibitors for 2C proteins were performed to uncover the inhibitor-binding characteristics of 2C. Among the potential inhibitors, the compound hydantoin exhibited the highest broad-spectrum antiviral activities with binding to 2C protein. The anti-enteroviral activity of GuaHCL, compound 19b, R523062, compound 12a, compound 12b, quinoline analogs 12a, compound 19d, N⁶-benzyladenosine, dibucaine derivatives 6i, TBZE-029, fluoxetine analogs 2b, dibucaine, 2-(α-hydroxybenzyl)-benzimidazole (HBB), metrifudil, pirlindole, MRL-1237, quinoline analogs 10a, zuclopenthixol, fluoxetine, fluoxetine HCl, and quinoline analogs 12c showed a trend of gradual decrease. In addition, the free energy with 22 compounds binding to EV 2C ranged from −0.35 to −88.18 kcal/mol. Our in silico studies will provide important information for the development of pan-enterovirus antiviral agents based on 2C.

Statistical potentials from the Gaussian scaling behaviour of chain fragments buried within protein globules

Knowledge-based approaches use the statistics collected from protein data-bank structures to estimate effective interaction potentials between amino acid pairs. Empirical relations are typically employed that are based on the crucial choice of a reference state associated to the null interaction case. Despite their significant effectiveness, the physical interpretation of knowledge-based potentials has been repeatedly questioned, with no consensus on the choice of the reference state. Here we use the fact that the Flory theorem, originally derived for chains in a dense polymer melt, holds also for chain fragments within the core of globular proteins, if the average over buried fragments collected from different non-redundant native structures is considered. After verifying that the ensuing Gaussian statistics, a hallmark of effectively non-interacting polymer chains, holds for a wide range of fragment lengths, although with significant deviations at short spatial scales, we use it to define a ‘bona fide’ reference state. Notably, despite the latter does depend on fragment length, deviations from it do not. This allows to estimate an effective interaction potential which is not biased by the presence of correlations due to the connectivity of the protein chain. We show how different sequence-independent effective statistical potentials can be derived using this approach by coarse-graining the protein representation at varying levels. The possibility of defining sequence-dependent potentials is explored.

Computational chemistry methods to investigate the effects caused by DNA variants linked with disease

Computational chemistry offers variety of tools to study properties of biological macromolecules. These tools vary in terms of levels of details from quantum mechanical treatment to numerous macroscopic approaches. Here, we provide a review of computational chemistry algorithms and tools for modeling the effects of genetic variations and their association with diseases. Particular emphasis is given on modeling the effects of missense mutations on stability, conformational dynamics, binding, hydrogen bond network, salt bridges, and pH-dependent properties of the corresponding macromolecules. It is outlined that the disease may be caused by alteration of one or several of above-mentioned biophysical characteristics, and a successful prediction of pathogenicity requires detailed analysis of how the alterations affect the function of involved macromolecules. The review provides a short list of most commonly used algorithms to predict the molecular effects of mutations as well.