Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

The roles of cell size and cell number in determining ovariole number in Drosophila

All insect ovaries are composed of functional units called ovarioles, which contain sequentially developing egg chambers. The number of ovarioles varies between and within species. Ovariole number is an important determinant of fecundity and thus affects individual fitness. Although Drosophila oogenesis has been intensively studied, the genetic and cellular basis for determination of ovariole number remains unknown. Ovariole formation begins during larval development with the morphogenesis of terminal filament cells (TFCs) into stacks called terminal filaments (TFs). We induced changes in ovariole number in Drosophila melanogaster by genetically altering cell size and cell number in the TFC population, and analyzed TF morphogenesis in these ovaries to understand the cellular basis for the changes in ovariole number. Increasing TFC size contributed to higher ovariole number by increasing TF number. Similarly, increasing total TFC number led to higher ovariole number via an increase in TF number. By analyzing ovarian morphogenesis in another Drosophila species we showed that TFC number regulation is a target of evolutionary change that affects ovariole number. In contrast, temperature-dependent plasticity in ovariole number was due to changes in cell-cell sorting during TF morphogenesis, rather than changes in cell size or cell number. We have thus identified two distinct developmental processes that regulate ovariole number: establishment of total TFC number, and TFC sorting during TF morphogenesis. Our data suggest that the genetic changes underlying species-specific ovariole number may alter the total number of TFCs available to contribute to TF formation. This work provides for the first time specific and quantitative developmental tools to investigate the evolution of a highly conserved reproductive structure.     

Utilization of Sialylated Glycans as Coreceptors Enhances the Neurovirulence of Serotype 3 Reovirus

Mammalian reoviruses display serotype-specific patterns of tropism and disease in the murine central nervous system (CNS) attributable to polymorphisms in viral attachment protein σ1. While all reovirus serotypes use junctional adhesion molecule-A as a cellular receptor, they differ in their utilization of carbohydrate coreceptors. This observation raises the possibility that carbohydrate binding by σ1 influences reovirus pathology in the CNS. In this study, we sought to define the function of carbohydrate binding in reovirus neuropathogenesis. Newborn mice were inoculated intramuscularly with wild-type strain type 3 Dearing (T3D) and T3D-σ1R202W, a point mutant T3D derivative that does not bind sialic acid (SA). Infected mice were monitored for survival, and viral loads at the sites of primary and secondary replication were quantified. Fewer mice inoculated with the wild-type virus survived in comparison to those inoculated with the mutant virus. The wild-type virus also produced higher titers in the spinal cord and brain at late times postinoculation but lower titers in the liver in comparison to those produced by the mutant virus. In addition, the wild-type virus was more virulent and produced higher titers in the brain than the mutant following intracranial inoculation. These animal infectivity studies suggest that T3D-σ1R202W harbors a defect in neural growth. Concordantly, compared with the wild-type virus, the mutant virus displayed a decreased capacity to infect and replicate in primary cultures of cortical neurons, a property dependent on cell surface SA. These results suggest that SA binding enhances the kinetics of reovirus replication in neural tissues and highlight a functional role for sialylated glycans as reovirus coreceptors in the CNS.     

Use of specific-pathogen-free (SPF) rhesus macaques to better model oral pediatric cytomegalovirus infection

Background Congenital human cytomegalovirus (HCMV) infection can result in lifelong neurological deficits. Seronegative pregnant woman often acquire primary HCMV from clinically asymptomatic, but HCMV‐shedding children. Methods Potential age‐related differences in viral and immune parameters of primary RhCMV infection were examined in an oral rhesus CMV infection model in specific pathogen free macaques. RhCMV shedding was measured by real time PCR in plasma, saliva and urine. Immune parameters, including neutralizing and binding antibodies and RhCMV‐specific T cell responses, were assessed in longitudinally collected blood samples. Results The oral RhCMV infection model in infant SPF rhesus macaques demonstrated that (i) the susceptibility to oral RhCMV infection declines with age, and (ii) infant macaques shed RhCMV more persistently and at higher titers compared to adult macaques. (iii) Conclusions The oral infant RhCMV infection model appears to reflect viral pathogenesis in human HCMV‐infected children. Larger studies are needed to define immune parameters associated with better control of RhCMV in adult compared to young animals.     

