Characterizing Rhodopsin-Arrestin Interactions with the Fragment Molecular Orbital (FMO) Method

Arrestin binding to G protein-coupled receptors (GPCRs) plays a vital role in receptor signaling. Recently, the crystal structure of rhodopsin bound to activated visual arrestin was resolved using XFEL (X-ray free electron laser). However, even with the crystal structure in hand, our ability to understand GPCR-arrestin binding is limited by the availability of accurate tools to explore receptor-arrestin interactions. We applied fragment molecular orbital (FMO) method to explore the interactions formed between the residues of rhodopsin and arrestin. FMO enables ab initio approaches to be applied to systems that conventional quantum mechanical (QM) methods would be too compute-expensive. The FMO calculations detected 35 significant interactions involved in rhodopsin-arrestin binding formed by 25 residues of rhodopsin and 28 residues of arrestin. Two major regions of interaction were identified: at the C-terminal tail of rhodopsin (D330-S343) and where the "finger loop" (G69-T79) of arrestin directly inserts into rhodopsin active core. Out of these 35 interactions, 23 were mainly electrostatic and 12 hydrophobic in nature

Characterizing Protein-Protein Interactions with the Fragment Molecular Orbital Method

Proteins are vital components of living systems, serving as building blocks, molecular machines, enzymes, receptors, ion channels, sensors, and transporters. Protein-protein interactions (PPIs) are a key part of their function. There are more than 645,000 reported disease-relevant PPIs in the human interactome, but drugs have been developed for only 2% of these targets. The advances in PPI-focused drug discovery are highly dependent on the availability of structural data and accurate computational tools for analysis of this data. Quantum mechanical approaches are often too expensive computationally, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. FMO provides essential information for PPI drug discovery, namely, identification of key interactions formed between residues of two proteins, including their strength (in kcal/mol) and their chemical nature (electrostatic or hydrophobic). In this chapter, we have demonstrated how three different FMO-based approaches (pair interaction energy analysis (PIE analysis), subsystem analysis (SA) and analysis of protein residue networks (PRNs)) have been applied to study PPI in three protein-protein complexes

Preimaginal Stages of the Emerald Ash Borer, Agrilus planipennis Fairmaire (Coleoptera: Buprestidae): An Invasive Pest on Ash Trees (Fraxinus)

This study provides the most detailed description of the immature stages of Agrilus planipennis Fairmaire to date and illustrates suites of larval characters useful in distinguishing among Agrilus Curtis species and instars. Immature stages of eight species of Agrilus were examined and imaged using light and scanning electron microscopy. For A. planipennis all preimaginal stages (egg, instars I-IV, prepupa and pupa) were described. A combination of 14 character states were identified that serve to identify larvae of A. planipennis. Our results support the segregation of Agrilus larvae into two informal assemblages based on characters of the mouthparts, prothorax, and abdomen: the A. viridis and A. ater assemblages, with A. planipennis being more similar to the former. Additional evidence is provided in favor of excluding A. planipennis from the subgenus Uragrilus

Analyzing GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

G-protein-coupled receptors (GPCRs) have enormous physiological and biomedical importance, and therefore it is not surprising that they are the targets of many prescribed drugs. Further progress in GPCR drug discovery is highly dependent on the availability of protein structural information. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanics (QM) approaches are often too computationally expensive to be of practical use in time-sensitive situations, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed, and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule toward ligand binding, including an analysis of their chemical nature. Such information is essential for an efficient structure-based drug design (SBDD) process. In this chapter, we describe how to use FMO in the characterization of GPCR-ligand interactions

AGTR2 and sprint/power performance: a case-control replication study for rs11091046 polymorphism in two ethnicities

We aimed to replicate, in a specific athletic event cohort (only track and field) and in two different ethnicities (Japanese and East European, i.e. Russian and Polish), original findings showing the association of the angiotensin-II receptor type-2 gene (AGTR2) rs11091046 A>C polymorphism with athlete status. We compared genotypic frequencies of the AGTR2 rs11091046 polymorphism among 282 track and field sprint/ power athletes (200 men and 82 women), including several national record holders and Olympic medallists (214 Japanese, 68 Russian and Polish), and 2024 control subjects (842 men and 1182 women) (804 Japanese, 1220 Russian and Polish). In men, a meta-analysis from the two combined cohorts showed a significantly higher frequency of the C allele in athletes than in controls (odds ratio: 1.62, P=0.008, heterogeneity index I 2 =0%). With regard to respective cohorts, C allele frequency was higher in Japanese male athletes than in controls (67.7% vs. 55.9%, P=0.022), but not in Russian/Polish male athletes (61.9% vs. 51.0%, P=0.172). In women, no significant results were obtained by meta-analysis for the two cohorts combination (P=0.850). The AC genotype frequency was significantly higher in Russian/Polish women athletes than in controls (69.2% vs. 42.1%, P=0.022), but not in Japanese women athletes (P=0.226). Our results, in contrast to previous findings, suggested by meta-analysis that the C allele of the AGTR2 rs11091046 polymorphism is associated with sprint/ power track and field athlete status in men, but not in women

