Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Development time and survival of Verrallina funerea (Theobald) (Diptera : Culicidae) immatures and other brackish water mosquito species in southeast Queensland, Australia

Verrallina funerea (Theobald) is a brackish water mosquito that is recognised as an important pest and vector in southeast Queensland, Australia. Immature development time and survival of Ve. funerea was defined in the laboratory in response to a range of temperatures (17-34 degrees C) and salinities (0-35 parts per thousand (p.p.t)). The expression of autogeny in this species was also assessed. Salinity only had a slight effect on mean development time from hatching to adult emergence (7.0-7.4 d at salinities of 0, 17.5 and 31.5 p.p.t) and survival was uniformly high (97.5-99.0%). Mean development times were shorter at 26, 29 and 32 degrees C (7.0, 6.8 and 6.8 d, respectively) and longest at 17 degrees C (12.2 d). The threshold temperature (t) was 5.8 degrees C and the thermal constant (K) was 142.9 degree-days above t. Survival to adulthood decreased from > 95% (at 17-29 degrees C) to 78% (at 32 degrees C) and 0% (at 34 degrees C). No expression of autogeny was observed. Immature development times of Ve. funerea, Ochlerotatus vigilax (Skuse) and Oc. procax (Skuse) were then determined under field conditions at Maroochy Shire. Following tide and rain inundation, cohorts of newly hatched larvae were monitored daily by dipping, and time until pupation was noted. Tidal inundation triggered hatching of Ve. funerea and Oc. vigilax larvae whereas Oc. procax larvae were found only after rain inundation. Estimates of Ve. funerea and Oc. vigilax field development times were similar (8-9 d) while Oc. procax development time was slightly longer (9-10 d). Based on these survey results, control activities targeting Ve. funerea must be initiated 4 d (if using Bacillus thuringiensis var. israelensis de Barjac) or 5 d (if using s-methoprene) after inundation. However, Casuarina glauca Sieber canopy and branchlets covering breeding habitats may present a problem for the penetration of such treatments.

Soil Carbon Turnover Measurement by Physical Fractionation at a Forest-to-Pasture Chronosequence in the Brazilian Amazon

The effect of conversion from forest-to-pasture upon soil carbon stocks has been intensively discussed, but few studies focus on how this land-use change affects carbon (C) distribution across soil fractions in the Amazon basin. We investigated this in the 20 cm depth along a chronosequence of sites from native forest to three successively older pastures. We performed a physicochemical fractionation of bulk soil samples to better understand the mechanisms by which soil C is stabilized and evaluate the contribution of each C fraction to total soil C. Additionally, we used a two-pool model to estimate the mean residence time (MRT) for the slow and active pool C in each fraction. Soil C increased with conversion from forest-to-pasture in the particulate organic matter (> 250 mu m), microaggregate (53-250 mu m), and d-clay (< 2 mu m) fractions. The microaggregate comprised the highest soil C content after the conversion from forest-to-pasture. The C content of the d-silt fraction decreased with time since conversion to pasture. Forest-derived C remained in all fractions with the highest concentration in the finest fractions, with the largest proportion of forest-derived soil C associated with clay minerals. Results from this work indicate that microaggregate formation is sensitive to changes in management and might serve as an indicator for management-induced soil carbon changes, and the soil C changes in the fractions are dependent on soil texture.

The crash involvement of older drivers is associated with their hazard perception latencies

Hazard perception in driving is the one of the few driving-specific skills associated with crash involvement. However, this relationship has only been examined in studies where the majority of individuals were younger than 65. We present the first data revealing an association between hazard perception and self-reported crash involvement in drivers aged 65 and over. In a sample of 271 drivers, we found that individuals whose mean response time to traffic hazards was slower than 6.68 seconds (the ROC-curve derived pass mark for the test) were 2.32 times (95% CI 1.46, 3.22) more likely to have been involved in a self-reported crash within the previous five years than those with faster response times. This likelihood ratio became 2.37 (95% CI 1.49, 3.28) when driving exposure was controlled for. As a comparison, individuals who failed a test of useful field of view were 2.70 (95% CI 1.44, 4.44) times more likely to crash than those who passed. The hazard perception test and the useful field of view measure accounted for separate variance in crash involvement. These findings indicate that hazard perception testing and training could be potentially useful for road safety interventions for this age group.

