A diabetic retinopathy detection method using an improved pillar K-means algorithm

The paper presents a new approach for medical image segmentation. Exudates are a visible sign of diabetic retinopathy that is the major reason of vision loss in patients with diabetes. If the exudates extend into the macular area, blindness may occur. Automated detection of exudates will assist ophthalmologists in early diagnosis. This segmentation process includes a new mechanism for clustering the elements of high-resolution images in order to improve precision and reduce computation time. The system applies K-means clustering to the image segmentation after getting optimized by Pillar algorithm; pillars are constructed in such a way that they can withstand the pressure. Improved pillar algorithm can optimize the K-means clustering for image segmentation in aspects of precision and computation time. This evaluates the proposed approach for image segmentation by comparing with Kmeans and Fuzzy C-means in a medical image. Using this method, identification of dark spot in the retina becomes easier and the proposed algorithm is applied on diabetic retinal images of all stages to identify hard and soft exudates, where the existing pillar K-means is more appropriate for brain MRI images. This proposed system help the doctors to identify the problem in the early stage and can suggest a better drug for preventing further retinal damage.     

Affinity comparison of different THCA synthase to CBGA using modeling computational approaches

The Δ^9-Tetrahydrocannabinol (THCA) is the primary psychoactive compound of Cannabis Sativa. It is produced by Δ^1- Tetrahydrocannabinolic acid synthase (THCA) which catalyzes the oxidative cyclization of cannabigerolic acid (CBGA) the precursor of the THCA. In this study, we were interested by the three dimensional structure of THCA synthase protein. Generation of models were done by MODELLER v9.11 and homology modeling with Δ1-tetrahydrocannabinolic acid (THCA) synthase X ray structure (PDB code 3VTE) on the basis of sequences retrieved from GenBank. Procheck, Errat, and Verify 3D tools were used to verify the reliability of the six 3D models obtained, the overall quality factor and the Prosa Z-score were also used to check the quality of the six modeled proteins. The RMSDs for C-alpha atoms, main-chain atoms, side-chain atoms and all atoms between the modeled structures and the corresponding template ranged between 0.290 Å-1.252 Å, reflecting the good quality of the obtained models. Our study of the CBGA-THCA synthase docking demonstrated that the active site pocket was successfully recognized using computational approach. The interaction energy of CBGA computed in ‘fiber types’ proteins ranged between -4.1 95 kcal/mol and -5.95 kcal/mol whereas in the ‘drug type’ was about -7.02 kcal/mol to -7.16 kcal/mol, which maybe indicate the important role played by the interaction energy of CBGA in the determination of the THCA level in Cannabis Sativa L. varieties. Finally, we have proposed an experimental design in order to explore the binding energy source of ligand-enzyme in Cannabis Sativa and the production level of the THCA in the absence of any information regarding the correlation between the enzyme affinity and THCA level production. This report opens the doors to more studies predicting the binding site pocket with accuracy from the perspective of the protein affinity and THCA level produced in Cannabis Sativa.     

In vitro alpha-amylase inhibition and in vivo antioxidant potential of Amaranthus spinosus in alloxan-induced oxidative stress in diabetic rats

Amaranthus spinosus Linn. (Amaranthaceae), commonly known as “Mulluharivesoppu” in Kannada, is used in the Indian traditional system of medicine for the treatment of diabetes. The present study deals with the scientific evaluation of alpha amylase and the antioxidant potential of methanol extract of A. spinosus (MEAS). The aim of this study was to investigate in vitro alpha-amylase enzyme inhibition by CNPG3 (2-chloro-4-nitrophenol α-d-maltotrioside) and in vivo antioxidant potential of malondialdehyde (MDA), glutathione (GSH), catalase (CAT) and total thiols (TT) in alloxan-induced diabetic rats of a methanolic extract of A. spinosus. Blood sugar was also determined in MEAS-treated alloxan-induced diabetic rats. MEAS showed significant inhibition of alpha-amylase activity and IC₅₀ 46.02 μg/ml. Oral administration of MEAS (200 and 400 mg/kg) for 15 days showed significant reduction in the elevated blood glucose, MDA and restores GSH, CAT and TT levels as compared with a diabetic control. The present study provides evidence that the methanolic extract of A. spinosus has potent alpha amylase, anti-diabetic and antioxidant activities.     

Atherogenic effect of Arecoline: A computational study

There are over 600 million people worldwide covering Asian and Oceanic countries including India have the habit of chewing areca nut as masticator in different forms. Arecoline (C(8)H(13)NO(2)) has been reported as one of the abundant constituents of areca nut. A good number of scientific publications have made Arecoline responsible for oral cancer. Based on observation from clinical situation in North East India, one of the most betel quid chewing region of the country, we suspected a link between consumption of areca nut and Cerebro Vascular Disease like stroke. Therefore, we considered Low Density Lipoprotein (LDL) receptor as target and Arecoline as ligand and studied ligand -target interaction using computational tools. Also we considered High Density Lipoprotein (HDL) receptor as another target to see if Arecoline has any binding potential with it over and above LDL receptor. Docking result indicated that Arecoline and Cholesterol both, have affinity towards extracellular domain of Human LDL receptor but affinity of Arecoline is much higher (-12.3560.) than that of Cholesterol(-0.1810). Docking of Arecoline and 1, 2-Hexyl-1- cyclopentanone thiosemicarbazone (thiosemicarbazone) with Bovine HDL receptor showed that Arecoline also has the potential (Score, -6.2690Kcal/Mol) to block HDL receptor though its potential is less than that (score, -10.0509 Kcal/Mol) of control (thiosemicarbazone). We, therefore, suggest that by inhibiting endocytosis of LDL cholesterol because of blocking LDL receptor function and also by preventing LDL cholesterol uptake by liver from blood because of interference with HDL receptor, Arecoline may contribute to atherosclerosis. The study therefore, indicates a positive correlation between chewing of betel quid and Cerebro Vascular Disease.     

Assessing computational tools for the discovery of small RNA genes in bacteria

Over the past decade, a number of biocomputational tools have been developed to predict small RNA (sRNA) genes in bacterial genomes. In this study, several of the leading biocomputational tools, which use different methodologies, were investigated. The performance of the tools, both individually and in combination, was evaluated on ten sets of benchmark data, including data from a novel RNA-seq experiment conducted in this study. The results of this study offer insight into the utility as well as the limitations of the leading biocomputational tools for sRNA identification and provide practical guidance for users of the tools.     

Structure modeling and antidiabetic activity of a seed protein of Momordica charantia in non-obese diabetic (NOD) mice

Momordica charantia is a well known medicinal plant used in the traditional medicinal system for the treatment of various diseases including diabetes mellitus. Recently, a novel protein termed as ADMc1 from the seed extract of M. charantia has been identified and isolated showing significant antihyperglycemic activity in type 1 diabetic rats in which diabetes was induced. However, the structure of this protein has not yet been analyzed. Homology modeling approach was used to generate a high quality protein 3D structure for the amino acid sequence of the ADMc1 protein in this study. The comparative assessment of secondary structures revealed ADMc1 as an all-alpha helix protein with random coils. Tertiary structure predicted on the template structure of Napin of B. Napus (PDB ID: 1SM7) with which the ADMc1 showed significant sequence similarity, was validated using protein structure validation tools like PROCHECK, WHAT_CHECK, VERIFY3D and ProSA. Arrangement of disulfide bridges formed by cysteine residues were predicted by the Dianna 1.1 server. The presence of multiple disulfide bond confers the stable nature of the ADMc1 protein. Further, the biological activity of the ADMc1 was assessed in non-obese diabetic (NOD) mice which are spontaneous model of type 1 diabetes. Significant reduction in the blood glucose levels of NOD mice was observed up to 8 h post administration of the rADMc1 protein. Overall, the structural characterizations with antihyperglycemic activity of this seed protein of Momordica charantia demonstrate its potential as an antidiabetic agent.     

A comparative analysis of leaf shape of wheat, barley and maize using an empirical shape model

Background and Aims

The phenotypes of grasses show differences depending on growth conditions and ontogenetic stage. Understanding these responses and finding suitable mathematical formalizations are an essential part of the development of plant and crop models. Usually, a marked change in architecture between juvenile and adult plants is observed, where dimension and shape of leaves are likely to change. In this paper, the plasticity of leaf shape is analysed according to growth conditions and ontogeny.

Methods

Leaf shape of Triticum aestivum, Hordeum vulgare and Zea mays cultivars grown under varying conditions was measured using digital image processing. An empirical leaf shape model was fitted to measured shape data of single leaves. Obtained values of model parameters were used to analyse the patterns in leaf shape.

Key Results

The model was able to delineate leaf shape of all studied species. The model error was small. Differences in leaf shape between juvenile and adult leaves in T. aestivum and H. vulgare were observed. Varying growth conditions impacted leaf dimensions but did not impact leaf shape of the respective species.

Conclusions

Leaf shape of the studied T. aestivum and H. vulgare cultivars was remarkably stable for a comparable ontogenetic stage (leaf rank), but differed between stages. Along with other aspects of grass architecture, leaf shape changed during the transition from juvenile to adult growth phase. Model-based analysis of leaf shape is a method to investigate these differences. Presented results can be integrated into architectural models of plant development to delineate leaf shape for different species, cultivars and environmental conditions.     

Characterization and evaluation of an arbitrary primed Polymerase Chain Reaction (PCR) product for the specific detection of Brucella species

Laboratory detection of Brucella is based largely on bacterial isolation and phenotypic characterization. These methods are lengthy and labor-intensive and have been associated with a heightened risk of laboratory-acquired infection. Antibody based indirect detection methods also suffer from limitations in proper diagnosis of the organism. To overcome these problems, nucleic acid amplification has been explored for rapid detection and confirmation of the presence of Brucella spp. PCR-based diagnostics is useful for screening large populations of livestock to identify infected individuals and confirms the presence of the pathogen. Random Amplification of Polymorphic DNA (RAPD) was performed and identified a 1.3 kb PCR fragment specifically amplifiable from DNA isolated from Brucella. A BLAST search revealed no significant homology with the reported sequences from species other than the members of Brucella. The isolated fragment seems to be a part of d-alanine–d-alanine ligase gene in Brucella sp. Translational BLAST revealed certain degree of homology of this sequence with orthologs of this gene reported from other microbial species at the deduced amino acid level. The sequence information was used to develop PCR based assays to detect Brucella sp. from various samples. The minimum detection limit of Brucella from blood and milk samples spiked with Brucella DNA was found to be 1 ng/ml and 10 ng/ml, respectively. In conclusion, we demonstrated that the PCR based detection protocol was successfully used for the detection of Brucella from various organs and spiked samples of diseased sheep. Diagnosis of Brucellosis by PCR based method reported in this study is relatively rapid, specific and simple.     

Saturation transfer difference NMR and computational modeling of a sialoadhesin-sialyl lactose complex

The siglecs are a family of I-type lectins binding to sialic acids on the cell surface. Sialoadhesin (siglec-1) is expressed at much higher levels in inflammatory macrophages and specifically binds to alpha-2,3-sialylated N-acetyl lactosamine residues of glycan chains. The terminal disaccharide alpha-D-Neu5Ac-(2-->3)-beta-D-Gal is thought to be the main epitope recognized by sialoadhesin. To understand the basis of this biological recognition reaction we combined NMR experiments with a molecular modeling study. We employed saturation transfer difference (STD) NMR experiments to characterize the binding epitope of alpha-2,3-sialylated lactose, alpha-D-Neu5Ac-(2-->3)-beta-D-Gal-(1-->4)-D-Glc 1 to sialoadhesin at atomic resolution. The experimental results were compared to a computational docking model and to X-ray data of a complex of sialyl lactose and sialoadhesin. The data reveal that sialoadhesin mainly recognizes the N-acetyl neuraminic acid and a small part of the galactose moiety of 1. The crystal structure of a complex of sialoadhesin with sialyl lactose 1 was used as a basis for a modeling study using the FlexiDock algorithm. The model generated was very similar to the original crystal structure. Therefore, the X-ray data were used to predict theoretical STD values utilizing the CORCEMA-STD protocol. The good agreement between experimental and theoretical STD values indicates that a combined modeling/STD NMR approach yields a reliable structural model for the complex of sialoadhesin with alpha-D-Neu5Ac-(2-->3)-beta-D-Gal-(1-->4)-D-Glc 1 in aqueous solution.     

Receptor for advanced glycation end-product (RAGE) gene polymorphism 2245G/A is associated with pro-inflammatory,oxidative-glycation markers and sRAGE in diabetic retinopathy

Background

Receptor for advanced glycation end-product (RAGE) gene polymorphism 2245G/A is associated with diabetic retinopathy (DR). However, the mechanism on how it affects the disease development is still unclear.

Aim

This study aims to investigate the relationship between 2245G/A RAGE gene polymorphism and selected pro-inflammatory, oxidative-glycation markers in DR patients.

Methods

A total of 371 unrelated type 2 diabetic patients [200 with retinopathy, 171 without retinopathy (DNR)] and 235 healthy subjects were recruited. Genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism method followed by DNA sequencing. The nuclear and cytosolic extracts from peripheral blood mononuclear cells were used for nuclear factor kappa B (NF-κB) p65 and superoxide dismutase activity measurement respectively. Plasma was used for glutathione peroxidase activity, advanced oxidation protein product (AOPP), monocyte chemoattractant protein (MCP)-1, pentosidine and soluble RAGE (sRAGE) measurements.

Results

DR patients with 2245GA genotype had significantly elevated levels of activated NF-κB p65, plasma MCP-1, AOPP and pentosidine but lower level of sRAGE when compared to DR patients with wild-type 2245GG.

Conclusion

The RAGE gene polymorphism 2245G/A is associated with pro-inflammatory, oxidative-glycation markers and circulating sRAGE in DR patients. Patients with 2245GA RAGE genotype could aggravate DR possibly via NF-κB mediated inflammatory pathway.     

Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments

The advent of rapid complete genome sequencing, and the potential to capture this information in genome-scale metabolic models, provide the possibility of comprehensively modeling microbial community interactions. For example, Rhodoferax and Geobacter species are acetate-oxidizing Fe(III)-reducers that compete in anoxic subsurface environments and this competition may have an influence on the in situ bioremediation of uranium-contaminated groundwater. Therefore, genome-scale models of Geobacter sulfurreducens and Rhodoferax ferrireducens were used to evaluate how Geobacter and Rhodoferax species might compete under diverse conditions found in a uranium-contaminated aquifer in Rifle, CO. The model predicted that at the low rates of acetate flux expected under natural conditions at the site, Rhodoferax will outcompete Geobacter as long as sufficient ammonium is available. The model also predicted that when high concentrations of acetate are added during in situ bioremediation, Geobacter species would predominate, consistent with field-scale observations. This can be attributed to the higher expected growth yields of Rhodoferax and the ability of Geobacter to fix nitrogen. The modeling predicted relative proportions of Geobacter and Rhodoferax in geochemically distinct zones of the Rifle site that were comparable to those that were previously documented with molecular techniques. The model also predicted that under nitrogen fixation, higher carbon and electron fluxes would be diverted toward respiration rather than biomass formation in Geobacter, providing a potential explanation for enhanced in situ U(VI) reduction in low-ammonium zones. These results show that genome-scale modeling can be a useful tool for predicting microbial interactions in subsurface environments and shows promise for designing bioremediation strategies.     

PCR detection of Anaplasma phagocytophilum in goat flocks in an area endemic for tick-borne fever in Switzerland

Central Switzerland is a highly endemic region for tick-borne fever (TBF) in cattle, however, little is known about A. phagocytophilum in goats. In the present study, 72 animals from six goat flocks (373 EDTA blood-samples) in Central Switzerland were analysed for A. phagocytophilum DNA. A real-time PCR targeting the msp2 gene of A. phagocytophilum was performed and in positive samples the partial 165 rRNA, groEL and msp4 gene were amplified for sequence analysis. Four DNA extracts were positive. Different sequence types on basis of the amplified genes were found. For comparison, sequences of A. phagocytophilum from 12 cattle (originating from Switzerland and Southern Germany) were analysed. The 165 rRNA gene sequences from cattle were all identical amongst each other, but the groEL and msp4 gene differed depending on the origin of the cattle samples and differed from the variants from goats. This study clearly provides molecular evidence for the presence of different types of A. phagocytophilum in goat flocks in Switzerland, a fact which deserves more thorough attention in clinical studies.     

