Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Heparin/heparan sulphate-based drugs

Glycosaminoglycans (GAGs) are an untapped source of novel chemical entities and, therefore, offer exciting new opportunities for the development of novel drug molecules because of their unique physical and biological properties. Advances in the functional understanding of GAG–protein interactions are enabling the development of GAG mimetics for use as anti-angiogenic, anti-metastatic, anti-inflammatory, anticoagulant and anti-thrombotic agents. Many anti-thrombotic molecules, such as Fondaparinux and Idraparinux, have been successful in clinical trials, and a new generation of heparin mimetic oligosaccharides and small molecules are currently in different stages of clinical development. In particular, the recent increased activity in the development of new mimetics by altering the composition of sulphated GAGs is very encouraging. This article reviews structurally defined heparin-mimetic oligosaccharides and small molecules currently in development or clinical trials...

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

The Transcriptional Regulation of Retrotransposon BARE

The purpose of this research project was to understand the steps of the retrotransposon BARE (BArley REtrotransposon) life cycle, from regulation of transcription to Virus-Like Particle (VLP) formation and ultimate integration back into the genome. Our study concentrates mainly on BARE1 transcriptional regulation because transcription is the crucial first step in the retrotransposon life cycle. The BARE element is a Class I LTR (Long Terminal Repeat) retrotransposon belonging to the Copia superfamily and was originally isolated in our research group. The LTR retrotransposons are transcribed from promoters in the LTRs and encode proteins for packaging of their transcripts, the reverse transcription of the transcripts into cDNA, and integration of the cDNA back into the genome. BARE1 is translated as a single polyprotein and cleaved into the capsid protein (GAG), integrase (IN), and reverse transcriptase-RNaseH (RT-RH) by the integral aspartic proteinase (AP). The BARE retrotransposon family comprises more than 104 copies in the barley (Hordeum vulgare) genome. The element is bound by long terminal repeats (LTRs, 1829 bp) containing promoters required for replication, signals for RNA processing, and motifs necessary for the integration of the cDNA. Members of the BARE1 subfamily are transcribed, translated, and form virus-like particles. Several basic questions concerning transcription are explored in the thesis: BARE1 transcription control, promoter choice in different barley tissues, start and termination sites for BARE transcripts, and BARE1 transcript polyadenylation (I). Polyadenylation is an important step during mRNA maturation, and determines its stability and translatability among other characteristics. Our work has found a novel way used by BARE1 to make extra GAG protein, which is critical for VLP formation. The discovery that BARE1 uses one RNA population for protein synthesis and another RNA population for making cDNA has established the most important step of the BARE1 life cycle (III). The relationship between BARE1 and BARE2 has been investigated. Besides BARE, we have examined the retrotransposon Cassandra (II), which uses a very different transcriptional mechanism and a fully parasitic life cycle. In general, this work is focused on BARE1 promoter activity, transcriptional regulation including differential promoter usage and RNA pools, extra GAG protein production and VLP formation. The results of this study give new insights into transcription regulation of LTR retrotransposons.

Molecular Epidemiology of HIV-1 Subtypes in India: Origin and Evolutionary History of the Predominant Subtype C

Background: India has the third largest HIV-1 epidemic with 2.4 million infected individuals. Molecular epidemiological analysis has identified the predominance of HIV-1 subtype C (HIV-1C). However, the previous reports have been limited by sample size, and uneven geographical distribution. The introduction of HIV-1C in India remains uncertain due to this lack of structured studies. To fill the gap, we characterised the distribution pattern of HIV-1 subtypes in India based on data collection from nationwide clinical cohorts between 2007 and 2011. We also reconstructed the time to the most recent common ancestor (tMRCA) of the predominant HIV-1C strains. Methodology/Principal Findings: Blood samples were collected from 168 HIV-1 seropositive subjects from 7 different states. HIV-1 subtypes were determined using two or three genes, gag, pol, and env using several methods. Bayesian coalescent-based approach was used to reconstruct the time of introduction and population growth patterns of the Indian HIV-1C. For the first time, a high prevalence (10%) of unique recombinant forms (BC and A1C) was observed when two or three genes were used instead of one gene (p<0.01; p = 0.02, respectively). The tMRCA of Indian HIV-1C was estimated using the three viral genes, ranged from 1967 (gag) to 1974 (env). Pol-gene analysis was considered to provide the most reliable estimate 1971, (95% CI: 1965-1976)]. The population growth pattern revealed an initial slow growth phase in the mid-1970s, an exponential phase through the 1980s, and a stationary phase since the early 1990s. Conclusions/Significance: The Indian HIV-1C epidemic originated around 40 years ago from a single or few genetically related African lineages, and since then largely evolved independently. The effective population size in the country has been broadly stable since the 1990s. The evolving viral epidemic, as indicated by the increase of recombinant strains, warrants a need for continued molecular surveillance to guide efficient disease intervention strategies.