Diurnal Variation of Retinal Thickness with Spectral Domain OCT

Purpose. To investigate whether diurnal variation occurs in retinal thickness measures derived from spectral domain optical coherence tomography (SD-OCT). Methods. Twelve healthy adult subjects had retinal thickness measured with SD-OCT every 2 h over a 10 h period. At each measurement session, three average B-scan images were derived from a series of multiple B-scans (each from a 5 mm horizontal raster scan along the fovea, containing 1500 A-scans/B-scan) and analyzed to determine the thickness of the total retina, as well as the thickness of the outer retinal layers. Average thickness values were calculated at the foveal center, at the 0.5 mm diameter foveal region, and for the temporal parafovea (1.5 mm from foveal center) and nasal parafovea (1.5 mm from foveal center). Results. Total retinal thickness did not exhibit significant diurnal variation in any of the considered retinal regions (p > 0.05). Evidence of significant diurnal variation was found in the thickness of the outer retinal layers (p < 0.05), with the most prominent changes observed in the photoreceptor layers at the foveal center. The photoreceptor inner and outer segment layer thickness exhibited mean amplitude (peak to trough) of daily change of 7 ± 3 μm at the foveal center. The peak in thickness was typically observed at the third measurement session (mean measurement time, 13:06). Conclusions. The total retinal thickness measured with SD-OCT does not exhibit evidence of significant variation over the course of the day. However, small but significant diurnal variation occurs in the thickness of the foveal outer retinal layers.

Strategy use of children with Down syndrome in a delay of gratification situation

Background: The capacity to delay gratification has been shown to be a very important developmental task for children who are developing typically. There is evidence that children with Down syndrome have more difficulty with a delay of gratification task than typically developing children of the same mental age. This study focused on the strategies children with Down syndrome use while in a delay of gratification situation to ascertain if these contribute to the differences in delay times from those of typically developing children. Method: Thirty-two children with Down syndrome (15 females) and 50 typically developing children participated in the study. Children with Down syndrome had a mental age, as measured by the Stanford-Binet IV, between 36 and 66 months (M = 45.66). The typically developing children had a mean chronological age of 45.76 months. Children participated in a delay of gratification task where they were offered two or one small treats and asked which they preferred. They were then told that they could have the two treats if they waited for the researcher to return (an undisclosed time of 15 min). If they did not want to wait any longer they could call the researcher back but then they could have only one treat. Twenty-two of the children with Down syndrome and 43 of the typically developing children demonstrated understanding of the task and their data are included here. Sessions were videotaped for later analysis. Results: There were significant differences in the mean waiting times of the two groups. The mean of the waiting times for children with Down syndrome was 181.32 s (SD = 347.62) and was 440.21 s (SD = 377.59) for the typically developing children. Eighteen percent of the group with Down syndrome waited for the researcher to return in comparison to 35% of the typically developing group. Sixty-four percent of children with Down syndrome called the researcher back and the remainder (18%) violated. In the typically developing group 37% called the researcher back and 28% violated. The mean waiting time for the group of children with Down syndrome who called the researcher back was 24 s. Examination of strategy use in this group was therefore very limited. There appeared to be quite similar strategy use across the groups who waited the full 15 min. Conclusions: These results confirm the difficulty children with Down syndrome have in delaying gratification. Teaching strategies for waiting, using information drawn from the behaviours of children who are developing typically may be a useful undertaking. Examination of other contributors to delay ability (e.g., language skills) is also likely to be helpful in understanding the difficulties demonstrated in delaying gratification.

Bayesian hierarchical models in statistical quality control methods to improve healthcare in hospitals