A molecular modeling based screening for potential inhibitors to alpha hemolysin from Staphylococcus aureus

Staphylococcus aureus, a Gram-positive bacterium is pathogenic in nature. It is known that secreted toxins remain active after antibiotic treatment. The alpha hemolysin or alpha toxin damages cell membrane and induces apoptosis and degradation of DNA. The titer of alphahemolysin increases and causes hemostasis disturbances, thrombocytopenia, and pulmonary lesions during staphylococcal infection. Therefore, it is of interest to inhibit alpha hemolysin using novel compounds. We used the structure of alpha hemolysin(PDB: 7AHL) to screen structures for 100,000 compounds from the ZINC database using molecular docking with AutoDock VINA. Nine (9) successive hits were then subjected for pharmacokinetic and toxicity properties by PROTOX (a webserver for the prediction of oral toxicities of small molecules) and FAFDrugs (a tool for prediction of ADME and Toxicity). This exercise further identified hit #1 ({[3a-(Dihydroxymethyl)-6-hydroxy-2,2-dimethyl-1,3,4-trioxatetrahydro-2H-pentalen-5- yl]methyl}amino(9H-fluoren-9-yl)acetate with binding affinity: -10.3 kcal/mol) and hit #2 (6-(Dihydroxymethyl)-2-{2-[3- (methylamino)propyl]-2-azatricyclo[9.4.0.03,8]pentadeca-1(11),3,5,7,12,14-hexaen-6-yloxy}tetrahydro-2H-pyran-3,4,5-triol with binding affinity: -9.6 kcal/mol) with acceptable toxicity and ADME properties for potential predicted hemolysin inhibition. These compounds should then be evaluated in vitro using inhibitory studies.     

Choke point analysis of metabolic pathways in E.histolytica: a computational approach for drug target identification

With the Entamoeba genome essentially complete, the organism can be studied from a whole genome standpoint. The understanding of cellular mechanisms and interactions between cellular components is instrumental to the development of new effective drugs and vaccines. Metabolic pathway analysis is becoming increasingly important for assessing inherent network properties in reconstructed biochemical reaction networks. Metabolic pathways illustrate how proteins work in concert to produce cellular compounds or to transmit information at different levels. Identification of drug targets in E. histolytica through metabolic pathway analysis promises to be a novel approach in this direction. This article focuses on the identification of drug targets by subjecting the Entamoeba genome to BLAST with the e-value inclusion threshold set to 0.005 and choke point analysis. A total of 86.9 percent of proposed drug targets with biological evidence are chokepoint reactions in Entamoeba genome database.     

Real-time PCR versus conventional PCR for malaria parasite detection in low-grade parasitemia

We have optimized a faster and cheaper real-time PCR and developed a conventional genus specific PCR based on 18S rRNA gene to detect malaria parasites in low-grade parasitemias. Additionally, we compared these PCRs to the OptiMAL-IT test. Since there is no consensus on choice of standard quantitative curve in real-time assays, we decided to investigate the performance of parasite DNA from three different sources: "genome", amplicon and plasmid. The amplicon curve showed the best efficiency in quantifying parasites. Both PCR assays detected 100% of the clinical samples tested; the sensitivity threshold was 0.5 parasite/mul and no PCR positive reaction occurred when malaria parasites were not present. Conversely, if OptiMAL-IT were employed for malaria diagnosis, 30% of false-negative results could be expected. We conclude that PCR assays have potential for detecting malaria parasites in asymptomatic infections, in evaluation of malaria vaccine molecule candidates, for screening blood donors, especially in endemic areas, or even in monitoring malaria therapy.     

Evaluation and optimisation of indel detection workflows for ion torrent sequencing of the BRCA1 and BRCA2 genes

Background

The Ion Torrent PGM is a popular benchtop sequencer that shows promise in replacing conventional Sanger sequencing as the gold standard for mutation detection. Despite the PGM’s reported high accuracy in calling single nucleotide variations, it tends to generate many false positive calls in detecting insertions and deletions (indels), which may hinder its utility for clinical genetic testing.

