Prediction of novel and analogous folds using fragment assembly and fold recognition

A number of new and newly improved methods for predicting protein structure developed by the Jones–University College London group were used to make predictions for the CASP6 experiment. Structures were predicted with a combination of fold recognition methods (mGenTHREADER, nFOLD, and THREADER) and a substantially enhanced version of FRAGFOLD, our fragment assembly method. Attempts at automatic domain parsing were made using DomPred and DomSSEA, which are based on a secondary structure parsing algorithm and additionally for DomPred, a simple local sequence alignment scoring function. Disorder prediction was carried out using a new SVM-based version of DISOPRED. Attempts were also made at domain docking and “microdomain” folding in order to build complete chain models for some targets.

High-strength, healable, supramolecular polymer nanocomposites

A supramolecular polymer blend, formed via π-π interactions between a π-electron rich pyrenyl endcapped oligomer and a chain-folding oligomer containing pairs of π-electron poor naphthalene-diimide (NDI) units, has been reinforced with cellulose nanocrystals (CNCs) to afford a healable nanocomposite material. Nanocomposites with varying weight percentage of CNCs (from 1.25 to 20.0 wt.%) within the healable supramolecular polymeric matrix have been prepared via solvent casting followed by compression molding, and their mechanical properties and healing behavior have been evaluated. It is found that homogeneously dispersed films can be formed with CNCs at less than 10 wt.%. Above 10 wt.% CNC heterogeneous nanocomposites were obtained. All the nanocomposites formed could be re-healed upon exposure to elevated temperatures although, for the homogeneous films, it was found that the healing rate was reduced with increasing CNC content. The best combination of healing efficiency and mechanical properties was obtained with the 7.5 wt.% CNC nanocomposite which exhibited a tensile modulus enhanced by as much as a factor of 20 over the matrix material alone and could be fully re-healed at 85 °C within 30 minutes. Thus it is demonstrated that supramolecular nanocomposites can afford greatly enhanced mechanical properties relative to the unreinforced polymer, while still allowing efficient thermal healing.

Large is fast, small is tight: determinants of speed and affinity in subunit capture by a periplasmic chaperone

The chaperone/usher pathway assembles surface virulence organelles of Gram-negative bacteria, consisting of fibers of linearly polymerized protein subunits. Fiber subunits are connected through 'donor strand complementation': each subunit completes the immunoglobulin (Ig)-like fold of the neighboring subunit by donating the seventh β-strand in trans. Whereas the folding of Ig domains is a fast first-order process, folding of Ig modules into the fiber conformation is a slow second-order process. Periplasmic chaperones separate this process in two parts by forming transient complexes with subunits. Interactions between chaperones and subunits are also based on the principle of donor strand complementation. In this study, we have performed mutagenesis of the binding motifs of the Caf1M chaperone and Caf1 capsular subunit from Yersinia pestis and analyzed the effect of the mutations on the structure, stability, and kinetics of Caf1M-Caf1 and Caf1-Caf1 interactions. The results suggest that a large hydrophobic effect combined with extensive main-chain hydrogen bonding enables Caf1M to rapidly bind an early folding intermediate of Caf1 and direct its partial folding. The switch from the Caf1M-Caf1 contact to the less hydrophobic, but considerably tighter and less dynamic Caf1-Caf1 contact occurs via the zip-out-zip-in donor strand exchange pathway with pocket 5 acting as the initiation site. Based on these findings, Caf1M was engineered to bind Caf1 faster, tighter, or both faster and tighter. To our knowledge, this is the first successful attempt to rationally design an assembly chaperone with improved chaperone function.

Conductance of amyloid β based peptide filaments: structure–function relations.