Initiation with Elongator tRNAs

In all domains of life, initiator tRNA functions exclusively at the first step of protein synthesis while elongator tRNAs extend the polypeptide chain. Unique features of initiator tRNA enable it to preferentially bind the ribosomal P site and initiate translation. Recently, we showed that the abundance of initiator tRNA also contributes to its specialized role. This motivates the question, can a cell also use elongator tRNA to initiate translation under certain conditions? To address this, we introduced non-AUG initiation codons CCC (Pro), GAG (Glu), GGU (Gly), UCU (Ser), UGU (Cys), ACG (Thr), AAU (Asn), and AGA (Arg) into the uracil DNA glycosylase gene (ung) used as a reporter gene. Enzyme assays from log-phase cells revealed initiation from non-AUG codons when intracellular initiator tRNA levels were reduced. The activity increased significantly in stationary phase. Further increases in initiation from non-AUG codons occurred in both growth phases upon introduction of plasmid-borne genes of cognate elongator tRNAs. Since purine-rich Shine-Dalgarno sequences occur frequently on mRNAs (in places other than the canonical AUG codon initiation contexts), initiation with elongator tRNAs from the alternate contexts may generate proteome diversity under stress without compromising genomic integrity. Thus, by changing the relative amounts of initiator and elongator tRNAs within the cell, we have blurred the distinction between the two classes of tRNAs thought to be frozen through years of evolution.

Electrostatic Interactions Between Glycosaminoglycan Molecules

The electrostatic interactions between nearest-neighbouring chondroitin sulfate glycosaminoglycan (CS-GAG) molecular chains are obtained on the bottle brush conformation of proteoglycan aggrecan based on an asymptotic solution of the Poisson-Boltzmann equation the CS-GAGs satisfy under the physiological conditions of articular cartilage. The present results show that the interactions are associated intimately with the minimum separation distance and mutual angle between the molecular chains themselves. Further analysis indicates that the electrostatic interactions are not only expressed to be purely exponential in separation distance and decrease with the increasing mutual angle but also dependent sensitively on the saline concentration in the electrolyte solution within the tissue, which is in agreement with the existed relevant conclusions.

A probability-based approach to setting annual catch levels.

The requirement of setting annual catch limits to prevent overfishing has been added to the Magnuson-Stevens Fishery Conservation and Management Reauthorization Act of 2006 (MSRA). Because this requirement is new, a body of applied scientific practice for deriving annual catch limits and accompanying targets does not yet exist. This article demonstrates an approach to setting levels of catch that is intended to keep the probability of future overfishing at a preset low level. The proposed framework is based on stochastic projection with uncertainty in population dynamics. The framework extends common projection methodology by including uncertainty in the limit reference point and in management implementation, and by making explicit the risk of overfishing that managers consider acceptable. The approach is illustrated with application to gag (Mycteroperca microlepis), a grouper that inhabits the waters off the southeastern United States. Although devised to satisfy new legislation of the MSRA, the framework has potential application to any fishery where the management goal is to limit the risk of overfishing by controlling catch.

Integration of technologies for understanding the functional relationship between reef habitat and fish growth and production

