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

Computational identification of ubiquitylation sites from protein sequences

<p>Abstract</p> <p>Background</p> <p>Ubiquitylation plays an important role in regulating protein functions. Recently, experimental methods were developed toward effective identification of ubiquitylation sites. To efficiently explore more undiscovered ubiquitylation sites, this study aims to develop an accurate sequence-based prediction method to identify promising ubiquitylation sites.</p> <p>Results</p> <p>We established an ubiquitylation dataset consisting of 157 ubiquitylation sites and 3676 putative non-ubiquitylation sites extracted from 105 proteins in the UbiProt database. This study first evaluates promising sequence-based features and classifiers for the prediction of ubiquitylation sites by assessing three kinds of features (amino acid identity, evolutionary information, and physicochemical property) and three classifiers (support vector machine, <it>k</it>-nearest neighbor, and NaïveBayes). Results show that the set of used 531 physicochemical properties and support vector machine (SVM) are the best kind of features and classifier respectively that their combination has a prediction accuracy of 72.19% using leave-one-out cross-validation.</p> <p>Consequently, an informative physicochemical property mining algorithm (IPMA) is proposed to select an informative subset of 531 physicochemical properties. A prediction system UbiPred was implemented by using an SVM with the feature set of 31 informative physicochemical properties selected by IPMA, which can improve the accuracy from 72.19% to 84.44%. To further analyze the informative physicochemical properties, a decision tree method C5.0 was used to acquire if-then rule-based knowledge of predicting ubiquitylation sites. UbiPred can screen promising ubiquitylation sites from putative non-ubiquitylation sites using prediction scores. By applying UbiPred, 23 promising ubiquitylation sites were identified from an independent dataset of 3424 putative non-ubiquitylation sites, which were also validated by using the obtained prediction rules.</p> <p>Conclusion</p> <p>We have proposed an algorithm IPMA for mining informative physicochemical properties from protein sequences to build an SVM-based prediction system UbiPred. UbiPred can predict ubiquitylation sites accompanied with a prediction score each to help biologists in identifying promising sites for experimental verification. UbiPred has been implemented as a web server and is available at <url>http://iclab.life.nctu.edu.tw/ubipred</url>.</p

Ocean and land forcing of the record-breaking Dust Bowl heat waves across central United States

International audienceThe severe drought of the 1930s Dust Bowl decade coincided with record-breaking summer heatwaves that contributed to the socioeconomic and ecological disaster over North America's Great Plains. It remains unresolved to what extent these exceptional heatwaves, hotter than in historically forced coupled climate model simulations, were forced by sea surface temperatures (SSTs) and exacerbated through human-induced deterioration of land cover. Here we show, using an atmospheric-only model, that anomalously warm North Atlantic SSTs enhance heatwave activity through an association with drier spring conditions resulting from weaker moisture transport. Model devegetation simulations, that represent the widespread exposure of bare soil in the 1930s, suggest human activity fueled stronger and more frequent heatwaves through greater evaporative drying in the warmer months. This study highlights the potential for the amplification of naturally occurring extreme events like droughts by vegetation feedbacks to create more extreme heatwaves in a warmer world

Multislice CT angiography in the selection of patients with ruptured intracranial aneurysms suitable for clipping or coiling

Introduction We sought to establish whether CT angiography (CTA) can be applied to the planning and performance of clipping or coiling in ruptured intracranial aneurysms without recourse to intraarterial digital subtraction angiography (IA-DSA). Methods Over the period April 2003 to January 2006 in all patients presenting with a subarachnoid haemorrhage CTA was performed primarily. If CTA demonstrated an aneurysm, coiling or clipping was undertaken. IA-DSA was limited to patients with negative or inconclusive CTA findings. We compared CTA images with findings at surgery or coiling in patients with positive CTA findings and in patients with negative and inconclusive findings in whom IA-DSA had been performed. Results In this study, 224 consecutive patients (mean age 52.7 years, 135 women) were included. In 133 patients (59%) CTA demonstrated an aneurysm, and CTA was followed directly by neurosurgical (n=55) or endovascular treatment (n=78). In 31 patients (14%) CTA findings were categorized as inconclusive, and in 60 (27%) CTA findings were negative. One patient received surgical treatment on the basis of false-positive CTA findings. In 17 patients in whom CTA findings were inconclusive, IA-DSA provided further diagnostic information required for correct patient selection for any therapy. Five ruptured aneurysms in patients with a nonperimesencephalic SAH were negative on CTA, and four of these were also false-negative on IA-DSA. On a patient basis the positive predictive value, negative predictive value, sensitivity, specificity and accuracy of CTA for symptomatic aneurysms were 99%, 90%, 96%, 98% and 96%, respectively. Conclusion CTA should be used as the first diagnostic modality in the selection of patients for surgical or endovascular treatment of ruptured intracranial aneurysms. If CTA renders inconclusive results, IA-DSA should be performed. With negative CTA results the complementary value of IA-DSA is marginal. IA-DSA is not needed in patients with negative CTA and classic perimesencephalic SAH. Repeat IA-DSA or CTA should still be performed in patients with a nonperimesencephalic SAH

Amerindian Helicobacter pylori Strains Go Extinct, as European Strains Expand Their Host Range