Controlling the morphology of self-assembled peptide nanostructures, particularly those based on amyloid peptides, has been the focus of intense research. In order to exploit these structures in electronic applications, further understanding of their electronic behavior is required. In this work, the role of peptide morphology in determining electronic conduction along self-assembled peptide nanofilament networks is demonstrated. The peptides used in this work were based on the sequence AAKLVFF, which is an extension of a core sequence from the amyloid b peptide. We show that the incorporation of a non-natural amino acid, 2-thienylalanine, instead of phenylalanine improves the obtained conductance with respect to that obtained for a similar structure based on the native sequence, which was not the case for the incorporation of 3-thienylalanine. Furthermore, we demonstrate that the morphology of the self-assembled structures, which can be controlled by the solvent used in the assembly process, strongly affects the conductance, with larger conduction obtained for a morphology of long, straight filaments. Our results demonstrate that, similar to natural systems, the assembly and folding of peptides could be of great importance for optimizing their function as components of electronic devices. Hence, sequence design and assembly conditions can be used to control the performance of peptide based structures in such electronic applications.

Radiative forcing due to aviation water vapour emissions

Three emissions inventories have been used with a fully Lagrangian trajectory model to calculate the stratospheric accumulation of water vapour emissions from aircraft, and the resulting radiative forcing. The annual and global mean radiative forcing due to present-day aviation water vapour emissions has been found to be 0.9 [0.3 to 1.4] mW m^2. This is around a factor of three smaller than the value given in recent assessments, and the upper bound is much lower than a recently suggested 20 mW m^2 upper bound. This forcing is sensitive to the vertical distribution of emissions, and, to a lesser extent, interannual variability in meteorology. Large differences in the vertical distribution of emissions within the inventories have been identified, which result in the choice of inventory being the largest source of differences in the calculation of the radiative forcing due to the emissions. Analysis of Northern Hemisphere trajectories demonstrates that the assumption of an e-folding time is not always appropriate for stratospheric emissions. A linear model is more representative for emissions that enter the stratosphere far above the tropopause.

Towards data warehousing and mining of protein unfolding simulation data

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Grid warehousing of molecular dynamics protein unfolding data

With the increasing awareness of protein folding disorders, the explosion of genomic information, and the need for efficient ways to predict protein structure, protein folding and unfolding has become a central issue in molecular sciences research. Molecular dynamics computer simulations are increasingly employed to understand the folding and unfolding of proteins. Running protein unfolding simulations is computationally expensive and finding ways to enhance performance is a grid issue on its own. However, more and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. This paper describes efforts to provide a grid-enabled data warehouse for protein unfolding data. We outline the challenge and present first results in the design and implementation of the data warehouse.

Realistic mean field forward predictions for the integration of co-registered EEG/fMRI

Brain activity can be measured non-invasively with functional imaging techniques. Each pixel in such an image represents a neural mass of about 105 to 107 neurons. Mean field models (MFMs) approximate their activity by averaging out neural variability while retaining salient underlying features, like neurotransmitter kinetics. However, MFMs incorporating the regional variability, realistic geometry and connectivity of cortex have so far appeared intractable. This lack of biological realism has led to a focus on gross temporal features of the EEG. We address these impediments and showcase a "proof of principle" forward prediction of co-registered EEG/fMRI for a full-size human cortex in a realistic head model with anatomical connectivity, see figure 1. MFMs usually assume homogeneous neural masses, isotropic long-range connectivity and simplistic signal expression to allow rapid computation with partial differential equations. But these approximations are insufficient in particular for the high spatial resolution obtained with fMRI, since different cortical areas vary in their architectonic and dynamical properties, have complex connectivity, and can contribute non-trivially to the measured signal. Our code instead supports the local variation of model parameters and freely chosen connectivity for many thousand triangulation nodes spanning a cortical surface extracted from structural MRI. This allows the introduction of realistic anatomical and physiological parameters for cortical areas and their connectivity, including both intra- and inter-area connections. Proper cortical folding and conduction through a realistic head model is then added to obtain accurate signal expression for a comparison to experimental data. To showcase the synergy of these computational developments, we predict simultaneously EEG and fMRI BOLD responses by adding an established model for neurovascular coupling and convolving "Balloon-Windkessel" hemodynamics. We also incorporate regional connectivity extracted from the CoCoMac database [1]. Importantly, these extensions can be easily adapted according to future insights and data. Furthermore, while our own simulation is based on one specific MFM [2], the computational framework is general and can be applied to models favored by the user. Finally, we provide a brief outlook on improving the integration of multi-modal imaging data through iterative fits of a single underlying MFM in this realistic simulation framework.

