Intelligent resolution: Integrating Cryo-EM with AI-driven multi-resolution simulations to observe the severe acute respiratory syndrome coronavirus-2 replication-transcription machinery in action

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) replication transcription complex (RTC) is a multi-domain protein responsible for replicating and transcribing the viral mRNA inside a human cell. Attacking RTC function with pharmaceutical compounds is a pathway to treating COVID-19. Conventional tools, e.g., cryo-electron microscopy and all-atom molecular dynamics (AAMD), do not provide sufficiently high resolution or timescale to capture important dynamics of this molecular machine. Consequently, we develop an innovative workflow that bridges the gap between these resolutions, using mesoscale fluctuating finite element analysis (FFEA) continuum simulations and a hierarchy of AI-methods that continually learn and infer features for maintaining consistency between AAMD and FFEA simulations. We leverage a multi-site distributed workflow manager to orchestrate AI, FFEA, and AAMD jobs, providing optimal resource utilization across HPC centers. Our study provides unprecedented access to study the SARS-CoV-2 RTC machinery, while providing general capability for AI-enabled multi-resolution simulations at scale.

Do Sideline Tests of Vestibular and Oculomotor Function Accurately Diagnose Sports-Related Concussion in Adults? A Systematic Review and Meta-analysis

BACKGROUND

Sports-related concussion (SRC) assessment tools are primarily based on subjective assessments of somatic, cognitive, and psychosocial/emotional symptoms. SRC symptoms remain underreported, and objective measures of SRC impairments would be valuable to assist diagnosis. Measurable impairments to vestibular and oculomotor processing have been shown to occur after SRC and may provide valid objective assessments.

PURPOSE

Determine the diagnostic accuracy of sideline tests of vestibular and oculomotor dysfunction to identify SRC in adults.

STUDY DESIGN

Systematic review; Level of evidence, 4.

METHODS

Electronic databases and gray literature were searched from inception until February 12, 2020. Physically active individuals (>16 years of age) who participated in sports were included. The reference standard for SRC was a combination of clinical signs and symptoms (eg, the Sport Concussion Assessment Tool [SCAT]), and index tests included any oculomotor assessment tool. The QUADAS tool was used to assess risk of bias, with the credibility of the evidence being rated according to GRADE.

RESULTS

A total of 8 studies were included in this review. All included studies used the King-Devick test, with no other measures being identified. Meta-analysis was performed on 4 studies with a summary sensitivity and specificity of 0.77 and 0.82, respectively. The overall credibility of the evidence was rated as very low.

CONCLUSION

Caution must be taken when interpreting these results given the very low credibility of the evidence, and the true summary sensitivity and specificity may substantially differ from the values calculated within this systematic review. Therefore, we recommend that clinicians using the King-Devick test to diagnose SRC in adults do so in conjunction with other tools such as the SCAT.

PROSPERO REGISTRATION

CRD42018106632.

The C-terminal head domain of Burkholderia pseudomallei BpaC has a striking hydrophilic core with an extensive solvent network

Gram‐negative pathogens like Burkholderia pseudomallei use trimeric autotransporter adhesins such as BpaC as key molecules in their pathogenicity. Our 1.4 Å crystal structure of the membrane‐proximal part of the BpaC head domain shows that the domain is exclusively made of left‐handed parallel β‐roll repeats. This, the largest such structure solved, has two unique features. First, the core, rather than being composed of the canonical hydrophobic Ile and Val, is made up primarily of the hydrophilic Thr and Asn, with two different solvent channels. Second, comparing BpaC to all other left‐handed parallel β‐roll structures showed that the position of the head domain in the protein correlates with the number and type of charged residues. In BpaC, only negatively charged residues face the solvent—in stark contrast to the primarily positive surface charge of the left‐handed parallel β‐roll “type” protein, YadA. We propose extending the definitions of these head domains to include the BpaC‐like head domain as a separate subtype, based on its unusual sequence, position, and charge. We speculate that the function of left‐handed parallel β‐roll structures may differ depending on their position in the structure.

Modelling the active SARS-CoV-2 helicase complex as a basis for structure-based inhibitor design

The RNA helicase (non-structural protein 13, NSP13) of SARS-CoV-2 is essential for viral replication, and it is highly conserved among the coronaviridae family, thus a prominent drug target to treat COVID-19. We present here structural models and dynamics of the helicase in complex with its native substrates based on thorough analysis of homologous sequences and existing experimental structures. We performed and analysed microseconds of molecular dynamics (MD) simulations, and our model provides valuable insights to the binding of the ATP and ssRNA at the atomic level. We identify the principal motions characterising the enzyme and highlight the effect of the natural substrates on this dynamics. Furthermore, allosteric binding sites are suggested by our pocket analysis. Our obtained structural and dynamical insights are important for subsequent studies of the catalytic function and for the development of specific inhibitors at our characterised binding pockets for this promising COVID-19 drug target.