Functional linkage between reef habitat quality and fish growth and production has remained elusive. Most current research is focused on correlative relationships between a general habitat type and presence/absence of a species, an index of species abundance, or species diversity. Such descriptive information largely ignores how reef attributes regulate reef fish abundance (density-dependent habitat selection), trophic interactions, and physiological performance (growth and condition). To determine the functional relationship between habitat quality, fish abundance, trophic interactions, and physiological performance, we are using an experimental reef system in the northeastern Gulf of Mexico where we apply advanced sensor and biochemical technologies. Our study site controls for reef attributes (size, cavity space, and reef mosaics) and focuses on the processes that regulate gag grouper (Mycteroperca microlepis) abundance, behavior and performance (growth and condition), and the availability of their pelagic prey. We combine mobile and fixed-active (fisheries) acoustics, passive acoustics, video cameras, and advanced biochemical techniques. Fisheries acoustics quantifies the abundance of pelagic prey fishes associated with the reefs and their behavior. Passive acoustics and video allow direct observation of gag and prey fish behavior and the acoustic environment, and provide a direct visual for the interpretation of fixed fisheries acoustics measurements. New application of biochemical techniques, such as Electron Transport System (ETS) assay, allow the in situ measurement of metabolic expenditure of gag and relates this back to reef attributes, gag behavior, and prey fish availability. Here, we provide an overview of our integrated technological approach for understanding and quantifying the functional relationship between reef habitat quality and one element of production – gag grouper growth on shallow coastal reefs.

Excreção urinária de glicosaminoglicanos em crianças com bexiga neurogênica secundária a mielomeningocele

A excreção urinária de glicosaminoglicanos (GAG) está alterada em várias patologias do trato urinário; o padrão de excreção pode estar associado com o estado da doença. A excreção urinária de GAG em crianças com bexiga neurogênica (BN) secundária a mielomeningocele (MMC) pode também estar alterada, mas até a presente data não há detalhamento epidemiológico dos pacientes e não se correlacionou o padrão de excreção com grau de disfunção vesical. Analisamos a excreção urinária de um grupo bem definido de crianças com MMC e correlacionamos os resultados com escore cistométrico. As amostras de urina de 17 pacientes com MMC, 10 meninos e 7 meninas (média de idade DP de 4,6 2,9 anos) foram obtidas durante o exame cistométrico. As amostras do grupo controle foram obtidas de 18 crianças normais, 13 meninos e 5 meninas (6,9 2,2 anos). Todas as crianças não estavam com infecção urinária, tinham função renal normal e não estavam sob tratamento farmacológico. A quantificação do GAG urinário total foi expressa em μg de ácido hexurônico / mg de creatinina e a proporção dos diferentes tipos de GAGs sulfatados foi obtida por eletroforese em gel de agarose. A avaliação cistométrica foi realizada utilizando aparelho de urodinâmica Dynapack modelo MPX816 (Dynamed, São Paulo, Brasil), a partir da qual o escore cistométrico foi calculado de acordo com procedimento recente publicado. [14]. Não observamos diferença significativa na excreção urinária de GAG total entre meninos e meninas tanto no grupo com MMC ( 0,913 0,528 vs 0,867 0,434, p>0,05) como no grupo controle (0,546 0,240 vs 0,699 0,296, p>0,05). Os resultados mostraram também que a excreção de GAG urinário não se correlacionou com a idade tanto no grupo com MMC ( r = -0,28, p> 0,05) como no grupo controle (r = -0,40, p> 0,05). Entretanto, a comparação dos dois grupos mostrou que o grupo com MMC excretava 52% a mais de GAG total que o grupo controle (0,894 0,477 vs 0,588 0,257, p <0,04). Nesses pacientes a excreção de GAG total não se correlacionou com a complacência vesical isoladamente (r = -0,18, p> 0,05) mas foi significativa e negativamente correlacionada ao escore cistométrico (r= -0,56, p<0,05). Em média, os pacientes com piores escores (<9) excretaram 81% a mais de GAG que os pacientes com melhor escore (>9) (1,157 0,467 vs 0,639 0,133, p<0,04). O sulfato de condroitin foi o GAG sulfatado predominante nos grupos neurogênico e controles (92,5 7,6% vs 96,4 4,8%, respectivamente, p> 0,05), enquanto o sulfato do heparan estava presente em quantidades marcadamente menores; o dermatam sulfato não foi detectado. A excreção urinária de GAG em pacientes com MMC é significativamente maior que a excreção das crianças normais e os altos valores encontrados estão correlacionados a um maior compromentimento da função vesical. Evidências em modelos animais com MMC induzida sugerem que alterações no detrusor estão associadas a um elevado turnover da matriz extra celular (MEC) vesical, o que pode explicar a elevada excreção de GAG nos pacientes com MMC. Além disso, esses resultados indicam que a excreção urinária de GAG pode ser usada como fator adjuvante para a caracterização da disfunção vesical em pacientes com MMC.