A Lagrangian model of air-mass photochemistry and mixing using a trajectory ensemble: the Cambridge Tropospheric Trajectory model of Chemistry And Transport (CiTTyCAT) version 4.2

A Lagrangian model of photochemistry and mixing is described (CiTTyCAT, stemming from the Cambridge Tropospheric Trajectory model of Chemistry And Transport), which is suitable for transport and chemistry studies throughout the troposphere. Over the last five years, the model has been developed in parallel at several different institutions and here those developments have been incorporated into one "community" model and documented for the first time. The key photochemical developments include a new scheme for biogenic volatile organic compounds and updated emissions schemes. The key physical development is to evolve composition following an ensemble of trajectories within neighbouring air-masses, including a simple scheme for mixing between them via an evolving "background profile", both within the boundary layer and free troposphere. The model runs along trajectories pre-calculated using winds and temperature from meteorological analyses. In addition, boundary layer height and precipitation rates, output from the analysis model, are interpolated to trajectory points and used as inputs to the mixing and wet deposition schemes. The model is most suitable in regimes when the effects of small-scale turbulent mixing are slow relative to advection by the resolved winds so that coherent air-masses form with distinct composition and strong gradients between them. Such air-masses can persist for many days while stretching, folding and thinning. Lagrangian models offer a useful framework for picking apart the processes of air-mass evolution over inter-continental distances, without being hindered by the numerical diffusion inherent to global Eulerian models. The model, including different box and trajectory modes, is described and some output for each of the modes is presented for evaluation. The model is available for download from a Subversion-controlled repository by contacting the corresponding authors.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

Suppression of gliadins results in altered protein body morphology in wheat

Wheat gluten proteins, gliadins and glutenins, are of great importance in determining the unique biomechanical properties of wheat. Studies have therefore been carried out to determine their pathways and mechanisms of synthesis, folding, and deposition in protein bodies. In the present work, a set of transgenic wheat lines has been studied with strongly suppressed levels of γ-gliadins and/or all groups of gliadins, using light and fluorescence microscopy combined with immunodetection using specific antibodies for γ-gliadins and HMW glutenin subunits. These lines represent a unique material to study the formation and fusion of protein bodies in developing seeds of wheat. Higher amounts of HMW subunits were present in most of the transgenic lines but only the lines with suppression of all gliadins showed differences in the formation and fusion of the protein bodies. Large rounded protein bodies were found in the wild-type lines and the transgenic lines with reduced levels of γ-gliadins, while the lines with all gliadins down-regulated had protein bodies of irregular shape and irregular formation. The size and number of inclusions, which have been reported to contain triticins, were also higher in the protein bodies in the lines with all the gliadins down-regulated. Changes in the protein composition and PB morphology reported in the transgenic lines with all gliadins down-regulated did not result in marked changes in the total protein content or instability of the different fractions.

Trafficking of storage proteins in developing grain of wheat

The processing properties of the wheat flour are largely determined by the structures and interactions of the grain storage proteins (also called gluten proteins) which form a continuous visco-elastic network in dough. Wheat gluten proteins are classically divided into two groups, the monomeric gliadins and the polymeric glutenins, with the latter being further classified into low molecular weight (LMW) and high molecular weight (HMW) subunits. The synthesis, folding and deposition of the gluten proteins take place within the endomembrane system of the plant cell. However, determination of the precise routes of trafficking and deposition of individual gluten proteins in developing wheat grain has been limited in the past by the difficulty of developing monospecific antibodies. To overcome this limitation, a single gluten protein (a LMW subunit) was expressed in transgenic wheat with a C-terminal epitope tag, allowing the protein to be located in the cells of the developing grain using highly specific antibodies. This approach was also combined with the use of wider specificity antibodies to compare the trafficking and deposition of different gluten protein groups within the same endosperm cells. These studies are in agreement with previous suggestions that two trafficking pathways occur in wheat, with the proteins either being transported via the Golgi apparatus into the vacuole or accumulating directly within the lumen of the ER. They also suggest that the same individual protein could be trafficked by either pathway, possibly depending on the stage of development, and that segregation of gluten proteins both between and within protein bodies may occur.

