Fr-TM-align: a new protein structural alignment method based on fragment alignments and the TM-score

©2008 Pandit and Skolnick; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This article is available from: http://www.biomedcentral.com/1471-2105/9/531doi:10.1186/1471-2105-9-531Background: Protein tertiary structure comparisons are employed in various fields of contemporary structural biology. Most structure comparison methods involve generation of an initial seed alignment, which is extended and/or refined to provide the best structural superposition between a pair of protein structures as assessed by a structure comparison metric. One such metric, the TM-score, was recently introduced to provide a combined structure quality measure of the coordinate root mean square deviation between a pair of structures and coverage. Using the TM-score, the TM-align structure alignment algorithm was developed that was often found to have better accuracy and coverage than the most commonly used structural alignment programs; however, there were a number of situations when this was not true. Results: To further improve structure alignment quality, the Fr-TM-align algorithm has been developed where aligned fragment pairs are used to generate the initial seed alignments that are then refined using dynamic programming to maximize the TM-score. For the assessment of the structural alignment quality from Fr-TM-align in comparison to other programs such as CE and TMalign, we examined various alignment quality assessment scores such as PSI and TM-score. The assessment showed that the structural alignment quality from Fr-TM-align is better in comparison to both CE and TM-align. On average, the structural alignments generated using Fr-TM-align have a higher TM-score (~9%) and coverage (~7%) in comparison to those generated by TM-align. Fr- TM-align uses an exhaustive procedure to generate initial seed alignments. Hence, the algorithm is computationally more expensive than TM-align. Conclusion: Fr-TM-align, a new algorithm that employs fragment alignment and assembly provides better structural alignments in comparison to TM-align. The source code and executables of Fr- TM-align are freely downloadable at: http://cssb.biology.gatech.edu/skolnick/files/FrTMalign/

A complex interaction between glycine/NMDA receptors and serotonergic/noradrenergic antidepressants in the forced swim test in mice

Both clinical and preclinical studies demonstrate the antidepressant activity of the functional NMDA receptor antagonists. In this study, we assessed the effects of two glycine/NMDA receptor ligands, namely L-701,324 (antagonist) and d-cycloserine (a partial agonist) on the action of antidepressant drugs with different pharmacological profiles in the forced swim test in mice. Swim sessions were conducted by placing mice individually in glass cylinders filled with warmed water for 6 min. The duration of behavioral immobility during the last 4 min of the test was evaluated. The locomotor activity of mice was measured with photoresistor actimeters. L-701,324 and d-cycloserine given with reboxetine (administered in subeffective doses) did not change the behavior of animals in the forced swim test. A potentiating effect was seen when both tested glycine site ligands were given concomitantly with imipramine or fluoxetine in this test. The lesion of noradrenaline nerve terminals produced by DSP-4 neither altered the baseline activity nor influenced the antidepressant-like action of L-701,324 or d-cycloserine. The depletion of serotonin by p-CPA did not alter baseline activity in the forced swim test. However, it completely antagonized the antidepressant-like action produced by L-701,324 and d-cycloserine. Moreover, the antidepressant-like effects of imipramine, fluoxetine and reboxetine were abolished by d-serine, a full agonist of glycine/NMDA receptors. The present study demonstrates that glycine/NMDA receptor functional antagonists enhance the antidepressant-like action of serotonin, but not noradrenaline-based antidepressants and such their activity seems to depend on serotonin rather than noradrenaline pathway

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

HMMerThread: Detecting Remote, Functional Conserved Domains in Entire Genomes by Combining Relaxed Sequence-Database Searches with Fold Recognition

Conserved domains in proteins are one of the major sources of functional information for experimental design and genome-level annotation. Though search tools for conserved domain databases such as Hidden Markov Models (HMMs) are sensitive in detecting conserved domains in proteins when they share sufficient sequence similarity, they tend to miss more divergent family members, as they lack a reliable statistical framework for the detection of low sequence similarity. We have developed a greatly improved HMMerThread algorithm that can detect remotely conserved domains in highly divergent sequences. HMMerThread combines relaxed conserved domain searches with fold recognition to eliminate false positive, sequence-based identifications. With an accuracy of 90%, our software is able to automatically predict highly divergent members of conserved domain families with an associated 3-dimensional structure. We give additional confidence to our predictions by validation across species. We have run HMMerThread searches on eight proteomes including human and present a rich resource of remotely conserved domains, which adds significantly to the functional annotation of entire proteomes. We find ∼4500 cross-species validated, remotely conserved domain predictions in the human proteome alone. As an example, we find a DNA-binding domain in the C-terminal part of the A-kinase anchor protein 10 (AKAP10), a PKA adaptor that has been implicated in cardiac arrhythmias and premature cardiac death, which upon stress likely translocates from mitochondria to the nucleus/nucleolus. Based on our prediction, we propose that with this HLH-domain, AKAP10 is involved in the transcriptional control of stress response. Further remotely conserved domains we discuss are examples from areas such as sporulation, chromosome segregation and signalling during immune response. The HMMerThread algorithm is able to automatically detect the presence of remotely conserved domains in proteins based on weak sequence similarity. Our predictions open up new avenues for biological and medical studies. Genome-wide HMMerThread domains are available at http://vm1-hmmerthread.age.mpg.de

Signaling pathways responsible for the rapid antidepressant-like effects of a GluN2A-preferring NMDA receptor antagonist

In a previous study we found that the preferring GluN2A receptor antagonist, NVP-AAM077, elicited rapid antidepressant-like effects in the forced swim test that was related to the release of glutamate and serotonin in the medial prefrontal cortex. In the present work we sought to examine the duration of this behavioral effect as well as the molecular readouts involved. Our results showed that NVP-AAM077 reduced the immobility in the forced swim test 30?min and 24?h after its administration. However, this effect waned 7 days later. The rapid antidepressant-like response seems to be associated with increases in the GluA1 subunit of ?-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors, mammalian target of rapamycin (mTOR) signaling, glia markers such as glial fibrillary acidic protein (GFAP) and excitatory amino acid transporter 1 (EAAT1), and a rapid mobilization of intracellular stores of brain-derived neurotrophic factor (BDNF) in the medial prefrontal cortex.Acknowledgements: M.G.-S. was recipient of a contract from the Sistema Nacional de Garantía Juvenil co-funded by the European Social Fund. We also thank Novartis for the generous gift of NVP-AAM077. This work was supported by the Instituto de Salud Carlos III, Subdirección General del Evaluación y Fomento de la Investigación (FIS Grants PI13/00038 and PI16/00217) that were co-funded by the European Regional Development Fund (‘A way to build Europe’). Funding from the Centro de Investigación Biomédica en Red de Salud Mental (CIBERSAM), Instituto de Salud Carlos III is also acknowledged

New statistical potential for quality assessment of protein models and a survey of energy functions

Abstract Background Scoring functions, such as molecular mechanic forcefields and statistical potentials are fundamentally important tools in protein structure modeling and quality assessment. Results The performances of a number of publicly available scoring functions are compared with a statistical rigor, with an emphasis on knowledge-based potentials. We explored the effect on accuracy of alternative choices for representing interaction center types and other features of scoring functions, such as using information on solvent accessibility, on torsion angles, accounting for secondary structure preferences and side chain orientation. Partially based on the observations made, we present a novel residue based statistical potential, which employs a shuffled reference state definition and takes into account the mutual orientation of residue side chains. Atom- and residue-level statistical potentials and Linux executables to calculate the energy of a given protein proposed in this work can be downloaded from http://www.fiserlab.org/potentials. Conclusions Among the most influential terms we observed a critical role of a proper reference state definition and the benefits of including information about the microenvironment of interaction centers. Molecular mechanical potentials were also tested and found to be over-sensitive to small local imperfections in a structure, requiring unfeasible long energy relaxation before energy scores started to correlate with model quality.</p

Neurogenic inflammation after traumatic brain injury and its potentiation of classical inflammation

Background: The neuroinflammatory response following traumatic brain injury (TBI) is known to be a key secondary injury factor that can drive ongoing neuronal injury. Despite this, treatments that have targeted aspects of the inflammatory pathway have not shown significant efficacy in clinical trials. Main body: We suggest that this may be because classical inflammation only represents part of the story, with activation of neurogenic inflammation potentially one of the key initiating inflammatory events following TBI. Indeed, evidence suggests that the transient receptor potential cation channels (TRP channels), TRPV1 and TRPA1, are polymodal receptors that are activated by a variety of stimuli associated with TBI, including mechanical shear stress, leading to the release of neuropeptides such as substance P (SP). SP augments many aspects of the classical inflammatory response via activation of microglia and astrocytes, degranulation of mast cells, and promoting leukocyte migration. Furthermore, SP may initiate the earliest changes seen in blood-brain barrier (BBB) permeability, namely the increased transcellular transport of plasma proteins via activation of caveolae. This is in line with reports that alterations in transcellular transport are seen first following TBI, prior to decreases in expression of tight-junction proteins such as claudin-5 and occludin. Indeed, the receptor for SP, the tachykinin NK1 receptor, is found in caveolae and its activation following TBI may allow influx of albumin and other plasma proteins which directly augment the inflammatory response by activating astrocytes and microglia. Conclusions: As such, the neurogenic inflammatory response can exacerbate classical inflammation via a positive feedback loop, with classical inflammatory mediators such as bradykinin and prostaglandins then further stimulating TRP receptors. Accordingly, complete inhibition of neuroinflammation following TBI may require the inhibition of both classical and neurogenic inflammatory pathways.Frances Corrigan, Kimberley A. Mander, Anna V. Leonard and Robert Vin