Rate of loop formation in peptides: a simulation study

Experimental techniques with high temporal and spatial resolution extend our knowledge of how biological macromolecules self-organise and function. Here, we provide an illustration of the convergence between simulation and experiment made possible by techniques such as triplet-triplet energy transfer and fluorescence quenching with long-lifetime and fast-quenching fluorescent probes. These techniques have recently been used to determine the average time needed for two residues in a peptide or protein segment to form a contact. The timescale of this process is accessible to computer simulation, providing a microscopic interpretation of the data and yielding new insight into the disordered state of proteins. Conversely, such experimental data also provide a test of the validity of alternative choices for the molecular models used in simulations, indicating their possible deficiencies. We carried out simulations of peptides of various composition and length using several models. End-to-end contact formation rates and their dependence on peptide length agree with experimental estimates for some sequences and some force fields but not for others. The deviations are due to artefactual structuring of some peptides, which is not observed when an atomistic model for the solvation water is used. Simulations show that the observed experimental rates are compatible with considerably different distributions of the end-to-end distance; for realistic models, these are never Gaussian but indicative of a rugged energy landscape.     

Conserved central domains control the quaternary structure of type I and type II Hsp40 molecular chaperones

Heat shock protein (Hsp)40s play an essential role in protein metabolism by regulating the polypeptide binding and release cycle of Hsp70. The Hsp40 family is large, and specialized family members direct Hsp70 to perform highly specific tasks. Type I and Type II Hsp40s, such as yeast Ydj1 and Sis1, are homodimers that dictate functions of cytosolic Hsp70, but how they do so is unclear. Type I Hsp40s contain a conserved, centrally located cysteine-rich domain that is replaced by a glycine- and methionine-rich region in Type II Hsp40s, but the mechanism by which these unique domains influence Hsp40 structure and function is unknown. This is the case because high-resolution structures of full-length forms of these Hsp40s have not been solved. To fill this void, we built low-resolution models of the quaternary structure of Ydj1 and Sis1 with information obtained from biophysical measurements of protein shape, small-angle X-ray scattering, and ab initio protein modeling. Low-resolution models were also calculated for the chimeric Hsp40s YSY and SYS, in which the central domains of Ydj1 and Sis1 were exchanged. Similar to their human homologs, Ydj1 and Sis1 each has a unique shape with major structural differences apparently being the orientation of the J domains relative to the long axis of the dimers. Central domain swapping in YSY and SYS correlates with the switched ability of YSY and SYS to perform unique functions of Sis1 and Ydj1, respectively. Models for the mechanism by which the conserved cysteine-rich domain and glycine- and methionine-rich region confer structural and functional specificity to Type I and Type II Hsp40s are discussed.     

Comparison of genetic diversity at microsatellite loci and quantitative traits in hatchery populations of Japanese flounder Paralichthys olivaceus

Relationships of genetic diversity at microsatellite loci and quantitative traits were examined in hatchery-produced populations of Japanese flounder using a relatively straightforward experiment. Five hatchery populations produced by wild-caught and domesticated broodstocks were used to examine the effects of different levels (one to three generations) of domestication on the genetic characteristics of hatchery populations. Allelic richness at seven microsatellite loci in all hatchery populations was lower than that in a wild population. Genetic variation measured by allelic richness and heterozygosity tended to decrease with an increase in generations of domestication. In addition, the degree of genetic differentiation from a wild population increased with an increase in generations of domestication. Significant differences in three morphometric traits (dorsal and anal fin ray counts and vertebral counts) and three physiological traits (high temperature, salinity and formalin tolerance) were observed among the hatchery populations. The degree of phenotypic difference among populations was larger in morphometric traits than in physiological traits. The divergence pattern of some quantitative traits was similar to that observed at microsatellite loci, suggesting that domestication causes the decrease of genetic variation and the increase of genetic differentiation for some quantitative traits concomitantly with those for microsatellite loci. Significant positive correlation was observed between F _ST and the degree of phenotypic difference in the three morphometric traits and formalin tolerance, indicating that genetic variation at microsatellite loci predicts the degree of phenotypic divergence in some quantitative traits.     

