Free energy basin-hopping

A global optimisation scheme is presented using basin-hopping with the acceptance criterion based on approximate free energy for the corresponding local minima of the potential energy. The method is illustrated for atomic and colloidal clusters and peptides to examine how the predicted global free energy minimum changes with temperature. Using estimates for the local free energies based on harmonic vibrational densities of states provides a computationally effective framework for predicting trends in structure at finite temperature. The resulting scheme represents a powerful tool for exploration of energy landscapes throughout molecular science.We are grateful to the EPSRC and the ERC for financial support under grants EP/1001352/1 and 267369, respectively.This is the final version of the article. It first appeared at http://www.sciencedirect.com/science/article/pii/S0009261415000974#

Yukawa potentials in systems with partial periodic boundary conditions I : Ewald sums for quasi-two dimensional systems

Yukawa potentials are often used as effective potentials for systems as colloids, plasmas, etc. When the Debye screening length is large, the Yukawa potential tends to the non-screened Coulomb potential ; in this small screening limit, or Coulomb limit, the potential is long ranged. As it is well known in computer simulation, a simple truncation of the long ranged potential and the minimum image convention are insufficient to obtain accurate numerical data on systems. The Ewald method for bulk systems, i.e. with periodic boundary conditions in all three directions of the space, has already been derived for Yukawa potential [cf. Y., Rosenfeld, {\it Mol. Phys.}, \bm{88}, 1357, (1996) and G., Salin and J.-M., Caillol, {\it J. Chem. Phys.}, \bm{113}, 10459, (2000)], but for systems with partial periodic boundary conditions, the Ewald sums have only recently been obtained [M., Mazars, {\it J. Chem. Phys.}, {\bf 126}, 056101 (2007)]. In this paper, we provide a closed derivation of the Ewald sums for Yukawa potentials in systems with periodic boundary conditions in only two directions and for any value of the Debye length. A special attention is paid to the Coulomb limit and its relation with the electroneutrality of systems.Comment: 40 pages, 5 figures and 4 table

An avian influenza A(H11N1) virus from a wild aquatic bird revealing a unique Eurasian-American genetic reassortment

Influenza surveillance in different wild bird populations is critical for understanding the persistence, transmission and evolution of these viruses. Avian influenza (AI) surveillance was undertaken in wild migratory and resident birds during the period 2007–2008, in view of the outbreaks of highly pathogenic AI (HPAI) H5N1 in poultry in India since 2006. In this study, we present the whole genome sequence data along with the genetic and virological characterization of an Influenza A(H11N1) virus isolated from wild aquatic bird for the first time from India. The virus was low pathogenicity and phylogenetic analysis revealed that it was distinct from reported H11N1 viruses. The hemagglutinin (HA) gene showed maximum similarity with A/semipalmatedsandpiper/Delaware/2109/2000 (H11N6) and A/shorebird/Delaware/236/2003(H11N9) while the neuraminidase (NA) gene showed maximum similarity with A/duck/Mongolia/540/2001(H1N1). The virus thus possessed an HA gene of the American lineage. The NA and other six genes were of the Eurasian lineage and showed closer relatedness to non-H11 viruses. Such a genetic reassortment is unique and interesting, though the pathways leading to its emergence and its future persistence in the avian reservoir is yet to be fully established

Structural and Functional Roles of Coevolved Sites in Proteins

Understanding the residue covariations between multiple positions in protein families is very crucial and can be helpful for designing protein engineering experiments. These simultaneous changes or residue coevolution allow protein to maintain its overall structural-functional integrity while enabling it to acquire specific functional modifications. Despite the significant efforts in the field there is still controversy in terms of the preferable locations of coevolved residues on different regions of protein molecules, the strength of coevolutionary signal and role of coevolution in functional diversification.In this paper we study the scale and nature of residue coevolution in maintaining the overall functionality and structural integrity of proteins. We employed a large scale study to investigate the structural and functional aspects of coevolved residues. We found that the networks representing the coevolutionary residue connections within our dataset are in general of 'small-world' type as they have clustering coefficient values higher than random networks and also show smaller mean shortest path lengths similar and/or lower than random and regular networks. We also found that altogether 11% of functionally important sites are coevolved with any other sites. Active sites are found more frequently to coevolve with any other sites (15%) compared to protein (11%) and ligand (9%) binding sites. Metal binding and active sites are also found to be more frequently coevolved with other metal binding and active sites, respectively. Analysis of the coupling between coevolutionary processes and the spatial distribution of coevolved sites reveals that a high fraction of coevolved sites are located close to each other. Moreover, approximately 80% of charge compensatory substitutions within coevolved sites are found at very close spatial proximity (<or= 5A), pointing to the possible preservation of salt bridges in evolution.Our findings show that a noticeable fraction of functionally important sites undergo coevolution and also point towards compensatory substitutions as a probable coevolutionary mechanism within spatially proximal coevolved functional sites

