A network model predicts the intensity of residue-protein thermal coupling.

MOTIVATION: The flow of vibrational energy in proteins has been shown not to obey expectations for isotropic media. The existence of preferential pathways for energy transport, with probable connections to allostery mechanisms, has been repeatedly demonstrated. Here, we investigate whether, by representing a set of protein structures as networks of interacting amino acid residues, we are able to model heat diffusion and predict residue-protein vibrational couplings, as measured by the Anisotropic Thermal Diffusion (ATD) computational protocol of modified molecular dynamics simulations. RESULTS: We revisit the structural rationales for the precise definition of a contact between amino acid residues. Using this definition to describe a set of proteins as contact networks where each node corresponds to a residue, we show that node centrality, particularly closeness centrality and eigenvector centrality , correlates to the strength of the vibrational coupling of each residue to the rest of the structure. We then construct an analytically solvable model of heat diffusion on a network, whose solution incorporates an explicit dependence on the connectivity of the heated node, as described by a perturbed graph Laplacian Matrix. AVAILABILITY AND IMPLEMENTATION: An implementation of the described model is available at http://leandro.iqm.unicamp.br/atd-scripts . CONTACT: leandro@iqm.unicamp.br.

Ligand- and receptor-based docking with LiBELa.

Methodologies on molecular docking are constantly improving. The problem consists on finding an optimal interplay between the computational cost and a satisfactory physical description of ligand-receptor interaction. In pursuit of an advance in current methods we developed a mixed docking approach combining ligand- and receptor-based strategies in a docking engine, where tridimensional descriptors for shape and charge distribution of a reference ligand guide the initial placement of the docking molecule and an interaction energy-based global minimization follows. This hybrid docking was evaluated with soft-core and force field potentials taking into account ligand pose and scoring. Our approach was found to be competitive to a purely receptor-based dock resulting in improved logAUC values when evaluated with DUD and DUD-E. Furthermore, the smoothed potential as evaluated here, was not advantageous when ligand binding poses were compared to experimentally determined conformations. In conclusion we show that a combination of ligand- and receptor-based strategy docking with a force field energy model results in good reproduction of binding poses and enrichment of active molecules against decoys. This strategy is implemented in our tool, LiBELa, available to the scientific community.