Ecology and life history affect different aspects of the population structure of 27 high-alpine plants.

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics.

Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima.

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.

EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Molecular characterization of the processing of arenavirus envelope glycoprotein precursors by subtilisin kexin isozyme-1/site-1 protease.

A crucial step in the life cycle of arenaviruses is the biosynthesis of the mature fusion-active viral envelope glycoprotein (GP) that is essential for virus-host cell attachment and entry. The maturation of the arenavirus GP precursor (GPC) critically depends on proteolytic processing by the cellular proprotein convertase (PC) subtilisin kexin isozyme-1 (SKI-1)/site-1 protease (S1P). Here we undertook a molecular characterization of the SKI-1/S1P processing of the GPCs of the prototypic arenavirus lymphocytic choriomeningitis virus (LCMV) and the pathogenic Lassa virus (LASV). Previous studies showed that the GPC of LASV undergoes processing in the endoplasmic reticulum (ER)/cis-Golgi compartment, whereas the LCMV GPC is cleaved in a late Golgi compartment. Herein we confirm these findings and provide evidence that the SKI-1/S1P recognition site RRLL, present in the SKI-1/S1P prodomain and LASV GPC, but not in the LCMV GPC, is crucial for the processing of the LASV GPC in the ER/cis-Golgi compartment. Our structure-function analysis revealed that the cleavage of arenavirus GPCs, but not cellular substrates, critically depends on the autoprocessing of SKI-1/S1P, suggesting differences in the processing of cellular and viral substrates. Deletion mutagenesis showed that the transmembrane and intracellular domains of SKI-1/S1P are dispensable for arenavirus GPC processing. The expression of a soluble form of the protease in SKI-I/S1P-deficient cells resulted in the efficient processing of arenavirus GPCs and rescued productive virus infection. However, exogenous soluble SKI-1/S1P was unable to process LCMV and LASV GPCs displayed at the surface of SKI-I/S1P-deficient cells, indicating that GPC processing occurs in an intracellular compartment. In sum, our study reveals important differences in the SKI-1/S1P processing of viral and cellular substrates.

Structure and function of RecA-DNA complexes.

While the E. coli RecA protein has been the most intensively studied enzyme of homologous recombination, the unusual RecA-DNA filament has stood alone until very recently. It now appears that this protein is part of a universal family that spans all of biology, and the filament that is formed by the protein on DNA is a universal structure. With RecA's role in recombination given new and greatly increased significance, we focus in this review on the energetics of the RecA-mediated strand exchange and the relation between the energetics and recombination spanning heterologous inserts.

Species-specific population structure in rock-specialized sympatric cichlid species in Lake Tanganyika, East Africa.

Species richness and geographical phenotypic variation in East African lacustrine cichlids are often correlated with ecological specializations and limited dispersal. This study compares mitochondrial and microsatellite genetic diversity and structure among three sympatric rock-dwelling cichlids of Lake Tanganyika, Eretmodus cyanostictus, Tropheus moorii, and Ophthalmotilapia ventralis. The species represent three endemic, phylogenetically distinct tribes (Eretmodini, Tropheini, and Ectodini), and display divergent ecomorphological and behavioral specialization. Sample locations span both continuous, rocky shoreline and a potential dispersal barrier in the form of a muddy bay. High genetic diversity and population differentiation were detected in T. moorii and E. cyanostictus, whereas much lower variation and structure were found in O. ventralis. In particular, while a 7-km-wide muddy bay curtails dispersal in all three species to a similar extent, gene flow along mostly continuous habitat appeared to be controlled by distance in E. cyanostictus, further restricted by site philopatry and/or minor habitat discontinuities in T. moorii, and unrestrained in O. ventralis. In contrast to the general pattern of high gene flow along continuous shorelines in rock-dwelling cichlids of Lake Malawi, our study identifies differences in population structure among stenotopic Lake Tanganyika species. The amount of genetic differentiation among populations was not related to the degree of geographical variation of body color, especially since more phenotypic variation is observed in O. ventralis than in the genetically highly structured E. cyanostictus.

Comparison of the organization and fine structure of a chicken and a Xenopus laevis vitellogenin gene.

The structural organization and the coding nucleotide sequence of the Xenopus laevis A2 and the chicken major vitellogenin genes have been compared. Both genes show the same exon-intron organization. However, the degree of homology between the nucleotide and derived amino acid sequences varies extensively along the genes. Several of the 35 exons are quite similar, and a unique cysteine motif in the lipovitellin II domain is conserved between the two genes. In contrast, one internal region is quite divergent. Part of this region encodes phosvitin, which appears to have evolved rapidly by both point mutations and duplications of serines or short other amino acid stretches. On the basis of these observations, we discuss the possible mechanism of evolution of phosvitin in vertebrates.

