SLTCAP: A Simple Method for Calculating the Number of Ions Needed for MD Simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Schmit JD
PROVIDER: S-EPMC6022740 | biostudies-literature | 2018 Apr
REPOSITORIES: biostudies-literature
Publications
SLTCAP: A Simple Method for Calculating the Number of Ions Needed for MD Simulation.
Schmit Jeremy D JD Kariyawasam Nilusha L NL Needham Vince V Smith Paul E PE
Journal of chemical theory and computation 20180316 4
An accurate depiction of electrostatic interactions in molecular dynamics requires the correct number of ions in the simulation box to capture screening effects. However, the number of ions that should be added to the box is seldom given by the bulk salt concentration because a charged biomolecule solute will perturb the local solvent environment. We present a simple method for calculating the number of ions that requires only the total solute charge, solvent volume, and bulk salt concentration ...[more]
Similar Datasets
Project description:Hunting records have proven useful for examining the historical status of wildlife populations. The number of animals harvested can provide information on past population sizes that would have been required to support harvest yields. Therefore, when statistical data on annual harvests are available, a minimum estimate of past population sizes can be calculated. A very simple method for estimating the sizes of historic wildlife populations using only annual hunting records and the maximum annual population increase rate is presented in this study. This method was applied to estimate past population sizes for Japanese sika deer (Cervus nippon yesoensis) in Hokkaido Island, Japan, using hunting records from 1873 to 1882, and assuming 15% and 35% population increase rates. The annual number of deer harvested during 1873 to 1882 ranged from 15,000 to 129,000. The minimum population size in 1873 was estimated as 349,000-473,000. This method was validated by applying it to the eastern population of Hokkaido Island in 1993 when the population size was approximately 260,000, and population sizes estimated by this method were 0.50-1.17 times the nominal population size. Thus, the population estimates from this method were approximately equal to or less than the expected population sizes, and this method can be used to obtain minimum estimates of wildlife populations. Because shorter durations of hunting records result in population size underestimates, it would be better to use hunting record of 10 years or longer in this method. In addition, the degree of underestimation may change with hunting pressure intensity on the populations, other causes of mortality, and maximum annual increase rates of the species. The method can be applied to any wildlife species for which records of annual harvest and maximum annual population increase rates of the species are available. The estimates obtained can provide benchmarks for the population size required for ecosystem conservation, and can be useful for wildlife management as they indicate the lowest limit to maintain the population.
| S-EPMC5997341 | biostudies-literature
Project description:Variations and alterations of copy numbers (CNVs and CNAs) carry disease susceptibility and drug responsiveness implications. Although there are many molecular methods to measure copy numbers, sensitivity, reproducibility, cost, and time issues remain. In the present study, we were able to solve those problems utilizing our modified real competitive PCR method with cloned competitors (mrcPCR). First, the mrcPCR for ERBB2 copy number was established, and the results were comparable to current standard methods but with a shorter assay time and a lower cost. Second, the mrcPCR assays for 24 drug-target genes were established, and the results in a panel of NCI-60 cells were comparable to those from real-time PCR and microarray. Third, the mrcPCR results for FCGR3A and the FCGR3B CNVs were comparable to those by the paralog ratio test (PRT), but without PRT's limitations. These results suggest that mrcPCR is comparable to the currently available standard or the most sensitive methods. In addition, mrcPCR would be invaluable for measurement of CNVs in genes with variants of similar structures, because combination of the other methods is not necessary, along with its other advantages such as short assay time, small sample amount requirement, and applicability to all sequences and genes.
| S-EPMC3728337 | biostudies-literature
Project description:BackgroundWe propose a new measure of treatment effect based on the expected reduction in the number of patients to treat (RNT) which is defined as the difference of the reciprocals of clinical measures of interest between two arms. Compared with the conventional number needed to treat (NNT), RNT shows superiority with both binary and time-to-event endpoints in randomized controlled trials (RCTs).MethodsFive real RCTs, two with binary endpoints and three with survival endpoints, are used to illustrate the concept of RNT and compare the performances between RNT and NNT. For survival endpoints, we propose two versions of RNT: one is based on the survival rate and the other is based on the restricted mean survival time (RMST). Hypothetical scenarios are also constructed to explore the advantages and disadvantages of RNT and NNT.ResultsBecause there is no baseline for computation of NNT, it fails to differentiate treatment effect in the absolute scale. In contrast, RNT conveys more information than NNT due to its reversed order of differencing and inverting. For survival endpoints, two versions of RNT calculated as the difference of the reciprocals of survival rates and RMSTs are complementary to each other. The RMST-based RNT can capture the entire follow-up profile and thus is clinically more intuitive and meaningful, as it inherits the time-to-event characteristics for survival endpoints instead of using truncated binary endpoints at a specific time point.ConclusionsThe RNT can serve as an alternative measure for quantifying treatment effect in RCTs, which complements NNT to help patients and clinicians better understand the magnitude of treatment benefit.
| S-EPMC7945324 | biostudies-literature
Project description:Most fish swim with body undulations that result from fluid-structure interactions between the fish's internal tissues and the surrounding water. Gaining insight into these complex fluid-structure interactions is essential to understand how fish swim. To this end, we developed a dedicated experimental-numerical inverse dynamics approach to calculate the lateral bending moment distributions for a large-amplitude undulatory swimmer that moves freely in three-dimensional space. We combined automated motion tracking from multiple synchronised high-speed video sequences, computation of fluid dynamic stresses on the swimmer's body from computational fluid dynamics, and bending moment calculations using these stresses as input for a novel beam model of the body. The bending moment, which represent the system's net actuation, varies over time and along the fish's central axis due to muscle actions, passive tissues, inertia, and fluid dynamics. Our three-dimensional analysis of 113 swimming events of zebrafish larvae ranging in age from 3 to 12 days after fertilisation shows that these bending moment patterns are not only relatively simple but also strikingly similar throughout early development and from fast starts to periodic swimming. This suggests that fish larvae may produce and adjust swimming movements relatively simply, yet effectively, while restructuring their neuromuscular control system throughout their rapid development.
| S-EPMC7481021 | biostudies-literature
Project description:Based on the original work of Rutherford (Radio-activity, 1905) and Bateman (Proc Camb Philos Soc 15:423-427, 1910), the authors designed two schemes consisting of explicit equations as simple methods for accurately obtaining activity of 90Sr and 90Y before they reach secular equilibrium. Application of the methods to the 90Sr/90Y system where 90Sr and 90Y are not in equilibrium will substantially reduce the time needed for determining activity of 90Sr because neither sequential measurements of 90Sr or 90Y (up to about 2 weeks) nor waiting for 90Sr and 90Y to reach equilibrium (more than 3 weeks) will be needed. We also implemented the explicit equations for decay/ingrowth correction of progeny's activity and applied them to the 95Zr/95Nb system. Using the equations, the authors corrected activity concentrations of 95Nb to a designated reference time from the activity concentrations measured from samples at different times, for instance, 2, 8, 15, and 29 days after a reference time. During those measurement times, 95Nb and 95Zr were not in equilibrium. The corrected 95Nb activity concentrations were within an accuracy of - 10%.
| S-EPMC6980259 | biostudies-literature
Project description:PURPOSE:We propose a new method to calculate proptosis by using the simple Heron's formula and analyze its feasibility. METHOD:It was a none-inferiority trial. The registration number was ChiCTR1900026490. The absolute value of proptosis in 120 eyes, 60 patients without eye injury or diseases, was measured by computed tomography (CT) and simple Heron's formula. We did regression analysis and analyzed the differences between the two methods with Medcalc software version 19.0.4. The result was showed by Passing-Bablok regression analysis diagram and Bland and Altman plot. RESULTS:The Passing-Bablok showed that the result of proptosis measured by CT and simple Heron's formula showed good positive correlation. A 95% limit of agreement in proptosis between CT and Heron's formula method was -0.46 to 0.54 mm in right eye and -0.45 to 0.46 mm in left eye. 1.66% (1/60) point was outside 95% LoA in both eyes. Moreover, a 95% limit of agreement between CT and Heron's formula method was -0.42 to 0.56 mm in difference of both eyes. 3.33% (2/60) points were outside 95% LoA. The points in all Bland and Altman plots were lower than 5%. It means that the results of comparison between the two methods had a good consistency in the measurement of proptosis. CONCLUSIONS:Heron's formula could be applied to calculate proptosis and has a good consistency compared with computed tomography (CT). This method is practical in proptosis assessment because of its accuracy, reliability and simplicity.
| S-EPMC7263846 | biostudies-literature
Project description:Dating of wood is a major task in historical research, archaeology and paleoclimatology. Currently, the most important dating techniques are dendrochronology and radiocarbon dating. Our approach is based on molecular decay over time under specific preservation conditions. In the models presented here, construction wood, cold soft waterlogged wood and wood from living trees are combined. Under these conditions, molecular decay as a usable clock for dating purposes takes place with comparable speed. Preservation conditions apart from those presented here are not covered by the model and cannot currently be dated with this method. For example, samples preserved in a clay matrix seem not to fit into the model. Other restrictions are discussed in the paper. One model presented covers 7,500 years with a root mean square error (RMSE) of 682 years for a single measurement. Another model reduced to the time period of the last 800 years results in a RMSE of 92 years. As multiple measurements can be performed on a single object, the total error for the whole object will be even lower.
| S-EPMC7347527 | biostudies-literature
Project description:The purpose of this study was to develop a quantitative framework to estimate the likelihood of multifocal epilepsy based on the number of unifocal seizures observed in the epilepsy monitoring unit (EMU).Patient records from the EMU at Massachusetts General Hospital (MGH) from 2012 to 2014 were assessed for the presence of multifocal seizures as well the presence of multifocal interictal discharges and multifocal structural imaging abnormalities during the course of the EMU admission. Risk factors for multifocal seizures were assessed using sensitivity and specificity analysis. A Kaplan-Meier survival analysis was used to estimate the risk of multifocal epilepsy for a given number of consecutive seizures. To overcome the limits of the Kaplan-Meier analysis, a parametric survival function was fit to the EMU subjects with multifocal seizures and this was used to develop a Bayesian model to estimate the risk of multifocal seizures during an EMU admission.Multifocal interictal discharges were a significant predictor of multifocal seizures within an EMU admission with a p < 0.01, albeit with only modest sensitivity 0.74 and specificity 0.69. Multifocal potentially epileptogenic lesions on MRI were not a significant predictor p = 0.44. Kaplan-Meier analysis was limited by wide confidence intervals secondary to significant patient dropout and concern for informative censoring. The Bayesian framework provided estimates for the number of unifocal seizures needed to predict absence of multifocal seizures. To achieve 90% confidence for the absence of multifocal seizure, three seizures are needed when the pretest probability for multifocal epilepsy is 20%, seven seizures for a pretest probability of 50%, and nine seizures for a pretest probability of 80%.These results provide a framework to assist clinicians in determining the utility of trying to capture a specific number of seizures in EMU evaluations of candidates for epilepsy surgery.
| S-EPMC4877132 | biostudies-literature
Project description:We consider a setting in which we have a treatment and a potentially large number of covariates for a set of observations, and wish to model their relationship with an outcome of interest. We propose a simple method for modeling interactions between the treatment and covariates. The idea is to modify the covariate in a simple way, and then fit a standard model using the modified covariates and no main effects. We show that coupled with an efficiency augmentation procedure, this method produces clinically meaningful estimators in a variety of settings. It can be useful for practicing personalized medicine: determining from a large set of biomarkers the subset of patients that can potentially benefit from a treatment. We apply the method to both simulated datasets and real trial data. The modified covariates idea can be used for other purposes, for example, large scale hypothesis testing for determining which of a set of covariates interact with a treatment variable.
| S-EPMC4338439 | biostudies-literature
Project description:Objective:Schizophrenia clinical trials commonly measure observed changes in Positive and Negative Syndrome Scale (PANSS) total score. However, it is more intuitive to think of response vs nonresponse, a binary outcome. Assessing binary outcomes enables calculation of number needed to treat (NNT) for therapeutic outcomes, number needed to harm (NNH) for adverse outcomes, and likelihood to be helped or harmed (LHH) to demonstrate benefit/risk tradeoffs. Here, NNT, NNH, and LHH were used to evaluate the clinical usefulness of aripiprazole lauroxil in patients with an acute schizophrenia exacerbation. Methods:Categorical efficacy and tolerability data were taken from the pivotal Phase 3 trial evaluating aripiprazole lauroxil for treatment of an acute exacerbation of schizophrenia. NNT and NNH values, with 95% CIs, were calculated in this post hoc analysis. Results:Using the intent-to-treat population for the pooled doses of aripiprazole lauroxil (441 mg [n=196] and 882 mg [n=204] q4w), responder rates (?30% improvement from baseline PANSS total score) were 35.3% for aripiprazole lauroxil arms vs 18.4% for placebo (n=196), yielding a NNT of 6 (95% CI: 5-11). Discontinuation rates due to adverse events (AEs) were higher among patients randomized to placebo than to either aripiprazole lauroxil dose. Akathisia was the only AE with an incidence ?5% in each aripiprazole lauroxil group and at least twice that of placebo (11.6%, 11.5%, and 4.3% of the patients receiving aripiprazole lauroxil 441 mg, 882 mg, and placebo, respectively), producing a NNH of 14 (95% CI: 9-33) for pooled aripiprazole lauroxil doses vs placebo. Calculating LHH for therapeutic response vs akathisia, aripiprazole lauroxil was 2.3 times more likely to result in a therapeutic response than an incident of akathisia. Conclusion:Using metrics of NNT, NNH, and LHH, aripiprazole lauroxil was an efficacious and well-tolerated intervention in a pivotal study in patients with an acute schizophrenia exacerbation.
| S-EPMC6751763 | biostudies-literature