Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer's, Parkinson's and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth

Crystal Structure of a Charge Engineered Human Lysozyme Having Enhanced Bactericidal Activity

Human lysozyme is a key component of the innate immune system, and recombinant forms of the enzyme represent promising leads in the search for therapeutic agents able to treat drug-resistant infections. The wild type protein, however, fails to participate effectively in clearance of certain infections due to inherent functional limitations. For example, wild type lysozymes are subject to electrostatic sequestration and inactivation by anionic biopolymers in the infected airway. A charge engineered variant of human lysozyme has recently been shown to possess improved antibacterial activity in the presence of disease associated inhibitory molecules. Here, the 2.04 Å crystal structure of this variant is presented along with an analysis that provides molecular level insights into the origins of the protein's enhanced performance. The charge engineered variant's two mutated amino acids exhibit stabilizing interactions with adjacent native residues, and from a global perspective, the mutations cause no gross structural perturbations or loss of stability. Importantly, the two substitutions dramatically expand the negative electrostatic potential that, in the wild type enzyme, is restricted to a small region near the catalytic residues. The net result is a reduction in the overall strength of the engineered enzyme's electrostatic potential field, and it appears that the specific nature of this remodeled field underlies the variant's reduced susceptibility to inhibition by anionic biopolymers

Production of phi mesons at mid-rapidity in sqrt(s_NN) = 200 GeV Au+Au collisions at RHIC

We present the first results of meson production in the K^+K^- decay channel from Au+Au collisions at sqrt(s_NN) = 200 GeV as measured at mid-rapidity by the PHENIX detector at RHIC. Precision resonance centroid and width values are extracted as a function of collision centrality. No significant variation from the PDG accepted values is observed. The transverse mass spectra are fitted with a linear exponential function for which the derived inverse slope parameter is seen to be constant as a function of centrality. These data are also fitted by a hydrodynamic model with the result that the freeze-out temperature and the expansion velocity values are consistent with the values previously derived from fitting single hadron inclusive data. As a function of transverse momentum the collisions scaled peripheral.to.central yield ratio RCP for the is comparable to that of pions rather than that of protons. This result lends support to theoretical models which distinguish between baryons and mesons instead of particle mass for explaining the anomalous proton yield.Comment: 326 authors, 24 pages text, 23 figures, 6 tables, RevTeX 4. To be submitted to Physical Review C as a regular article. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Local conservation scores without a priori assumptions on neutral substitution rates

<p>Abstract</p> <p>Background</p> <p>Comparative genomics aims to detect signals of evolutionary conservation as an indicator of functional constraint. Surprisingly, results of the ENCODE project revealed that about half of the experimentally verified functional elements found in non-coding DNA were classified as unconstrained by computational predictions. Following this observation, it has been hypothesized that this may be partly explained by biased estimates on neutral evolutionary rates used by existing sequence conservation metrics. All methods we are aware of rely on a comparison with the neutral rate and conservation is estimated by measuring the deviation of a particular genomic region from this rate. Consequently, it is a reasonable assumption that inaccurate neutral rate estimates may lead to biased conservation and constraint estimates.</p> <p>Results</p> <p>We propose a conservation signal that is produced by local Maximum Likelihood estimation of evolutionary parameters using an optimized sliding window and present a Kullback-Leibler projection that allows multiple different estimated parameters to be transformed into a conservation measure. This conservation measure does not rely on assumptions about neutral evolutionary substitution rates and little a priori assumptions on the properties of the conserved regions are imposed. We show the accuracy of our approach (KuLCons) on synthetic data and compare it to the scores generated by state-of-the-art methods (phastCons, GERP, SCONE) in an ENCODE region. We find that KuLCons is most often in agreement with the conservation/constraint signatures detected by GERP and SCONE while qualitatively very different patterns from phastCons are observed. Opposed to standard methods KuLCons can be extended to more complex evolutionary models, e.g. taking insertion and deletion events into account and corresponding results show that scores obtained under this model can diverge significantly from scores using the simpler model.</p> <p>Conclusion</p> <p>Our results suggest that discriminating among the different degrees of conservation is possible without making assumptions about neutral rates. We find, however, that it cannot be expected to discover considerably different constraint regions than GERP and SCONE. Consequently, we conclude that the reported discrepancies between experimentally verified functional and computationally identified constraint elements are likely not to be explained by biased neutral rate estimates.</p

How to do an evaluation: pitfalls and traps

The recent literature is replete with papers evaluating computational tools (often those operating on 3D structures) for their performance in a certain set of tasks. Most commonly these papers compare a number of docking tools for their performance in cognate re-docking (pose prediction) and/or virtual screening. Related papers have been published on ligand-based tools: pose prediction by conformer generators and virtual screening using a variety of ligand-based approaches. The reliability of these comparisons is critically affected by a number of factors usually ignored by the authors, including bias in the datasets used in virtual screening, the metrics used to assess performance in virtual screening and pose prediction and errors in crystal structures used

