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

Analysis of Gene Regulatory Networks in the Mammalian Circadian Rhythm

Circadian rhythm is fundamental in regulating a wide range of cellular, metabolic, physiological, and behavioral activities in mammals. Although a small number of key circadian genes have been identified through extensive molecular and genetic studies in the past, the existence of other key circadian genes and how they drive the genomewide circadian oscillation of gene expression in different tissues still remains unknown. Here we try to address these questions by integrating all available circadian microarray data in mammals. We identified 41 common circadian genes that showed circadian oscillation in a wide range of mouse tissues with a remarkable consistency of circadian phases across tissues. Comparisons across mouse, rat, rhesus macaque, and human showed that the circadian phases of known key circadian genes were delayed for 4–5 hours in rat compared to mouse and 8–12 hours in macaque and human compared to mouse. A systematic gene regulatory network for the mouse circadian rhythm was constructed after incorporating promoter analysis and transcription factor knockout or mutant microarray data. We observed the significant association of cis-regulatory elements: EBOX, DBOX, RRE, and HSE with the different phases of circadian oscillating genes. The analysis of the network structure revealed the paths through which light, food, and heat can entrain the circadian clock and identified that NR3C1 and FKBP/HSP90 complexes are central to the control of circadian genes through diverse environmental signals. Our study improves our understanding of the structure, design principle, and evolution of gene regulatory networks involved in the mammalian circadian rhythm

Modulation of gene expression in heart and liver of hibernating black bears (Ursus americanus)

<p>Abstract</p> <p>Background</p> <p>Hibernation is an adaptive strategy to survive in highly seasonal or unpredictable environments. The molecular and genetic basis of hibernation physiology in mammals has only recently been studied using large scale genomic approaches. We analyzed gene expression in the American black bear, <it>Ursus americanus</it>, using a custom 12,800 cDNA probe microarray to detect differences in expression that occur in heart and liver during winter hibernation in comparison to summer active animals.</p> <p>Results</p> <p>We identified 245 genes in heart and 319 genes in liver that were differentially expressed between winter and summer. The expression of 24 genes was significantly elevated during hibernation in both heart and liver. These genes are mostly involved in lipid catabolism and protein biosynthesis and include RNA binding protein motif 3 (<it>Rbm3</it>), which enhances protein synthesis at mildly hypothermic temperatures. Elevated expression of protein biosynthesis genes suggests induction of translation that may be related to adaptive mechanisms reducing cardiac and muscle atrophies over extended periods of low metabolism and immobility during hibernation in bears. Coordinated reduction of transcription of genes involved in amino acid catabolism suggests redirection of amino acids from catabolic pathways to protein biosynthesis. We identify common for black bears and small mammalian hibernators transcriptional changes in the liver that include induction of genes responsible for fatty acid β oxidation and carbohydrate synthesis and depression of genes involved in lipid biosynthesis, carbohydrate catabolism, cellular respiration and detoxification pathways.</p> <p>Conclusions</p> <p>Our findings show that modulation of gene expression during winter hibernation represents molecular mechanism of adaptation to extreme environments.</p

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl

From computational discovery to experimental characterization of a high hole mobility organic crystal

For organic semiconductors to find ubiquitous electronics applications, the development of new materials with high mobility and air stability is critical. Despite the versatility of carbon, exploratory chemical synthesis in the vast chemical space can be hindered by synthetic and characterization difficulties. Here we show that in silico screening of novel derivatives of the dinaphtho[2,3-b:2′,3′-f]thieno[3,2-b]thiophene semiconductor with high hole mobility and air stability can lead to the discovery of a new high-performance semiconductor. On the basis of estimates from the Marcus theory of charge transfer rates, we identified a novel compound expected to demonstrate a theoretic twofold improvement in mobility over the parent molecule. Synthetic and electrical characterization of the compound is reported with single-crystal field-effect transistors, showing a remarkable saturation and linear mobility of 12.3 and 16 cm2 V−1 s−1, respectively. This is one of the very few organic semiconductors with mobility greater than 10 cm2 V−1 s−1 reported to date

Identification of a novel angiogenic peptide from periostin

Angiogenic peptides have therapeutic potential for the treatment of chronic ischemic diseases. Periostin, an extracellular matrix protein expressed in injured tissues, promotes angiogenesis and tissue repair. We previously reported that in vivo administration of both recombinant full-length protein and the first FAS I domain of periostin alleviated peripheral artery occlusive disease by stimulating the migration of humane endothelial colony forming cells (ECFCs) and subsequent angiogenesis. In the present study, we ascertained the peptide sequence responsible for the periostin-induced angiogenesis. By serial deletion mapping of the first FAS I domain, we identified a peptide sequence (amino acids 142-151) of periostin for stimulation of chemotactic migration, adhesion, proliferation and endothelial tube formation of human ECFCs in vitro. Chemotactic migration of ECFCs induced by the periostin peptide was blocked by pre-incubation with an anti-??5 integrin neutralizing antibody. Treatment of ECFCs with the periostin peptide led to phosphorylation of both AKT and ERK, and pretreatment of ECFCs with the MEK-ERK pathway inhibitor U0126 or the PI3K-AKT pathway inhibitors, LY294002 or Wortmannin, blocked the periostin peptide-stimulated migration of ECFCs. These results suggest that the synthetic periostin peptide can be applied for stimulating angiogenic and therapeutic potentials of ECFCs