Spectrum of pathogenic variants and founder effects in amelogenesis imperfecta associated with MMP20

Amelogenesis imperfecta (AI) describes a heterogeneous group of developmental enamel defects that typically have Mendelian inheritance. Exome sequencing of 10 families with recessive hypomaturation AI revealed four novel and one known variants in the matrix metallopeptidase 20 (MMP20) gene that were predicted to be pathogenic. MMP20 encodes a protease that cleaves the developing extracellular enamel matrix and is necessary for normal enamel crystal growth during amelogenesis. New homozygous missense changes were shared between four families of Pakistani heritage (c.625G>C; p.(Glu209Gln)) and two of Omani origin (c.710C>A; p.(Ser237Tyr)). In two families of UK origin and one from Costa Rica, affected individuals were homozygous for the previously reported c.954-2A>T; p.(Ile319Phefs*19) variant. For each of these variants, microsatellite haplotypes appeared to exclude a recent founder effect, but elements of haplotype were conserved, suggesting more distant founding ancestors. New compound heterozygous changes were identified in one family of the European heritage: c.809_811+12delinsCCAG; p.(?) and c.1122A>C; p.(Gln374His). This report further elucidates the mutation spectrum of MMP20 and the probable impact on protein function, confirms a consistent hypomaturation phenotype and shows that mutations in MMP20 are a common cause of autosomal recessive AI in some communities.

Base-pair resolution analysis of the effect of supercoiling on DNA flexibility and major groove recognition by triplex-forming oligonucleotides

In the cell, DNA is arranged into highly-organised and topologically-constrained (supercoiled) structures. It remains unclear how this supercoiling affects the detailed double-helical structure of DNA, largely because of limitations in spatial resolution of the available biophysical tools. Here, we overcome these limitations, by a combination of atomic force microscopy (AFM) and atomistic molecular dynamics (MD) simulations, to resolve structures of negatively-supercoiled DNA minicircles at base-pair resolution. We observe that negative superhelical stress induces local variation in the canonical B-form DNA structure by introducing kinks and defects that affect global minicircle structure and flexibility. We probe how these local and global conformational changes affect DNA interactions through the binding of triplex-forming oligonucleotides to DNA minicircles. We show that the energetics of triplex formation is governed by a delicate balance between electrostatics and bonding interactions. Our results provide mechanistic insight into how DNA supercoiling can affect molecular recognition, that may have broader implications for DNA interactions with other molecular species.

New variants and in silico analyses in GRK1 associated Oguchi disease

Biallelic mutations in G-Protein coupled receptor kinase 1 (GRK1) cause Oguchi disease, a rare subtype of congenital stationary night blindness (CSNB). The purpose of this study was to identify disease causing GRK1 variants and use in-depth bioinformatic analyses to evaluate how their impact on protein structure could lead to pathogenicity. Patients' genomic DNA was sequenced by whole genome, whole exome or focused exome sequencing. Disease associated variants, published and novel, were compared to nondisease associated missense variants. The impact of GRK1 missense variants at the protein level were then predicted using a series of computational tools. We identified twelve previously unpublished cases with biallelic disease associated GRK1 variants, including eight novel variants, and reviewed all GRK1 disease associated variants. Further structure-based scoring revealed a hotspot for missense variants in the kinase domain. In addition, to aid future clinical interpretation, we identified the bioinformatics tools best able to differentiate disease associated from nondisease associated variants. We identified GRK1 variants in Oguchi disease patients and investigated how disease-causing variants may impede protein function in-silico.

Combined Force-Torque Spectroscopy of Proteins by Means of Multiscale Molecular Simulation

Assessing the structural properties of large proteins is important to gain an understanding of their function in, e.g., biological systems or biomedical applications. We propose a method to examine the mechanical properties of proteins subject to applied forces by means of multiscale simulation. Both stretching and torsional forces are considered, and these may be applied independently of each other. As a proof of principle, we apply torsional forces to a coarse-grained continuum model of the antibody protein immunoglobulin G using fluctuating finite element analysis and use it to identify the area of strongest deformation. This region is essential to the torsional properties of the molecule as a whole because it represents the softest, most deformable domain. Zooming in, this part of the molecule is subjected to torques and stretching forces using molecular dynamics simulations on an atomistically resolved level to investigate its torsional properties. We calculate the torsional resistance as a function of the rotation of the domain while subjecting it to various stretching forces. From this, we assess how the measured twist-torque profiles develop with increasing stretching force and show that they exhibit torsion stiffening, in qualitative agreement with experimental findings. We argue that combining the twist-torque profiles for various stretching forces effectively results in a combined force-torque spectroscopy analysis, which may serve as a mechanical signature for a biological macromolecule.

KOBRA: a fluctuating elastic rod model for slender biological macromolecules

