Accuracy limit of rigid 3-point water models.
Ontology highlight
ABSTRACT:
SUBMITTER: Izadi S
PROVIDER: S-EPMC4991989 | biostudies-other | 2016 Aug
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Rigid n-point water models are widely used in atomistic simulations, but have known accuracy drawbacks. Increasing the number of point charges, as well as adding electronic polarizability, are two common strategies for accuracy improvements. Both strategies come at considerable computational cost, which weighs heavily against modest possible accuracy improvements in practical simulations. In an effort to provide guidance for model development, here we have explored the limiting accuracy of "electrostatically globally optimal" n-point water models in terms of their ability to reproduce properties of water dimer-a mimic of the condensed state of water. For a given n, each model is built upon a set of reference multipole moments (e.g. ab initio) and then optimized to reproduce water dimer total dipole moment. The models are then evaluated with respect to the accuracy of reproducing the geometry of the water dimer. We find that global optimization of the charge distribution alone can deliver high accuracy of the water model: for n = 4 or n = 5, the geometry of the resulting water dimer can be almost within 50 of the ab initio reference, which is half that of the experimental error margin. Thus, global optimization of the charge distribution of classical n-point water models can lead to high accuracy models. We also find that while the accuracy improvement in going from n = 3 to n = 4 is substantial, the additional accuracy increase in going from n = 4 to n = 5 is marginal. Next, we have explored accuracy limitations of the standard practice of adding electronic polarizability (via a Drude particle) to a "rigid base"-pre-optimization rigid n-point water model. The resulting model (n = 3) shows a relatively small improvement in accuracy, suggesting that the strategy of merely adding the polarizability to an inferior accuracy water model used as the base cannot fix the defects of the latter. An alternative strategy in which the parameters of the rigid base model are globally optimized along with the polarizability parameter is much more promising: the resulting 3-point polarizable model out-performs even the 5-point optimal rigid model by a large margin. We suggest that future development efforts consider 3- and 4-point polarizable models where global optimization of the "rigid base" is coupled to optimization of the polarizability to deliver globally optimal solutions.
| S-EPMC6855648 | biostudies-literature
Project description:Mutational optimization of two long-range interactions first observed in Ac-WINGKWT-NH2, (a) bifurcated H-bonding involving the threonine amide H(N) and side chain OH and the N-terminal acetyl carbonyl and (b) an H-bond between the entgegen-H(N) of the C-terminal amide and the indole ring of Trp6 that stabilizes a face-to-edge indole/indole interaction between Trp1 and Trp6, has afforded < or = 10 residue systems that yield a remarkably stable fold in water. Optimization was achieved by designing a hydrophobic cluster that sequesters these H-bonds from solvent exposure. The structures and extent of amide H/D exchange protection for CH3CH2CO-WI pGXWTGPS (p = D-Pro, X = Leu or Ile) were determined. These two systems are greater than 94% folded at 298 K (97.5% at 280 K) with melting temperatures > 75 degrees C. The fold appears to display minimal fluxionality; a well-converged NMR structure rationalizes all of the large structuring shifts observed, and we suggest that these designed constructs can be viewed as microproteins.
| S-EPMC2677142 | biostudies-literature
Project description:A new biradical polarizing agent, bTbtk-py, for dynamic nuclear polarization (DNP) experiments in aqueous media is reported. The synthesis is discussed in light of the requirements of the optimum, theoretical, biradical system. To date, the DNP NMR signal enhancement resulting from bTbtk-py is the largest of any biradical in the ideal glycerol/water solvent matrix, ? = 230. EPR and X-ray crystallography are used to characterize the molecule and suggest approaches for further optimizing the biradical distance and relative orientation.
| S-EPMC3304014 | biostudies-literature
Project description:PurposeRapid antigen-detecting tests (Ag-RDTs) for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can transform pandemic control. Thus far, sensitivity (≤ 85%) of lateral-flow assays has limited scale-up. Conceivably, microfluidic immunofluorescence Ag-RDTs could increase sensitivity for SARS-CoV-2 detection.MethodsThis multi-centre diagnostic accuracy study investigated performance of the microfluidic immunofluorescence LumiraDx™ assay, enrolling symptomatic and asymptomatic participants with suspected SARS-CoV-2 infection. Participants collected a supervised nasal mid-turbinate (NMT) self-swab for Ag-RDT testing, in addition to a professionally collected nasopharyngeal (NP) swab for routine testing with reverse transcriptase polymerase chain reaction (RT-PCR). Results were compared to calculate sensitivity and specificity. Sub-analyses investigated the results by viral load, symptom presence and duration. An analytical study assessed exclusivity and limit-of-detection (LOD). In addition, we evaluated ease-of-use.ResultsThe study was conducted between November 2nd 2020 and 4th of December 2020. 761 participants were enrolled, with 486 participants reporting symptoms on testing day. 120 out of 146 RT-PCR positive cases were detected positive by LumiraDx™, resulting in a sensitivity of 82.2% (95% CI 75.2-87.5%). Specificity was 99.3% (CI 98.3-99.7%). Sensitivity was increased in individuals with viral load ≥ 7 log10 SARS-CoV2 RNA copies/ml (93.8%; CI 86.2-97.3%). Testing against common respiratory commensals and pathogens showed no cross-reactivity and LOD was estimated to be 2-56 PFU/mL. The ease-of-use-assessment was favourable for lower throughput settings.ConclusionThe LumiraDx™ assay showed excellent analytical sensitivity, exclusivity and clinical specificity with good clinical sensitivity using supervised NMT self-sampling.Trial registration number and registration dateDRKS00021220 and 01.04.2020.
| S-EPMC8358901 | biostudies-literature
Project description:Vertebral fractures are common in the elderly, but efforts to reduce their incidence have been hampered by incomplete understanding of the failure processes that are involved. This study's goal was to elucidate failure processes in the lumbar vertebra and to assess the accuracy of quantitative computed tomography (QCT)-based finite element (FE) simulations of these processes. Following QCT scanning, spine segments (n?=?27) consisting of L1 with adjacent intervertebral disks and neighboring endplates of T12 and L2 were compressed axially in a stepwise manner. A microcomputed tomography scan was performed at each loading step. The resulting time-lapse series of images was analyzed using digital volume correlation (DVC) to quantify deformations throughout the vertebral body. While some diversity among vertebrae was observed on how these deformations progressed, common features were large strains that developed progressively in the superior third and, concomitantly, in the midtransverse plane, in a manner that was associated with spatial variations in microstructural parameters such as connectivity density. Results of FE simulations corresponded qualitatively to the measured failure patterns when boundary conditions were derived from DVC displacements at the endplate. However, quantitative correspondence was often poor, particularly when boundary conditions were simplified to uniform compressive loading. These findings suggest that variations in trabecular microstructure are one cause of the differences in failure patterns among vertebrae and that both the lack of incorporation of these variations into QCT-based FE models and the oversimplification of boundary conditions limit the accuracy of these models in simulating vertebral failure.
| S-EPMC5816243 | biostudies-literature
Project description:Data classification, the process of analyzing data and organizing it into categories or clusters, is a fundamental computing task of natural and artificial information processing systems. Both supervised classification and unsupervised clustering work best when the input vectors are distributed over the data space in a highly non-uniform way. These tasks become however challenging in weakly structured data sets, where a significant fraction of data points is located in between the regions of high point density. We derive the theoretical limit for classification accuracy that arises from this overlap of data categories. By using a surrogate data generation model with adjustable statistical properties, we show that sufficiently powerful classifiers based on completely different principles, such as perceptrons and Bayesian models, all perform at this universal accuracy limit under ideal training conditions. Remarkably, the accuracy limit is not affected by certain non-linear transformations of the data, even if these transformations are non-reversible and drastically reduce the information content of the input data. We further compare the data embeddings that emerge by supervised and unsupervised training, using the MNIST data set and human EEG recordings during sleep. We find for MNIST that categories are significantly separated not only after supervised training with back-propagation, but also after unsupervised dimensionality reduction. A qualitatively similar cluster enhancement by unsupervised compression is observed for the EEG sleep data, but with a very small overall degree of cluster separation. We conclude that the handwritten letters in MNIST can be considered as 'natural kinds', whereas EEG sleep recordings are a relatively weakly structured data set, so that unsupervised clustering will not necessarily re-cover the human-defined sleep stages.
| S-EPMC9772417 | biostudies-literature
Project description:Background:Free water (FW) in neuroimaging is non-flowing extracellular water in the cranium and brain tissue, and includes both cerebral spinal fluid (CSF) and fluid in intercellular space or edema. For a region such as a voxel (spatial unit of measurement in neuroimaging), the FW fraction is defined as the volume fraction of FW within that volume. Quantifying the FW fraction allows estimating contamination by fluid of neuroimaging or magnetic resonance spectroscopy measurements within a voxel. New method:An upper limit to the fraction of FW within a voxel, based on any diffusion tensor imaging (DTI) sequence including a standard single shell at one b-value, can be derived from the standard diffusion tensor by scaling the third eigenvalue of the diffusion tensor. Assuming a two-compartment model, the diffusivity of a voxel is a combination of tissue and FW diffusivity. FW fraction is FW volume divided by voxel volume. Assuming FW diffuses equally in all directions, the diffusivity component is representable by a single, non-tensor diffusivity value. Since the diffusivity of water is known for a given temperature, and brain temperature is relatively constant, the FW diffusivity value can be assumed constant. The third eigenvector of the voxel diffusion tensor is the direction of least diffusivity and since the FW component of diffusivity is equal in all directions, we show that FW diffusivity cannot be lower than the third eigenvalue. Assuming FW contributes proportionally to voxel diffusivity, we show that the third eigenvalue divided by water diffusivity (as a constant based on known water diffusivity at 36.7 °C) forms an upper limit on the FW-fraction (fUL ). Results:We calculated fUL for 384 subjects from the IXI dataset. Values mostly ranged from 0.1 to 0.6, and were closely related to radial diffusivity.Comparison with Existing Methods:fUL is easily calculated from any DTI data, but is not a true estimate of FW-fraction. Conclusions:The fUL measure offers a starting point in calculating the true FW-fraction of a voxel, or an easy-to-calculate voxel characteristic.
| S-EPMC6072896 | biostudies-literature
Project description:We present a rigid-body-based technique (called rigid-cluster elastic network interpolation) to generate feasible transition pathways between two distinct conformations of a macromolecular assembly. Many biological molecules and assemblies consist of domains which act more or less as rigid bodies during large conformational changes. These collective motions are thought to be strongly related with the functions of a system. This fact encourages us to simply model a macromolecule or assembly as a set of rigid bodies which are interconnected with distance constraints. In previous articles, we developed coarse-grained elastic network interpolation (ENI) in which, for example, only Calpha atoms are selected as representatives in each residue of a protein. We interpolate distance differences of two conformations in ENI by using a simple quadratic cost function, and the feasible conformations are generated without steric conflicts. Rigid-cluster interpolation is an extension of the ENI method with rigid-clusters replacing point masses. Now the intermediate conformations in an anharmonic pathway can be determined by the translational and rotational displacements of large clusters in such a way that distance constraints are observed. We present the derivation of the rigid-cluster model and apply it to a variety of macromolecular assemblies. Rigid-cluster ENI is then modified for a hybrid model represented by a mixture of rigid clusters and point masses. Simulation results show that both rigid-cluster and hybrid ENI methods generate sterically feasible pathways of large systems in a very short time. For example, the HK97 virus capsid is an icosahedral symmetric assembly composed of 60 identical asymmetric units. Its original Hessian matrix size for a Calpha coarse-grained model is >(300,000)(2). However, it reduces to (84)(2) when we apply the rigid-cluster model with icosahedral symmetry constraints. The computational cost of the interpolation no longer scales heavily with the size of structures; instead, it depends strongly on the minimal number of rigid clusters into which the system can be decomposed.
| S-EPMC1366543 | biostudies-literature
Project description:ObjectivesTo apply contact measurement and reference point system (RPS) alignment techniques to establish a method for 3D reconstruction of the edentulous jaw models with centric relation and to quantitatively evaluate its accuracy.MethodsUpper and lower edentulous jaw models were clinically prepared, 10 pairs of resin cylinders with same size were adhered to axial surfaces of upper and lower models. The occlusal bases and the upper and lower jaw models were installed in the centric relation position. Faro Edge 1.8m was used to directly obtain center points of the base surface of the cylinders (contact method). Activity 880 dental scanner was used to obtain 3D data of the cylinders and the center points were fitted (fitting method). 3 pairs of center points were used to align the virtual model to centric relation. An observation coordinate system was interactively established. The straight-line distances in the X (horizontal left/right), Y (horizontal anterior/posterior), and Z (vertical) between the remaining 7 pairs of center points derived from contact method and fitting method were measured respectively and analyzed using a paired t-test.ResultsThe differences of the straight-line distances of the remaining 7 pairs of center points between the two methods were X: 0.074 ± 0.107 mm, Y: 0.168 ± 0.176 mm, and Z: -0.003± 0.155 mm. The results of paired t-test were X and Z: p >0.05, Y: p <0.05.ConclusionBy using contact measurement and the reference point system alignment technique, highly accurate reconstruction of the vertical distance and centric relation of a digital edentulous jaw model can be achieved, which meets the design and manufacturing requirements of the complete dentures. The error of horizontal anterior/posterior jaw relation was relatively large.
| S-EPMC4319956 | biostudies-literature
Project description:BackgroundThe accuracy of point-of-care blood glucose (BG) meters is important for the detection of dysglycemia, calculation of insulin doses, and the calibration of continuous glucose monitors. The objective of this study was to compare the accuracy of commercially available glucose meters in a challenging laboratory study using samples with a wide range of reference BG and hemoglobin values.MethodsFresh, discarded blood samples from a hospital STAT laboratory were either used without modification, spiked with a glucose solution, or incubated at 37°C to produce 347 samples with an even distribution across reference BG levels from 20 to 440 mg/dl and hemoglobin values from 9 to 16 g/dl. We measured the BG of each sample with 17 different commercially available glucose meters and the reference method (YSI 2300) at the same time. We determined the mean absolute relative difference (MARD) for each glucose meter, overall and stratified by reference BG and by hemoglobin level.ResultsThe accuracy of different meters widely, exhibiting a range of MARDs from 5.6% to 20.8%. Accuracy was lower in the hypoglycemic range, but was not consistently lower in samples with anemic blood hemoglobin levels.ConclusionsThe accuracy of commercially available glucose meters varies widely. Although the sample mix in this study was much more challenging than those that would be collected under most use conditions, some meters were robust to these challenges and exhibited high accuracy in this setting. These data on relative accuracy and robustness to challenging samples may be useful in informing the choice of a glucose meter.
| S-EPMC5505415 | biostudies-literature