Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

On the statistical analysis of the GS-NS0 cell proteome: Imputation, clustering and variability testing

We have undertaken two-dimensional gel electrophoresis proteomic profiling on a series of cell lines with different recombinant antibody production rates. Due to the nature of gel-based experiments not all protein spots are detected across all samples in an experiment, and hence datasets are invariably incomplete. New approaches are therefore required for the analysis of such graduated datasets. We approached this problem in two ways. Firstly, we applied a missing value imputation technique to calculate missing data points. Secondly, we combined a singular value decomposition based hierarchical clustering with the expression variability test to identify protein spots whose expression correlates with increased antibody production. The results have shown that while imputation of missing data was a useful method to improve the statistical analysis of such data sets, this was of limited use in differentiating between the samples investigated, and highlighted a small number of candidate proteins for further investigation. (c) 2006 Elsevier B.V. All rights reserved.

Phytotoxicity of surface waters of the Thames and Brisbane River Estuaries: A combined chemical analysis and bioassay approach for the comparison of two systems

The Thames Estuary, UK, and the Brisbane River, Australia, are comparable in size and catchment area. Both are representative of the large and growing number of the world's estuaries associated with major cities. Principle differences between the two systems relate to climate and human population pressures. In order to assess the potential phytotoxic impact of herbicide residues in the estuaries, surface waters were analysed with a PAM fluorometry-based bioassay that employs the photosynthetic efficiency (photosystem II quantum yield) of laboratory cultured microalgae, as an endpoint measure of phytotoxicity. In addition, surface waters were chemically analysed for a limited number of herbicides. Diuron atrazine and simazine were detected in both systems at comparable concentrations. In contrast, bioassay results revealed that whilst detected herbicides accounted for the observed phytotoxicity of Brisbane River extracts with great accuracy, they consistently explained only around 50% of the phytotoxicity induced by Thames Estuary extracts. Unaccounted for phytotoxicity in Thames surface waters is indicative of unidentified phytotoxins. The greatest phytotoxic response was measured at Charing Cross, Thames Estuary, and corresponded to a diuron equivalent concentration of 180 ng L-1. The study employs relative potencies (REP) of PSII impacting herbicides and demonstrates that chemical analysis alone is prone to omission of valuable information. Results of the study provide support for the incorporation of bioassays into routine monitoring programs where bioassay data may be used to predict and verify chemical contamination data, alert to unidentified compounds and provide the user with information regarding cumulative toxicity of complex mixtures. (c) 2005 Elsevier B.V. All rights reserved.

Preliminary crystallographic characterization of an RNA helicase from Kunjin virus

Kunjin virus is a member of the Flavivirus genus and is an Australian variant of West Nile virus. The C-terminal domain of the Kunjin virus NS3 protein displays helicase activity. The protein is thought to separate daughter and template RNA strands, assisting the initiation of replication by unwinding RNA secondary structure in the 3' nontranslated region. Expression, purification and preliminary crystallographic characterization of the NS3 helicase domain are reported. It is shown that Kunjin virus helicase may adopt a dimeric assembly in absence of nucleic acids, oligomerization being a means to provide the helicases with multiple nucleic acid-binding capability, facilitating translocation along the RNA strands. Kunjin virus NS3 helicase domain is an attractive model for studying the molecular mechanisms of flavivirus replication, while simultaneously providing a new basis for the rational development of anti-flaviviral compounds.

Parasitic castration by the digenian trematode Allopodocotyle sp alters gene expression in the brain of the host mollusc Haliotis asinina

Infection of molluscs by digenean trematode parasites typically results in the repression of reproduction - the so-called parasitic castration. This is known to occur by altering the expression of a range of host neuropeptide genes. Here we analyse the expression levels of 10 members of POU, Pax, Sox and Hox transcription factor gene families, along with genes encoding FNIRFamide, prohormone convertase and P-tubulin, in the brain ganglia of actively reproducing (summer), non-reproducing (winter) and infected Haliotis asinina (a vetigastropod mollusc). A number of the regulatory genes are differentially expressed in parasitised H. asinina, but in only a few cases do expression patterns in infected animals match those occurring in animals where reproduction is normally repressed. (c) 2006 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

Predicting residue-wise contact orders in proteins by support vector regression

Background: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

Reference values for whole body and cerebral multi-frequency bio-impedance data in neonates less than 12 h postpartum