KOBRA (KirchOff Biological Rod Algorithm) is an algorithm and software package designed to perform dynamical simulations of elongated biomolecules such as those containing alpha-helices and coiled-coils. It represents these as coarsely-discretised Kirchoff rods, with linear elements that can stretch, bend and twist independently. These rods can have anisotropic and inhomogeneous parameters and bent or twisted equilibrium structures, allowing for a coarse-grained parameterisation of complex biological structures. Each element is non-inertial and subject to thermal fluctuations. The speed and simplicity of the algorithm allows KOBRA rods to easily access timescales from nanoseconds to seconds. To demonstrate this functionality, a KOBRA rod was parameterised using data from all-atom simulations of the Ndc80 protein complex, and compared against these simulations and negative-stain EM images. The distribution of bend angles and principal components were highly correlated between KOBRA, all-atom molecular dynamics, and experimental data. The properties of a hinge region, thought to be found at an unstructured loop, were studied. A C++ implementation of KOBRA is available under the GNU GPLv3 free software licence.

Factors affecting the species of Campylobacter colonizing chickens reared for meat

AIM

To investigate factors influencing Campylobacter spp. colonization of broiler chickens.

METHODS AND RESULTS

Campylobacters were isolated from caeca from 319 flocks of two different breeds (199 Cobb and 120 Hubbard), reared as standard (199), Freedom Food/corn fed (57), free-range (47) or organic (16). The standard category exclusively used Cobb birds slaughtered at 38-41 days. The Freedom Food/corn-fed and free-range Hubbard birds were slaughtered at 49-56 days and the organic flocks at 70 days. Campylobacters were picked at random from direct plates. Both breed of chicken (Hubbard) and age at slaughter were independently associated with increased likelihood of colonization by Campylobacter coli rather than Campylobacter jejuni, but breed could not be separated from other aspects of husbandry with the data available.

CONCLUSIONS

Chickens are frequently colonized by C. jejuni and C. coli and most human infections originate from poultry. In most developed countries approximately 90% of human infections are caused by C. jejuni, but fewer than 10% by C. coli. This might be due to C. coli being less pathogenic than C. jejuni to humans, and/or to chicken meat carrying fewer C. coli than C. jejuni. More investigations are needed into these aspects before it can be concluded that slaughtering older birds from slower-growing breeds would reduce the risk of human Campylobacter disease.

SIGNIFICANCE AND IMPACT OF THE STUDY

Meat from certain breeds of poultry are predominantly colonized by C. coli rather than C. jejuni. More research is needed to understand the impact this may have on the number and severity of human campylobacter infections.

Continuum mechanical parameterisation of cytoplasmic dynein from atomistic simulation

Cytoplasmic dynein is responsible for intra-cellular transport in eukaryotic cells. Using Fluctuating Finite Element Analysis (FFEA), a novel algorithm that represents proteins as continuum viscoelastic solids subject to thermal noise, we are building computational tools to study the mechanics of these molecular machines. Here we present a methodology for obtaining the material parameters required to represent the flexibility of cytoplasmic dynein within FFEA from atomistic molecular dynamics (MD) simulations, and show that this continuum representation is sufficient to capture the principal dynamic properties of the motor.

Nano-encapsulated Escherichia coli Divisome Anchor ZipA, and in Complex with FtsZ

The E. coli membrane protein ZipA, binds to the tubulin homologue FtsZ, in the early stage of cell division. We isolated ZipA in a Styrene Maleic Acid lipid particle (SMALP) preserving its position and integrity with native E. coli membrane lipids. Direct binding of ZipA to FtsZ is demonstrated, including FtsZ fibre bundles decorated with ZipA. Using Cryo-Electron Microscopy, small-angle X-ray and neutron scattering, we determine the encapsulated-ZipA structure in isolation, and in complex with FtsZ to a resolution of 1.6 nm. Three regions can be identified from the structure which correspond to, SMALP encapsulated membrane and ZipA transmembrane helix, a separate short compact tether, and ZipA globular head which binds FtsZ. The complex extends 12 nm from the membrane in a compact structure, supported by mesoscale modelling techniques, measuring the movement and stiffness of the regions within ZipA provides molecular scale analysis and visualisation of the early divisome.

Securing the future of research computing in the biosciences

Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data.

Fluctuating Finite Element Analysis (FFEA): A continuum mechanics software tool for mesoscale simulation of biomolecules

Fluctuating Finite Element Analysis (FFEA) is a software package designed to perform continuum mechanics simulations of proteins and other globular macromolecules. It combines conventional finite element methods with stochastic thermal noise, and is appropriate for simulations of large proteins and protein complexes at the mesoscale (length-scales in the range of 5 nm to 1 μm), where there is currently a paucity of modelling tools. It requires 3D volumetric information as input, which can be low resolution structural information such as cryo-electron tomography (cryo-ET) maps or much higher resolution atomistic co-ordinates from which volumetric information can be extracted. In this article we introduce our open source software package for performing FFEA simulations which we have released under a GPLv3 license. The software package includes a C ++ implementation of FFEA, together with tools to assist the user to set up the system from Electron Microscopy Data Bank (EMDB) or Protein Data Bank (PDB) data files. We also provide a PyMOL plugin to perform basic visualisation and additional Python tools for the analysis of FFEA simulation trajectories. This manuscript provides a basic background to the FFEA method, describing the implementation of the core mechanical model and how intermolecular interactions and the solvent environment are included within this framework. We provide prospective FFEA users with a practical overview of how to set up an FFEA simulation with reference to our publicly available online tutorials and manuals that accompany this first release of the package.