No evidence of a common DNA variant profile specific to world class endurance athletes

There are strong genetic components to cardiorespiratory fitness and its response to exercise training. It would be useful to understand the differences in the genomic profile of highly trained endurance athletes of world class caliber and sedentary controls. An international consortium (GAMES) was established in order to compare elite endurance athletes and ethnicity-matched controls in a case-control study design. Genome-wide association studies were undertaken on two cohorts of elite endurance athletes and controls (GENATHLETE and Japanese endurance runners), from which a panel of 45 promising markers was identified. These markers were tested for replication in seven additional cohorts of endurance athletes and controls: from Australia, Ethiopia, Japan, Kenya, Poland, Russia and Spain. The study is based on a total of 1520 endurance athletes (835 who took part in endurance events in World Championships and/or Olympic Games) and 2760 controls. We hypothesized that world-class athletes are likely to be characterized by an even higher concentration of endurance performance alleles and we performed separate analyses on this subsample. The meta-analysis of all available studies revealed one statistically significant marker (rs558129 at GALNTL6 locus, p = 0.0002), even after correcting for multiple testing. As shown by the low heterogeneity index (I2 = 0), all eight cohorts showed the same direction of association with rs558129, even though p-values varied across the individual studies. In summary, this study did not identify a panel of genomic variants common to these elite endurance athlete groups. Since GAMES was underpowered to identify alleles with small effect sizes, some of the suggestive leads identified should be explored in expanded comparisons of world-class endurance athletes and sedentary controls and in tightly controlled exercise training studies. Such studies have the potential to illuminate the biology not only of world class endurance performance but also of compromised cardiac functions and cardiometabolic diseases

The Earliest Evidence of Holometabolan Insect Pupation in Conifer Wood

Background: The pre-Jurassic record of terrestrial wood borings is poorly resolved, despite body fossil evidence of insect diversification among xylophilic clades starting in the late Paleozoic. Detailed analysis of borings in petrified wood provides direct evidence of wood utilization by invertebrate animals, which typically comprises feeding behaviors.\ud \ud Methodology/Principal Findings: We describe a U-shaped boring in petrified wood from the Late Triassic Chinle Formation of southern Utah that demonstrates a strong linkage between insect ontogeny and conifer wood resources. Xylokrypta durossi new ichnogenus and ichnospecies is a large excavation in wood that is backfilled with partially digested xylem, creating a secluded chamber. The tracemaker exited the chamber by way of a small vertical shaft. This sequence of behaviors is most consistent with the entrance of a larva followed by pupal quiescence and adult emergence — hallmarks of holometabolous insect ontogeny. Among the known body fossil record of Triassic insects, cupedid beetles (Coleoptera: Archostemata) are deemed the most plausible tracemakers of Xylokrypta, based on their body size and modern xylobiotic lifestyle.\ud \ud Conclusions/Significance: This oldest record of pupation in fossil wood provides an alternative interpretation to borings once regarded as evidence for Triassic bees. Instead Xylokrypta suggests that early archostematan beetles were leaders in exploiting wood substrates well before modern clades of xylophages arose in the late Mesozoic

Genome-wide association study identifies three novel genetic markers associated with elite endurance performance

