Introduction to protein folding for physicists

The prediction of the three-dimensional native structure of proteins from the knowledge of their amino acid sequence, known as the protein folding problem, is one of the most important yet unsolved issues of modern science. Since the conformational behaviour of flexible molecules is nothing more than a complex physical problem, increasingly more physicists are moving into the study of protein systems, bringing with them powerful mathematical and computational tools, as well as the sharp intuition and deep images inherent to the physics discipline. This work attempts to facilitate the first steps of such a transition. In order to achieve this goal, we provide an exhaustive account of the reasons underlying the protein folding problem enormous relevance and summarize the present-day status of the methods aimed to solving it. We also provide an introduction to the particular structure of these biological heteropolymers, and we physically define the problem stating the assumptions behind this (commonly implicit) definition. Finally, we review the 'special flavor' of statistical mechanics that is typically used to study the astronomically large phase spaces of macromolecules. Throughout the whole work, much material that is found scattered in the literature has been put together here to improve comprehension and to serve as a handy reference.Comment: 53 pages, 18 figures, the figures are at a low resolution due to arXiv restrictions, for high-res figures, go to http://www.pabloechenique.co

A Bayesian test for the appropriateness of a model in the biomagnetic inverse problem

This paper extends the work of Clarke [1] on the Bayesian foundations of the biomagnetic inverse problem. It derives expressions for the expectation and variance of the a posteriori source current probability distribution given a prior source current probability distribution, a source space weight function and a data set. The calculation of the variance enables the construction of a Bayesian test for the appropriateness of any source model that is chosen as the a priori infomation. The test is illustrated using both simulated (multi-dipole) data and the results of a study of early latency processing of images of human faces. [1] C.J.S. Clarke. Error estimates in the biomagnetic inverse problem. Inverse Problems, 10:77--86, 1994.Comment: 13 pages, 16 figures. Submitted to Inverse Problem

Angular Forces Around Transition Metals in Biomolecules

Quantum-mechanical analysis based on an exact sum rule is used to extract an semiclassical angle-dependent energy function for transition metal ions in biomolecules. The angular dependence is simple but different from existing classical potentials. Comparison of predicted energies with a computer-generated database shows that the semiclassical energy function is remarkably accurate, and that its angular dependence is optimal.Comment: Tex file plus 4 postscript figure

Interpretation of the Precision Matrix and Its Application in Estimating Sparse Brain Connectivity during Sleep Spindles from Human Electrocorticography Recordings

The correlation method from brain imaging has been used to estimate functional connectivity in the human brain. However, brain regions might show very high correlation even when the two regions are not directly connected due to the strong interaction of the two regions with common input from a third region. One previously proposed solution to this problem is to use a sparse regularized inverse covariance matrix or precision matrix (SRPM) assuming that the connectivity structure is sparse. This method yields partial correlations to measure strong direct interactions between pairs of regions while simultaneously removing the influence of the rest of the regions, thus identifying regions that are conditionally independent. To test our methods, we first demonstrated conditions under which the SRPM method could indeed find the true physical connection between a pair of nodes for a spring-mass example and an RC circuit example. The recovery of the connectivity structure using the SRPM method can be explained by energy models using the Boltzmann distribution. We then demonstrated the application of the SRPM method for estimating brain connectivity during stage 2 sleep spindles from human electrocorticography (ECoG) recordings using an 8 x 8 electrode array. The ECoG recordings that we analyzed were from a 32-year-old male patient with long-standing pharmaco-resistant left temporal lobe complex partial epilepsy. Sleep spindles were automatically detected using delay differential analysis and then analyzed with SRPM and the Louvain method for community detection. We found spatially localized brain networks within and between neighboring cortical areas during spindles, in contrast to the case when sleep spindles were not present

Direction and magnitude of nicotine effects on the fMRI BOLD response are related to nicotine effects on behavioral performance