Defining Dynamics of Membrane-Bound Pyrophosphatases by Experimental and Computational Single-Molecule FRET

Membrane-bound pyrophosphatases couple the hydrolysis of inorganic pyrophosphate to the pumping of ions (sodium or protons) across a membrane in order to generate an electrochemical gradient. This class of membrane protein is widely conserved across plants, fungi, archaea, and bacteria, but absent in multicellular animals, making them a viable target for drug design against protozoan parasites such as Plasmodium falciparum. An excellent understanding of many of the catalytic states throughout the enzymatic cycle has already been afforded by crystallography. However, the dynamics and kinetics of the catalytic cycle between these static snapshots remain to be elucidated. Here, we employ single-molecule Förster resonance energy transfer (FRET) measurements to determine the dynamic range and frequency of conformations available to the enzyme in a lipid bilayer during the catalytic cycle. First, we explore issues related to the introduction of fluorescent dyes by cysteine mutagenesis; we discuss the importance of residue selection for dye attachment, and the balance between mutating areas of the protein that will provide useful dynamics while not altering highly conserved residues that could disrupt protein function. To complement and guide the experiments, we used all-atom molecular dynamics simulations and computational methods to estimate FRET efficiency distributions for dye pairs at different sites in different protein conformational states. We present preliminary single-molecule FRET data that points to insights about the binding modes of different membrane-bound pyrophosphatase substrates and inhibitors.

Peptide strand length controls the energetics of self-assembly and morphology of β-sheet fibrils

Self-assembling peptides can be used as versatile, natural, and multifunctional building blocks to produce a variety of well-defined nanostructures, materials and devices for applications in medicine and nanotechnology. Here, we concentrate on the 1D self-assembly of de novo designed Px-2 peptide β-strands into anti-parallel β-sheet tapes and higher order aggregates. We study six members of the Px-2 family, ranging from 3 amino acids (aa) to 13 aa in length, using a range of complementary experimental techniques, computer simulation and theoretical statistical mechanics. The critical concentration for self-assembly (c*) is found to increase systematically with decreasing peptide length. The shortest peptide found to self-assemble into soluble β-tapes in water is a 5 amino acid residue peptide. These investigations help decipher the role of the peptide length in controlling self-assembly, aggregate morphology, and material properties. By extracting free energies from these data using a statistical mechanical analysis and combining the results with computer simulations at the atomistic level, we can extract the entropy of association for individual β-strands.

Interference between Triplex and Protein Binding to Distal Sites on Supercoiled DNA

We have explored the interdependence of the binding of a DNA triplex and a repressor protein to distal recognition sites on supercoiled DNA minicircles using MD simulations. We observe that the interaction between the two ligands through their influence on their DNA template is determined by a subtle interplay of DNA mechanics and electrostatics, that the changes in flexibility induced by ligand binding play an important role and that supercoiling can instigate additional ligand-DNA contacts that would not be possible in simple linear DNA sequences.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks determine the levels of transcription of a given gene at a given time. Using three complementary computational techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are constrained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any requirement for external signalling factors. Our models combine atomistic molecular dynamics (MD) with coarse-grained simulations and statistical mechanical calculations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denaturation in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informational transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles.

Skp is a multivalent chaperone of outer-membrane proteins

The trimeric chaperone Skp sequesters outer-membrane proteins (OMPs) within a hydrophobic cage, thereby preventing their aggregation during transport across the periplasm in Gram-negative bacteria. Here, we studied the interaction between Escherichia coli Skp and five OMPs of varying size. Investigations of the kinetics of OMP folding revealed that higher Skp/OMP ratios are required to prevent the folding of 16-stranded OMPs compared with their 8-stranded counterparts. Ion mobility spectrometry-mass spectrometry (IMS-MS) data, computer modeling and molecular dynamics simulations provided evidence that 10- to 16-stranded OMPs are encapsulated within an expanded Skp substrate cage. For OMPs that cannot be fully accommodated in the expanded cavity, sequestration is achieved by binding of an additional Skp trimer. The results suggest a new mechanism for Skp chaperone activity involving the coordination of multiple copies of Skp in protecting a single substrate from aggregation.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Flexibility of KorA, a plasmid-encoded, global transcription regulator, in the presence and the absence of its operator

The IncP (Incompatibility group P) plasmids are important carriers in the spread of antibiotic resistance across Gram-negative bacteria. Gene expression in the IncP-1 plasmids is stringently controlled by a network of four global repressors, KorA, KorB, TrbA and KorC interacting cooperatively. Intriguingly, KorA and KorB can act as co-repressors at varying distances between their operators, even when they are moved to be on opposite sides of the DNA. KorA is a homodimer with the 101-amino acid subunits, folding into an N-terminal DNA-binding domain and a C-terminal dimerization domain. In this study, we have determined the structures of the free KorA repressor and two complexes each bound to a 20-bp palindromic DNA duplex containing its consensus operator sequence. Using a combination of X-ray crystallography, nuclear magnetic resonance spectroscopy, SAXS and molecular dynamics calculations, we show that the linker between the two domains is very flexible and the protein remains highly mobile in the presence of DNA. This flexibility allows the DNA-binding domains of the dimer to straddle the operator DNA on binding and is likely to be important in cooperative binding to KorB. Unexpectedly, the C-terminal domain of KorA is structurally similar to the dimerization domain of the tumour suppressor p53.

Contribution of fluorophore dynamics and solvation to resonant energy transfer in protein-DNA complexes: a molecular-dynamics study

In Förster resonance energy transfer (FRET) experiments, extracting accurate structural information about macromolecules depends on knowing the positions and orientations of donor and acceptor fluorophores. Several approaches have been employed to reduce uncertainties in quantitative FRET distance measurements. Fluorophore-position distributions can be estimated by surface accessibility (SA) calculations, which compute the region of space explored by the fluorophore within a static macromolecular structure. However, SA models generally do not take fluorophore shape, dye transition-moment orientation, or dye-specific chemical interactions into account. We present a detailed molecular-dynamics (MD) treatment of fluorophore dynamics for an ATTO donor/acceptor dye pair and specifically consider as case studies dye-labeled protein-DNA intermediates in Cre site-specific recombination. We carried out MD simulations in both an aqueous solution and glycerol/water mixtures to assess the effects of experimental solvent systems on dye dynamics. Our results unequivocally show that MD simulations capture solvent effects and dye-dye interactions that can dramatically affect energy transfer efficiency. We also show that results from SA models and MD simulations strongly diverge in cases where donor and acceptor fluorophores are in close proximity. Although atomistic simulations are computationally more expensive than SA models, explicit MD studies are likely to give more realistic results in both homogeneous and mixed solvents. Our study underscores the model-dependent nature of FRET analyses, but also provides a starting point to develop more realistic in silico approaches for obtaining experimental ensemble and single-molecule FRET data.

Remnant Pachira quinata pasture trees have greater opportunities to self and suffer reduced reproductive success due to inbreeding depression

Habitat fragmentation is extensive throughout the world, converting natural ecosystems into fragments of varying size, density and connectivity. The potential value of remnant trees in agricultural landscapes as seed sources and in connecting fragments has formed a fertile area of debate. This study contrasted the mating patterns of bat-pollinated Pachira quinata trees in a continuous forest to those in pasture through microsatellite-based paternity analysis of progeny. The breeding system was determined by analysis of pollen tube growth and seed production from controlled pollinations. Fitness of selfed and outcrossed seed was compared by germination and seedling growth. There was more inbreeding within pasture trees (outcrossing=0.828±0.015) compared with forest trees (0.926±0.005). Pasture trees had fewer sires contributing to mating events, but pollen dispersal distances were greater than those in the forest. Paternity analysis showed variation in outcrossing rates among pasture trees with high proportions of external and self pollen sources detected. A leaky self-incompatibility system was found, with self pollen having reduced germination on stigmas and slower growth rate through the style. Controlled pollinations also showed a varied ability to self among trees, which was reflected in the selfing rates among pasture trees shown by the paternity analysis (0-80% selfing). Self pollination resulted in lower seed set, germination and seedling growth compared with outcrossing. While remnant trees in agricultural landscapes are involved in broader mating patterns, they show increased but varied levels of inbreeding, which result in reduced fitness.

orchard: Paternity program for autotetraploid species

Advances in molecular marker technology have provided new opportunities to study the population genetics of polyploid taxa. Paternity analysis using microsatellite markers can be used in detection of gene flow between individuals and populations, in mating system analysis, to identify factors that influence fecundity and fertility, to identify behaviour of parent-offspring relationships and in the analysis of the reproductive success of different ecological groups. As there is no specific program for carrying out paternity analysis in tetraploid species, specialized software was designed for the assignment of paternity for autotetraploid species. orchard is a novel implementation of exclusion and likelihood statistics for carrying out paternity analysis of autotetraploids. First, the program performs an exclusion method, and then, a likelihood statistic is used with nonexcluded candidate fathers. Optional features include estimation of allele dosage of known mother trees and the estimation of pollen flow distances. orchard was tested using a data set of microsatellite data of Dipteryx odorata, a tetraploid Amazonian tree species.

In pursuit of an accurate spatial and temporal model of biomolecules at the atomistic level: a perspective on computer simulation

