T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

How curvature-generating proteins build scaffolds on membrane nanotubes

Bin/Amphiphysin/Rvs (BAR) domain proteins control the curvature of lipid membranes in endocytosis, trafficking, cell motility, the formation of complex sub-cellular structures, and many other cellular phenomena. They form three-dimensional assemblies, which act as molecular scaffolds to reshape the membrane and alter its mechanical properties. It is unknown, however, how a protein scaffold forms and how BAR domains interact in these assemblies at protein densities relevant for a cell. In this work, we employ various experimental, theoretical and simulation approaches to explore how BAR proteins organize to form a scaffold on a membrane nanotube. By combining quantitative microscopy with analytical modeling, we demonstrate that a highly curving BAR protein endophilin nucleates its scaffolds at the ends of a membrane tube, contrary to a weaker curving protein centaurin, which binds evenly along the tube’s length. Our work implies that the nature of local protein-membrane interactions can affect the specific localization of proteins on membrane-remodeling sites. Furthermore, we show that amphipathic helices are dispensable in forming protein scaffolds. Finally, we explore a possible molecular structure of a BAR-domain scaffold using coarse-grained molecular dynamics simulations. Together with fluorescence microscopy, the simulations show that proteins need only to cover 30–40% of a tube’s surface to form a rigid assembly. Our work provides mechanical and structural insights into the way BAR proteins may sculpt the membrane as a high-order cooperative assembly in important biological processes. 

Groundwater flow model of the Logan river alluvial aquifer system Josephville, South East Queensland

The study focuses on an alluvial plain situated within a large meander of the Logan River at Josephville near Beaudesert which supports a factory that processes gelatine. The plant draws water from on site bores, as well as the Logan River, for its production processes and produces approximately 1.5 ML per day (Douglas Partners, 2004) of waste water containing high levels of dissolved ions. At present a series of treatment ponds are used to aerate the waste water reducing the level of organic matter; the water is then used to irrigate grazing land around the site. Within the study the hydrogeology is investigated, a conceptual groundwater model is produced and a numerical groundwater flow model is developed from this. On the site are several bores that access groundwater, plus a network of monitoring bores. Assessment of drilling logs shows the area is formed from a mixture of poorly sorted Quaternary alluvial sediments with a laterally continuous aquifer comprised of coarse sands and fine gravels that is in contact with the river. This aquifer occurs at a depth of between 11 and 15 metres and is overlain by a heterogeneous mixture of silts, sands and clays. The study investigates the degree of interaction between the river and the groundwater within the fluvially derived sediments for reasons of both environmental monitoring and sustainability of the potential local groundwater resource. A conceptual hydrogeological model of the site proposes two hydrostratigraphic units, a basal aquifer of coarse-grained materials overlain by a thick semi-confining unit of finer materials. From this, a two-layer groundwater flow model and hydraulic conductivity distribution was developed based on bore monitoring and rainfall data using MODFLOW (McDonald and Harbaugh, 1988) and PEST (Doherty, 2004) based on GMS 6.5 software (EMSI, 2008). A second model was also considered with the alluvium represented as a single hydrogeological unit. Both models were calibrated to steady state conditions and sensitivity analyses of the parameters has demonstrated that both models are very stable for changes in the range of ± 10% for all parameters and still reasonably stable for changes up to ± 20% with RMS errors in the model always less that 10%. The preferred two-layer model was found to give the more realistic representation of the site, where water level variations and the numerical modeling showed that the basal layer of coarse sands and fine gravels is hydraulically connected to the river and the upper layer comprising a poorly sorted mixture of silt-rich clays and sands of very low permeability limits infiltration from the surface to the lower layer. The paucity of historical data has limited the numerical modelling to a steady state one based on groundwater levels during a drought period and forecasts for varying hydrological conditions (e.g. short term as well as prolonged dry and wet conditions) cannot reasonably be made from such a model. If future modelling is to be undertaken it is necessary to establish a regular program of groundwater monitoring and maintain a long term database of water levels to enable a transient model to be developed at a later stage. This will require a valid monitoring network to be designed with additional bores required for adequate coverage of the hydrogeological conditions at the Josephville site. Further investigations would also be enhanced by undertaking pump testing to investigate hydrogeological properties in the aquifer.

Modeling and exploiting behavior patterns in dynamic environments

