Protein Dynamics and Contact Topology Reveal Protein–DNA Binding Orientation

Theoretical study of the adiponectin receptors: binding site characterization and molecular dynamics of possible ligands for drug design

The two adiponectin receptors (AdipoR1 and AdipoR2) have been implicated in glucose and lipid regulation involved in several metabolic pathologies including type II diabetes. Their exact biochemical functions and mechanisms remain poorly understood. Moreover, these receptors do not yet have data on possible co-crystallized active ligands. In this study, we applied different computational methodologies to address three main unanswered questions: first, the localization and validation of possible binding sites; second, the generation of novel ligands with amenable characteristics to target the receptors; and third, the determination of important chemical interactions between the ligands and the receptors. Computational analysis of the binding site reveals that the residues triad R267, F271, and Y310 could be responsible for changes in the spatial arrangement and geometry of the binding pocket in AdipoR1. Molecular docking results in high docking scores of − 13.6 and − 16.5 kcal/mol for the top best ligands in AdipoR1 and AdipoR2 respectively. Finally, molecular dynamics suggests that hydrolytic activity may be possible with these compounds and that this reaction could be mediated by aspartic acid residues. The two adiponectin receptors have an endogenous protein ligand, adiponectin. However the synthesis is expensive and technically challenging. Although some debatable agonists have been proposed investigations of suitable synthetic ligands are indeed, very much needed for targeting these receptors and their associate pathologies and metabolic pathways. Furthermore, these findings provide a framework for further biochemical investigations of amenable compounds for drug discovery in order to target these receptors and their associated pathologies.

Molecular dynamics simulations and linear response theories jointly describe biphasic responses of myoglobin relaxation and reveal evolutionarily conserved frequent communicators

In this study, we provide a time-dependent mechanical model, taking advantage of molecular dynamics simulations, quasiharmonic analysis of molecular dynamics trajectories, and time-dependent linear response theories to describe vibrational energy redistribution within the protein matrix. The theoretical description explained the observed biphasic responses of specific residues in myoglobin to CO-photolysis and photoexcitation on heme. The fast responses were found to be triggered by impulsive forces and propagated mainly by principal modes <40 cm⁻¹. The predicted fast responses for individual atoms were then used to study signal propagation within the protein matrix and signals were found to propagate ~8 times faster across helices (4076 m/s) than within the helices, suggesting the importance of tertiary packing in the sensitivity of proteins to external perturbations. We further developed a method to integrate multiple intramolecular signal pathways and discover frequent “communicators”. These communicators were found to be evolutionarily conserved including those distant from the heme.

An Efficient Timer and Sizer of Biomacromolecular Motions

Life ticks as fast as how proteins move. Computationally expensive molecular dynamics simulation has been the only theoretical tool to gauge the time and sizes of these motions, though barely to their slowest ends. Here, we convert a computationally cheap elastic network model (ENM) into a molecular timer and sizer to gauge the slowest functional motions of structured biomolecules. Quasi-harmonic analysis, fluctuation profile matching, and the Wiener-Khintchine theorem are used to define the "time periods," t, for anharmonic principal components (PCs), which are validated by nuclear magnetic resonance (NMR) order parameters. The PCs with their respective "time periods" are mapped to the eigenvalues (λ_ENM) of the corresponding ENM modes. Thus, the power laws t(ns) = 56.1λ_ENM^-1.6 and σ²(Ų) = 32.7λ_ENM^-3.0 can be established allowing the characterization of the timescales of NMR-resolved conformers, crystallographic anisotropic displacement parameters, and important ribosomal motions, as well as motional sizes of the latter.

DynOmics: dynamics of structural proteome and beyond

DynOmics (dynomics.pitt.edu) is a portal developed to leverage rapidly growing structural proteomics data by efficiently and accurately evaluating the dynamics of structurally resolved systems, from individual molecules to large complexes and assemblies, in the context of their physiological environment. At the core of the portal is a newly developed server, ENM 1.0, which permits users to efficiently generate information on the collective dynamics of any structure in PDB format, user-uploaded or database-retrieved. ENM 1.0 integrates two widely used elastic network models (ENMs)—the Gaussian Network Model (GNM) and the Anisotropic Network Model (ANM), extended to take account of molecular environment. It enables users to assess potentially functional sites, signal transduction or allosteric communication mechanisms, and protein–protein and protein–DNA interaction poses, in addition to delivering ensembles of accessible conformers reconstructed at atomic details based on the global modes of motions predicted by the ANM. The ‘environment’ is defined in a flexible manner, from lipid bilayer and crystal contacts, to substrate or ligands bound to a protein, or surrounding subunits in a multimeric structure or assembly. User-friendly interactive features permit users to easily visualize how the environment alter the intrinsic dynamics of the query systems. ENM 1.0 can be accessed at http://enm.pitt.edu/ or http://dyn.life.nthu.edu.tw/oENM/.

An SVM-based method for assessment of transcription factor-DNA complex models

Background

Atomic details of protein-DNA complexes can provide insightful information for better understanding of the function and binding specificity of DNA binding proteins. In addition to experimental methods for solving protein-DNA complex structures, protein-DNA docking can be used to predict native or near-native complex models. A docking program typically generates a large number of complex conformations and predicts the complex model(s) based on interaction energies between protein and DNA. However, the prediction accuracy is hampered by current approaches to model assessment, especially when docking simulations fail to produce any near-native models.

Results

We present here a Support Vector Machine (SVM)-based approach for quality assessment of the predicted transcription factor (TF)-DNA complex models. Besides a knowledge-based protein-DNA interaction potential DDNA3, we applied several structural features that have been shown to play important roles in binding specificity between transcription factors and DNA molecules to quality assessment of complex models. To address the issue of unbalanced positive and negative cases in the training dataset, we applied hard-negative mining, an iterative training process that selects an initial training dataset by combining all of the positive cases and a random sample from the negative cases. Results show that the SVM model greatly improves prediction accuracy (84.2%) over two knowledge-based protein-DNA interaction potentials, orientation potential (60.8%) and DDNA3 (68.4%). The improvement is achieved through reducing the number of false positive predictions, especially for the hard docking cases, in which a docking algorithm fails to produce any near-native complex models.

Conclusions

A learning-based SVM scoring model with structural features for specific protein-DNA binding and an atomic-level protein-DNA interaction potential DDNA3 significantly improves prediction accuracy of complex models by successfully identifying cases without near-native structural models.

Electronic supplementary material

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