To investigate the association between multiple single-nucleotide polymorphisms (SNPs), aerobic performance and elite endurance athlete status in Russians. By using GWAS approach, we examined the association between 1,140,419 SNPs and relative maximal oxygen consumption rate (VO2max) in 80 international-level Russian endurance athletes (46 males and 34 females). To validate obtained results, we further performed case-control studies by comparing the frequencies of the most significant SNPs (with P<10-5-10-8) between 218 endurance athletes and opposite cohorts (192 Russian controls, 1367 European controls, and 230 Russian power athletes). Initially, six 'endurance alleles' were identified showing discrete associations with •VO2maxboth in males and females. Next, case-control studies resulted in remaining three SNPs (NFIA-AS2 rs1572312, TSHR rs7144481, RBFOX1 rs7191721) associated with endurance athlete status. The C allele of the most significant SNP, rs1572312, was associated with high values of •VO2max(males: P=0.0051; females: P=0.0005). Furthermore, the frequency of the rs1572312 C allele was significantly higher in elite endurance athletes (95.5%) in comparison with non-elite endurance athletes (89.8%, P=0.0257), Russian (88.8%, P=0.007) and European (90.6%, P=0.0197) controls and power athletes (86.2%, P=0.0005). The rs1572312 SNP is located on the nuclear factor I A antisense RNA 2 (NFIA-AS2) gene which is supposed to regulate the expression of the NFIA gene (encodes transcription factor involved in activation of erythropoiesis and repression of the granulopoiesis). Our data show that the NFIA-AS2 rs1572312, TSHR rs7144481 and RBFOX1 rs7191721 polymorphisms are associated with aerobic performance and elite endurance athlete status

Expression profile of genes regulated by activity of the Na-H exchanger NHE1

BACKGROUND: In mammalian cells changes in intracellular pH (pH(i)), which are predominantly controlled by activity of plasma membrane ion exchangers, regulate a diverse range of normal and pathological cellular processes. How changes in pH(i )affect distinct cellular processes has primarily been determined by evaluating protein activities and we know little about how pH(i )regulates gene expression. RESULTS: A global profile of genes regulated in mammalian fibroblasts by decreased pH(i )induced by impaired activity of the plasma membrane Na-H exchanger NHE1 was characterized by using cDNA microarrays. Analysis of selected genes by quantitative RT-PCR, TaqMan, and immunoblot analyses confirmed results obtained from cDNA arrays. Consistent with established roles of pH(i )and NHE1 activity in cell proliferation and oncogenic transformation, grouping regulated genes into functional categories and biological pathways indicated a predominant number of genes with altered expression were associated with growth factor signaling, oncogenesis, and cell cycle progression. CONCLUSION: A comprehensive analysis of genes selectively regulated by pH(i )provides insight on candidate targets that might mediate established effects of pH(i )on a number of normal and pathological cell functions

Homopolymer tract length dependent enrichments in functional regions of 27 eukaryotes and their novel dependence on the organism DNA (G+C)% composition

BACKGROUND: DNA homopolymer tracts, poly(dA).poly(dT) and poly(dG).poly(dC), are the simplest of simple sequence repeats. Homopolymer tracts have been systematically examined in the coding, intron and flanking regions of a limited number of eukaryotes. As the number of DNA sequences publicly available increases, the representation (over and under) of homopolymer tracts of different lengths in these regions of different genomes can be compared. RESULTS: We carried out a survey of the extent of homopolymer tract over-representation (enrichment) and over-proportional length distribution (above expected length) primarily in the single gene documents, but including some whole chromosomes of 27 eukaryotics across the (G+C)% composition range from 20 – 60%. A total of 5.2 × 10(7 )bases from 15,560 cleaned (redundancy removed) sequence documents were analyzed. Calculated frequencies of non-overlapping long homopolymer tracts were found over-represented in non-coding sequences of eukaryotes. Long poly(dA).poly(dT) tracts demonstrated an exponential increase with tract length compared to predicted frequencies. A novel negative slope was observed for all eukaryotes between their (G+C)% composition and the threshold length N where poly(dA).poly(dT) tracts exhibited over-representation and a corresponding positive slope was observed for poly(dG).poly(dC) tracts. Tract size thresholds where over-representation of tracts in different eukaryotes began to occur was between 4 – 11 bp depending upon the organism (G+C)% composition. The higher the GC%, the lower the threshold N value was for poly(dA).poly(dT) tracts, meaning that the over-representation happens at relatively lower tract length in more GC-rich surrounding sequence. We also observed a novel relationship between the highest over-representations, as well as lengths of homopolymer tracts in excess of their random occurrence expected maximum lengths. CONCLUSIONS: We discuss how our novel tract over-representation observations can be accounted for by a few models. A likely model for poly(dA).poly(dT) tract over-representation involves the known insertion into genomes of DNA synthesized from retroviral mRNAs containing 3' polyA tails. A proposed model that can account for a number of our observed results, concerns the origin of the isochore nature of eukaryotic genomes via a non-equilibrium GC% dependent mutation rate mechanism. Our data also suggest that tract lengthening via slip strand replication is not governed by a simple thermodynamic loop energy model