Calcein as a fluorescent iron chemosensor for the determination of low molecular weight iron in biological fluids

The fluorescence quenching of calcein (CA) is not iron specific and results in a negative calibration curve. In the present study, deferoxamine (DFO), a strong iron chelator, was used to regenerate the fluorescence quenched by iron. Therefore, the differences in fluorescence reading of the same sample with or without addition of DFO are positively and specifically proportional to the amounts of iron. We found that the same iron species but different anions (e.g. ferric sulfate or ferric citrate) differed in CA fluorescence quenching, so did the same anions but different iron (e.g. ferrous or ferric sulfates). Excessive amounts of citrate competed with CA for iron and citrate could be removed by barium precipitation. After optimizing the experimental conditions, the sensitivity of the fluorescent CA assay is 0.02 M of iron, at least 10 times more sensitive than the colorimetric assays. Sera from 6 healthy subjects were tested for low molecular weight (LMW) chelator bound iron in the filtrates of 10 kDa nominal molecular weight limit (NMWL). The LMW iron was marginally detectable in the normal sera. However, increased levels of LMW iron were obtained at higher transferrin (Tf) saturation (1.64–2.54 M range at 80% Tf saturation, 2.77–3.15 M range at 100% Tf saturation and 3.09–3.39 M range at 120% Tf saturation). The application of the assay was further demonstrated in the filtrates of human liver HepG2 and human lung epithelial A549 cells treated with iron or iron-containing dusts.     

Weak phylogenetic effects on ecological niches of Sylvia warblers

To understand the evolution of ecological niches it is important to know whether niche evolution is constrained by phylogeny. We approached this question for Sylvia warblers by testing if closely related species are more similar in 20 ecologically relevant morphological traits than distantly related species. Phylogenetic relatedness was quantified using a molecular phylogeny based on the mitochondrial cytochrome b gene. By Principal Component Analysis (PCA) two major niche axes were extracted. We tested the individual ecomorphological traits and the positions of the species on the PCA axes for phylogenetic effects using Mantel tests. The results demonstrated small but significant phylogenetic effects only for the length of the middle toe, a trait probably correlated with locomotion. In general, however, phylogenetic effects were very weak. This suggests that ecological niches in passerine birds have the potential to evolve rapidly and are not subject to major phylogenetic constraints.     

Testing the link between the latitudinal gradient in species richness and rates of molecular evolution

Numerous hypotheses have been proposed to explain latitudinal gradients in species richness, but all are subject to ongoing debate. Here we examine Rohde's (1978, 1992) hypothesis, which proposes that climatic conditions at low latitudes lead to elevated rates of speciation. This hypothesis predicts that rates of molecular evolution should increase towards lower latitudes, but this prediction has never been tested. We discuss potential links between rates of molecular evolution and latitudinal diversity gradients, and present the first test of latitudinal variation in rates of molecular evolution. Using 45 phylogenetically independent, latitudinally separated pairs of bird species and higher taxa, we compare rates of evolution of two mitochondrial genes and DNA-DNA hybridization distances. We find no support for an effect of latitude on rate of molecular evolution. This result casts doubt on the generality of a key component of Rohde's hypothesis linking climate and speciation.     

Molecular dynamics simulation of dark-adapted rhodopsin in an explicit membrane bilayer: coupling between local retinal and larger scale conformational change

The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature of this change, involving a coupling between a local process and larger changes throughout the protein, may be important for many membrane proteins. In addition, functional work has shown that this coupling occurs with different efficiency in different lipid settings. To begin to understand the nature of the efficiency of this coupling in different lipid settings, we present a molecular dynamics study of rhodopsin in an explicit dioleoyl-phosphatidylcholine bilayer. Our system was simulated for 40 ns and provides insights into the very early events of the visual cascade, before the full transition and activation have occurred. In particular, we see an event near 10 ns that begins with a change in hydrogen bonding near the retinal and that leads through a series of coupled changes to a shift in helical tilt. This type of event, though rare on the molecular dynamics time-scale, could be an important clue to the types of coupling that occur between local and large-scale conformational change in many membrane proteins.     

Torque-speed relationship of the Na+-driven flagellar motor of Vibrio alginolyticus

The torque-speed relationship of the Na(+)-driven flagellar motor of Vibrio alginolyticus was investigated. The rotation rate of the motor was measured by following the position of a bead, attached to a flagellar filament, using optical nanometry. In the presence of 50mM NaCl, the generated torque was relatively constant ( approximately 3800pNnm) at lower speeds (speeds up to approximately 300Hz) and then decreased steeply, similar to the H(+)-driven flagellar motor of Escherichia coli. When the external NaCl concentration was varied, the generated torque of the flagellar motor was changed over a wide range of speeds. This result could be reproduced using a simple kinetic model, which takes into consideration the association and dissociation of Na(+) onto the motor. These results imply that for a complete understanding of the mechanism of flagellar rotation it is essential to consider both the electrochemical gradient and the absolute concentration of the coupling ion.     

Drug recognition and stabilisation of the parallel-stranded DNA quadruplex d(TTAGGGT)4 containing the human telomeric repeat