Markedly discordant mitochondrial DNA and allozyme phylogenies of tube-nosed fruit bats, Nyctimene, at the Australian–Oriental biogeographical interface

Australo-Papuan tube-nosed fruit bats of the genus Nyctimene reach their most westerly distribution on the island archipelagos of eastern Indonesia. A recent morphological examination indicates three species occur on Moluccan islands. Both allozyme electrophoresis and mitochondrial DNA (mtDNA) sequencing reveal there are only two species present, Nyctimene cephalotes and Nyctimene albiventer , but there is considerable disparity between the results obtained from the two genetic data sets. Allozyme data indicate N. albiventer occurs on Wokam, which sits on the Australian continental shelf and was joined to Australo-Papua during the last glacial maximum, and N. cephalotes on the other Moluccan islands, all of which are off the continental shelf. Divergence of these two species is dated at approximately 1.5 Mya. By contrast, the mtDNA gene tree shows two deep clades, one containing all specimens from Wokam and Yamdena, and the second all the specimens from the other islands. This especially marked incongruence between the two genetic data sets is ascribed to either a single introgression event of N. albiventer mtDNA into N. cephalotes on Yamdena at around 100 000 ya or lineage sorting of very old mtDNA lineages that coalesce a rather long time before the speciation event. These results highlight that caution should be exercised when relying on mtDNA as barcodes in species taxonomy. The continental-associated N. albiventer on Wokam has higher levels of allozyme heterozygosity and mtDNA nucleotide diversity than the N. cephalotes populations occurring on the more remote islands, indicating that this colonizing species has experienced bottlenecks and/or a low effective population size since speciation.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 589–602.     

Limits to genetic bottlenecks and founder events imposed by the Allee effect

The Allee effect can result in a negative population growth rate at low population density. Consequently, populations below a minimum (critical) density are unlikely to persist. A lower limit on population size should constrain the loss of genetic variability due to genetic drift during population bottlenecks or founder events. We explored this phenomenon by modeling changes in genetic variability and differentiation during simulated bottlenecks of the alpine copepod, Hesperodiaptomus shoshone. Lake surveys, whole-lake re-introduction experiments and model calculations all indicate that H. shoshone should be unlikely to establish or persist at densities less than 0.5–5 individuals m⁻³. We estimated the corresponding range in minimum effective population size using the distribution of habitat (lake) sizes in nature and used these values to model the expected heterozygosity, allelic richness and genetic differentiation resulting from population bottlenecks. We found that during realistic bottlenecks or founder events, >90% of H. shoshone populations in the Sierra Nevada may be resistant to significant changes in heterozygosity or genetic distance, and 70–75% of populations may lose <10% of allelic richness. We suggest that ecological constraints on minimum population size be considered when using genetic markers to estimate historical population dynamics.     

Effects of habitat fragmentation, population size and demographic history on genetic diversity: the Cross River gorilla in a comparative context

In small and fragmented populations, genetic diversity may be reduced owing to increased levels of drift and inbreeding. This reduced diversity is often associated with decreased fitness and a higher threat of extinction. However, it is difficult to determine when a population has low diversity except in a comparative context. We assessed genetic variability in the critically endangered Cross River gorilla (Gorilla gorilla diehli), a small and fragmented population, using 11 autosomal microsatellite loci. We show that levels of diversity in the Cross River population are not evenly distributed across the three genetically identified subpopulations, and that one centrally located subpopulation has higher levels of variability than the others. All measures of genetic variability in the Cross River population were comparable to those of the similarly small mountain gorilla (G. beringei beringei) populations (Bwindi and Virunga). However, for some measures both the Cross River and mountain gorilla populations show lower levels of diversity than a sample from a large, continuous western gorilla population (Mondika, G. gorilla gorilla). Finally, we tested for the genetic signature of a bottleneck in each of the four populations. Only Cross River showed strong evidence of a reduction in population size, suggesting that the reduction in size of this population was more recent or abrupt than in the two mountain gorilla populations. These results emphasize the need for maintaining connectivity in fragmented populations and highlight the importance of allowing small populations to expand.     