Comparison of Human Memory CD8 T Cell Responses to Adenoviral Early and Late Proteins in Peripheral Blood and Lymphoid Tissue

Treatment of invasive adenovirus (Ad) disease in hematopoietic stem cell transplant (SCT) recipients with capsid protein hexon-specific donor T cells is under investigation. We propose that cytotoxic T cells (CTLs) targeted to the late protein hexon may be inefficient in vivo because the early Ad protein E3-19K downregulates HLA class I antigens in infected cells. In this study, CD8+ T cells targeted to highly conserved HLA A2-restricted epitopes from the early regulatory protein DNA polymerase (P-977) and late protein hexon (H-892) were compared in peripheral blood (PB) and tonsils of naturally infected adults. In tonsils, epitope-specific pentamers detected a significantly higher frequency of P-977+CD8+ T cells compared to H-892+CD8+ T cells; this trend was reversed in PB. Tonsil epitope-specific CD8+ T cells expressed IFN-γ and IL-2 but not perforin or TNF-α, whereas PB T cells were positive for IFN-γ, TNF-α, and perforin. Tonsil epitope-specific T cells expressed lymphoid homing marker CCR7 and exhibited lower levels of the activation marker CD25 but higher proliferative potential than PB T cells. Finally, in parallel with the kinetics of mRNA expression, P-977-specific CTLs lysed targets as early as 8 hrs post infection. In contrast, H-892-specific CTLs did not kill unless infected fibroblasts were pretreated with IFN-γ to up regulate HLA class I antigens, and cytotoxicity was delayed until 16–24 hours. These data show that, in contrast to hexon CTLs, central memory type DNA polymerase CTLs dominate the lymphoid compartment and kill fibroblasts earlier after infection without requiring exogenous IFN-γ. Thus, use of CTLs targeted to both early and late Ad proteins may improve the efficacy of immunotherapy for life-threatening Ad disease in SCT recipients

Neural Computation via Neural Geometry: A Place Code for Inter-whisker Timing in the Barrel Cortex?

The place theory proposed by Jeffress (1948) is still the dominant model of how the brain represents the movement of sensory stimuli between sensory receptors. According to the place theory, delays in signalling between neurons, dependent on the distances between them, compensate for time differences in the stimulation of sensory receptors. Hence the location of neurons, activated by the coincident arrival of multiple signals, reports the stimulus movement velocity. Despite its generality, most evidence for the place theory has been provided by studies of the auditory system of auditory specialists like the barn owl, but in the study of mammalian auditory systems the evidence is inconclusive. We ask to what extent the somatosensory systems of tactile specialists like rats and mice use distance dependent delays between neurons to compute the motion of tactile stimuli between the facial whiskers (or ‘vibrissae’). We present a model in which synaptic inputs evoked by whisker deflections arrive at neurons in layer 2/3 (L2/3) somatosensory ‘barrel’ cortex at different times. The timing of synaptic inputs to each neuron depends on its location relative to sources of input in layer 4 (L4) that represent stimulation of each whisker. Constrained by the geometry and timing of projections from L4 to L2/3, the model can account for a range of experimentally measured responses to two-whisker stimuli. Consistent with that data, responses of model neurons located between the barrels to paired stimulation of two whiskers are greater than the sum of the responses to either whisker input alone. The model predicts that for neurons located closer to either barrel these supralinear responses are tuned for longer inter-whisker stimulation intervals, yielding a topographic map for the inter-whisker deflection interval across the surface of L2/3. This map constitutes a neural place code for the relative timing of sensory stimuli

Autism as a disorder of neural information processing: directions for research and targets for therapy

The broad variation in phenotypes and severities within autism spectrum disorders suggests the involvement of multiple predisposing factors, interacting in complex ways with normal developmental courses and gradients. Identification of these factors, and the common developmental path into which theyfeed, is hampered bythe large degrees of convergence from causal factors to altered brain development, and divergence from abnormal brain development into altered cognition and behaviour. Genetic, neurochemical, neuroimaging and behavioural findings on autism, as well as studies of normal development and of genetic syndromes that share symptoms with autism, offer hypotheses as to the nature of causal factors and their possible effects on the structure and dynamics of neural systems. Such alterations in neural properties may in turn perturb activity-dependent development, giving rise to a complex behavioural syndrome many steps removed from the root causes. Animal models based on genetic, neurochemical, neurophysiological, and behavioural manipulations offer the possibility of exploring these developmental processes in detail, as do human studies addressing endophenotypes beyond the diagnosis itself