Analysis of IL28B variants in an Egyptian population defines the 20 kilobases minimal region involved in spontaneous clearance of hepatitis C virus

Spontaneous clearance of hepatitis C virus (HCV) occurs in ~30% of acute infections. Host genetics play a major role in HCV clearance, with a strong effect of single nucleotide polymorphisms (SNPs) of the IL28B gene already found in different populations, mostly infected with viral genotypes 1 and 3. Egypt has the highest prevalence of HCV infection in the world, which is mostly due to viral genotype 4. We investigated the role of several IL28B SNPs in HCV spontaneous clearance in an Egyptian population. We selected nine SNPs within the IL28B genomic region covering the linkage disequilibrium (LD) block known to be associated with HCV clearance in European populations. These SNPs were genotyped in 261 HCV-infected Egyptian subjects (130 with spontaneous clearance and 131 with chronic infection). The most associated SNPs were rs12979860 (P = 1.6 × 10(-7)) and the non-synonymous IL28B SNP, rs8103142 (P = 1.6 × 10(-7)). Interestingly, three SNPs at the two bounds of the region were monomorphic, reducing the size of the LD block in which the causal variants are potentially located to ～20 kilobases. HCV clearance in Egypt was associated with a region of IL28B smaller than that identified in European populations, and involved the non-synonymous IL28B SNP, rs8103142.     

Mobilizing the genome of Lepidoptera through novel sequence gains and end creation by non-autonomous Lep1 Helitrons

Transposable elements (TEs) can affect the structure of genomes through their acquisition and transposition of novel DNA sequences. The 134-bp repetitive elements, Lep1, are conserved non-autonomous Helitrons in lepidopteran genomes that have characteristic 5'-CT and 3'-CTAY nucleotide termini, a 3'-terminal hairpin structure, a 5'- and 3'-subterminal inverted repeat (SIR), and integrations that occur between AT or TT nucleotides. Lep1 Helitrons have acquired and propagated sequences downstream of their 3'-CTAY termini that are 57-344-bp in length and have termini composed of a 3'-CTRR preceded by a 3'-hairpin structure and a region complementary to the 5'-SIR (3'-SIRb). Features of both the Lep1 Helitron and multiple acquired sequences indicate that secondary structures at the 3'-terminus may have a role in rolling circle replication or genome integration mechanisms, and are a prerequisite for novel end creation by Helitron-like TEs. The preferential integration of Lep1 Helitrons in proximity to gene-coding regions results in the creation of genetic novelty that is shown to impact gene structure and function through the introduction of novel exon sequence (exon shuffling). These findings are important in understanding the structural requirements of genomic DNA sequences that are acquired and transposed by Helitron-like TEs.     

Effects of S-propargyl-cysteine (SPRC) in caerulein-induced acute pancreatitis in mice

Hydrogen sulfide (H(2)S), a novel gaseous messenger, is synthesized endogenously from L-cysteine by two pyridoxal-5'-phosphate-dependent enzymes, cystathionine β-synthase (CBS) and cystathionine γ-lyase (CSE). S-propargyl-cysteine (SPRC) is a slow H(2)S releasing drug that provides cysteine, a substrate of CSE. The present study was aimed to investigate the effects of SPRC in an in vivo model of acute pancreatitis (AP) in mice. AP was induced in mice by hourly caerulein injections (50 μg/kg) for 10 hours. Mice were treated with SPRC (10 mg/kg) or vehicle (distilled water). SPRC was administered either 12 h before or 3 h before the induction of pancreatitis. Mice were sacrificed 1 h after the last caerulein injection. Blood, pancreas and lung tissues were collected and processed to measure the plasma amylase, plasma H(2)S, myeloperoxidase (MPO) activities and cytokine levels in pancreas and lung. The results revealed that significant reduction of inflammation, both in pancreas and lung was associated with SPRC given 3 h prior to the induction of AP. Furthermore, the beneficial effects of SPRC were associated with reduction of pancreatic and pulmonary pro-inflammatory cytokines and increase of anti-inflammatory cytokine. SPRC administered 12 h before AP induction did not cause significant improvement in pancreatic and lung inflammation. Plasma H(2)S concentration showed significant difference in H(2)S levels between control, vehicle and SPRC (administered 3 h before AP) treatment groups. In conclusion, these data provide evidence for protective effects of SPRC in AP possibly by virtue of its slow release of endogenous H(2)S.