Evidence That Interspecies Polymorphism in the Human and Rat Cholecystokinin Receptor-2 Affects Structure of the Binding Site for the Endogenous Agonist Cholecystokinin

The cholecystokinin (CCK) receptor-2 exerts very important central and peripheral functions by binding the neuropeptides cholecystokinin or gastrin. Because this receptor is a potential therapeutic target, great interest has been devoted to the identification of efficient antagonists. However, interspecies genetic polymorphism that does not alter cholecystokinin-induced signaling was shown to markedly affect activity of synthetic ligands. In this context, precise structural study of the agonist binding site on the human cholecystokinin receptor-2 is a prerequisite to elucidating the molecular basis for its activation and to optimizing properties of synthetic ligands. In this study, using site-directed mutagenesis and molecular modeling, we delineated the binding site for CCK on the human cholecystokinin receptor-2 by mutating amino acids corresponding to that of the rat homolog. By doing so, we demonstrated that, although resembling that of rat homolog, the human cholecystokinin receptor-2 binding site also displays important distinct structural features that were demonstrated by susceptibility to several point mutations (F120A, Y189A, H207A). Furthermore, docking of CCK in the human and rat cholecystokinin receptor-2, followed by dynamic simulations, allowed us to propose a plausible structural explanation of the experimentally observed difference between rat and human cholecystokinin-2 receptors.

Rational Design of Acridine-Based Ligands with Selectivity for Human Telomeric Quadruplexes

Structure-based modeling methods have been used to design a series of disubstituted triazole-linked acridine compounds with selectivity for human telomeric quadruplex DNAs. A focused library of these compounds was prepared using click chemistry and the selectivity concept was validated against two promoter quadruplexes from the c-kit gene with known molecular structures, as well as with duplex DNA using a FRET-based melting method. Lead compounds were found to have reduced effects on the thermal stability of the c-kit quadruplexes and duplex DNA structures. These effects were further explored with a series of competition experiments, which confirmed that binding to duplex DNA is very low even at high duplex:telomeric quadruplex ratios. Selectivity to the c-kit quadruplexes is more complex, with some evidence of their stabilization at increasing excess over human telomeric quadruplex DNA. Selectivity is a result of the dimensions of the triazole-acridine compounds; and in particular the separation of the two alkyl-amino terminal groups. Both lead compounds also have selective inhibitory effects on the proliferation of cancer cell lines compared to a normal cell line, and one has been shown to inhibit the activity of the telomerase enzyme, which is selectively expressed in tumor cells, where it plays a role in maintaining telomere integrity and cellular immortalization.

Targeting Human Gastrointestinal Stromal Tumor Cells with a Quadruplex-Binding Small Molecule

Most of human gastrointestinal stromal tumors (GIST) are driven by activating mutations in the protooncogene KIT, a tyrosine kinase receptor. Clinical treatment with imatinib targets the kinase domain of KIT, but tumor regrowth occurs as a result of them development of resistant mutations in the kinase active site. An alternative small-molecule approach to GIST therapy is described, in which the KIT gene is directly targeted, and thus, kinase resistance may be circumvented. A naphthalene diimide derivative has been used to demonstrate the concept of dual quadruplex targeting. This compound strongly stabilizes both telomeric quadruplex DNA and quadruplex sites in the KIT promoter in vitro. It is shown here that the compound is a potent inducer of growth arrest in a patient-derived GIST cell line at a concentration (similar to 1 mu M) that also results in effective inhibition of telomerase activity and almost complete suppression of KIT mRNA and KIT protein expression. Molecular modeling studies with a telomeric quadruplex have been used to rationalize aspects of the experimental quadruplex melting data.

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers.

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150–200 nucleotides at the 3' end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.

P2Y(1) antagonists: Combining receptor-based modeling and QSAR for a quantitative prediction of the biological activity based on consensus scoring

P2Y(1) is an ADP-activated G protein-coupled receptor (GPCR). Its antagonists impede platelet aggregation in vivo and are potential antithrombotic agents. Combining ligand and structure-based modeling we generated a consensus model (LIST-CM) correlating antagonist structures with their potencies. We docked 45 antagonists into our rhodopsin-based human P2Y(1) homology model and calculated docking scores and free binding energies with the Linear Interaction Energy (LIE) method in continuum-solvent. The resulting alignment was also used to build QSAR based on CoMFA, CoMSIA, and molecular descriptors. To benefit from the strength of each technique and compensate for their limitations, we generated our LIST-CM with a PLS regression based on the predictions of each methodology. A test set featuring untested substituents was synthesized and assayed in inhibition of 2-MeSADP-stimulated PLC activity and in radioligand binding. LIST-CM outperformed internal and external predictivity of any individual model to predict accurately the potency of 75% of the test set.

Market Basket Survey Shows Elevated Levels of As in South Central U.S. Processed Rice Compared to California: Consequences for Human Dietary Exposure:consequences for human dietary exposure