The Chromatin Remodeler SPLAYED Regulates Specific Stress Signaling Pathways

Organisms are continuously exposed to a myriad of environmental stresses. Central to an organism's survival is the ability to mount a robust transcriptional response to the imposed stress. An emerging mechanism of transcriptional control involves dynamic changes in chromatin structure. Alterations in chromatin structure are brought about by a number of different mechanisms, including chromatin modifications, which covalently modify histone proteins; incorporation of histone variants; and chromatin remodeling, which utilizes ATP hydrolysis to alter histone-DNA contacts. While considerable insight into the mechanisms of chromatin remodeling has been gained, the biological role of chromatin remodeling complexes beyond their function as regulators of cellular differentiation and development has remained poorly understood. Here, we provide genetic, biochemical, and biological evidence for the critical role of chromatin remodeling in mediating plant defense against specific biotic stresses. We found that the Arabidopsis SWI/SNF class chromatin remodeling ATPase SPLAYED (SYD) is required for the expression of selected genes downstream of the jasmonate (JA) and ethylene (ET) signaling pathways. SYD is also directly recruited to the promoters of several of these genes. Furthermore, we show that SYD is required for resistance against the necrotrophic pathogen Botrytis cinerea but not the biotrophic pathogen Pseudomonas syringae. These findings demonstrate not only that chromatin remodeling is required for selective pathogen resistance, but also that chromatin remodelers such as SYD can regulate specific pathways within biotic stress signaling networks

Deletion of the Chd6 exon 12 affects motor coordination

Members of the CHD protein family play key roles in gene regulation through ATP-dependent chromatin remodeling. This is facilitated by chromodomains that bind histone tails, and by the SWI2/SNF2-like ATPase/helicase domain that remodels chromatin by moving histones. Chd6 is ubiquitously expressed in both mouse and human, with the highest levels of expression in the brain. The Chd6 gene contains 37 exons, of which exons 12-19 encode the highly conserved ATPase domain. To determine the biological role of Chd6, we generated mouse lines with a deletion of exon 12. Chd6 without exon 12 is expressed at normal levels in mice, and Chd6 Exon 12 −/− mice are viable, fertile, and exhibit no obvious morphological or pathological phenotype. Chd6 Exon 12 −/− mice lack coordination as revealed by sensorimotor analysis. Further behavioral testing revealed that the coordination impairment was not due to muscle weakness or bradykinesia. Histological analysis of brain morphology revealed no differences between Chd6 Exon 12 −/− mice and wild-type (WT) controls. The location of CHD6 on human chromosome 20q12 is overlapped by the linkage map regions of several human ataxias, including autosomal recessive infantile cerebellar ataxia (SCAR6), a nonprogressive cerebrospinal ataxia. The genomic location, expression pattern, and ataxic phenotype of Chd6 Exon 12 −/− mice indicate that mutations within CHD6 may be responsible for one of these ataxias

Erytrocyte membrane anionic charge in type 2 diabetic patients with retinopathy

BACKGROUND: The Steno hypothesis states that changes in basement membrane anionic charge leads to diabetic microvascular complications. In diabetic nephropathy, loss of basement membrane glycosaminoglycans and the association between glomerular basement membrane heparan sulphate and proteinuria has been documented. A correlation between erythrocyte surface and the glomerular capillary wall charges has also been observed. The aim of this study is to evaluate the relationship between retinopathy and erythrocyte anionic charge and urinary glycosaminoglycan excretion in type 2 diabetic patients. METHODS: 49 subjects (58 ± 7 yrs, M/F 27/22) with type 2 diabetes with proliferative retinopathy (n = 13), nonproliferative retinopathy (n = 13) and without retinopathy (n = 23) were included in the study. 38 healthy subjects were selected as control group (57 ± 5 yrs, M/F 19/19). Erythrocyte anionic charge (EAC) was determined by the binding of the cationic dye, alcian blue. Urinary glycosaminoglycan and microalbumin excretion were measured. RESULTS: EAC was significantly decreased in diabetic patients with retinopathy (255 ± 30 ng alcian blue/10(6 )RBC, 312 ± 30 ng alcian blue/10(6 )RBC for diabetic and control groups respectively, p < 0.001). We did not observe an association between urinary GAG and microalbumin excretion and diabetic retinopathy. EAC is found to be negatively corralated with microalbuminuria in all groups. CONCLUSIONS: We conclude that type 2 diabetic patients with low erythrocyte anionic charge are associated with diabetic retinopathy. Reduction of negative charge of basement membranes may indicate general changes in microvasculature rather than retinopathy. More prospective and large studies needs to clarify the role of glycosaminoglycans on progression of retinopathy in type 2 diabetic patients