Structure of the extracellular domains of human and Xenopus Fn14: implications in the evolution of TWEAK and Fn14 interactions.

TWEAK (TNF homologue with weak apoptosis-inducing activity) and Fn14 (fibroblast growth factor-inducible protein 14) are members of the tumor necrosis factor (TNF) ligand and receptor super-families. Having observed that Xenopus Fn14 cross-reacts with human TWEAK, despite its relatively low sequence homology to human Fn14, we examined the conservation in tertiary fold and binding interfaces between the two species. Our results, combining NMR solution structure determination, binding assays, extensive site-directed mutagenesis and molecular modeling, reveal that, in addition to the known and previously characterized β-hairpin motif, the helix-loop-helix motif makes an essential contribution to the receptor/ligand binding interface. We further discuss the insight provided by the structural analyses regarding how the cysteine-rich domains of the TNF receptor super-family may have evolved over time. DATABASE: Structural data are available in the Protein Data Bank/BioMagResBank databases under the accession codes 2KMZ, 2KN0 and 2KN1 and 17237, 17247 and 17252. STRUCTURED DIGITAL ABSTRACT: TWEAK binds to hFn14 by surface plasmon resonance (View interaction) xeFn14 binds to TWEAK by enzyme linked immunosorbent assay (View interaction) TWEAK binds to xeFn14 by surface plasmon resonance (View interaction) hFn14 binds to TWEAK by enzyme linked immunosorbent assay (View interaction).

Comparative population genetic structure in a plant-pollinator/seed predator system.

Comparative analyses of spatial genetic structure of populations of plants and the insects they interact with provide an indication of how gene flow, natural selection and genetic drift may jointly influence the distribution of genetic variation and potential for local co-adaptation for interacting species. Here, we analysed the spatial scale of genetic structure within and among nine populations of an interacting species pair, the white campion Silene latifolia and the moth Hadena bicruris, along a latitudinal gradient across Northern/Central Europe. This dioecious, short-lived perennial plant inhabits patchy, often disturbed environments. The moth H. bicruris acts both as its pollinator and specialist seed predator that reproduces by laying eggs in S. latifolia flowers. We used nine microsatellite markers for S. latifolia and eight newly developed markers for H. bicruris. We found high levels of inbreeding in most populations of both plant and pollinator/seed predator. Among populations, significant genetic structure was observed for S. latifolia but not for its pollinator/seed predator, suggesting that despite migration among populations of H. bicruris, pollen is not, or only rarely, carried over between populations, thus maintaining genetic structure among plant populations. There was a weak positive correlation between genetic distances of S. latifolia and H. bicruris. These results indicate that while significant structure of S. latifolia populations creates the potential for differentiation at traits relevant for the interaction with the pollinator/seed predator, substantial gene flow in H. bicruris may counteract this process in at least some populations.

Morphological and molecular diversity among cassava genotypes

The objective of this work was to characterize morphologically and molecularly the genetic diversity of cassava accessions, collected from different regions in Brazil. A descriptive analysis was made for 12 morphological traits in 419 accessions. Data was transformed into binary data for cluster analysis and analysis of molecular variance. A higher proportion of white or cream (71%) root cortex color was found, while flesh colors were predominantly white (49%) and cream (42%). Four accession groups were classified by the cluster analysis, but they were not grouped according to their origin, which indicates that diversity is not structured in space. The variation was greater within regions (95.6%). Sixty genotypes were also evaluated using 14 polymorphic microsatellite markers. Molecular results corroborated the morphological ones, showing the same random distribution of genotypes, with no grouping according to origin. Diversity indices were high for each region, and a greater diversity was found within regions, with: a mean number of alleles per locus of 3.530; observed and expected heterozygosity of 0.499 and 0.642, respectively; and Shannon index of 1.03. The absence of spatial structure among cassava genotypes according to their origins shows the anthropic influence in the distribution and movement of germplasm, both within and among regions.

Structure of the glycosyl-phosphatidylinositol membrane anchor of the Leishmania major promastigote surface protease.

