Systematics within Gyps vultures: a clade at risk

Background  

Populations of the Oriental White-backed Vulture (Gyps bengalensis) have declined by over 95% within the past decade. This decline is largely due to incidental consumption of the non-steroidal anti-inflammatory veterinary pharmaceutical diclofenac, commonly used to treat domestic livestock. The conservation status of other Gyps vultures in southern Asia is also of immediate concern, given the lack of knowledge regarding status of their populations and the continuing existence of taxonomic uncertainties. In this study, we assess phylogenetic relationships for all recognized species and the majority of subspecies within the genus Gyps. The continuing veterinary use of diclofenac is an unknown but potential risk to related species with similar feeding habits to Gyps bengalensis. Therefore, an accurate assessment of the phylogenetic relationships among Gyps vultures should aid in their conservation by clarifying taxonomic uncertainties, and enabling inference of their respective relatedness to susceptible G. bengalensis.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Distinct patterns in homooligomer tract sequence context in prokaryotic and eukaryotic DNA

The distributions of the junction sequences of homooligomer tracts of various lengths have been examined in prokaryotic DNA sequences and compared with those of eukaryotes. The general trends in the nearest and next to nearest neighbors to the tracts are similar for both groups. In both prokaryotes and eukaryotes A/T runs are preferentially flanked on either the 5' or the 3' ends by A and/or T. G/C runs are preferentially flanked by G and/or C. There is discrimination against A/T runs flanked by G or C and G/C runs flanked by A or T. However, whereas the distribution of prokaryotic homooligomer tract junction sequences was quite homogeneous, large variations were observed in the 5-fold larger eukaryotic database, increasing in magnitude from tracts of length 2 to 3 to 4 base pairs long. Possible DNA conformational implications and in particular DNA curvature and packaging aspects of prokaryotes and eukaryotes are discussed.     

Patterns of divergence during evolution of alpha 1-proteinase inhibitors in mammals

alpha 1-Proteinase inhibitor (alpha 1-PI), a member of the serine proteinase inhibitor superfamily, has a primary role in controlling neutrophil elastase activity within the mammalian circulation. Several studies have indicated that the reactive center region of alpha 1-PI, the amino acid sequence of which is critical to recognition of and binding to target proteinases, is highly divergent within and among species. This appears to be a consequence of accelerated rates of evolution that may have been driven by positive Darwinian selection. In order to examine this and other features of alpha 1-PI evolution in more detail, we have isolated and sequenced cDNAs representing alpha 1- PI mRNAs of the mouse species Mus saxicola and Mus minutoides and have compared these with a number of other mammalian alpha 1-PI mRNAs. Relative to other mammalian mRNAs, the extent of nonsynonymous substitution is generally high throughout the alpha 1-PI mRNA molecule, indicating greater overall rates of amino acid substitution. Within and among mouse species, the 5'-half of the mRNA, but not the 3'-half, has been homogenized by concerted evolution. Finally, the reactive center is under diversifying or positive Darwinian selection in murid rodents (rats, mice) and guinea pigs yet is under purifying selection in primates and artiodactyls. The significance of these findings to alpha 1-PI function and the possible selective forces driving evolution of serpins in general are discussed.      

Mitochondrial DNA polymorphisms in subterranean mole-rats of the Spalax ehrenbergi superspecies in Israel, and its peripheral isolates

Patterns of mitochondrial DNA (mtDNA) variation were examined in 133 mole-rats constituting all four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of the Spalax ehrenbergi superspecies in Israel, as well as the peripheral isolates of 2n = 60. In the main range of the complex, a total of 28 mtDNA haplotypes were found in 64 mole-rats, with most haplotypes being unique to either a single chromosomal species or population. mtDNA divergence increased from low to high diploid number in a north-to-south direction in Israel. Overall levels of mtDNA diversity were unexpectedly the highest in the 2n = 60, the youngest species of the complex. The mtDNA haplotypes can be separated into two major groups, 2n = 52-54 and 2n = 58-60, and a phylogenetic analysis for each group revealed evidence of a few haplotypes not sorted by diploid number. The overall patterns of mtDNA divergence seen within and among the four chromosomal species are consistent with the parapatric mode of speciation as suggested from previous studies of allozyme and DNA hybridization. In a separate data set the patterns of mtDNA variation were examined across the main geographic range and across peripheral semi-isolates and isolates of the 2n = 60 chromosomal species. Fifteen haplotypes were found in 69 mole-rats. High levels of mtDNA diversity characterized the main range, semi-isolated, and even some desert isolated populations. The peripheral isolates contain much mtDNA diversity, including novel haplotypes.      

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

Conformations of folded proteins in restricted spaces

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.     

Solvent effect on binding thermodynamics of biopolymers

The indirect solvent-induced effect on the free energy of binding of biopolymers is examined within the framework of classical statistical mechanics. We focus specifically on the role of the solute-solvent hydrogen bonding. In particular, we have estimated the first order solvent effect on the indirect interaction between two biopolymers. We find that the solvent-induced interactions between two hydrophilic groups through water-bridged hydrogen bonds could significantly enhance the binding free energy. Some preliminary estimates indicate that this effect is significant and perhaps could be crucial in molecular recognition processes. Furthermore, we have calculated, from crystal structure data, the distance distribution between all the oxygens and nitrogens on the surface of some proteins that do not belong to the binding domain. In most cases we found an enhanced peak in the range of 4-5 A, which is where we expect to find strong solvent-induced interactions.