EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an 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 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, 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 could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Fitness, inclusive fitness, and optimization

Individual-as-maximizing agent analogies result in a simple understanding of the functioning of the biological world. Identifying the conditions under which individuals can be regarded as fitness maximizing agents is thus of considerable interest to biologists. Here, we compare different concepts of fitness maximization, and discuss within a single framework the relationship between Hamilton's (J Theor Biol 7: 1-16, 1964) model of social interactions, Grafen's (J Evol Biol 20: 1243-1254, 2007a) formal Darwinism project, and the idea of evolutionary stable strategies. We distinguish cases where phenotypic effects are additive separable or not, the latter not being covered by Grafen's analysis. In both cases it is possible to define a maximand, in the form of an objective function phi(z), whose argument is the phenotype of an individual and whose derivative is proportional to Hamilton's inclusive fitness effect. However, this maximand can be identified with the expression for fecundity or fitness only in the case of additive separable phenotypic effects, making individual-as-maximizing agent analogies unattractive (although formally correct) under general situations of social interactions. We also feel that there is an inconsistency in Grafen's characterization of the solution of his maximization program by use of inclusive fitness arguments. His results are in conflict with those on evolutionary stable strategies obtained by applying inclusive fitness theory, and can be repaired only by changing the definition of the problem.

Large-scale evolutionary trends of Acrochordiceratidae Arthaber, 1911 (Ammonoidea, Middle Triassic) and Cope's rule

Directed evolution of life through millions of years, such as increasing adult body size, is one of the most intriguing patterns displayed by fossil lineages. Processes and causes of such evolutionary trends are still poorly understood. Ammonoids (externally shelled marine cephalopods) are well known to have experienced repetitive morphological evolutionary trends of their adult size, shell geometry and ornamentation. This study analyses the evolutionary trends of the family Acrochordiceratidae Arthaber, 1911 from the Early to Middle Triassic (251228 Ma). Exceptionally large and bed-rock-controlled collections of this ammonoid family were obtained from strata of Anisian age (Middle Triassic) in north-west Nevada and north-east British Columbia. They enable quantitative and statistical analyses of its morphological evolutionary trends. This study demonstrates that the monophyletic clade Acrochordiceratidae underwent the classical evolute to involute evolutionary trend (i.e. increasing coiling of the shell), an increase in its shell adult size (conch diameter) and an increase in the indentation of its shell suture shape. These evolutionary trends are statistically robust and seem more or less gradual. Furthermore, they are nonrandom with the sustained shift in the mean, the minimum and the maximum of studied shell characters. These results can be classically interpreted as being constrained by the persistence and common selection pressure on this mostly anagenetic lineage characterized by relatively moderate evolutionary rates. Increasing involution of ammonites is traditionally interpreted by increasing adaptation mostly in terms of improved hydrodynamics. However, this trend in ammonoid geometry can also be explained as a case of Copes rule (increasing adult body size) instead of functional explanation of coiling, because both shell diameter and shell involution are two possible paths for ammonoids to accommodate size increase.

Evolutionary insights on C4 photosynthetic subtypes in grasses from genomics and phylogenetics

In plants, an oligogene family encodes NADP-malic enzymes (NADP-me), which are responsible for various functions and exhibit different kinetics and expression patterns. In particular, a chloroplast isoform of NADP-me plays a key role in one of the three biochemical subtypes of C4 photosynthesis, an adaptation to warm environments that evolved several times independently during angiosperm diversification. By combining genomic and phylogenetic approaches, this study aimed at identifying the molecular mechanisms linked to the recurrent evolutions of C4-specific NADP-me in grasses (Poaceae). Genes encoding NADP-me (nadpme) were retrieved from genomes of model grasses and isolated from a large sample of C3 and C4 grasses. Genomic and phylogenetic analyses showed that 1) the grass nadpme gene family is composed of four main lineages, one of which is expressed in plastids (nadpme-IV), 2) C4-specific NADP-me evolved at least five times independently from nadpme-IV, and 3) some codons driven by positive selection underwent parallel changes during the multiple C4 origins. The C4 NADP-me being expressed in chloroplasts probably constrained its recurrent evolutions from the only plastid nadpme lineage and this common starting point limited the number of evolutionary paths toward a C4 optimized enzyme, resulting in genetic convergence. In light of the history of nadpme genes, an evolutionary scenario of the C4 phenotype using NADP-me is discussed.

Correcting for the bias due to expression specificity improves the estimation of constrained evolution of expression between mouse and human.

