Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

Y chromosome polymorphism is a strong determinant of male fitness in Drosophila melanogaster

In many species, the Y (or W) chromosome carries relatively few functional genes. This observation motivates the null hypothesis that the Y will be a minor contributor to genetic variation for fitness. Previous data and theory supported the null hypothesis, but evidence presented here shows that the Y of Drosophila melanogaster is a major determinant of a male's total fitness, with standing genetic variation estimated to be 68% of that of an entire X/autosome genomic haplotype. Most Y-linked genes are expressed during spermatogenesis, and correspondingly, we found that the Y influences fitness primarily through its effect on a male's reproductive success (sperm competition and/or mating success) rather than his egg-to-adult viability. But the fitness of a Y highly depended on the genetic makeup of its bearer, reverting from high to low in different genetic backgrounds. This pattern leads to large epistatic (inconsistent among backgrounds) but no additive (consistent among backgrounds) Y-linked genetic variance for fitness. On a microevolutionary scale, the observed large epistatic variation on the Y substantially reduces heritable variation for fitness among males, and on a macroevolutionary scale, the Y produces strong selection for genomic rearrangements that move interacting genes onto the nonrecombining region of the Y.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Symmetry and the energy landscapes of biomolecules

The role of symmetry in the folding of proteins is discussed using energy landscape theory. An analytical argument shows it is much easier to find sequences with funneled energy landscape capable of fast folding if the structure is symmetric. The analogy with phase transitions of small clusters with magic numbers is discussed.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

A chromosomal locus required for copper resistance, competitive fitness, and cytochrome c biogenesis in Pseudomonas fluorescens.

A chromosomal locus required for copper resistance and competitive fitness was cloned from a strain of Pseudomonas fluorescens isolated from copper-contaminated agricultural soil. Sequence analysis of this locus revealed six open reading frames with homology to genes involved in cytochrome c biogenesis in other bacteria, helC, cycJ, cycK, tipB, cycL, and cycH, with the closest similarity being to the aeg-46.5(yej) region of the Escherichia coli chromosome. The proposed functions of these genes in other bacteria include the binding, transport, and coupling of heme to apocytochrome c in the periplasm of these Gram-negative bacteria. Putative heme-binding motifs were present in the predicted products of cycK and cycL, and TipB contained a putative disulfide oxidoreductase active site proposed to maintain the heme-binding site of the apocytochrome in a reduced state for ligation of heme. Tn3-gus mutagenesis showed that expression of the genes was constitutive but enhanced by copper, and confirmed that the genes function both in copper resistance and production of active cytochrome c. However, two mutants in cycH were copper-sensitive and oxidase-positive, suggesting that the functions of these genes, rather than cytochrome c oxidase itself, were required for resistance to copper.

Mapping leaf surface landscapes.

Leaf surfaces provide the ecologically relevant landscapes to those organisms that encounter or colonize the leaf surface. Leaf surface topography directly affects microhabitat availability for colonizing microbes, microhabitat quality and acceptability for insects, and the efficacy of agricultural spray applications. Prior detailed mechanistic studies that examined particular fungi-plant and pollinator-plant interactions have demonstrated the importance of plant surface topography or roughness in determining the outcome of the interactions. Until now, however, it has not been possible to measure accurately the topography--i.e., the three-dimensional structure--of such leaf surfaces or to record precise changes in patterns of leaf surface elevation over time. Using contact mode atomic force microscopy, we measured three-dimensional coordinates of upper leaf surfaces of Vaccinium macrocarpon (cranberry), a perennial plant, on leaves of two age classes. We then produced topographic maps of these leaf surfaces, which revealed striking differences between age classes of leaves: old leaves have much rougher surfaces than those of young leaves. Atomic force microscope measurements were analyzed by lag (1) autocorrelation estimates of leaf surfaces by age class. We suggest that the changes in topography result from removal of epicuticular lipids and that the changes in leaf surface topography influence phylloplane ecology. Visualizing and mapping leaf surfaces permit detailed investigations into leaf surface-mediated phenomena, improving our understanding of phylloplane interactions.

An RNA structure involved in feedback regulation of splicing and of translation is critical for biological fitness.