Composição da matriz extracelular do ligamento cardinal de mulheres na pós-menopausa com e sem prolapso uterino

O prolapso uterino tem sua incidência aumentada na pós-menopausa. O objetivo deste estudo é identificar as alterações na matriz extracelular do ligamento cardinal associadas à menopausa e ao prolapso uterino. Ligamento cardinal de três diferentes grupos de mulheres, pré-menopausa, prolapso uterino e pós-menopausa, foram identificados e biopsiados durante 57 histerectomias abdominais ou vaginais. As amostras foram processadas por métodos bioquímicos para caracterização e quantificação de glicosaminoglicanos sulfatados e colágeno. As concentrações relativas de glicosaminoglicanos foram obtidas por eletroforese. Procedimentos histológicos foram feitos para identificar fibras elásticas (Weigert), distribuição de colágeno (Picro Sirius) e decorin (imunohistoquímica). Nossos resultados mostraram aumento na concentração de GAG de 72,2%, redução na concentração de colágeno de 37% e diminuição de 22% de fibras elásticas no grupo de prolapso uterino quando comparado ao grupo da pós-menopausa (p<0,05, p<0,04 e p<0,05 respectivamente). As concentrações relativas de glicosaminoglicanos sulfatados para condroitin sulfato, heparan sulfato e dermatan sulfato não mostraram diferenças entre os três grupos. A organização do colágeno foi similar entre os três grupos e a marcação do decorin pareceu estar diminuída no grupo de prolapso uterino. Nossos resultados indicam alterações no metabolismo do tecido conjuntivo. O ligamento cardinal da mulher na pós-menopausa possui uma matriz extracelular mais densa. Esta alteração não ocorre na mulher com prolapso uterino.

Análise da composição da matriz extracelular das camadas urotelial e muscular da bexiga de mulheres em diferentes idades

A complacência da bexiga depende de músculos lisos, fibras colágenas, fibras elásiticas e suas relações. O objetivo deste trabalho é determinar a composição da matriz extracelular em amostras de bexigas normais através de análise bioquímica de colágeno e glicosaminoglicanos em amostras obtidas de mulheres em diferentes grupos de idade, analisando separadamente as camadas urotelial e muscular. Avaliamos 17 amostras de bexiga divididas em três grupos: infância (N=5), menacme (N=6) e pós-menopausa (N=6). As bexigas foram analisadas para concentração de GAG total e colágeno e para análise qualitativa de GAG por eletroforese em gel de agarose. Na camada muscular, não houve diferença entre os grupos tanto para GAG quanto para colágeno. Na camada urotelial, a análise da concentração de colágeno não mostrou diferença entre os grupos, mas a concentração de GAG no grupo da pós-menopausa (0.21 0.12 μg de ácido hexurônico/mg de tecido seco) apresentou diferença em relação aos grupos do menacme (1.78 1.62 μg de ácido hexurônico/mg de tecido seco) e da infância ( 2.29 1.32 μg de ácido hexurônico/mg de tecido seco).Nosso trabalho concluiu que a concentração de GAG está substancialmente diminuída na camada urotelial da bexiga de mulheres na pós-menopausa.

Estudo proteômico da interação do Aspergillus fumigatus com células endoteliais da veia umbilical humana (HUVECs)