MOTIVATION: Comparative analyses of gene expression data from different species have become an important component of the study of molecular evolution. Thus methods are needed to estimate evolutionary distances between expression profiles, as well as a neutral reference to estimate selective pressure. Divergence between expression profiles of homologous genes is often calculated with Pearson's or Euclidean distance. Neutral divergence is usually inferred from randomized data. Despite being widely used, neither of these two steps has been well studied. Here, we analyze these methods formally and on real data, highlight their limitations and propose improvements. RESULTS: It has been demonstrated that Pearson's distance, in contrast to Euclidean distance, leads to underestimation of the expression similarity between homologous genes with a conserved uniform pattern of expression. Here, we first extend this study to genes with conserved, but specific pattern of expression. Surprisingly, we find that both Pearson's and Euclidean distances used as a measure of expression similarity between genes depend on the expression specificity of those genes. We also show that the Euclidean distance depends strongly on data normalization. Next, we show that the randomization procedure that is widely used to estimate the rate of neutral evolution is biased when broadly expressed genes are abundant in the data. To overcome this problem, we propose a novel randomization procedure that is unbiased with respect to expression profiles present in the datasets. Applying our method to the mouse and human gene expression data suggests significant gene expression conservation between these species. CONTACT: marc.robinson-rechavi@unil.ch; sven.bergmann@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis.

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Discordances between phylogenetic and morphological patterns in alpine leaf beetles attest to an intricate biogeographic history of lineages in postglacial Europe.

Pleistocene glacial and interglacial periods have moulded the evolutionary history of European cold-adapted organisms. The role of the different mountain massifs has, however, not been accurately investigated in the case of high-altitude insect species. Here, we focus on three closely related species of non-flying leaf beetles of the genus Oreina (Coleoptera, Chrysomelidae), which are often found in sympatry within the mountain ranges of Europe. After showing that the species concept as currently applied does not match barcoding results, we show, based on more than 700 sequences from one nuclear and three mitochondrial genes, the role of biogeography in shaping the phylogenetic hypothesis. Dating the phylogeny using an insect molecular clock, we show that the earliest lineages diverged more than 1 Mya and that the main shift in diversification rate occurred between 0.36 and 0.18 Mya. By using a probabilistic approach on the parsimony-based dispersal/vicariance framework (MP-DIVA) as well as a direct likelihood method of state change optimization, we show that the Alps acted as a cross-roads with multiple events of dispersal to and reinvasion from neighbouring mountains. However, the relative importance of vicariance vs. dispersal events on the process of rapid diversification remains difficult to evaluate because of a bias towards overestimation of vicariance in the DIVA algorithm. Parallels are drawn with recent studies of cold-adapted species, although our study reveals novel patterns in diversity and genetic links between European mountains, and highlights the importance of neglected regions, such as the Jura and the Balkanic range.

The origins of C4 grasslands: integrating evolutionary and ecosystem science.

The evolution of grasses using C4 photosynthesis and their sudden rise to ecological dominance 3 to 8 million years ago is among the most dramatic examples of biome assembly in the geological record. A growing body of work suggests that the patterns and drivers of C4 grassland expansion were considerably more complex than originally assumed. Previous research has benefited substantially from dialog between geologists and ecologists, but current research must now integrate fully with phylogenetics. A synthesis of grass evolutionary biology with grassland ecosystem science will further our knowledge of the evolution of traits that promote dominance in grassland systems and will provide a new context in which to evaluate the relative importance of C4 photosynthesis in transforming ecosystems across large regions of Earth.

Two's company, three's a crowd: experimental evaluation of the evolutionary maintenance of trioecy in Mercurialis annua (Euphorbiaceae).

Trioecy is an uncommon sexual system in which males, females, and hermaphrodites co-occur as three clearly different gender classes. The evolutionary stability of trioecy is unclear, but would depend on factors such as hermaphroditic sex allocation and rates of outcrossing vs. selfing. Here, trioecious populations of Mercurialis annua are described for the first time. We examined the frequencies of females, males and hermaphrodites across ten natural populations and evaluated the association between the frequency of females and plant densities. Previous studies have shown that selfing rates in this species are density-dependent and are reduced in the presence of males, which produce substantially more pollen than hermaphrodites. Accordingly, we examined the evolutionary stability of trioecy using an experiment in which we (a) indirectly manipulated selfing rates by altering plant densities and the frequency of males in a fully factorial manner across 20 experimental plots and (b) examined the effect of these manipulations on the frequency of the three sex phenotypes in the next generation of plants. In the parental generation, we measured the seed and pollen allocations of hermaphrodites and compared them with allocations by unisexual plants. In natural populations, females occurred at higher frequencies in denser patches, a finding consistent with our expectations. Under our experimental conditions, however, no combination of plant densities and male frequencies was associated with increased frequencies of females. Our results suggest that the factors that regulate female frequencies in trioecious populations of M. annua are independent of those regulating male frequencies (density), and that the stable co-existence of all three sex phenotypes within populations is unlikely.