Results

Recently, the proprietary analytical workflow for the Ion Torrent sequencer, Torrent Suite (TS), underwent a series of upgrades. We evaluated three major upgrades of TS by calling indels in the BRCA1 and BRCA2 genes. Our analysis revealed that false negative indels could be generated by TS under both default calling parameters and parameters adjusted for maximum sensitivity. However, indel calling with the same data using the open source variant callers, GATK and SAMtools showed that false negatives could be minimised with the use of appropriate bioinformatics analysis. Furthermore, we identified two variant calling measures, Quality-by-Depth (QD) and VARiation of the Width of gaps and inserts (VARW), which substantially reduced false positive indels, including non-homopolymer associated errors without compromising sensitivity. In our best case scenario that involved the TMAP aligner and SAMtools, we achieved 100% sensitivity, 99.99% specificity and 29% False Discovery Rate (FDR) in indel calling from all 23 samples, which is a good performance for mutation screening using PGM.

Conclusions

New versions of TS, BWA and GATK have shown improvements in indel calling sensitivity and specificity over their older counterpart. However, the variant caller of TS exhibits a lower sensitivity than GATK and SAMtools. Our findings demonstrate that although indel calling from PGM sequences may appear to be noisy at first glance, proper computational indel calling analysis is able to maximize both the sensitivity and specificity at the single base level, paving the way for the usage of this technology for future clinical genetic testing.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-516) contains supplementary material, which is available to authorized users.     

A quantitative Real Time PCR based method for the detection of Phytophthora infestans causing Late blight of potato,in infested soil

A fast and simple polymerase chain reaction method has been developed for detection of Phytophthora infestans oospores, the causal agent of Late blight of Potato in soil. The method involves the disruption of oospores by grinding dry soil, using abrasive properties, in the presence of glass powder and skimmed milk powder within a short time. The latter prevents loss of DNA by adsorption to soil particles or by degradation and reduces the co-extraction of PCR inhibitors with the DNA. After phenol/chloroform extraction; the DNA is suitable for direct PCR amplification without a precipitation step. This amplification leads to detection of pathogen in infested soils before planting of crop. The real-time PCR assay we describe is highly sensitive and specific, and has several advantages over conventional PCR assays used for P. infestans detection to confirm positive inoculum level in potato seeds and elsewhere. With increasing amounts of standard DNA templates, the respective threshold cycle (Ct) values were determined and a linear relationship was established between these Ct values and the logarithm of initial template amounts. The method is rapid, cost efficient, and when combined with suitable internal controls can be applied to the detection and quantification of P. infestans oospores on a large-scale basis.     

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K_d values that vary between 10^− 4 and 10^− 14 M, characteristics that are explained by a ‘dual-recognition’ mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential ΔΔGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental ΔΔG data, particularly for interactions not mediated by interfacial water molecules. ΔΔG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.     

Nocturnal insects use optic flow for flight control

To avoid collisions when navigating through cluttered environments, flying insects must control their flight so that their sensory systems have time to detect obstacles and avoid them. To do this, day-active insects rely primarily on the pattern of apparent motion generated on the retina during flight (optic flow). However, many flying insects are active at night, when obtaining reliable visual information for flight control presents much more of a challenge. To assess whether nocturnal flying insects also rely on optic flow cues to control flight in dim light, we recorded flights of the nocturnal neotropical sweat bee, Megalopta genalis, flying along an experimental tunnel when: (i) the visual texture on each wall generated strong horizontal (front-to-back) optic flow cues, (ii) the texture on only one wall generated these cues, and (iii) horizontal optic flow cues were removed from both walls. We find that Megalopta increase their groundspeed when horizontal motion cues in the tunnel are reduced (conditions (ii) and (iii)). However, differences in the amount of horizontal optic flow on each wall of the tunnel (condition (ii)) do not affect the centred position of the bee within the flight tunnel. To better understand the behavioural response of Megalopta, we repeated the experiments on day-active bumble-bees (Bombus terrestris). Overall, our findings demonstrate that despite the limitations imposed by dim light, Megalopta-like their day-active relatives-rely heavily on vision to control flight, but that they use visual cues in a different manner from diurnal insects.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.