We report the largest market basket survey of arsenic (As) in U.S. rice to date. Our findings show differences in transitional-metal levels between polished and unpolished rice and geographical variation in As and selenium (Se) between rice processed in California and the South Central U.S. The mean and median As grain levels for the South Central U.S. were 0.30 and 0.27 µg As g-1, respectively, for 107 samples. Levels for California were 41% lower than the South Central U.S., with a mean of 0.17 µg As g-1 and a median of 0.16 µg As g-1 for 27 samples. The mean and median Se grain levels for the South Central U.S. were 0.19 µg Se g-1. Californian rice levels were lower, averaging only 0.08 and 0.06 µg Se g-1 for mean and median values, respectively. The difference between the two regions was found to be significant for As and Se (General Linear Model (GLM):? As p < 0.001; Se p < 0.001). No statistically significant differences were observed in As or Se levels between polished and unpolished rice (GLM:? As p = 0.213; Se p = 0.113). No significant differences in grain levels of manganese (Mn), cobalt (Co), copper (Cu), or zinc (Zn) were observed between California and the South Central U.S. Modeling arsenic intake for the U.S. population based on this survey shows that for certain groups (namely Hispanics, Asians, sufferers of Celiac disease, and infants) dietary exposure to inorganic As from elevated levels in rice potentially exceeds the maximum intake of As from drinking water (based on consumption of 1 L of 0.01 mg L-1 In. As) and Californian state exposure limits. Further studies on the transformation of As in soil, grain As bioavailability in the human gastrointestinal tract, and grain elemental speciation trends are critical.

Genome-wide linkage analysis for human longevity: Genetics of healthy aging study

Clear evidence exists for heritability of humanlongevity, and much interest is focused on identifying genes associated with longer lives. To identify such longevity alleles, we performed the largest genome-wide linkage scan thus far reported. Linkage analyses included 2118nonagenarian Caucasian sibling pairs that have been enrolled in 15 study centers of 11 European countries as part of the Genetics of Healthy Aging (GEHA) project. In the joint linkage analyses, we observed four regions that show linkage with longevity; chromosome 14q11.2 (LOD = 3.47), chromosome 17q12-q22 (LOD = 2.95), chromosome 19p13.3-p13.11 (LOD = 3.76), and chromosome 19q13.11-q13.32 (LOD = 3.57). To fine map these regions linked to longevity, we performed association analysis using GWAS data in a subgroup of 1228 unrelated nonagenarian and 1907 geographically matched controls. Using a fixed-effect meta-analysis approach, rs4420638 at the TOMM40/ APOE/APOC1 gene locus showed significant association with longevity (P-value = 9.6 × 10). By combined modeling of linkage and association, we showed that association of longevity with APOEe4 and APOEe2 alleles explain the linkage at 19q13.11-q13.32 with P-value = 0.02 and P-value = 1.0 × 10, respectively. In the largest linkage scan thus far performed for human familial longevity, we confirm that the APOE locus is a longevity gene and that additional longevity loci may be identified at 14q11.2, 17q12-q22, and 19p13.3-p13.11. As the latter linkage results are not explained by common variants, we suggest that rare variants play an important role in human familial longevity.

Fluid Shear Induces Conformation Change in Human Blood Protein von Willebrand Factor in Solution

Many of the physiological functions of von Willebrand Factor (VWF), including its binding interaction with blood platelets, are regulated by the magnitude of applied fluid/hydrodynamic stress. We applied two complementary strategies to study the effect of fluid forces on the solution structure of VWF. First, small-angle neutron scattering was used to measure protein conformation changes in response to laminar shear rates (G) up to 3000/s. Here, purified VWF was sheared in a quartz Couette cell and protein conformation was measured in real time over length scales from 2-140 nm. Second, changes in VWF structure up to 9600/s were quantified by measuring the binding of a fluorescent probe 1,1'-bis(anilino)-4-,4'-bis(naphtalene)-8,8'-disulfonate (bis-ANS) to hydrophobic pockets exposed in the sheared protein. Small angle neutron scattering studies, coupled with quantitative modeling, showed that VWF undergoes structural changes at G < 3000/s. These changes were most prominent at length scales <10 nm (scattering vector (q) range >0.6/nm). A mathematical model attributes these changes to the rearrangement of domain level features within the globular section of the protein. Studies with bis-ANS demonstrated marked increase in bis-ANS binding at G > 2300/s. Together, the data suggest that local rearrangements at the domain level may precede changes at larger-length scales that accompany exposure of protein hydrophobic pockets. Changes in VWF conformation reported here likely regulate protein function in response to fluid shear.

Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment:Meeting new challenges

Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

Status of GPCR modeling and docking as reflected by community-wide GPCR Dock 2010 assessment

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe.

Automatic extraction of Reduced order models from CFD simulations for building energy modeling