Multiple frequency bio-electrical impedance analysis (MFBIA) may be useful for monitoring fluid balance in newborn infants or to provide early prediction of the outcome following perinatal asphyxia. A reference range of data is needed for identification of babies with abnormal impedance values. This was a cross-sectional observational study in 84 term and near-term healthy neonates less than 12 h postpartum. Whole body and cerebral MFBIA measurements were performed at the bedside in the post-natal ward. Gestational age, post-natal age, gender, birthweight, head circumference and foot length measures were recorded. Reference values for impedance at the characteristic frequency (Z(C)) and resistance at zero frequency (R-0) are reported for whole body and cerebral impedance. Significant correlations (p < 0.05) were observed between whole body impedance and birthweight, footlength and head circumference. Females had a significantly higher whole body R0 than males. Cerebral impedance did not correlate significantly with any of the demographic measures and therewere no gender differences observed for cerebral impedance. The reference range for whole body multi-frequency bio-impedance values in term and near-term infants within the first 12 h postpartum can be calculated from the footlength (FL) using the following equations: Z(C) = (942.9 - 4.818* FL) +/- 124.6 Omega; R-0 = (1042 - 4.520(*)FL) +/- 135.5 Omega. For cerebral impedance the reference range is 29.5-48.7 Omega for Z(C) and 33.7-58.0 Omega for R-0.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Are corals colorful?

Using in situ spectrometry data and visual system modeling, we investigate whether the colors conferred to the reef-building corals by GFP-like proteins would look colorful not only to humans, but also to fish occupying different ecological niches on the reef. Some GFP-like proteins, most notably fluorescent greens and nonfluorescent chromoproteins, indeed generate intense color signals. An unexpected finding was that fluorescent proteins might also make corals appear less colorful to fish, counterbalancing the effect of absorption by the photosynthetic pigments of the endosymbiotic algae, which might be a form of protection against herbivores. We conclude that GFP-determined coloration of corals may be an important factor in visual ecology of the reef fishes.

The cyclotide family of circular miniproteins: Nature's combinatorial peptide template

The cyclotides are a recently discovered family of miniproteins that contain a head-to-tail cyclized backbone and a knotted arrangement of disulfide bonds. They are approximately 30 amino acids in size and are present in high abundance in plants from the Violaceae, Rubiaceae, and Cucurbitaceae families, with individual plants containing a suite of up to 100 cyclotides. They have a diverse range of biological activities, including uterotonic, anti-HIV, antitumor, and antimicrobial activities, although their natural function is likely that of defending their host plants from pathogens and pests. This review focuses on the structural aspects of cyclotides, which may be thought of as a natural combinatorial peptide template in which a wide range of amino acids is displayed on a compact molecular core made up of the cyclic cystine knot structural motif. Cyclotides are exceptionally stable and are resistant to denaturation via thermal, chemical, or enzymatic treatments. The struclural features that contribute to their remarkable stability are described ill this review. (c) 2006 Wiley Periodicals, Inc.

The discriminative bilateral filter: An enhanced denoising filter for electron microscopy data

Advances in three-dimensional (313) electron microscopy (EM) and image processing are providing considerable improvements in the resolution of subcellular volumes, macromolecular assemblies and individual proteins. However, the recovery of high-frequency information from biological samples is hindered by specimen sensitivity to beam damage. Low dose electron cryo-microscopy conditions afford reduced beam damage but typically yield images with reduced contrast and low signal-to-noise ratios (SNRs). Here, we describe the properties of a new discriminative bilateral (DBL) filter that is based upon the bilateral filter implementation of Jiang et al. (Jiang, W., Baker, M.L., Wu, Q., Bajaj, C., Chin, W., 2003. Applications of a bilateral denoising filter in biological electron microscopy. J. Struc. Biol. 128, 82-97.). In contrast to the latter, the DBL filter can distinguish between object edges and high-frequency noise pixels through the use of an additional photometric exclusion function. As a result, high frequency noise pixels are smoothed, yet object edge detail is preserved. In the present study, we show that the DBL filter effectively reduces noise in low SNR single particle data as well as cellular tomograms of stained plastic sections. The properties of the DBL filter are discussed in terms of its usefulness for single particle analysis and for pre-processing cellular tomograms ahead of image segmentation. (c) 2006 Elsevier Inc. All rights reserved.

The role of histidine residues in low-pH-mediated viral membrane fusion

A central event in the invasion of a host cell by an enveloped virus is the fusion of viral and cell membranes. For many viruses, membrane fusion is driven by specific viral surface proteins that undergo large-scale conformational rearrangements, triggered by exposure to low pH in the endosome upon internalization. Here, we present evidence suggesting that in both class I (helical hairpin proteins) and class 11 (beta-structure-rich proteins) pH-dependent fusion proteins the protonation of specific histidine residues triggers fusion via an analogous molecular mechanism. These histidines are located in the vicinity of positively charged residues in the prefusion conformation, and they subsequently form salt bridges with negatively charged residues in the postfusion conformation. The molecular surfaces involved in the corresponding structural rearrangements leading to fusion are highly conserved and thus might provide a suitable common target for the design of antivirals, which could be active against a diverse range of pathogenic viruses.