Considerable variability across individuals has been reported in both the behavioral and fMRI blood oxygen level-dependent (BOLD) response to nicotine. We aimed to investigate (1) whether there is a heterogeneous effect of nicotine on behavioral and BOLD responses across participants and (2) if heterogeneous BOLD responses are associated with behavioral performance measures. In this double-blind, placebo-controlled, cross-over study, 41 healthy participants (19 smokers)—drawn from a larger population-based sample—performed a visual oddball task after acute challenge with 1 mg nasal nicotine. fMRI data and reaction time were recorded during performance of the task. Across the entire group of subjects, we found increased activation in the anterior cingulate cortex, middle frontal gyrus, superior temporal gyrus, post-central gyrus, planum temporal and frontal pole in the nicotine condition compared with the placebo condition. However, follow-up analyses of this difference in activation between the placebo and nicotine conditions revealed that some participants showed an increase in activation while others showed a decrease in BOLD activation from the placebo to the nicotine condition. A reduction of BOLD activation from placebo to nicotine was associated with a decrease in reaction time and reaction time variability and vice versa, suggesting that it is the direction of BOLD response to nicotine which is related to task performance. We conclude that the BOLD response to nicotine is heterogeneous and that the direction of response to nicotine should be taken into account in future pharmaco-fMRI research on the central action of nicotine

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation

A Chemical Genomics Approach to Drug Reprofiling in Oncology: Antipsychotic Drug Risperidone as a Potential Adenocarcinoma Treatment

Drug reprofiling is emerging as an effective paradigm for discovery of cancer treatments. Herein, an antipsychotic drug is immobilised using the Magic Tag® chemical genomics tool and screened against a T7 bacteriophage displayed library of polypeptides from Drosophila melanogaster, as a whole genome model, to uncover an interaction with a section of 17-β-HSD10, a proposed prostate cancer target. A computational study and enzyme inhibition assay with full length human 17-β-HSD10 identifies risperidone as a drug reprofiling candidate. When formulated with rumenic acid, risperidone slows proliferation of PC3 prostate cancer cells in vitro and retards PC3 prostate cancer tumour growth in vivo in xenografts in mice, presenting an opportunity to reprofile risperidone as a cancer treatment

Open Babel: An open chemical toolbox

Background: A frequent problem in computational modeling is the interconversion of chemical structures between different formats. While standard interchange formats exist (for example, Chemical Markup Language) and de facto standards have arisen (for example, SMILES format), the need to interconvert formats is a continuing problem due to the multitude of different application areas for chemistry data, differences in the data stored by different formats (0D versus 3D, for example), and competition between software along with a lack of vendorneutral formats. Results: We discuss, for the first time, Open Babel, an open-source chemical toolbox that speaks the many languages of chemical data. Open Babel version 2.3 interconverts over 110 formats. The need to represent such a wide variety of chemical and molecular data requires a library that implements a wide range of cheminformatics algorithms, from partial charge assignment and aromaticity detection, to bond order perception and canonicalization. We detail the implementation of Open Babel, describe key advances in the 2.3 release, and outline a variety of uses both in terms of software products and scientific research, including applications far beyond simple format interconversion. Conclusions: Open Babel presents a solution to the proliferation of multiple chemical file formats. In addition, it provides a variety of useful utilities from conformer searching and 2D depiction, to filtering, batch conversion, and substructure and similarity searching. For developers, it can be used as a programming library to handle chemical data in areas such as organic chemistry, drug design, materials science, and computational chemistry. It is freely available under an open-source license fro

Right hemisphere has the last laugh: neural dynamics of joke appreciation

Understanding a joke relies on semantic, mnemonic, inferential, and emotional contributions from multiple brain areas. Anatomically constrained magnetoencephalography (aMEG) combining high-density whole-head MEG with anatomical magnetic resonance imaging allowed us to estimate where the humor-specific brain activations occur and to understand their temporal sequence. Punch lines provided either funny, not funny (semantically congruent), or nonsensical (incongruent) replies to joke questions. Healthy subjects rated them as being funny or not funny. As expected, incongruous endings evoke the largest N400m in left-dominant temporo-prefrontal areas, due to integration difficulty. In contrast, funny punch lines evoke the smallest N400m during this initial lexical–semantic stage, consistent with their primed “surface congruity” with the setup question. In line with its sensitivity to ambiguity, the anteromedial prefrontal cortex may contribute to the subsequent “second take” processing, which, for jokes, presumably reflects detection of a clever “twist” contained in the funny punch lines. Joke-selective activity simultaneously emerges in the right prefrontal cortex, which may lead an extended bilateral temporo-frontal network in establishing the distant unexpected creative coherence between the punch line and the setup. This progression from an initially promising but misleading integration from left frontotemporal associations, to medial prefrontal ambiguity evaluation and right prefrontal reprocessing, may reflect the essential tension and resolution underlying humor