Despite huge advances in the computational techniques available for simulating biomolecules at the quantum-mechanical, atomistic and coarse-grained levels, there is still a widespread perception amongst the experimental community that these calculations are highly specialist and are not generally applicable by researchers outside the theoretical community. In this article, the successes and limitations of biomolecular simulation and the further developments that are likely in the near future are discussed. A brief overview is also provided of the experimental biophysical methods that are commonly used to probe biomolecular structure and dynamics, and the accuracy of the information that can be obtained from each is compared with that from modelling. It is concluded that progress towards an accurate spatial and temporal model of biomacromolecules requires a combination of all of these biophysical techniques, both experimental and computational.

Understanding the apparent stator-rotor connections in the rotary ATPase family using coarse-grained computer modeling

Advances in structural biology, such as cryo-electron microscopy (cryo-EM) have allowed for a number of sophisticated protein complexes to be characterized. However, often only a static snapshot of a protein complex is visualized despite the fact that conformational change is frequently inherent to biological function, as is the case for molecular motors. Computer simulations provide valuable insights into the different conformations available to a particular system that are not accessible using conventional structural techniques. For larger proteins and protein complexes, where a fully atomistic description would be computationally prohibitive, coarse-grained simulation techniques such as Elastic Network Modeling (ENM) are often employed, whereby each atom or group of atoms is linked by a set of springs whose properties can be customized according to the system of interest. Here we compare ENM with a recently proposed continuum model known as Fluctuating Finite Element Analysis (FFEA), which represents the biomolecule as a viscoelastic solid subject to thermal fluctuations. These two complementary computational techniques are used to answer a critical question in the rotary ATPase family; implicit within these motors is the need for a rotor axle and proton pump to rotate freely of the motor domain and stator structures. However, current single particle cryo-EM reconstructions have shown an apparent connection between the stators and rotor axle or pump region, hindering rotation. Both modeling approaches show a possible role for this connection and how it would significantly constrain the mobility of the rotary ATPase family.

Long-term impacts of selective logging on two Amazonian tree species with contrasting ecological and reproductive characteristics: inferences from Eco-gene model simulations

The impact of logging and subsequent recovery after logging is predicted to vary depending on specific life history traits of the logged species. The Eco-gene simulation model was used to evaluate the long-term impacts of selective logging over 300 years on two contrasting Brazilian Amazon tree species, Dipteryx odorata and Jacaranda copaia. D. odorata (Leguminosae), a slow growing climax tree, occurs at very low densities, whereas J. copaia (Bignoniaceae) is a fast growing pioneer tree that occurs at high densities. Microsatellite multilocus genotypes of the pre-logging populations were used as data inputs for the Eco-gene model and post-logging genetic data was used to verify the output from the simulations. Overall, under current Brazilian forest management regulations, there were neither short nor long-term impacts on J. copaia. By contrast, D. odorata cannot be sustainably logged under current regulations, a sustainable scenario was achieved by increasing the minimum cutting diameter at breast height from 50 to 100 cm over 30-year logging cycles. Genetic parameters were only slightly affected by selective logging, with reductions in the numbers of alleles and single genotypes. In the short term, the loss of alleles seen in J. copaia simulations was the same as in real data, whereas fewer alleles were lost in D. odorata simulations than in the field. The different impacts and periods of recovery for each species support the idea that ecological and genetic information are essential at species, ecological guild or reproductive group levels to help derive sustainable management scenarios for tropical forests.

Insecticide resistance in house flies from the United States: resistance levels and frequency of pyrethroid resistance alleles

Although insecticide resistance is a widespread problem for most insect pests, frequently the assessment of resistance occurs over a limited geographic range. Herein, we report the first widespread survey of insecticide resistance in the USA ever undertaken for the house fly, Musca domestica, a major pest in animal production facilities. The levels of resistance to six different insecticides were determined (using discriminating concentration bioassays) in 10 collections of house flies from dairies in nine different states. In addition, the frequencies of Vssc and CYP6D1 alleles that confer resistance to pyrethroid insecticides were determined for each fly population. Levels of resistance to the six insecticides varied among states and insecticides. Resistance to permethrin was highest overall and most consistent across the states. Resistance to methomyl was relatively consistent, with 65-91% survival in nine of the ten collections. In contrast, resistance to cyfluthrin and pyrethrins + piperonyl butoxide varied considerably (2.9-76% survival). Resistance to imidacloprid was overall modest and showed no signs of increasing relative to collections made in 2004, despite increasing use of this insecticide. The frequency of Vssc alleles that confer pyrethroid resistance was variable between locations. The highest frequencies of kdr, kdr-his and super-kdr were found in Minnesota, North Carolina and Kansas, respectively. In contrast, the New Mexico population had the highest frequency (0.67) of the susceptible allele. The implications of these results to resistance management and to the understanding of the evolution of insecticide resistance are discussed.

A comparison of the electromechanical properties of structurally diverse proteins by molecular dynamics simulation