This paper presents a new approach to improving the effectiveness of autonomous systems that deal with dynamic environments. The basis of the approach is to find repeating patterns of behavior in the dynamic elements of the system, and then to use predictions of the repeating elements to better plan goal directed behavior. It is a layered approach involving classifying, modeling, predicting and exploiting. Classifying involves using observations to place the moving elements into previously defined classes. Modeling involves recording features of the behavior on a coarse grained grid. Exploitation is achieved by integrating predictions from the model into the behavior selection module to improve the utility of the robot's actions. This is in contrast to typical approaches that use the model to select between different strategies or plays. Three methods of adaptation to the dynamic features of the environment are explored. The effectiveness of each method is determined using statistical tests over a number of repeated experiments. The work is presented in the context of predicting opponent behavior in the highly dynamic and multi-agent robot soccer domain (RoboCup).

Evidence for carbonate platform failure during rapid sea-level rise; ca 14 000 year old bioclastic flow deposits in the Lesser Antilles

Bioclastic flow deposits offshore from the Soufrie`re Hills volcano on Montserrat in the Lesser Antilles were deposited by the largest volume sediment flows near this active volcano in the last 26 kyr. The volume of these deposits exceeds that of the largest historic volcanic dome collapse in the world, which occurred on Montserrat in 2003. These flows were most probably generated by a large submarine slope failure of the carbonate shelf comprising the south west flank of Antigua or the east flank of Redonda; adjacent islands that are not volcanically active. The bioclastic flow deposits are relatively coarse-grained and either ungraded or poorly graded, and were deposited by non cohesive debris flow and high density turbidity currents. The bioclastic deposit often comprises multiple sub-units that cannot be correlated between core sites; some located just 2 km apart. Multiple sub-units in the bioclastic deposit result from either flow reflection, stacking of multiple debris flow lobes, and/or multi-stage collapse of the initial landslide. This study provides unusually precise constraints on the age of this mass flow event that occurred at ca 14 ka. Few large submarine landslides have been well dated, but the slope failures that have been dated are commonly associated with periods of rapid sea-level change.

Using static analysis for automatic assessment and mitigation of unwanted and malicious activities within Android applications

In the last decade, smartphones have gained widespread usage. Since the advent of online application stores, hundreds of thousands of applications have become instantly available to millions of smart-phone users. Within the Android ecosystem, application security is governed by digital signatures and a list of coarse-grained permissions. However, this mechanism is not fine-grained enough to provide the user with a sufficient means of control of the applications' activities. Abuse of highly sensible private information such as phone numbers without users' notice is the result. We show that there is a high frequency of privacy leaks even among widely popular applications. Together with the fact that the majority of the users are not proficient in computer security, this presents a challenge to the engineers developing security solutions for the platform. Our contribution is twofold: first, we propose a service which is able to assess Android Market applications via static analysis and provide detailed, but readable reports to the user. Second, we describe a means to mitigate security and privacy threats by automated reverse-engineering and refactoring binary application packages according to the users' security preferences.

Hierarchical multiscale model for biomechanics analysis of microfilament networks

The mechanisms of force generation and transference via microfilament networks are crucial to the understandings of mechanobiology of cellular processes in living cells. However, there exists an enormous challenge for all-atom physics simulation of real size microfilament networks due to scale limitation of molecular simulation techniques. Following biophysical investigations of constitutive relations between adjacent globular actin monomers on filamentous actin, a hierarchical multiscale model was developed to investigate the biomechanical properties of microfilament networks. This model was validated by previous experimental studies of axial tension and transverse vibration of single F-actin. The biomechanics of microfilament networks can be investigated at the scale of real eukaryotic cell size (10 μm). This multiscale approach provides a powerful modeling tool which can contribute to the understandings of actin-related cellular processes in living cells.

Two cycles of voluminous pyroclastic volcanism and sedimentation related to episodic granite emplacement during the late Archean : Eastern Yilgarn Craton, Western Australia

The thick piles of late-Archean volcaniclastic sedimentary successions that overlie the voluminous greenstone units of the eastern Yilgarn Craton, Western Australia, record the important transition from the cessation in mafic-ultramafic volcanism to cratonisation between about 2690 and 2655 Ma. Unfortunately, an inability to clearly subdivide the superficially similar sedimentary successions and correlate them between the various geological terranes and domains of the eastern Yilgarn Craton has led to uncertainty about the timing and nature of the region's palaeogeographic and palaeotectonic evolution. Here, we present the results of some 2025 U–Pb laser-ablation-ICP-MS analyses and 323 Sensitive High-Resolution Ion Microprobe (SHRIMP) analyses of detrital zircons from 14 late-Archean felsic clastic successions of the eastern Yilgarn Craton, which have enabled correlation of clastic successions. The results of our data, together with those compiled from previous studies, show that the post-greenstone sedimentary successions include two major cycles that both commenced with voluminous pyroclastic volcanism and ended with widespread exhumation and erosion associated with granite emplacement. Cycle One commences with an influx of rapidly reworked feldspar-rich pyroclastic debris. These units, here-named the Early Black Flag Group, are dominated by a single population of detrital zircons with an average age of 2690–2680 Ma. Thick (up to 2 km) dolerite bodies, such as the Golden Mile Dolerite, intrude the upper parts of the Early Black Flag Group at about 2680 Ma. Incipient development of large granite domes during Cycle One created extensional basins predominantly near their southeastern and northwestern margins (e.g., St Ives, Wallaby, Kanowna Belle and Agnew), into which the Early Black Flag Group and overlying coarse mafic conglomerate facies of the Late Black Flag Group were deposited. The clast compositions and detrital-zircon ages of the late Black Flag Group detritus match closely the nearby and/or stratigraphically underlying successions, thus suggesting relatively local provenance. Cycle Two involved a similar progression to that observed in Cycle One, but the age and composition of the detritus were notably different. Deposition of rapidly reworked quartz-rich pyroclastic deposits dominated by a single detrital-zircon age population of 2670–2660 Ma heralded the beginning of Cycle Two. These coarse-grained quartz-rich units, are name here the Early Merougil Group. The mean ages of the detrital zircons from the Early Merougil Group match closely the age of the peak in high-Ca (quartz-rich) granite magmatism in the Yilgarn Craton and thus probably represent the surface expression of the same event. Successions of the Late Merougil Group are dominated by coarse felsic conglomerate with abundant volcanic quartz. Although the detrital zircons in these successions have a broad spread of age, the principal sub-populations have ages of about 2665 Ma and thus match closely those of the Early Merougil Group. These successions occur most commonly at the northwestern and southeastern margins of the granite batholiths and thus are interpreted to represent resedimented units dominted by the stratigraphically underlying packages of the Early Merougil Group. The Kurrawang Group is the youngest sedimentary units identified in this study and is dominated by polymictic conglomerate with clasts of banded iron formation (BIF), granite and quartzite near the base and quartz-rich sandstone units containing detrital zircons aged up to 3500 Ma near the top. These units record provenance from deeper and/or more-distal sources. We suggest here that the principal driver for the major episodes of volcanism, sedimentation and deformation associated with basin development was the progressive emplacement of large granite batholiths. This interpretation has important implication for palaeogeographic and palaeotectonic evolution of all late-Archean terranes around the world.

Business process model abstraction

In order to execute, study, or improve operating procedures, companies document them as business process models. Often, business process analysts capture every single exception handling or alternative task handling scenario within a model. Such a tendency results in large process specifications. The core process logic becomes hidden in numerous modeling constructs. To fulfill different tasks, companies develop several model variants of the same business process at different abstraction levels. Afterwards, maintenance of such model groups involves a lot of synchronization effort and is erroneous. We propose an abstraction technique that allows generalization of process models. Business process model abstraction assumes a detailed model of a process to be available and derives coarse-grained models from it. The task of abstraction is to tell significant model elements from insignificant ones and to reduce the latter. We propose to learn insignificant process elements from supplementary model information, e.g., task execution time or frequency of task occurrence. Finally, we discuss a mechanism for user control of the model abstraction level – an abstraction slider.

Process model abstraction : a slider approach

Process models provide companies efficient means for managing their business processes. Tasks where process models are employed are different by nature and require models of various abstraction levels. However, maintaining several models of one business process involves a lot of synchronization effort and is erroneous. Business process model abstraction assumes a detailed model of a process to be available and derives coarse grained models from it. The task of abstraction is to tell significant model elements from insignificant ones and to reduce the latter. In this paper we argue that process model abstraction can be driven by different abstraction criteria. Criterion choice depends on a task which abstraction facilitates. We propose an abstraction slider - a mechanism that allows user control of the model abstraction level. We discuss examples of combining the slider with different abstraction criteria and sets of process model transformation rules.

Mechanics of microfilaments networks : a cross-scales study

The mechanical properties of microfilament networks are systematically summarized at different special scales in this paper. We have presented the mechanical models of single microfilaments and microfilament networks at microscale. By adopting a coarse-grained simulation strategy, the mechanical stability of microfilaments related cellular structures are analysed. Structural analysis is conducted to microfilament networks to understand the stress relaxation under compression. The nanoscale molecular mechanisms of the microfilaments deformation is also summarized from the viewpoint of molecular dynamics simulation. This paper provides the fundaments of multiscale modelling framework for the mechanical behaviours simulation of hierarchical microfilament networks.

Multivariate statistical analysis of hydrochemical data to assess alluvial aquifer–stream connectivity during drought and flood : Cressbrook Creek, southeast Queensland, Australia

A catchment-scale multivariate statistical analysis of hydrochemistry enabled assessment of interactions between alluvial groundwater and Cressbrook Creek, an intermittent drainage system in southeast Queensland, Australia. Hierarchical cluster analyses and principal component analysis were applied to time-series data to evaluate the hydrochemical evolution of groundwater during periods of extreme drought and severe flooding. A simple three-dimensional geological model was developed to conceptualise the catchment morphology and the stratigraphic framework of the alluvium. The alluvium forms a two-layer system with a basal coarse-grained layer overlain by a clay-rich low-permeability unit. In the upper and middle catchment, alluvial groundwater is chemically similar to streamwater, particularly near the creek (reflected by high HCO3/Cl and K/Na ratios and low salinities), indicating a high degree of connectivity. In the lower catchment, groundwater is more saline with lower HCO3/Cl and K/Na ratios, notably during dry periods. Groundwater salinity substantially decreased following severe flooding in 2011, notably in the lower catchment, confirming that flooding is an important mechanism for both recharge and maintaining groundwater quality. The integrated approach used in this study enabled effective interpretation of hydrological processes and can be applied to a variety of hydrological settings to synthesise and evaluate large hydrochemical datasets.

Physical mechanism of the compressive response of F-actin networks : significance of crosslinker unbinding events

A model of crosslinker unbinding is implemented in a highly coarsegrained granular model of F-actin cytoskeleton. We employ this specific granular model to study the mechanisms of the compressive responses of F-actin networks. It is found that the compressive response of F-actin cytoskeleton has dependency on the strain rate. The evolution of deformation energy in the network indicates that crosslinker unbinding events can induce the remodelling of F-actin cytoskeleton in response to external loadings. The internal stress in F-actin cytoskeleton can efficiently dissipate with the help of crosslinker unbinding, which could lead to the spontaneous relaxation of living cells.

Trans-scale granular modelling of cytoskeleton: A mini-review

Living cells are the functional unit of organs that controls reactions to their exterior. However, the mechanics of living cells can be difficult to characterize due to the crypticity of their microscale structures and associated dynamic cellular processes. Fortunately, multiscale modelling provides a powerful simulation tool that can be used to study the mechanical properties of these soft hierarchical, biological systems. This paper reviews recent developments in hierarchical multiscale modeling technique that aimed at understanding cytoskeleton mechanics. Discussions are expanded with respects to cytoskeletal components including: intermediate filaments, microtubules and microfilament networks. The mechanical performance of difference cytoskeleton components are discussed with respect to their structural and material properties. Explicit granular simulation methods are adopted with different coarse-grained strategies for these cytoskeleton components and the simulation details are introduced in this review.

Screened and unscreened phases in sedimenting suspensions

A coarse-grained stochastic hydrodynamical description of velocity and concentration fluctuations in steadily sedimenting suspensions is constructed and analyzed using self-consistent and renormalization-group methods. We find a nonequilibrium phase transition from an "unscreened" phase in which we recover the Caflisch-Luke [Phys. Fluids 28, 759 (1985)] divergence of the velocity variance to a "screened" phase where the fluctuations have a finite correlation length depending on the volume fraction phi as phi(-1/3), in agreement with Segre et al. [Phys. Rev. Lett. 79, 2574 (1997)] (if their observation of a phi-independent diffusivity is used), and the velocity variance is independent of system size.