Improving the performance of the PLB index for ligand-binding site prediction using dihedral angles and the solvent-accessible surface area

Unraveling Hydration Shell Dynamics and Viscosity Effects Around Cyanamide Probes via 2D IR Spectroscopy

Hydration dynamics and solvent viscosity play critical roles in the structure and function of biomolecules. An overwhelming body of evidence suggests that protein and membrane fluctuations are closely linked to solvent fluctuations. While extensive research exists on the use of vibrational probes to detect local interactions and solvent dynamics, fewer studies have explored how the behavior of these reporters changes in response to bulk viscosity. To address this gap, two-dimensional infrared spectroscopy (2D IR) was employed in this study to investigate the ultrafast hydration dynamics around a cyanamide (NCN) probe attached to a nucleoside, deoxycytidine, in aqueous solutions with varying glycerol content. The use of a small vibrational probe on a targeted nucleic acid offers the potential to capture more localized hydration dynamics than alternative methods. The time scales for the frequency correlation decays were found to increase linearly with bulk viscosity, ranging from 0.9 to 11.4 ps over viscosities of 0.96-49.1 cP. Additionally, molecular dynamics (MD) simulations were performed to model the local hydration dynamics around the NCN probe. Interestingly, increasing the glycerol content did not significantly alter the hydration of the deoxycytidine. The MD simulations further suggested that the NCN probe's frequency fluctuations were primarily influenced by the dynamics of water in the second solvation shell. Cage correlation functions, which measure the movement of water molecules in and out of the second solvation shell, exhibited decays that closely matched those of the frequency-fluctuation correlation function (FFCF). These findings offer new insights into hydration dynamics and the impact of viscosity on biological systems.

A unique approach for protein secondary structure comparison under TOPS representation

To unravel the intricate connection between protein function and protein structure, it is imperative to comprehensively evaluate protein secondary structure similarity from various perspectives. While numerous techniques have been suggested for comparing protein secondary structure elements (SSE), there continues to be a substantial need for finding alternative ways of comparing the same. In this paper, Topology of Protein Structure (TOPS) representations of protein secondary structures are considered to offer a new alignment-free method for evaluating similarities/dissimilarities of protein secondary structures. Initially, a two-dimensional numerical representation of the SSE is created, associating each point with a mass reflecting its frequency of occurrence. Then the means of coordinate values are determined by averaging weighted sums, and these mean values are subsequently used to calculate moments-of-inertia. Next, a four-component descriptor is generated out of the eigenvalues of the matrix and the mean values of the represented coordinates. Thereafter, Manhattan distance measure is used to obtain the distance matrix. This is finally applied to obtain the phylogenetic trees under the use of NJ method. SSE considered in the proposed method comprises 36-elements from the Chew-Kedem database giving five different taxa: globin, alpha-beta, tim-barrel, beta, and alpha. Phylogenetic trees were created for these SSE through the application of various methods: Clustal-Omega, LZ-Complexity, SED, TOPS + and TOC, to facilitate comparative analysis. Phylogenetic tree of the proposed method outperformed results of the previous methods when applied to the same SSE. Therefore, the method effectively constructs phylogenetic tree for analyzing protein secondary structure comparison.Communicated by Ramaswamy H. Sarma.

Design of rigid protein-protein interaction inhibitors enables targeting of undruggable Mcl-1

The structure-based design of small-molecule inhibitors targeting protein-protein interactions (PPIs) remains a huge challenge as the drug must bind typically wide and shallow protein sites. A PPI target of high interest for hematological cancer therapy is myeloid cell leukemia 1 (Mcl-1), a prosurvival guardian protein from the Bcl-2 family. Despite being previously considered undruggable, seven small-molecule Mcl-1 inhibitors have recently entered clinical trials. Here, we report the crystal structure of the clinical-stage inhibitor AMG-176 bound to Mcl-1 and analyze its interaction along with clinical inhibitors AZD5991 and S64315. Our X-ray data reveal high plasticity of Mcl-1 and a remarkable ligand-induced pocket deepening. Nuclear Magnetic Resonance (NMR)-based free ligand conformer analysis demonstrates that such unprecedented induced fit is uniquely achieved by designing highly rigid inhibitors, preorganized in their bioactive conformation. By elucidating key chemistry design principles, this work provides a roadmap for targeting the largely untapped PPI class more successfully.