Proteins are subjected to electric fields both within the cell and during routine biochemical analysis. We have used atomistic molecular dynamics simulations to study conformational changes within three structurally diverse proteins subjected to high electric fields. At electric fields in excess of .5 V/nm, major structural changes were observed in all three proteins due to charge redistribution within the biomolecule. However, the electromechanical resilience was found to be highly dependent on the protein secondary structure, with α-helices showing a particularly high susceptibility to deformation by the applied electric field.

Molecular dynamics simulations of mechanical failure in polymorphic arrangements of amyloid fibrils containing structural defects

We examine how the different steric packing arrangements found in amyloid fibril polymorphs can modulate their mechanical properties using steered molecular dynamics simulations. Our calculations demonstrate that for fibrils containing structural defects, their ability to resist force in a particular direction can be dominated by both the number and molecular details of the defects that are present. The simulations thereby suggest a hierarchy of factors that govern the mechanical resilience of fibrils, and illustrate the general principles that must be considered when quantifying the mechanical properties of amyloid fibres containing defects.

Thermodynamics of writhe in DNA minicircles from molecular dynamics simulations

DNA supercoiling plays a role in genetic control by imposing torsional stress. This can induce writhe, which changes the global shape of the DNA. We have used atomistic molecular dynamics simulations to partition the free energy changes driving the writhing and unwrithing transitions in supercoiled minicircle DNA. The calculations show that while writhing is energetically driven, the unwrithing transition occurs because the circular state has a higher configurational entropy than the plectoneme. Writhing improves the van der Waals interactions between stacked bases, but can be suppressed by electrostatic repulsion within the negatively charged backbone strands in low salt conditions where electrostatic screening is poor. The free energy difference between circular and plectonemic DNA is determined by such a delicate balance of opposing thermodynamic terms that any perturbation in the environment, such as a change in salt concentration, can be sufficient to convert between these two states. This switchable behavior provides a mechanism for supercoiled DNA to store and communicate biological information physically as well as chemically.

Multiple origins of kdr-type resistance in the house fly, Musca domestica

Insecticide resistance is a model phenotype that can be used to investigate evolutionary processes underlying the spread of alleles across a global landscape, while offering valuable insights into solving the problems that resistant pests present to human health and agriculture. Pyrethroids are one of the most widely used classes of insecticides world-wide and they exert their toxic effects through interactions with the voltage-sensitive sodium channel (Vssc). Specific mutations in Vssc (kdr, kdr-his and super-kdr) are known to cause resistance to pyrethroid insecticides in house flies. In order to determine the number of evolutionary origins of kdr, kdr-his and super-kdr, we sequenced a region of Vssc from house flies collected in the USA, Turkey and China. Our phylogenetic analysis of Vssc unequivocally supports the hypothesis of multiple independent origins of kdr, super-kdr and kdr-his on an unprecedented geographic scale. The implications of these evolutionary processes on pest management are discussed.

Behavior of 3-isobutyl-2-hydroxypyrazine (IBHP), a key intermediate in 3-isobutyl-2-methoxypyrazine (IBMP) metabolism, in ripening wine grapes

3-Isobutyl-2-hydroxypyrazine (IBHP) is thought to be a key intermediate in both the biosynthesis and degradation of the herbaceous smelling 3-isobutyl-2-methoxypyrazine (IBMP), but its behavior during the growing season is not well understood. First, an improved method for IBHP quantification was developed. A deuterated IBHP standard was added to samples prior to isolation by mixed-mode cation exchange solid phase extraction. Extracts were silylated prior to quantification by GC-MS. A limit of detection of ca. 20 ng/L could be achieved for a 100 mL juice sample. This method was used to quantify IBHP during the 2010 growing season in berries of two clones of Cabernet franc in the Finger Lakes region of New York and of Merlot grown in the California Central Valley. For all three sources, IBHP was detectable at the earliest sampling point, and its concentration per berry increased to a maximum around veraison, 208-477 pg/berry. On a per berry basis, IBHP peaked and began to decline 1-2 weeks after IBMP, indicating that previous studies that sampled preveraison fruit have missed the true maximum value of IBHP. The highest per berry concentration of IBHP observed was in the California Merlot. However, after veraison, IBHP declined more rapidly in the California Merlot than in the New York Cabernet franc, such that the Merlot had the lowest IBHP concentration at harvest. Thus, IBHP at harvest cannot be used as a proxy for IBMP at veraison as was previously suggested.

Phosphorylation as a tool to modulate aggregation propensity and to predict fibril architecture

