Onset of quark-hadron duality in pion electroproduction

A large data set of charged-pion (pi+/-) electroproduction from both hydrogen and deuterium targets has been obtained spanning the low-energy residual-mass region. These data conclusively show the onset of the quark-hadron duality phenomenon, as predicted for high-energy hadron electroproduction. We construct several ratios from these data to exhibit the relation of this phenomenon to the high-energy factorization ansatz of electron-quark scattering and subsequent quark-->pion production mechanisms.

Determination of the pion charge form factor at Q2=1.60 and 2.45 (GeV/c)2

The 1H(e,e'pi+)n cross section was measured at four-momentum transfers of Q2=1.60 and 2.45 GeV2 at an invariant mass of the photon nucleon system of W=2.22 GeV. The charged pion form factor (F(pi)) was extracted from the data by comparing the separated longitudinal pion electroproduction cross section to a Regge model prediction in which F(pi) is a free parameter. The results indicate that the pion form factor deviates from the charge-radius constrained monopole form at these values of Q2 by one sigma, but is still far from its perturbative quantum chromodynamics prediction.

Yellow leaf of sugarcane is caused by at least three different genotypes of sugarcane yellow leaf virus, one of which predominates on the Island of Réunion

The genetic diversity of sugarcane yellow leaf virus (SCYLV) was analyzed with 43 virus isolates from Réunion Island and 17 isolates from world-wide locations. We attempted to amplify by reverse-transcription polymerase chain reaction (RT-PCR), clone, and sequence four different fragments covering 72% of the genome of these virus isolates. The number of amplified isolates and useful sequence information varied according to each fragment, whereas an amplicon was obtained with diagnostic primers for 59 out of 60 isolates (98%). Phylogenetic analyses of the sequences determined here and additional sequences of 11 other SCYLV isolates available from GenBank showed that SCYLV isolates were distributed in different phylogenetic groups or belonged to single genotypes. The majority of isolates from Réunion Island were grouped in phylogenetic clusters that did not contain any isolates from other origins. The complete six ORFs (5612 bp) of five SCYLV isolates (two from Réunion Island, one from Brazil, one from China, and one from Peru) were amplified, cloned, and sequenced. The existence of at least three distinct genotypes of SCYLV was shown by phylogenetic analysis of the sequences of these isolates and additional published sequences of three SCYLV isolates (GenBank accessions). The biological significance of these genotypes and of the origin of the distinct lineage of SCYLV in Réunion Island remains to be determined.

Effects of nitrate on the tadpoles of two Ranids (Rana catesbeiana and R. clamitans)

Nitrate is a potential stressor of amphibian larvae. Previous research has shown variation in the effects of nitrate among species and even populations of amphibians. However, relatively few species and populations of amphibians have been examined for nitrate tolerance, especially from populations in the agricultural regions of the United States. Using a 15-day laboratory experiment, we investigated the effects of nitrate on the larvae of two species of anurans (Rana catesbeiana and R. clamitans) that are common throughout the agricultural Midwest. Survival of R. catesbeiana and R. clamitans tadpoles was negatively affected by higher concentrations of nitrate (20 mg L(-1)), with R. clamitans being more affected than R. catesbeiana. The final mass of R. catesbeiana tadpoles was affected by nitrate concentration, with tadpoles exposed to intermediate concentrations of nitrate (5 mg L(-1)) being the heaviest. The final mass of R. clamitans tadpoles was not affected by nitrate concentration. These experiments suggest that these two species differ in their responses to the toxic effects of nitrate and, along with previous results on Ranids, suggest there is a great deal of variation in nitrate tolerance in this family.

Strange-quark contributions to parity-violating asymmetries in the forward g0 electron-proton scattering experiment

We have measured parity-violating asymmetries in elastic electron-proton scattering over the range of momentum transfers 0.12 < or =Q2 < or =1.0 GeV2. These asymmetries, arising from interference of the electromagnetic and neutral weak interactions, are sensitive to strange-quark contributions to the currents of the proton. The measurements were made at Jefferson Laboratory using a toroidal spectrometer to detect the recoiling protons from a liquid hydrogen target. The results indicate nonzero, Q2 dependent, strange-quark contributions and provide new information beyond that obtained in previous experiments.

Incorporation of flexibility into rigid-body docking: Applications in rounds 3-5 of CAPRI

We have submitted models for all 9 targets in Rounds 3-5 of CAPRI and have predicted at least 30% of the correct contacts for 4 of the targets and at least 10% of the correct contacts for another 4 targets. We have employed a variety of techniques but have had the greatest success by combining established rigid-body docking with a variety of initial conformations generated by molecular dynamics.

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.

A link between sequence conservation and domain motion within the AAA+ family

The AAA+ family of proteins play fundamental roles in all three kingdoms of life. It is thought that they act as molecular chaperones in aiding the assembly or disassembly of proteins or protein complexes. Recent structural studies on a number of AAA+ family proteins have revealed that they share similar structural elements. These structures provide a possible link between nucleotide binding/hydrolysis and the conformational changes which are then amplified to generate mechanical forces for their specific functions. However, from these individual studies it is far from clear whether AAA+ proteins in general share properties in terms of nucleotide induced conformational changes. In this study, we analyze sequence conservation within the AAA+ family and identify two subfamilies, each with a distinct conserved linker sequence that may transfer conformational changes upon ATP binding/release to movements between subdomains and attached domains. To investigate the relation of these linker sequences to conformational changes, molecular dynamics (MD) simulations on X-ray structures of AAA+ proteins from each subfamily have been performed. These simulations show differences in both the N-linker peptide, subdomain motion, and cooperativity between elements of quaternary structure. Extrapolation of subdomain movements from one MD simulation enables us to produce a structure in close agreement with cryo-EM experiments.

A variable region of the sugarcane bacilliform virus (SCBV) genome can be used to generate promoters for transgene expression in sugarcane

Four promoters derived from sugarcane bacilliform virus (SCBV) were compared and characterised. Three were obtained by PCR amplification of purified virion DNA extracted from three sugarcane cultivars. The fourth promoter was obtained by subcloning from an almost genome-length clone of SCBV. All promoters were able to drive stable expression of beta-glucuronidase in sugarcane. The PCR-derived promoter sequences shared more DNA homology with banana streak virus than to the subcloned SCBV. The subcloned promoter was the strongest expressing and was able to drive reporter gene expression in vitro and in the leaves, meristems and roots of glasshouse-grown sugarcane. Expression levels were at least equal to or higher than those measured for the maize polyubiquitin promoter.

Etanercept treatment of renal amyloidosis complicating rheumatoid arthritis

Rheumatoid arthritis and juvenile arthritis represent the commonest diseases complicated by AA amyloidosis in developed countries. Up to 5% of patients with rheumatoid arthritis will develop AA amyloidosis, with renal failure being the commonest cause of mortality. To date, treatment of this condition has focused on suppressing the underlying inflammatory condition with drugs such as cyclophosphamide and chlorambucil, but both these drugs are associated with myelotoxicity, leukaemia and sterility. Tumour necrosis factor-alpha (TNF-alpha) is thought to be involved in amyloid deposition. The efficacy of anti-TNF-alpha therapy (etanercept) in the treatment of renal amyloidosis complicating rheumatoid arthritis is demonstrated here and the current scientific data on this subject are presented.

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.

Molecular analysis of Fiji disease virus genome segments 5, 6, 8 and 10

The complete sequences of Fiji disease virus (FDV) genome segments 5 (S5), S6, S8 and S10 were obtained and comprised 3150 nt, 2831 nt, 1959 nt and 1819 nt, respectively. Each segment contained a single ORF which encoded putative proteins of 115 kDa, 97 kDa, 69 kDa and 63.0 kDa, respectively. The putative amino acid sequences encoded by S5 and S6 contained putative leucine zipper motifs while FDV S5 and S8 each contained an ATP-GTP-binding motif. At the amino acid level, FDV S5, S6, S8 and S10 showed most similarity to the corresponding segments of Rice black-streaked dwarf virus. Based on sequence similarities, it is predicted that FDV S8 encodes a minor core protein, while FDV S10 encodes an outer capsid protein. The evolutionary relationships of FDV to other reoviruses are discussed.

Sugarcane bacilliform virus encapsidates genome concatamers and does not appear to integrate into the Saccharum officinarum genome

Sugarcane bacilliform virus (SCBV) DNA molecules larger than the complete genome length of 7.6 kbp were detected in infected plants and in virions. We have confirmed that these high molecular weight nucleic acids were open circular DNA and viral in origin. Due to their open circular conformation, accurate size determination of the DNA molecules was not possible using conventional electrophoresis. Using field inversion gel electrophoresis (FIGE), however, the DNA appeared to increase in genome size increments, with sizes ranging from 1 to 4 genomes (31 kbp) detected. The DNA was packaged into virions, which may explain the observation of purified virions with lengths corresponding to one, two or three times the modal length of 130 nm. The DNA products were possibly concatamers formed during replication as a result of a terminal overlap on the sense strand, and were shown to be overlapped individual genome-length molecules and not covalently-bonded continuous DNA strands. Southern analysis indicated that SCBV sequences are not integrated into the sugarcane genome and that the high molecular weight DNA observed in the sugarcane accessions analysed represents SCBV concatamers.

Novel use of a genetic algorithm for protein structure prediction: Searching template and sequence alignment space

A novel genetic algorithm was applied to all CASP5 targets. The algorithm simultaneously searches template and alignment space. Results show that the current implementation of the method is perhaps most useful in recognizing and refining remote homology targets. This new method is briefly described and results are analyzed. Strengths and weaknesses of the current implementation of the algorithm are discussed.

Functional Recycling of C2 Domains Throughout Evolution: A Comparative Study of Synaptotagmin, Protein Kinase C and Phospholipase C by Sequence, Structural and Modelling Approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level, the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various sub-families. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces, our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics.

Molecular dynamics simulations of the bacterial outer membrane protein FhuA: a comparative study of the ferrichrome-free and bound states

FhuA is one of the more complex members of the superfamily of bacterial outer membrane proteins. Its primary function is to provide a binding site on the outer membrane surface for siderophores, such as ferrichrome, and subsequently to facilitate their energy-dependent transport across the membrane, presumably powered by the TonB-ExbBD protein complex that resides in the cytoplasmic membrane. Crystal structures of FhuA with and without a bound ferrichrome molecule have provided some clues as to the initial stages of the siderophore transport mechanism. In the current study, we have employed 10-ns duration molecular dynamics simulations of FhuA and of the FhuA-ferrichrome complex, both embedded in a phospholipid bilayer, to probe the short timescale dynamics of this integral membrane protein, and to explore possible mechanistic implications of this dynamic behavior. Analysis of the dynamics of the protein suggests that the extracellular loops move as relatively rigid entities relative to the transmembrane beta-barrel. Comparison of the two simulations (with and without bound ferrichrome) revealed some ligand-induced changes in loop mobility. Specifically, loop L8 appears to be involved in a mechanism whereby the binding site is gated closed upon ligand binding. Analysis of the dynamics of water molecules within the core of the FhuA protein provided no evidence for a water-permeable protopore through which the ferrichrome might pass without a major perturbation of the FhuA protein. Overall, these simulations support the proposal that binding of ferrichrome initiates a signaling mechanism that ultimately leads to the TonB-mediated partial or total removal of the core domain from the beta-barrel, thus opening up a permeable pore. These simulations are among the longest that have been performed on a complex membrane protein. However, a simple analysis of sampling reveals that the description of protein motions is far from complete.

Evaluation of the 3D-Dock protein docking suite in rounds 1 and 2 of the CAPRI blind trial

The 3D-Dock suite of programs has been used to make predictions for the seven targets in rounds 1 and 2 of the CAPRI method evaluation exercise. Some correct contacts were obtained in at least one prediction for four of seven targets. Target 06 was predicted very well, with an RMSD of the ligand after superimposition of the receptor of only 0.77 A. We investigate the performance of the various stages of the method, with the aim of finding where improvements need to be made, and in particular whether the manual interventions that were made were essential, and whether results of the level of accuracy obtained for target 06 may be expected with confidence.

Sequence analysis of an Australian isolate of sugarcane bacilliform badnavirus

The genome of an Australian isolate of Sugarcane bacilliform virus (SCBV-IM) was cloned, sequenced and analysed. The genome consisted of 7687 nucleotides and contained three open reading frames which were similar in size and organisation to those of other badnaviruses. SCBV-IM was found to be most similar to the SCBV-Morocco isolate with amino acid sequence similarity of 91.4 %, 83.8 % and 85.3 % in the ORF I, II and III coding regions, respectively. Phylogenetic analysis of the SCBV-IM ORF III deduced amino acid sequence showed that SCBV isolates were more closely related to each other than to other badnaviruses. Amplification of SCBV sequences from three different sugarcane varieties revealed considerable variability in the viral populations, both within single infected plants as well as between infected plants, suggesting that the SCBV isolates sequenced to date may not be representative of the range of virus variability.

Free energy of a potassium ion in a model of the channel formed by an amphipathic leucine-serine peptide

We use molecular dynamics simulations to investigate the position-dependent free energy of a potassium ion in a model of an ion channel formed by the synthetic amphipathic leucine-serine peptide, LS3. The channel model is a parallel bundle of six LS3 helices around which are packed 146 methane-like spheres in order to mimic a membrane. At either end of and within the channel are 1051 water molecules, plus four ions (two potassium and two chloride). The free energy of a potassium ion in the channel was estimated using the weighted histogram analysis (WHAM) method. This is the first time to our knowledge that such a calculation has been carried out as a function of the position of an ion in three dimensions within a channel. The results indicate that for this channel, which is lined by hydrophilic serine sidechains, there is a relatively weak dependence of the free energy on the axial/off-axial position of the ion. There are some off-axis local minima, especially in the C-terminal half of the channel. Using the free energy results, a single channel current-voltage curve was estimated using a one-dimensional Nernst-Planck equation. Although reasonable agreement with experiment is achieved for K(+) ions flowing from the N-terminal to the C-terminal mouth, in the opposite direction the current is underestimated. This underestimation may be a consequence of under-sampling of the conformational dynamics of the channel. We suggest that our simulations may have captured, for example, a sub-conductance level (i.e. an incompletely open state) of the LS3 channel.

Setting up and optimization of membrane protein simulations

In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins.

Novel alpha7-like nicotinic acetylcholine receptor subunits in the nematode Caenorhabditis elegans

We have used reverse-transcription-polymerase chain reaction (RT-PCR) and DNA sequencing techniques to confirm the transcription of seven (six alpha and one non-alpha) novel candidate nicotinic acetylcholine receptor (nAChR) subunit-encoding genes identified in the genome sequence of the nematode Caenorhabditis elegans. Compared to vertebrate nAChR subunits, they most closely resemble the homomer-forming, neuronal alpha7 subunit. Comparison of the predicted amino acid sequences of the new nAChR subunits with those described previously in C. elegans reveals five subunits (four alpha and one non-alpha) which resemble the DEG-3-like group of subunits. To date, this highly divergent nAChR subunit group is unique to C. elegans. ACR-22 is the first non-alpha member of the DEG-3-like group of subunits to be identified. Two new members of the related ACR-16-like nAChR group of subunits have also been shown to be transcribed, making the ACR-16-like subunit group the largest in C. elegans. Residues in the alpha subunit second transmembrane region (M2) which contribute to the channel lining show variations with implications for channel function. For example, in ACR-22, the highly conserved 0' lysine of M2 is replaced by histidine. Restrained molecular dynamics simulations have been used to generate molecular models of homo-pentameric M2 helix bundles for the novel subunits, enabling identification and display of pore-lining and protein interface residues. The number and diversity of genes encoding C. elegans nAChR subunits with similarities to the homomer-forming vertebrate alpha7 subunits and the identification of related non-alpha subunits, only found in C. elegans to date, suggest that at least some of these subunits may contribute to heteromers in vivo.

Treatment of depression in rural Arkansas: policy implications for improving care

Policy-makers have long suspected that greater barriers to care result in depressed rural residents being less likely to receive high-quality treatment. This study recruited 470 depressed community residents in a 1992 telephone survey, followed 95 percent of them through one year, and abstracted additional data on their health care utilization from insurance claims, medical and pharmacy records. Bivariate and multivariate models demonstrated that during the year following the baseline, there were no significant rural-urban differences in the rate (probability of any outpatient depression treatment), type (probability of receiving general medical depression care only), or quality (completion of guideline-concordant acute-stage care) of outpatient depression treatment. Annual expenditures for outpatient depression treatment were lower for rural subjects compared with their urban counterparts. Rural subjects had 3.05 times the odds of being admitted to a hospital for physical problems and 3.06 times the odds of being admitted to a hospital for mental health problems during the year following baseline compared with urban subjects. Cost-offset analyses demonstrate that every dollar invested in depression treatment was associated with a $2.61 decrease in the cost of treating physical problems in depressed rural residents. Limited insurance coverage and limited availability of services were the most significant barriers to specialty and general medical outpatient treatment for depression in both rural and urban residents. More than 80 percent of depressed residents in both rural and urban areas visited a primary care provider during the year following baseline. The potential cost offset of depression treatment in rural populations plus the improvement in productivity observed in both rural and urban populations indicate that it may be economically possible to improve quality of care for depression without bankrupting an already strained health care budget.

Prediction of protein-protein interactions by docking methods

Recently, developments have been made in predicting the structure of docked complexes when the coordinates of the components are known. The process generally consists of a stage during which the components are combined rigidly and then a refinement stage. Several rapid new algorithms have been introduced in the rigid docking problem and promising refinement techniques have been developed, based on modified molecular mechanics force fields and empirical measures of desolvation, combined with minimisations that switch on the short-range interactions gradually. There has also been progress in developing a benchmark set of targets for docking and a blind trial, similar to the trials of protein structure prediction, has taken place.