While studies of the regulation of gene expression have generally concerned qualitative changes in the selection or the level of expression of a gene, much of the regulation that occurs within a cell involves the continuous subtle optimization of the levels of proteins used in macromolecular complexes. An example is the biosynthesis of the ribosome, in which equimolar amounts of nearly 80 ribosomal proteins must be supplied by the cytoplasm to the nucleolus. We have found that the transcript of one of the ribosomal protein genes of Saccharomyces cerevisiae, RPL32, participates in such fine tuning. Sequences from exon I of the RPL32 transcript interact with nucleotides from the intron to form a structure that binds L32 to regulate splicing. In the spliced transcript, the same sequences interact with nucleotides from exon II to form a structure that binds L32 to regulate translation, thus providing two levels of autoregulation. We now show, by using a sensitive cocultivation assay, that these RNA structures and their interaction with L32 play a role in the fitness of the cell. The change of a single nucleotide within the 5' leader of the RPL32 transcript, which abolishes the site for L32 binding, leads to detectably slower growth and to eventual loss of the mutant strain from the culture. Experiments designed to assess independently the regulation of splicing and the regulation of translation are presented. These observations demonstrate that, in evolutionary terms, subtle regulatory compensations can be critical. The change in structure of an RNA, due to alteration of just one noncoding nucleotide, can spell the difference between biological success and failure.

Muller's ratchet decreases fitness of a DNA-based microbe.

Muller proposed that an asexual organism will inevitably accumulate deleterious mutations, resulting in an increase of the mutational load and an inexorable, ratchet-like, loss of the least mutated class [Muller, H.J. (1964) Mutat. Res. 1, 2-9]. The operation of Muller's ratchet on real populations has been experimentally demonstrated only in RNA viruses. However, these cases are exceptional in that the mutation rates of the RNA viruses are extremely high. We have examined whether Muller's ratchet operates in Salmonella typhimurium, a DNA-based organism with a more typical genomic mutation rate. Cells were grown asexually under conditions expected to result in high genetic drift, and the increase in mutational load was determined. S. typhimurium accumulated mutations under these conditions such that after 1700 generations, 1% of the 444 lineages tested had suffered an obvious loss of fitness, as determined by decreased growth rate. These results suggest that in the absence of sex and with high genetic drift, genetic mechanisms, such as back or compensatory mutations, cannot compensate for the accumulation of deleterious mutations. In addition, we measured the appearance of auxotrophs, which allowed us to calculate an average spontaneous mutation rate of approximately 0.3-1.5 x 10(-9) mutations per base pair per generation. This rate is measured for the largest genetic target studied so far, a collection of about 200 genes.

Exponential increases of RNA virus fitness during large population transmissions.

The great adaptability shown by RNA viruses is a consequence of their high mutation rates. Here we investigate the kinetics of virus fitness gains during repeated transfers of large virus populations in cell culture. Results always show that fitness increases exponentially. Low fitness clones exhibit regular increases observed as biphasic periods of exponential evolutionary improvement, while neutral clones show monophasic kinetics. These results are significant for RNA virus epidemiology, optimal handling of attenuated live virus vaccines, and routine laboratory procedures.

Organismal duplication, inclusive fitness theory, and altruism: understanding the evolution of endosperm and the angiosperm reproductive syndrome.

For almost a century, events relating to the evolutionary origin of endosperm, a unique embryo-nourishing tissue that is essential to the reproductive process in flowering plants, have remained a mystery. Integration of recent advances in phylogenetic reconstruction, comparative reproductive biology, and genetic theory can be used to elucidate the evolutionary events and forces associated with the establishment of endosperm. Endosperm is shown to be derived from one of two embryos formed during a rudimentary process of "double fertilization" that evolved in the ancestors of angiosperms. Acquisition of embryo-nourishing behavior (with accompanying loss of individual fitness) by this supernumerary fertilization product was dependent upon compensatory gains in the inclusive fitness of related embryos. The result of the loss of individual fitness by one of the two original products of double fertilization was the establishment of endosperm, a highly modified embryo/organism that reproduces cryptically through behavior that enhances the fitness of its associated embryo within a seed. Finally, although triploid endosperm remains a synapomorphy of angiosperms, inclusive fitness analysis demonstrates that the embryo-nourishing properties of endosperm initially evolved in a diploid condition.