The NMR structure of the parallel-stranded DNA quadruplex d(TTAGGGT)(4), containing the human telomeric repeat, has been determined in solution in complex with a fluorinated pentacyclic quino[4,3,2-kl]acridinium cation (RHPS4). RHPS4 has been identified as a potent inhibitor of telomerase at submicromolar levels (IC(50) value of 0.33(+/-0.13)microM), exhibiting a wide differential between telomerase inhibition and acute cellular toxicity. All of the data point to RHPS4 exerting its chemotherapeutic potency through interaction with, and stabilisation of, four-stranded G-quadruplex structures. RHPS4 forms a dynamic interaction with d(TTAGGGT)(4), as evident from 1H and 19F linewidths, with fast exchange between binding sites induced at 318 K. Perturbations to DNA chemical shifts and 24 intermolecular nuclear Overhauser effects (NOEs) identify the 5'-ApG and 5'-GpT steps as the principle intercalation sites; a structural model has been refined using NOE-restrained molecular dynamics. The central G-tetrad core remains intact, with drug molecules stacking at the ends of the G-quadruplex. The partial positive charge on position 13-N of the acridine ring appears to act as a "pseudo" potassium ion and is positioned above the centre of the G-tetrad in the region of high negative charge density. In both ApG and GpT intercalation sites, the drug is seen to converge to the same orientation in which the pi-system of the drug overlaps primarily with two bases of each G-tetrad. The drug is held in place by stacking interactions with the G-tetrads; however, there is some evidence for a more dynamic, weakly stabilised A-tetrad that stacks partially on top of the drug at the 5'-end of the sequence. Together, the interactions of RHPS4 increase the t(m) of the quadruplex by approximately 20 degrees C. There is no evidence for drug intercalation within the G-quadruplex; however, the structural model strongly supports end-stacking interactions with the terminal G-tetrads.     

A structure-based anatomy of the E.coli metabolome

The Escherichia coli metabolome has been characterised using the two-dimensional structures of 745 metabolites, obtained from the EcoCyc and KEGG databases. Physicochemical properties of the metabolome have been calculated to provide an overview of this set of cognate ligands. A library of fragments commonly found among these molecules has been employed to reveal the main constituents of metabolites, and to assist a broad classification of the metabolome into biochemically relevant classes. Fragment-based fingerprints reveal the metabolome as a continuum in the two-dimensional structural space, where clusters of molecules sharing similar scaffolds can be identified, but are generally overlapping. Nucleotide, carbohydrate and amino acid-like molecules are the most prominent, but at high levels of similarity, a more detailed classification is possible. Classification schemes for the metabolome are a promising tool for understanding the chemical diversity of the metabolome. When used in conjunction with existing classifications of the proteome, they can help to elucidate the binding preferences and promiscuity of proteins and their cognate substrates.     

Folding of Cu, Zn superoxide dismutase and familial amyotrophic lateral sclerosis

Cu, Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity.     

Protein dynamics in a family of laboratory evolved thermophilic enzymes

Molecular dynamics simulations were employed to study how protein solution structure and dynamics are affected by adaptation to high temperature. Simulations were carried out on a para-nitrobenzyl esterase (484 residues) and two thermostable variants that were generated by laboratory evolution. Although these variants display much higher melting temperatures than wild-type (up to 18 degrees C higher) they are both >97% identical in sequence to the wild-type. In simulations at 300 K the thermostable variants remain closer to their crystal structures than wild-type. However, they also display increased fluctuations about their time-averaged structures. Additionally, both variants show a small but significant increase in radius of gyration relative to wild-type. The vibrational density of states was calculated for each of the esterases. While the density of states profiles are similar overall, both thermostable mutants show increased populations of the very lowest frequency modes (<10 cm(-1)), with the more stable mutant showing the larger increase. This indicates that the thermally stable variants experience increased concerted motions relative to wild-type. Taken together, these data suggest that adaptation for high temperature stability has resulted in a restriction of large deviations from the native state and a corresponding increase in smaller scale fluctuations about the native state. These fluctuations contribute to entropy and hence to the stability of the native state. The largest changes in localized dynamics occur in surface loops, while other regions, particularly the active site residues, remain essentially unchanged. Several mutations, most notably L313F and H322Y in variant 8G8, are in the region showing the largest increase in fluctuations, suggesting that these mutations confer more flexibility to the loops. As a validation of our simulations, the fluctuations of Trp102 were examined in detail, and compared with Trp102 phosphorescence lifetimes that were previously measured. Consistent with expectations from the theory of phosphorescence, an inverse correlation between out-of-plane fluctuations on the picosecond time scale and phosphorescence lifetime was observed.     

Structure of a G-quadruplex-ligand complex

Stabilisation of G-quadruplex structures formed from telomeric DNA, by means of quadruplex-selective ligands, is a means of inhibiting the telomerase enzyme from catalysing the synthesis of telomeric DNA repeats. In order to understand the molecular basis of ligand-quadruplex recognition, the crystal structure has been determined of such a complex, at 1.75A resolution. This complex is between a dimeric antiparallel G-quadruplex formed from the Oxytricha nova telomeric DNA sequence d(GGGGTTTTGGGG), and a di-substituted aminoalkylamido acridine compound. The structure shows that the acridine moiety is bound at one end of the stack of G-quartets, within one of the thymine loops. It is held in place by a combination of stacking interactions and specific hydrogen bonds with thymine bases. The stability of the ligand in this binding site has been confirmed by a 2ns molecular dynamics simulation.