Social evolution in the sweat bee Halictus scabiosae

The evolution of eusociality is one of the major evolutionary transitions of life on earth. For investigating the conditions and processes that are central to the origin of such integrated social organization, it is best to study organisms in which individuals have retained some flexibility in their reproductive strategies. Halictid bees are especially well suited as model organisms, because they show huge variation in social systems, both within and between species. In this thesis, I investigated female reproductive strategies in the primitively eusocial bee Halictus scabiosae, with a focus on the role of helpers, in order to get insight into the mechanisms governing the evolution and maintenance of eusociality. This species produces two broods per year. The females from the first brood can stay in the natal nest to help raise a second brood of males and gynes that become the next-generation foundresses in spring. We first compared the morphology of females from the two broods, as well as the nutrition they receive as larvae. Then we conducted a helper- removal experiment in the field to quantify the effects of the presence of helpers on colony survival and productivity. Finally, we reconstructed pedigree relationships of individuals using microsatellite markers in order to detect who reproduces in the nest and how much individuals drift between nests. We found that first brood females had a uniformly small size and low fat reserves, which may be caused by the restricted pollen and nectar provisions on which they develop. Colony survival and productivity was increased by the presence of a single helper, but the effect was small and mostly limited to small colonies. By inferring parentage within and across colonies, we could determine that females from the first brood rarely reproduce in their natal nests. However, foundresses are frequently replaced, and foundresses and females from the first brood occasionally move to and reproduce in foreign colonies. As a result, colonies often contain offspring from unrelated individuals, and the relatedness of females to the brood they rear is low. Overall, this thesis shows that the reproductive system of H. scabiosae is highly flexible. The production of helpers in the first brood is important for colony success and productivity, but there is a high colony failure rate and part of the first brood females drift and reproduce in foreign nests. Both foundresses and helpers appear to be constrained by harsh environmental conditions or social factors limiting reproduction and independent colony founding. - L'origine des insectes sociaux est un domaine fascinant pour la recherche. Pour comprendre les mécanismes et les conditions qui sont nécessaires pour l'évolution et le maintien de la vie en société, il est judicieux d'étudier des sociétés primitives d'insectes, où toutes les femelles ont conservé la capacité de se reproduire, même si leur rôle comportemental dans la colonie est d'aider sans se reproduire. Une des familles d'abeilles, les halictes, est idéale pour cette sorte de recherche, en raison de la grande variabilité dans leur comportement social. Dans cette thèse, j'ai étudié les stratégies reproductives des femelles de Halictus scabiosae pour mieux comprendre les mécanismes qui influencent l'évolution de la vie en société. Cette espèce produit deux cohortes de couvain par année. Les femelles du premier couvain restent souvent dans leur nid natal pour aider à élever le deuxième couvain, tandis que les femelles du deuxième couvain s'accouplent et hibernent pour devenir les nouvelles fondatrices au printemps suivant. Nous avons d'abord comparé la morphologie des femelles issues des deux couvains ainsi que leur nutrition au stade de larve. Puis, dans une expérience sur le terrain, nous avons quantifié l'apport d'une ouvrière pour la survie et la productivité de la colonie. Finalement, nous avons reconstruit des pedigrees en utilisant des marqueurs génétiques, pour savoir qui se reproduit dans la colonie et combien d'individus migrent entre colonies. Les résultats montrent que les femelles du premier couvain sont uniformément plus petites et plus maigres, ce qui indique que les fondatrices réduisent les provisions de nourriture pour leur premier couvain afin de les inciter à aider dans le nid au lieu de se reproduire indépendamment. Dans l'expérience sur le terrain, la survie et la productivité de la colonie augmentaient avec la présence d'une ouvrière additionnelle, mais l'effet était petit et limité aux petites colonies. Par la reconstruction de pedigrees, nous pouvions constater que les femelles du premier couvain pondent rarement dans leurs nids natals. Les fondatrices cependant sont souvent remplacées en cours de saison, et migrent fréquemment entre nids, tandis que les femelles du premier couvain pondent parfois des oeufs dans des nids étrangers. De ce fait, les colonies contiennent souvent des descendants d'individus étrangers, et la parenté génétique entre les femelles et le deuxième couvain est basse. Cette thèse démontre que le système reproductif de H. scabiosae est très flexible. La production d'ouvrières est importante pour la survie de la colonie et sa productivité, mais le taux d'échec est élevé et une partie des femelles du premier couvain migrent et pondent dans une colonie étrangère. Autant les fondatrices que les ouvrières semblent être contraintes par des conditions environnementales ou sociales qui limitent la reproduction et les nouvelles fondations de colonie. - Die Entstehung von sozialen Lebensformen ist eines der wichtigsten Entwicklungen in der Geschichte des Lebens. Um die Bedingungen oder Prozesse zu verstehen, welche bei der Entstehung und dem Erhalt von sozialen Merkmalen wichtig sind, sollte man Lebewesen untersuchen, welche je nach Umwelteinflüßen ihr soziales Verhalten flexibel ändern können. Furchenbienen (Halictidae) gehören dazu. Diese weisen nämlich ein breites Spektrum verschiedener sozialer Organisationsformen auf, oftmals sogar innerhalb der einzelnen Arten. In meiner Doktorarbeit befasste ich mich mit den Fortpflanzungsstrategien der Weibchen der Skabiosen-Furchenbiene Halictus scabiosae. Diese Art produziert zwei Brüten pro Jahr. Die Weibchen der ersten Brut bleiben dabei meist als Arbeiterinnen in ihrem Geburtsnest, wohingegen die Weibchen der zweiten Brut nach der Paarung überwintern, um im nächsten Frühling neue Kolonien zu gründen. In einem ersten Schritt verglichen wir die beiden Brüten bezüglich der Grösse und der Fettreserven der Weibchen sowie der Pollen-Nektar-Vorräte für die Larven. Dann bestimmten wir in einem Feldexperiment, wieviel eine zusätzliche Arbeiterin zum Überleben und zur Produktiviät der Kolonie beiträgt. Schliesslich ermittelten wir durch genetische Tests die Verwandtschaftsbeziehungen zwischen den Bienen, um herauszufinden, wer in den Kolonien tatsächlich die Eier legt und ob und wieviel die Bienen zwischen verschiedenen Nestern wandern. Wir stellten fest, dass die Weibchen von der ersten Brut einheitlich kleiner sind und weniger Fettreserven besitzen. Das weist daraufhin, dass die Nestgründerin die erste Brut unterernährt, um die Wahrscheinlichkeit zu erhöhen, dass diese Weibchen als Arbeiterinnen im Nest bleiben anstatt sich unabhängig fortzupflanzen. Schon eine einzelne zusätzliche Arbeiterin verbesserte die Überlebenschancen und Produktivität der Kolonie, der Effekt war allerdings klein und auf kleine Kolonien beschränkt. Die Verwandtschaftsanalysen zeigten, dass die Arbeiterinnen nur sehr selten ein Ei in ihr Geburtsnest legen. Erstaunlicherweise wanderten die Nestgründerinnen oft zwischen verschiedenen Nestern. Einige Weibchen der ersten Brut wanderten auch in ein fremdes Nest und produzierten dort Nachkommen. Diese Doktorarbeit zeigt, dass die Fortpflanzungsstrategien der Skabiosen-Furchenbiene tatsächlich sehr flexibel sind. Die Anwesenheit von Arbeiterinnen ist wichtig für das Überleben und die Produktivität der Kolonie. Die Misserfolgsraten bleiben jedoch hoch, und ein Teil der Weibchen der ersten Brut pflanzt sich in fremden Nestern fort. Sowohl die Nestgründerinnen als auch die Weibchen der ersten Brut scheinen durch Umweltsbedingungen oder durch soziale Faktoren in der Wahl ihrer Fortpflanzungs¬strategie eingeschränkt zu sein.

Evolutionary genomics and HIV restriction factors.

PURPOSE OF REVIEW: To provide updated insights into innate antiviral immunity and highlight prototypical evolutionary features of well characterized HIV restriction factors. RECENT FINDINGS: Recently, a new HIV restriction factor, Myxovirus resistance 2, has been discovered and the region/residue responsible for its activity identified using an evolutionary approach. Furthermore, IFI16, an innate immunity protein known to sense several viruses, has been shown to contribute to the defense to HIV-1 by causing cell death upon sensing HIV-1 DNA. SUMMARY: Restriction factors against HIV show characteristic signatures of positive selection. Different patterns of accelerated sequence evolution can distinguish antiviral strategies--offense or defence--as well as the level of specificity of the antiviral properties. Sequence analysis of primate orthologs of restriction factors serves to localize functional domains and sites responsible for antiviral action. We use recent discoveries to illustrate how evolutionary genomic analyses help identify new antiviral genes and their mechanisms of action.