Quality oriented management systems and methods have become the dominant business and governance paradigm. From this perspective, satisfying customers’ expectations by supplying reliable, good quality products and services is the key factor for an organization and even government. During recent decades, Statistical Quality Control (SQC) methods have been developed as the technical core of quality management and continuous improvement philosophy and now are being applied widely to improve the quality of products and services in industrial and business sectors. Recently SQC tools, in particular quality control charts, have been used in healthcare surveillance. In some cases, these tools have been modified and developed to better suit the health sector characteristics and needs. It seems that some of the work in the healthcare area has evolved independently of the development of industrial statistical process control methods. Therefore analysing and comparing paradigms and the characteristics of quality control charts and techniques across the different sectors presents some opportunities for transferring knowledge and future development in each sectors. Meanwhile considering capabilities of Bayesian approach particularly Bayesian hierarchical models and computational techniques in which all uncertainty are expressed as a structure of probability, facilitates decision making and cost-effectiveness analyses. Therefore, this research investigates the use of quality improvement cycle in a health vii setting using clinical data from a hospital. The need of clinical data for monitoring purposes is investigated in two aspects. A framework and appropriate tools from the industrial context are proposed and applied to evaluate and improve data quality in available datasets and data flow; then a data capturing algorithm using Bayesian decision making methods is developed to determine economical sample size for statistical analyses within the quality improvement cycle. Following ensuring clinical data quality, some characteristics of control charts in the health context including the necessity of monitoring attribute data and correlated quality characteristics are considered. To this end, multivariate control charts from an industrial context are adapted to monitor radiation delivered to patients undergoing diagnostic coronary angiogram and various risk-adjusted control charts are constructed and investigated in monitoring binary outcomes of clinical interventions as well as postintervention survival time. Meanwhile, adoption of a Bayesian approach is proposed as a new framework in estimation of change point following control chart’s signal. This estimate aims to facilitate root causes efforts in quality improvement cycle since it cuts the search for the potential causes of detected changes to a tighter time-frame prior to the signal. This approach enables us to obtain highly informative estimates for change point parameters since probability distribution based results are obtained. Using Bayesian hierarchical models and Markov chain Monte Carlo computational methods, Bayesian estimators of the time and the magnitude of various change scenarios including step change, linear trend and multiple change in a Poisson process are developed and investigated. The benefits of change point investigation is revisited and promoted in monitoring hospital outcomes where the developed Bayesian estimator reports the true time of the shifts, compared to priori known causes, detected by control charts in monitoring rate of excess usage of blood products and major adverse events during and after cardiac surgery in a local hospital. The development of the Bayesian change point estimators are then followed in a healthcare surveillances for processes in which pre-intervention characteristics of patients are viii affecting the outcomes. In this setting, at first, the Bayesian estimator is extended to capture the patient mix, covariates, through risk models underlying risk-adjusted control charts. Variations of the estimator are developed to estimate the true time of step changes and linear trends in odds ratio of intensive care unit outcomes in a local hospital. Secondly, the Bayesian estimator is extended to identify the time of a shift in mean survival time after a clinical intervention which is being monitored by riskadjusted survival time control charts. In this context, the survival time after a clinical intervention is also affected by patient mix and the survival function is constructed using survival prediction model. The simulation study undertaken in each research component and obtained results highly recommend the developed Bayesian estimators as a strong alternative in change point estimation within quality improvement cycle in healthcare surveillances as well as industrial and business contexts. The superiority of the proposed Bayesian framework and estimators are enhanced when probability quantification, flexibility and generalizability of the developed model are also considered. The empirical results and simulations indicate that the Bayesian estimators are a strong alternative in change point estimation within quality improvement cycle in healthcare surveillances. The superiority of the proposed Bayesian framework and estimators are enhanced when probability quantification, flexibility and generalizability of the developed model are also considered. The advantages of the Bayesian approach seen in general context of quality control may also be extended in the industrial and business domains where quality monitoring was initially developed.

Moments of action provide insight into critical times for advection–diffusion–reaction processes

In 2010 Berezhkovskii and coworkers introduced the concept of local accumulation time (LAT) as a finite measure of the time required for the transient solution of a reaction diffusion equation to effectively reach steady state(Biophys J. 99, L59 (2010); Phys Rev E. 83, 051906 (2011)). Berezhkovskii’s approach is a particular application of the concept of mean action time (MAT) that was introduced previously by McNabb (IMA J Appl Math. 47, 193 (1991)). Here, we generalize these previous results by presenting a framework to calculate the MAT, as well as the higher moments, which we call the moments of action. The second moment is the variance of action time; the third moment is related to the skew of action time, and so on. We consider a general transition from some initial condition to an associated steady state for a one–dimensional linear advection–diffusion–reaction partial differential equation(PDE). Our results indicate that it is possible to solve for the moments of action exactly without requiring the transient solution of the PDE. We present specific examples that highlight potential weaknesses of previous studies that have considered the MAT alone without considering higher moments. Finally, we also provide a meaningful interpretation of the moments of action by presenting simulation results from a discrete random walk model together with some analysis of the particle lifetime distribution. This work shows that the moments of action are identical to the moments of the particle lifetime distribution for certain transitions.