We studied the diversity of bacteria and host in the H. pylori-human model. The human indigenous bacterium H. pylori diverged along with humans, into African, European, Asian and Amerindian groups. Of these, Amerindians have the least genetic diversity. Since niche diversity widens the sets of resources for colonizing species, we predicted that the Amerindian H. pylori strains would be the least diverse. We analyzed the multilocus sequence (7 housekeeping genes) of 131 strains: 19 cultured from Africans, 36 from Spanish, 11 from Koreans, 43 from Amerindians and 22 from South American Mestizos. We found that all strains that had been cultured from Africans were African strains (hpAfrica1), all from Spanish were European (hpEurope) and all from Koreans were hspEAsia but that Amerindians and Mestizos carried mixed strains: hspAmerind and hpEurope strains had been cultured from Amerindians and hpEurope and hpAfrica1 were cultured from Mestizos. The least genetically diverse H. pylori strains were hspAmerind. Strains hpEurope were the most diverse and showed remarkable multilocus sequence mosaicism (indicating recombination). The lower genetic structure in hpEurope strains is consistent with colonization of a diversity of hosts. If diversity is important for the success of H. pylori, then the low diversity of Amerindian strains might be linked to their apparent tendency to disappear. This suggests that Amerindian strains may lack the needed diversity to survive the diversity brought by non-Amerindian hosts

Gene expression clines reveal local adaptation and associated trade-offs at a continental scale

Local adaptation, where fitness in one environment comes at a cost in another, should lead to spatial variation in trade-offs between life history traits and may be critical for population persistence. Recent studies have sought genomic signals of local adaptation, but often have been limited to laboratory populations representing two environmentally different locations of a species' distribution. We measured gene expression, as a proxy for fitness, in males of Drosophila subobscura, occupying a 20° latitudinal and 11 °C thermal range. Uniquely, we sampled six populations and studied both common garden and semi-natural responses to identify signals of local adaptation. We found contrasting patterns of investment: transcripts with expression positively correlated to latitude were enriched for metabolic processes, expressed across all tissues whereas negatively correlated transcripts were enriched for reproductive processes, expressed primarily in testes. When using only the end populations, to compare our results to previous studies, we found that locally adaptive patterns were obscured. While phenotypic trade-offs between metabolic and reproductive functions across widespread species are well-known, our results identify underlying genetic and tissue responses at a continental scale that may be responsible for this. This may contribute to understanding population persistence under environmental change

Bacteria Modulate the CD8+ T Cell Epitope Repertoire of Host Cytosol-Exposed Proteins to Manipulate the Host Immune Response

The main adaptive immune response to bacteria is mediated by B cells and CD4+ T-cells. However, some bacterial proteins reach the cytosol of host cells and are exposed to the host CD8+ T-cells response. Both gram-negative and gram-positive bacteria can translocate proteins to the cytosol through type III and IV secretion and ESX-1 systems, respectively. The translocated proteins are often essential for the bacterium survival. Once injected, these proteins can be degraded and presented on MHC-I molecules to CD8+ T-cells. The CD8+ T-cells, in turn, can induce cell death and destroy the bacteria's habitat. In viruses, escape mutations arise to avoid this detection. The accumulation of escape mutations in bacteria has never been systematically studied. We show for the first time that such mutations are systematically present in most bacteria tested. We combine multiple bioinformatic algorithms to compute CD8+ T-cell epitope libraries of bacteria with secretion systems that translocate proteins to the host cytosol. In all bacteria tested, proteins not translocated to the cytosol show no escape mutations in their CD8+ T-cell epitopes. However, proteins translocated to the cytosol show clear escape mutations and have low epitope densities for most tested HLA alleles. The low epitope densities suggest that bacteria, like viruses, are evolutionarily selected to ensure their survival in the presence of CD8+ T-cells. In contrast with most other translocated proteins examined, Pseudomonas aeruginosa's ExoU, which ultimately induces host cell death, was found to have high epitope density. This finding suggests a novel mechanism for the manipulation of CD8+ T-cells by pathogens. The ExoU effector may have evolved to maintain high epitope density enabling it to efficiently induce CD8+ T-cell mediated cell death. These results were tested using multiple epitope prediction algorithms, and were found to be consistent for most proteins tested

Current provision of simulation in the UK and Republic of Ireland trauma and orthopaedic specialist training: a national survey.

Aims: The primary aim of the survey was to map the current provision of simulation training within UK and Republic of Ireland (RoI) trauma and orthopaedic (T&O) specialist training programmes to inform future design of a simulation based-curriculum. The secondary aims were to characterize; the types of simulation offered to trainees by stage of training, the sources of funding for simulation, the barriers to providing simulation in training, and to measure current research activity assessing the educational impact of simulation. Methods: The development of the survey was a collaborative effort between the authors and the British Orthopaedic Association Simulation Group. The survey items were embedded in the Performance and Opportunity Dashboard, which annually audits quality in training across several domains on behalf of the Speciality Advisory Committee (SAC). The survey was sent via email to the 30 training programme directors in March 2019. Data were retrieved and analyzed at the Warwick Clinical Trials Unit, UK. Results: Overall, 28 of 30 programme directors completed the survey (93%). 82% of programmes had access to high-fidelity simulation facilities such as cadaveric laboratories. More than half (54%) had access to a non-technical skills simulation training. Less than half (43%) received centralized funding for simulation, a third relied on local funding such as the departmental budget, and there was a heavy reliance on industry sponsorship to partly or wholly fund simulation training (64%). Provision was higher in the mid-stages (ST3-5) compared to late-stages (ST6-8) of training, and was formally timetabled in 68% of prostgrammes. There was no assessment of the impact of simulation training using objective behavioural measures or real-world clinical outcomes. Conclusion: There is currently widespread, but variable, provision of simulation in T&O training in the UK and RoI, which is likely to expand further with the new curriculum. It is important that research activity into the impact of simulation training continues, to develop an evidence base to support investment in facilities and provision