O Aspergillus fumigatus é o principal agente etiológico da aspergilose invasiva, uma infecção fúngica oportunista que acomete, principalmente, pacientes de Unidades Hematológicas, como aqueles com neutropenia profunda e prolongada. Após a filamentação este fungo angioinvasivo é capaz de ativar e causar danos em células endoteliais de veia umbilical humana (HUVEC) que passam a expressar um fenótipo pró-trombótico. A ativação destas células, dependente de contato célulacélula, é mediada por TNF-α e caracterizada pela expressão de moléculas próinflamatórias, como citocinas, quimiocinas e moléculas de adesão. Recentemente, nosso grupo comparou a ativação endotelial de HUVECs desafiadas com cepas selvagens e uma cepa mutante para o gene UGM1. Nestes experimentos a cepa mutante Δugm1, que apresenta um fenótipo de maior produção de galactosaminogalactana (GAG) na parede celular, mostrou um fenótipo hiperadesivo e uma capacidade maior de ativar células endoteliais. Entretanto, os receptores e as vias de sinalização envolvidos nesta ativação permanecem desconhecidos. Assim, o objetivo deste trabalho foi verificar as proteínas envolvidas nestes processos através do estudo das proteínas diferencialmente expressas nas HUVECs após a interação com A. fumigatus, usando a técnica proteômica 2D-DIGE. Brevemente, as HUVECs foram infectadas com tubos germinativos da cepa selvagem (AF293) e da cepa Δugm1 de A. fumigatus. Em seguida, as proteínas foram marcadas com diferentes fluorocromos e separadas por eletroforese bidimensional. A análise quantitativa foi realizada utilizando o software DeCyder. Foram identificadas por MS/MS cinco proteínas diferencialmente expressas, incluindo a galectina-1 e a anexina A2, ambas mais expressas após a interação, sendo a primeira ~25% mais expressa após a interação com a mutante Δugm1. Este trabalho propõe que a galectina-1 poderia ser o receptor endotelial para polímeros de galactose presentes na parede celular do A. fumigatus, e que a Anexina A2 poderia estar envolvida na sinalização intracelular em resposta a este patógeno. No entanto, experimentos complementares, em curso, são necessários para comprovar esta hipótese.

On-going generation of multiple forms of HIV-1 intersubtype recombinants in the Yunnan Province of China

Objectives: To investigate the molecular epidemiology of HIV in China's Yunnan Province, where the initial HIV-1 outbreak among injecting drug users (IDU) occurred in 1989, and to analyse the genesis and interrelationship of the epidemic with that in surrounding areas. Design: A molecular epidemiological investigation was conducted among IDU in three prefectures in Yunnan Province, including Wenshan (east), Honghe (southeast) and Dehong (west). Methods: Thirty-nine specimens were collected from consenting IDU in 2000-2001. The nucleotide sequences of 2.6 kb gag-RT and 340 base pair (bp) env (C2/V3) regions were determined. Phylogenetic tree and recombination breakpoint analyses were performed. Results: The circulating recombinant form (CRF), CRF08_BC, predominated in east Yunnan near Guangxi Province (89% in Wenshan and 81% in Honghe), whereas it was not detected in Dehong(0/14) in the west. In contrast, 71% (10/14) of the Dehong isolates were unique recombinant forms (URF), mostly between subtypes B' (Thailand variant of subtype B) and C, with distinct profiles of recombination breakpoints. The subtype B' accounts for the remaining 29% (4/14) of Dehong isolates. Interestingly, two Honghe isolates (2/16) shared some of the precise B'/C recombination breakpoints with CRF07_BC. Conclusion: New recombinant strains are arising continually in west Yunnan near the Myanmar border. Some appeared to be secondary recombinants derived from CRF07_BC that had further recombined with other strains. The uneven distribution of subtypes, CRF and URF, suggests the presence of independent transmission networks and clusters among IDU in Yunnan. (C) 2002 Lippincott Williams Wilkins.

Identification and characterization of a new class of human immunodeficiency virus type 1 recombinants comprised of two circulating recombinant forms, CRF07_BC and CRF08_BC, in China

We identified a new class of human immunodeficiency virus type 1 (HIV-1) recombinants (00CN-HH069 and 00CN-HH086) in which further recombination occurred between two established circulating recombinant forms (CRFs). These two isolates were found among 57 HIV-1 samples from a cohort of injecting drug users in eastern Yunnan Province of China. Informative-site analysis in conjunction with bootscanning plots and exploratory tree analysis revealed that these two strains were closely related mosaics comprised of CRF07_BC and CRF08_BC, which are found in China. The genotype screening based on gag-reverse transcriptase sequences if 57 samples from eastern Yunnan identified 47 CRF08_BC specimens (82.5%), 5 CRF07_BC specimens (8.8%), and 3 additional specimens with the novel recombinant structure. These new "second-generation" recombinants thus constitute a substantial proportion (5 of 57; 8.8%) of HIV-1 strains in this population and may belong to a new but yet-undefined class of CRF. This might be the first example of CRFs recombining with each other, leading to the evolution of second-generation inter-CRF recombinants.