Improved surgical safety after laparoscopic compared to open surgery for apparent early stage endometrial cancer : results from a randomised controlled trial

AIM: To compare Total Laparoscopic Hysterectomy (TLH) and Total Abdominal Hysterectomy (TAH) with regard to surgical safety. METHODS: Between October 2005 and June 2010, 760 patients with apparent early stage endometrial cancer were enroled in a multicentre, randomised clinical trial (LACE) comparing outcomes following TLH or TAH. The main study end points for this analysis were surgical adverse events (AE), hospital length of stay, conversion from laparoscopy to laparotomy, including 753 patients who completed at least 6 weeks of follow-up. Postoperative AEs were graded according to Common Toxicity Criteria (V3), and those immediately life-threatening, requiring inpatient hospitalisation or prolonged hospitalisation, or resulting in persistent or significant disability/incapacity were regarded as serious AEs. RESULTS: The incidence of intra-operative AEs was comparable in either group. The incidence of post-operative AE CTC grade 3+ (18.6% in TAH, 12.9% in TLH, p 0.03) and serious AE (14.3% in TAH, 8.2% in TLH, p 0.007) was significantly higher in the TAH group compared to the TLH group. Mean operating time was 132 and 107 min, and median length of hospital stay was 2 and 5 days in the TLH and TAH group, respectively (p<0.0001). The decline of haemoglobin from baseline to day 1 postoperatively was 2g/L less in the TLH group (p 0.006). CONCLUSIONS: Compared to TAH, TLH is associated with a significantly decreased risk of major surgical AEs. A laparoscopic surgical approach to early stage endometrial cancer is safe.

Comment on "Local accumulation times for source, diffusion, and degradation models in two and three dimensions" [J. Chem. Phys. 138, 104121 (2013)]

In a recent paper, Gordon, Muratov, and Shvartsman studied a partial differential equation (PDE) model describing radially symmetric diffusion and degradation in two and three dimensions. They paid particular attention to the local accumulation time (LAT), also known in the literature as the mean action time, which is a spatially dependent timescale that can be used to provide an estimate of the time required for the transient solution to effectively reach steady state. They presented exact results for three-dimensional applications and gave approximate results for the two-dimensional analogue. Here we make two generalizations of Gordon, Muratov, and Shvartsman’s work: (i) we present an exact expression for the LAT in any dimension and (ii) we present an exact expression for the variance of the distribution. The variance provides useful information regarding the spread about the mean that is not captured by the LAT. We conclude by describing further extensions of the model that were not considered by Gordon,Muratov, and Shvartsman. We have found that exact expressions for the LAT can also be derived for these important extensions...

How long does it take for aquifer recharge or aquifer discharge processes to reach steady state?

Groundwater flow models are usually characterized as being either transient flow models or steady state flow models. Given that steady state groundwater flow conditions arise as a long time asymptotic limit of a particular transient response, it is natural for us to seek a finite estimate of the amount of time required for a particular transient flow problem to effectively reach steady state. Here, we introduce the concept of mean action time (MAT) to address a fundamental question: How long does it take for a groundwater recharge process or discharge processes to effectively reach steady state? This concept relies on identifying a cumulative distribution function, $F(t;x)$, which varies from $F(0;x)=0$ to $F(t;x) \to \infty$ as $t\to \infty$, thereby providing us with a measurement of the progress of the system towards steady state. The MAT corresponds to the mean of the associated probability density function $f(t;x) = \dfrac{dF}{dt}$, and we demonstrate that this framework provides useful analytical insight by explicitly showing how the MAT depends on the parameters in the model and the geometry of the problem. Additional theoretical results relating to the variance of $f(t;x)$, known as the variance of action time (VAT), are also presented. To test our theoretical predictions we include measurements from a laboratory–scale experiment describing flow through a homogeneous porous medium. The laboratory data confirms that the theoretical MAT predictions are in good agreement with measurements from the physical model.