Chemical Probes that Competitively and Selectively Inhibit Stat3 Activation

Signal transducer and activator of transcription (Stat) 3 is an oncogene constitutively activated in many cancer systems where it contributes to carcinogenesis. To develop chemical probes that selectively target Stat3, we virtually screened 920,000 small drug-like compounds by docking each into the peptide-binding pocket of the Stat3 SH2 domain, which consists of three sites—the pY-residue binding site, the +3 residue-binding site and a hydrophobic binding site, which served as a selectivity filter. Three compounds satisfied criteria of interaction analysis, competitively inhibited recombinant Stat3 binding to its immobilized pY-peptide ligand and inhibited IL-6-mediated tyrosine phosphorylation of Stat3. These compounds were used in a similarity screen of 2.47 million compounds, which identified 3 more compounds with similar activities. Examination of the 6 active compounds for the ability to inhibit IFN-γ-mediated Stat1 phosphorylation revealed that 5 of 6 were selective for Stat3. Molecular modeling of the SH2 domains of Stat3 and Stat1 bound to compound revealed that compound interaction with the hydrophobic binding site was the basis for selectivity. All 5 selective compounds inhibited nuclear-to-cytoplasmic translocation of Stat3, while 3 of 5 compounds induced apoptosis preferentially of breast cancer cell lines with constitutive Stat3 activation. Thus, virtual ligand screening of compound libraries that targeted the Stat3 pY-peptide binding pocket identified for the first time 3 lead compounds that competitively inhibited Stat3 binding to its pY-peptide ligand; these compounds were selective for Stat3 vs. Stat1 and induced apoptosis preferentially of breast cancer cells lines with constitutively activated Stat3

The complete mitochondrial genome of the citrus red mite Panonychus citri (Acari: Tetranychidae): high genome rearrangement and extremely truncated tRNAs

<p>Abstract</p> <p>Background</p> <p>The family Tetranychidae (Chelicerata: Acari) includes ~1200 species, many of which are of agronomic importance. To date, mitochondrial genomes of only two Tetranychidae species have been sequenced, and it has been found that these two mitochondrial genomes are characterized by many unusual features in genome organization and structure such as gene order and nucleotide frequency. The scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). Information on Tetranychidae mitochondrial genomes is quite important for phylogenetic evaluation and population genetics, as well as the molecular evolution of functional genes such as acaricide-resistance genes. In this study, we sequenced the complete mitochondrial genome of <it>Panonychus citri </it>(Family Tetranychidae), a worldwide citrus pest, and provide a comparison to other Acari.</p> <p>Results</p> <p>The mitochondrial genome of <it>P. citri </it>is a typical circular molecule of 13,077 bp, and contains the complete set of 37 genes that are usually found in metazoans. This is the smallest mitochondrial genome within all sequenced Acari and other Chelicerata, primarily due to the significant size reduction of protein coding genes (PCGs), a large rRNA gene, and the A + T-rich region. The mitochondrial gene order for <it>P. citri </it>is the same as those for <it>P. ulmi </it>and <it>Tetranychus urticae</it>, but distinctly different from other Acari by a series of gene translocations and/or inversions. The majority of the <it>P. citri </it>mitochondrial genome has a high A + T content (85.28%), which is also reflected by AT-rich codons being used more frequently, but exhibits a positive GC-skew (0.03). The Acari mitochondrial <it>nad1 </it>exhibits a faster amino acid substitution rate than other genes, and the variation of nucleotide substitution patterns of PCGs is significantly correlated with the G + C content. Most tRNA genes of <it>P. citri </it>are extremely truncated and atypical (44-65, 54.1 ± 4.1 bp), lacking either the T- or D-arm, as found in <it>P. ulmi</it>, <it>T. urticae</it>, and other Acariform mites.</p> <p>Conclusions</p> <p>The <it>P. citri </it>mitochondrial gene order is markedly different from those of other chelicerates, but is conserved within the family Tetranychidae indicating that high rearrangements have occurred after Tetranychidae diverged from other Acari. Comparative analyses suggest that the genome size, gene order, gene content, codon usage, and base composition are strongly variable among Acari mitochondrial genomes. While extremely small and unusual tRNA genes seem to be common for Acariform mites, further experimental evidence is needed.</p

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

First narrow-band search for continuous gravitational waves from known pulsars in advanced detector data

Spinning neutron stars asymmetric with respect to their rotation axis are potential sources of continuous gravitational waves for ground-based interferometric detectors. In the case of known pulsars a fully coherent search, based on matched filtering, which uses the position and rotational parameters obtained from electromagnetic observations, can be carried out. Matched filtering maximizes the signalto- noise (SNR) ratio, but a large sensitivity loss is expected in case of even a very small mismatch between the assumed and the true signal parameters. For this reason, narrow-band analysis methods have been developed, allowing a fully coherent search for gravitational waves from known pulsars over a fraction of a hertz and several spin-down values. In this paper we describe a narrow-band search of 11 pulsars using data from Advanced LIGO’s first observing run. Although we have found several initial outliers, further studies show no significant evidence for the presence of a gravitational wave signal. Finally, we have placed upper limits on the signal strain amplitude lower than the spin-down limit for 5 of the 11 targets over the bands searched; in the case of J1813-1749 the spin-down limit has been beaten for the first time. For an additional 3 targets, the median upper limit across the search bands is below the spin-down limit. This is the most sensitive narrow-band search for continuous gravitational waves carried out so far