Despite the importance of post-translational modifications in controlling the solubility and conformational properties of proteins and peptides, precisely how the aggregation propensity of different peptide sequences is modulated by chemical modification remains unclear. Here we have investigated the effect of phosphorylation on the aggregation propensity of a 13-residue synthetic peptide incorporating one or more phosphate groups at seven different sites at various pH values. Fibril formation was shown to be inhibited when a single phosphate group was introduced at all seven locations in the peptide sequence at pH 7.5, when the phosphate group is fully charged. By contrast, when the same peptides were analysed at pH 1.1, when the phosphate is fully protonated, fibrils from all seven peptide sequences form rapidly. At intermediate pH values (pH 3.6) when the phosphate group is mono-anionic, the aggregation propensity of the peptides was found to be highly dependent on the position of the phosphate group in the peptide sequence. Using this information, combined with molecular dynamics (MD) simulations of the peptide sequence, we provide evidence consistent with the peptide forming amyloid fibrils with a class 7 architecture. The results highlight the potential utility of phosphorylation as a method of reversibly controlling the aggregation kinetics of peptide sequences both during and after synthesis. Moreover, by exploiting the ability of the phosphate group to adopt different charge states as a function of pH, and combining experimental insights with atomistic information calculated from MD simulations as pH is varied, we show how the resulting information can be used to predict fibril structures consistent with both datasets, and use these to rationalise their sensitivity of fibrillation kinetics both to the location of the phosphate group and its charge state.

Priority issues in occupational cancer research: Ontario stakeholder perspectives

INTRODUCTION

Workers are potentially exposed to known and suspected carcinogens in the workplace, many of which have not been fully evaluated. Despite persistent need, research on occupational cancer appears to have declined in recent decades. The formation of the Occupational Cancer Research Centre (OCRC) is an effort to counter this downward trend in Ontario. The OCRC conducted a survey of the broad stakeholder community to learn about priority issues on occupational cancer research.

METHODS

The OCRC received 177 responses to its survey from academic, health care, policy, industry, and labour-affiliated stakeholders. Responses were analyzed based on workplace exposures, at-risk occupations and cancers by organ system, stratified by respondents' occupational role.

DISCUSSION

Priority issues identified included workplace exposures such as chemicals, respirable dusts and fibres (e.g. asbestos), radiation (e.g. electromagnetic fields), pesticides, and shift work; and occupations such as miners, construction workers, and health care workers. Insufficient funding and a lack of exposure data were identified as the central barriers to conducting occupational cancer research.

CONCLUSION

The results of this survey underscore the great need for occupational cancer research in Ontario and beyond. They will be very useful as the OCRC develops its research agenda.

Systematic examination of polymorphism in amyloid fibrils by molecular-dynamics simulation

Amyloid fibrils often exhibit polymorphism. Polymorphs are formed when proteins or peptides with identical sequences self-assemble into fibrils containing substantially different arrangements of the β-strands. We used atomistic molecular-dynamics simulation to examine the thermodynamic stability of a amyloid fibrils in different polymorphic forms by performing a systematic investigation of sequence and symmetry space for a series of peptides with a range of physicochemical properties. We show that the stability of fibrils depends on both sequence and the symmetry because these factors determine the availability of favorable interactions between the peptide strands within a sheet and in intersheet packing. By performing a detailed analysis of these interactions as a function of symmetry, we obtained a series of simple design rules that can be used to determine which polymorphs of a given sequence are most likely to form thermodynamically stable fibrils. These rules can potentially be employed to design peptide sequences that aggregate into a preferred polymorphic form for nanotechnological purposes.

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

Sporophytic self-incompatibility in Senecio squalidus (Asteraceae): S allele dominance interactions and modifiers of cross-compatibility and selfing rates

Understanding genetic mechanisms of self-incompatibility (SI) and how they evolve is central to understanding the mating behaviour of most outbreeding angiosperms. Sporophytic SI (SSI) is controlled by a single multi-allelic locus, S, which is expressed in the diploid (sporophyte) plant to determine the SI phenotype of its haploid (gametophyte) pollen. This allows complex patterns of independent S allele dominance interactions in male (pollen) and female (pistil) reproductive tissues. Senecio squalidus is a useful model for studying the genetic regulation and evolution of SSI because of its population history as an alien invasive species in the UK. S. squalidus maintains a small number of S alleles (7-11) with a high frequency of dominance interactions. Some S. squalidus individuals also show partial selfing and/or greater levels of cross-compatibility than expected under SSI. We previously speculated that these might be adaptations to invasiveness. Here we describe a detailed characterization of the regulation of SSI in S. squalidus. Controlled crosses were used to determine the S allele dominance hierarchy of six S alleles and effects of modifiers on cross-compatibility and partial selfing. Complex dominance interactions among S alleles were found with at least three levels of dominance and tissue-specific codominance. Evidence for S gene modifiers that increase selfing and/or cross-compatibility was also found. These empirical findings are discussed in the context of theoretical predictions for maintenance of S allele dominance interactions, and the role of modifier loci in the evolution of SI.

An in silico design for a DNA nanomechanical switch

We have calculated how the charge transfer properties of DNA change in response to the application of an external stretching force. Since charge transfer occurs through the DNA π-stack, any disruption to this stacking causes dramatic changes in the transport properties of the biomolecule, as our calculations demonstrate. We therefore propose that the mechanical response of DNA to an applied stretching force might be used in the design of a nanomechanical switch.