Pre-crash traffic flow trend analysis on motorways

Crashes on motorway contribute to a significant proportion (40-50%) of non-recurrent motorway congestions. Hence reduce crashes will help address congestion issues (Meyer, 2008). Crash likelihood estimation studies commonly focus on traffic conditions in a Short time window around the time of crash while longer-term pre-crash traffic flow trends are neglected. In this paper we will show, through data mining techniques, that a relationship between pre-crash traffic flow patterns and crash occurrence on motorways exists, and that this knowledge has the potential to improve the accuracy of existing models and opens the path for new development approaches. The data for the analysis was extracted from records collected between 2007 and 2009 on the Shibuya and Shinjuku lines of the Tokyo Metropolitan Expressway in Japan. The dataset includes a total of 824 rear-end and sideswipe crashes that have been matched with traffic flow data of one hour prior to the crash using an incident detection algorithm. Traffic flow trends (traffic speed/occupancy time series) revealed that crashes could be clustered with regards of the dominant traffic flow pattern prior to the crash. Using the k-means clustering method allowed the crashes to be clustered based on their flow trends rather than their distance. Four major trends have been found in the clustering results. Based on these findings, crash likelihood estimation algorithms can be fine-tuned based on the monitored traffic flow conditions with a sliding window of 60 minutes to increase accuracy of the results and minimize false alarms.

Hyperopic defocus and diurnal changes in human choroid and axial length

Purpose To investigate the influence of monocular hyperopic defocus on the normal diurnal rhythms in axial length and choroidal thickness of young adults. Methods A series of axial length and choroidal thickness measurements (collected at ~3 hourly intervals, with the first measurement at ~9 am and the final measurement at ~9 pm) were obtained for 15 emmetropic young adults over three consecutive days. The natural diurnal rhythms (Day 1, no defocus), diurnal rhythms with monocular hyperopic defocus (Day 2, – 2.00 DS spectacle lens over the right eye), and the recovery from any defocus induced changes (Day 3, no defocus) in diurnal rhythms were examined. Results Both axial length and choroidal thickness underwent significant diurnal changes on each of the three measurement days (p<0.0001). The introduction of monocular hyperopic defocus resulted in significant changes in the diurnal variations observed in both parameters (p<0.05). A significant (p<0.001) increase in the mean amplitude (peak to trough) of change in axial length (mean increase, 0.016 ± 0.005 mm) and choroidal thickness (mean increase, 0.011 ± 0.003 mm) was observed on day 2 with hyperopic defocus compared to the two ‘no defocus’ days (days 1 and 3). At the second measurement (mean time 12:10 pm) on the day with hyperopic defocus, the eye was significantly longer by 0.012 ± 0.002 mm compared to the other two days (p<0.05). No significant difference was observed in the average timing of the daily peaks in axial length (mean peak time 12:12 pm) and choroidal thickness (21:02 pm) over the three days. Conclusions The introduction of monocular hyperopic defocus resulted in a significant increase in the amplitude of the diurnal change in axial length and choroidal thickness that returned to normal the following day after removal of the blur stimulus.

A simplified approach for calculating the moments of action for linear reaction-diffusion equations

The mean action time is the mean of a probability density function that can be interpreted as a critical time, which is a finite estimate of the time taken for the transient solution of a reaction-diffusion equation to effectively reach steady state. For high-variance distributions, the mean action time under-approximates the critical time since it neglects to account for the spread about the mean. We can improve our estimate of the critical time by calculating the higher moments of the probability density function, called the moments of action, which provide additional information regarding the spread about the mean. Existing methods for calculating the nth moment of action require the solution of n nonhomogeneous boundary value problems which can be difficult and tedious to solve exactly. Here we present a simplified approach using Laplace transforms which allows us to calculate the nth moment of action without solving this family of boundary value problems and also without solving for the transient solution of the underlying reaction-diffusion problem. We demonstrate the generality of our method by calculating exact expressions for the moments of action for three problems from the biophysics literature. While the first problem we consider can be solved using existing methods, the second problem, which is readily solved using our approach, is intractable using previous techniques. The third problem illustrates how the Laplace transform approach can be used to study coupled linear reaction-diffusion equations.