In common with many other plasma membrane glycoproteins of eukaryotic origin, the promastigote surface protease (PSP) of the protozoan parasite Leishmania contains a glycosyl-phosphatidylinositol (GPI) membrane anchor. The GPI anchor of Leishmania major PSP was purified following proteolysis of the PSP and analyzed by two-dimensional 1H-1H NMR, compositional and methylation linkage analyses, chemical and enzymatic modifications, and amino acid sequencing. From these results, the structure of the GPI-containing peptide was found to be Asp-Gly-Gly-Asn-ethanolamine-PO4-6Man alpha 1-6Man alpha 1-4GlcN alpha 1-6myo-inositol-1-PO4-(1-alkyl-2-acyl-glycerol). The glycan structure is identical to the conserved glycan core regions of the GPI anchor of Trypanosoma brucei variant surface glycoprotein and rat brain Thy-1 antigen, supporting the notion that this portion of GPIs are highly conserved. The phosphatidylinositol moiety of the PSP anchor is unusual, containing a fully saturated, unbranched 1-O-alkyl chain (mainly C24:0) and a mixture of fully saturated unbranched 2-O-acyl chains (C12:0, C14:0, C16:0, and C18:0). This lipid composition differs significantly from those of the GPIs of T. brucei variant surface glycoprotein and mammalian erythrocyte acetylcholinesterase but is similar to that of a family of glycosylated phosphoinositides found uniquely in Leishmania.

The use of external mesh reinforcement to reduce intimal hyperplasia and preserve the structure of human saphenous veins.

The saphenous vein is the conduit of choice in bypass graft procedures. Haemodynamic factors play a major role in the development of intimal hyperplasia (IH), and subsequent bypass failure. To evaluate the potential protective effect of external reinforcement on such a failure, we developed an ex vivo model for the perfusion of segments of human saphenous veins under arterial shear stress. In veins submitted to pulsatile high pressure (mean pressure at 100 mmHg) for 3 or 7 days, the use of an external macroporous polyester mesh 1) prevented the dilatation of the vessel, 2) decreased the development of IH, 3) reduced the apoptosis of smooth muscle cells, and the subsequent fibrosis of the media layer, 4) prevented the remodelling of extracellular matrix through the up-regulation of matrix metalloproteinases (MMP-2, MMP-9) and plasminogen activator type I. The data show that, in an experimental ex vivo setting, an external scaffold decreases IH and maintains the integrity of veins exposed to arterial pressure, via increase in shear stress and decrease wall tension, that likely contribute to trigger selective molecular and cellular changes.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

How a haemosporidian parasite of bats gets around: the genetic structure of a parasite, vector and host compared.

Parasite population structure is often thought to be largely shaped by that of its host. In the case of a parasite with a complex life cycle, two host species, each with their own patterns of demography and migration, spread the parasite. However, the population structure of the parasite is predicted to resemble only that of the most vagile host species. In this study, we tested this prediction in the context of a vector-transmitted parasite. We sampled the haemosporidian parasite Polychromophilus melanipherus across its European range, together with its bat fly vector Nycteribia schmidlii and its host, the bent-winged bat Miniopterus schreibersii. Based on microsatellite analyses, the wingless vector, and not the bat host, was identified as the least structured population and should therefore be considered the most vagile host. Genetic distance matrices were compared for all three species based on a mitochondrial DNA fragment. Both host and vector populations followed an isolation-by-distance pattern across the Mediterranean, but not the parasite. Mantel tests found no correlation between the parasite and either the host or vector populations. We therefore found no support for our hypothesis; the parasite population structure matched neither vector nor host. Instead, we propose a model where the parasite's gene flow is represented by the added effects of host and vector dispersal patterns.

Localization and fine structure of a vaccinia virus gene encoding an envelope antigen.

The major antigen on the envelope of extracellular vaccinia virus particles is a polypeptide with an apparent molecular weight of 37,000 (p37K; G. Hiller and K. Weber, J. Virol. 55:651-659, 1985). The gene encoding p37K was mapped in the vaccinia virus genome by hybrid selection of RNA followed by in vitro translation. p37K was then identified among the in vitro translation products by immunoprecipitation with a monoclonal antibody. The gene is located close to the right-hand end of the HindIII F fragment. The corresponding region of the DNA was sequenced, and an open reading frame encoding a polypeptide of 41,748 daltons was observed. The 5' end of the mRNA, as defined by nuclease S1 analysis, maps within only a few nucleotides of the translation initiation codon. Examination of the DNA sequence around the putative initiation site of transcription revealed a characteristic sequence, TAAATG, which includes the ATG translation initiation codon and which is conserved in all but one late gene so far analyzed. It is therefore likely that this sequence is an important regulatory signal for late gene expression in vaccinia virus.