Thermal comfort is defined as “that condition of mind which expresses satisfaction with the thermal environment’ [1] [2]. Field studies have been completed in order to establish the governing conditions for thermal comfort [3]. These studies showed that the internal climate of a room was the strongest factor in establishing thermal comfort. Direct manipulation of the internal climate is necessary to retain an acceptable level of thermal comfort. In order for Building Energy Management Systems (BEMS) strategies to be efficiently utilised it is necessary to have the ability to predict the effect that activating a heating/cooling source (radiators, windows and doors) will have on the room. The numerical modelling of the domain can be challenging due to necessity to capture temperature stratification and/or different heat sources (radiators, computers and human beings). Computational Fluid Dynamic (CFD) models are usually utilised for this function because they provide the level of details required. Although they provide the necessary level of accuracy these models tend to be highly computationally expensive especially when transient behaviour needs to be analysed. Consequently they cannot be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. The test case used in this work is a room of the Environmental Research Institute (ERI) Building at the University College Cork (UCC). ROMs have shown that they are sufficiently accurate with a total error of less than 1% and successfully retain a satisfactory representation of the phenomena modelled. The number of zones in a ROM defines the size and complexity of that ROM. It has been observed that ROMs with a higher number of zones produce more accurate results. As each ROM has a time to solution of less than 20 seconds they can be integrated into the BEMS of a building which opens the potential to real time physics based building energy modelling.

Modeling Klebsiella pneumoniae pathogenesis by infection of the wax moth Galleria mellonella

The implementation of infection models that approximate human disease is essential for understanding pathogenesis at the molecular level and for testing new therapies before they are entered into clinical stages. Insects are increasingly being used as surrogate hosts because they share, with mammals, essential aspects of the innate immune response to infections. We examined whether the larva of the wax moth Galleria mellonella could be used as a host model to conceptually approximate Klebsiella pneumoniae-triggered pneumonia. We report that the G. mellonella model is capable of distinguishing between pathogenic and nonpathogenic Klebsiella strains. Moreover, K. pneumoniae infection of G. mellonella models some of the known features of Klebsiella-induced pneumonia, i.e., cell death associated with bacterial replication, avoidance of phagocytosis by phagocytes, and the attenuation of host defense responses, chiefly the production of antimicrobial factors. Similar to the case for the mouse pneumonia model, activation of innate responses improved G. mellonella survival against subsequent Klebsiella challenge. Virulence factors necessary in the mouse pneumonia model were also implicated in the Galleria model. We found that mutants lacking capsule polysaccharide, lipid A decorations, or the outer membrane proteins OmpA and OmpK36 were attenuated in Galleria. All mutants activated G. mellonella defensive responses. The Galleria model also allowed us to monitor Klebsiella gene expression. The expression levels of cps and the loci implicated in lipid A remodeling peaked during the first hours postinfection, in a PhoPQ- and PmrAB-governed process. Taken together, these results support the utility of G. mellonella as a surrogate host for assessing infections with K. pneumoniae.

Activation of human TLR4/MD-2 by hypoacylated lipopolysaccharide from a clinical isolate of Burkholderia cenocepacia

Lung infection by Burkholderia species, in particular B. cenocepacia, accelerates tissue damage and increase post-lung transplant mortality in cystic fibrosis patients. Host- microbes interplay largely depends on interactions between pathogen specific molecules and innate immune receptors such as the Toll-like receptor 4 (TLR4), which recognizes the lipid A moiety of the bacterial lipopolysaccharide (LPS). The human TLR4/MD-2 LPS receptor complex is strongly activated by hexa-acylated lipid A and poorly activated by underacylated lipid A. Here, we report that B. cenocepacia LPS strongly activates human TLR4/MD-2 despite its lipid A having only five acyl chains. Further, we show that aminoarabinose residues in lipid A contribute to TLR4-lipid A interactions, and experiments in a mouse model of LPS-induced endotoxic shock confirmed the pro- inflammatory potential of B. cenocepacia penta-acylated lipid A. Molecular modeling, combined with mutagenesis of TLR4-MD2 interactive surfaces, suggests that longer acyl chains and the aminoarabinose residues in the B. cenocepacia lipid A allow exposure of the fifth acyl chain on the surface of MD-2 enabling interactions with TLR4 and its dimerization. Our results provide a molecular model for activation of the human TLR4/MD- 2 complex by penta-acylated lipid A, explaining the ability of hypoacylated B. cenocepacia LPS to promote pro- inflammatory responses associated to the severe pathogenicity of this opportunistic bacterium.

Microphysiological modeling of the reproductive tract: A fertile endeavor

Preclinical toxicity testing in animal models is a cornerstone of the drug development process, yet it is often unable to predict adverse effects and tolerability issues in human subjects. Species-specific responses to investigational drugs have led researchers to utilize human tissues and cells to better estimate human toxicity. Unfortunately, human cell-derived models are imperfect because toxicity is assessed in isolation, removed from the normal physiologic microenvironment. Microphysiological modeling often referred to as 'organ-on-a-chip' or 'human-on-a-chip' places human tissue into a microfluidic system that mimics the complexity of human in vivo physiology, thereby allowing for toxicity testing on several cell types, tissues, and organs within a more biologically relevant environment. Here we describe important concepts when developing a repro-on-a-chip model. The development of female and male reproductive microfluidic systems is critical to sex-based in vitro toxicity and drug testing. This review addresses the biological and physiological aspects of the male and female reproductive systems in vivo and what should be considered when designing a microphysiological human-on-a-chip model. Additionally, interactions between the reproductive tract and other systems are explored, focusing on the impact of factors and hormones produced by the reproductive tract and disease pathophysiology.