Expression of epitope-tagged LMW glutenin subunits in the starchy endosperm of transgenic wheat and their incorporation into glutenin polymers

The low-molecular-weight (LMW) glutenin subunits are components of the highly cross-linked glutenin polymers that confer viscoelastic properties to gluten and dough. They have both quantitative and qualitative effects on dough quality that may relate to differences in their ability to form the inter-chain disulphide bonds that stabilise the polymers. In order to determine the relationship between dough quality and the amounts and properties of the LMW subunits, we have transformed the pasta wheat cultivars Svevo and Ofanto with three genes encoding proteins, which differ in their numbers or positions of cysteine residues. The transgenes were delivered under control of the high-molecular-weight (HMW) subunit 1Dx5 gene promoter and terminator regions, and the encoded proteins were C-terminally tagged by the introduction of the c-myc epitope. Stable transformants were obtained with both cultivars, and the use of a specific antibody to the c-myc epitope tag allowed the transgene products to be readily detected in the complex mixture of LMW subunits. A range of transgene expression levels was observed. The addition of the epitope tag did not compromise the correct folding of the trangenic subunits and their incorporation into the glutenin polymers. Our results demonstrate that the ability to specifically epitope-tag LMW glutenin transgenes can greatly assist in the elucidation of their individual contributions to the functionality of the complex gluten system.

Spiral jetty, thirty seven years later: the cinematic time of James Benning

This article considers cinematic time in James Benning’s film, casting a glance (2007), in relation to its subject, Robert Smithson’s 1970 earthwork Spiral Jetty, and his film of the same name. The radicalism of Smithson’s thinking on time has been widely acknowledged, and his influence continues to pervade contemporary artistic practice. The relationship of Benning’s films with this legacy may appear somewhat oblique, given their apparent phenomenological rendition of ‘real time’. However, closer examination of Benning’s formal strategies reveals a more complex temporal construction, characterized by uncertain intervals that interrupt the folding of cinematic time into the flow of consciousness. Smithson’s film uses cinematic analogy to gesture towards vast reaches of geological time; Benning’s film creates a simulated timescale to evoke the short history of the earthwork itself. Smithson’s embrace of the entropic was a counter-cultural stance at the end of the1960s, but under the shadow of ecological disaster, this orientation has come to appear melancholy and romantic rather than radical. Benning’s film returns the jetty to anthropic time, but raises questions about the ways we inhabit time. His practice of working with ‘borrowed time’ is particularly suited to the cultural and historical moment of his later work.

Multivalency in healable supramolecular polymers: the effect of supramolecular cross-link density on the mechanical properties and healing of noncovalent polymer networks

Polymers with the ability to heal themselves could provide access to materials with extended lifetimes in a wide range of applications such as surface coatings, automotive components and aerospace composites. Here we describe the synthesis and characterisation of two novel, stimuli-responsive, supramolecular polymer blends based on π-electron-rich pyrenyl residues and π-electron-deficient, chain-folding aromatic diimides that interact through complementary π–π stacking interactions. Different degrees of supramolecular “cross-linking” were achieved by use of divalent or trivalent poly(ethylene glycol)-based polymers featuring pyrenyl end-groups, blended with a known diimide–ether copolymer. The mechanical properties of the resulting polymer blends revealed that higher degrees of supramolecular “cross-link density” yield materials with enhanced mechanical properties, such as increased tensile modulus, modulus of toughness, elasticity and yield point. After a number of break/heal cycles, these materials were found to retain the characteristics of the pristine polymer blend, and this new approach thus offers a simple route to mechanically robust yet healable materials.