Recognizing Ion Ligand-Binding Residues by Random Forest Algorithm Based on Optimized Dihedral Angle

The prediction of ion ligand–binding residues in protein sequences is a challenging work that contributes to understand the specific functions of proteins in life processes. In this article, we selected binding residues of 14 ion ligands as research objects, including four acid radical ion ligands and 10 metal ion ligands. Based on the amino acid sequence information, we selected the composition and position conservation information of amino acids, the predicted structural information, and physicochemical properties of amino acids as basic feature parameters. We then performed a statistical analysis and reclassification for dihedral angle and proposed new methods on the extraction of feature parameters. The methods mainly included applying information entropy on the extraction of polarization charge and hydrophilic–hydrophobic information of amino acids and using position weight matrices on the extraction of position conservation information. In the prediction model, we used the random forest algorithm and obtained better prediction results than previous works. With the independent test, the Matthew's correlation coefficient and accuracy of 10 metal ion ligand–binding residues were larger than 0.07 and 52%, respectively; the corresponding evaluation values of four acid radical ion ligand–binding residues were larger than 0.15 and 86%, respectively. Further, we classified and combined the phi and psi angles and optimized prediction model for each ion ligand–binding residue.

Fragment-Based Ligand-Protein Contact Statistics: Application to Docking Simulations

In this work, the information contained in the contacts between fragments of small-molecule ligands and protein residues has been collected and its exploitability has been verified by using the scoring of docking simulations as a test case for bringing about a proof of concept. Contact statistics between small-molecule fragments and binding site residues were collected and analyzed using a dataset composed of 200,000+ binding sites and associated ligands, derived from the database of the LIBRA ligand binding site recognition software, as a starting point. The fragments were generated by applying the decomposition algorithm implemented in BRICS. A simple "potential" based on the contact frequencies was tested against the CASF-2013 benchmark; its performance was then evaluated through the rescoring of docking poses generated for the DUD-E dataset. The results obtained indicate that this approach, its simplicity notwithstanding, yields promising results that are comparable, and in some cases, superior, to those obtained with other, more complex scoring functions.

RaptorX-Angle: real-value prediction of protein backbone dihedral angles through a hybrid method of clustering and deep learning

Background

Protein dihedral angles provide a detailed description of protein local conformation. Predicted dihedral angles can be used to narrow down the conformational space of the whole polypeptide chain significantly, thus aiding protein tertiary structure prediction. However, direct angle prediction from sequence alone is challenging.

Results

In this article, we present a novel method (named RaptorX-Angle) to predict real-valued angles by combining clustering and deep learning. Tested on a subset of PDB25 and the targets in the latest two Critical Assessment of protein Structure Prediction (CASP), our method outperforms the existing state-of-art method SPIDER2 in terms of Pearson Correlation Coefficient (PCC) and Mean Absolute Error (MAE). Our result also shows approximately linear relationship between the real prediction errors and our estimated bounds. That is, the real prediction error can be well approximated by our estimated bounds.

Conclusions

Our study provides an alternative and more accurate prediction of dihedral angles, which may facilitate protein structure prediction and functional study.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2065-x) contains supplementary material, which is available to authorized users.

Systematic Identification of Machine-Learning Models Aimed to Classify Critical Residues for Protein Function from Protein Structure

Protein structure and protein function should be related, yet the nature of this relationship remains unsolved. Mapping the critical residues for protein function with protein structure features represents an opportunity to explore this relationship, yet two important limitations have precluded a proper analysis of the structure-function relationship of proteins: (i) the lack of a formal definition of what critical residues are and (ii) the lack of a systematic evaluation of methods and protein structure features. To address this problem, here we introduce an index to quantify the protein-function criticality of a residue based on experimental data and a strategy aimed to optimize both, descriptors of protein structure (physicochemical and centrality descriptors) and machine learning algorithms, to minimize the error in the classification of critical residues. We observed that both physicochemical and centrality descriptors of residues effectively relate protein structure and protein function, and that physicochemical descriptors better describe critical residues. We also show that critical residues are better classified when residue criticality is considered as a binary attribute (i.e., residues are considered critical or not critical). Using this binary annotation for critical residues 8 models rendered accurate and non-overlapping classification of critical residues, confirming the multi-factorial character of the structure-function relationship of proteins.