Case study and application of process analytical technology (PAT) towards bioprocessing: II. Use of ultra-performance liquid chromatography (UPLC) for making real-time pooling decisions for process chromatography

Process analytical technology (PAT) has been gaining a lot of momentum in the biopharmaceutical community due to the potential for continuous real-time quality assurance resulting in improved operational control and compliance. Two of the key goals that have been outlined for PAT are "variability is managed by the process" and "product quality attributes can be accurately and reliably predicted over the design space established for materials used, process parameters, manufacturing, environmental, and other conditions". Recently, we have been examining the feasibility of applying different analytical tools for designing PAT applications for bioprocessing. We have previously shown that a commercially available online high performance liquid chromatography (HPLC) system can be used for analysis that can facilitate real-time decisions for column pooling based on product quality attributes (Rathore et al., 2008). In this article we test the feasibility of using a commercially available ultra- performance liquid chromatography (UPLC) system for real-time pooling of process chromatography columns. It is demonstrated that the UPLC system offers a feasible approach and meets the requirements of a PAT application. While the application presented here is of a reversed phase assay, the approach and the hardware can be easily applied to other modes of liquid chromatography.     

On-line near infrared monitoring of glycerol-boosted anaerobic digestion processes: evaluation of process analytical technologies

A study of NIR as a tool for process monitoring of thermophilic anaerobic digestion boosted by glycerol has been carried out, aiming at developing simple and robust Process Analytical Technology modalities for on-line surveillance in full scale biogas plants. Three 5 L laboratory fermenters equipped with on-line NIR sensor and special sampling stations were used as a basis for chemometric multivariate calibration. NIR characterisation using Transflexive Embedded Near Infra-Red Sensor (TENIRS) equipment integrated into an external recurrent loop on the fermentation reactors, allows for representative sampling, of the highly heterogeneous fermentation bio slurries. Glycerol is an important by-product from the increasing European bio-diesel production. Glycerol addition can boost biogas yields, if not exceeding a limiting 5-7 g L(-1) concentration inside the fermenter-further increase can cause strong imbalance in the anaerobic digestion process. A secondary objective was to evaluate the effect of addition of glycerol, in a spiking experiment which introduced increasing organic overloading as monitored by volatile fatty acids (VFA) levels. High correlation between on-line NIR determinations of glycerol and VFA contents has been documented. Chemometric regression models (PLS) between glycerol and NIR spectra needed no outlier removals and only one PLS-component was required. Test set validation resulted in excellent measures of prediction performance, precision: r(2) = 0.96 and accuracy = 1.04, slope of predicted versus reference fitting. Similar prediction statistics for acetic acid, iso-butanoic acid and total VFA proves that process NIR spectroscopy is able to quantify all pertinent levels of both volatile fatty acids and glycerol.     

Effect of Calcium Ions on Hydrodynamic Properties of Pentameric and Decameric C-Reactive Protein in Solution

The molecular mass and sedimentation coefficient of native C-reactive protein in solution were determined by analytical ultracentrifugation in the presence and absence of calcium ions. Pentameric C-reactive protein was shown to be the major macroscopic form of this protein in solution. The removal of calcium ions from solution caused decompaction of the protein accompanied by changes in its hydrodynamic parameters. The sedimentation coefficient s ⁰ _20,w of pentameric C-reactive protein in solution containing 2 mM Ca²⁺ (6.6S) exceeded that for C-reactive protein in solution containing 2 mM EDTA (6.4S). Analysis of average molecular masses M _w and M _z obtained from sedimentation data demonstrated that the solution of highly purified protein was not homogeneous. As shown by intermolecular crosslinking, the solution also contained the 241-kDa decamer of C-reactive protein (9.5S) as a separate macroscopic form, whose share hardly reached 10% in the presence of 2 mM Ca²⁺ and increased after removal of calcium ions. The decamers were shown to result from intermolecular association of the pentamers.