How to do an evaluation: pitfalls and traps

The recent literature is replete with papers evaluating computational tools (often those operating on 3D structures) for their performance in a certain set of tasks. Most commonly these papers compare a number of docking tools for their performance in cognate re-docking (pose prediction) and/or virtual screening. Related papers have been published on ligand-based tools: pose prediction by conformer generators and virtual screening using a variety of ligand-based approaches. The reliability of these comparisons is critically affected by a number of factors usually ignored by the authors, including bias in the datasets used in virtual screening, the metrics used to assess performance in virtual screening and pose prediction and errors in crystal structures used

Comparative analyses of vertebrate posterior HoxD clusters reveal atypical cluster architecture in the caecilian Typhlonectes natans

<p>Abstract</p> <p>Background</p> <p>The posterior genes of the <it>HoxD </it>cluster play a crucial role in the patterning of the tetrapod limb. This region is under the control of a global, long-range enhancer that is present in all vertebrates. Variation in limb types, as is the case in amphibians, can probably not only be attributed to variation in <it>Hox </it>genes, but is likely to be the product of differences in gene regulation. With a collection of vertebrate genome sequences available today, we used a comparative genomics approach to study the posterior <it>HoxD </it>cluster of amphibians. A frog and a caecilian were included in the study to compare coding sequences as well as to determine the gain and loss of putative regulatory sequences.</p> <p>Results</p> <p>We sequenced the posterior end of the <it>HoxD </it>cluster of a caecilian and performed comparative analyses of this region using <it>HoxD </it>clusters of other vertebrates. We determined the presence of conserved non-coding sequences and traced gains and losses of these footprints during vertebrate evolution, with particular focus on amphibians. We found that the caecilian <it>HoxD </it>cluster is almost three times larger than its mammalian counterpart. This enlargement is accompanied with the loss of one gene and the accumulation of repeats in that area. A similar phenomenon was observed in the coelacanth, where a different gene was lost and expansion of the area where the gene was lost has occurred. At least one phylogenetic footprint present in all vertebrates was lost in amphibians. This conserved region is a known regulatory element and functions as a boundary element in neural tissue to prevent expression of <it>Hoxd </it>genes.</p> <p>Conclusion</p> <p>The posterior part of the <it>HoxD </it>cluster of <it>Typhlonectes natans </it>is among the largest known today. The loss of <it>Hoxd-12 </it>and the expansion of the intergenic region may exert an influence on the limb enhancer, by having to bypass a distance seven times that of regular <it>HoxD </it>clusters. Whether or not there is a correlation with the loss of limbs remains to be investigated. These results, together with data on other vertebrates show that the tetrapod <it>Hox </it>clusters are more variable than previously thought.</p

Dissecting the Transcriptional Regulatory Properties of Human Chromosome 16 Highly Conserved Non-Coding Regions

Non-coding DNA conservation across species has been often used as a predictor for transcriptional enhancer activity. However, only a few systematic analyses of the function of these highly conserved non-coding regions (HCNRs) have been performed. Here we use zebrafish transgenic assays to perform a systematic study of 113 HCNRs from human chromosome 16. By comparing transient and stable transgenesis, we show that the first method is highly inefficient, leading to 40% of false positives and 20% of false negatives. When analyzed in stable transgenic lines, a great majority of HCNRs were active in the central nervous system, although some of them drove expression in other organs such as the eye and the excretory system. Finally, by testing a fraction of the HCNRs lacking enhancer activity for in vivo insulator activity, we find that 20% of them may contain enhancer-blocking function. Altogether our data indicate that HCNRs may contain different types of cis-regulatory activity, including enhancer, insulators as well as other not yet discovered functions

Distinct Functional Constraints Partition Sequence Conservation in a cis-Regulatory Element

Different functional constraints contribute to different evolutionary rates across genomes. To understand why some sequences evolve faster than others in a single cis-regulatory locus, we investigated function and evolutionary dynamics of the promoter of the Caenorhabditis elegans unc-47 gene. We found that this promoter consists of two distinct domains. The proximal promoter is conserved and is largely sufficient to direct appropriate spatial expression. The distal promoter displays little if any conservation between several closely related nematodes. Despite this divergence, sequences from all species confer robustness of expression, arguing that this function does not require substantial sequence conservation. We showed that even unrelated sequences have the ability to promote robust expression. A prominent feature shared by all of these robustness-promoting sequences is an AT-enriched nucleotide composition consistent with nucleosome depletion. Because general sequence composition can be maintained despite sequence turnover, our results explain how different functional constraints can lead to vastly disparate rates of sequence divergence within a promoter

Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: conformational analysis and binding mode of multisite inhibitors

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer’s disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challeng- ing. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments—huprine and rhein— that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8–14, 154–169 and 307–318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable “floppy” pocket would enable the bind- ing of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligo- methylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications intro- duced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors