An Efficient and Accurate Formalism for the Treatment of Large Amplitude Intramolecular Motion.
Ontology highlight
ABSTRACT:
SUBMITTER: Reinisch G
PROVIDER: S-EPMC3419465 | biostudies-other | 2012 Aug
REPOSITORIES: biostudies-other
Similar Datasets
Project description:A new hybrid-model fitting program for methylamine-like molecules has been developed, on the basis of an effective Hamiltonian in which the ammonia-like inversion motion is treated using a tunneling formalism, whereas the internal-rotation motion is treated using an explicit kinetic energy operator and potential energy function. The Hamiltonian in the computer program is set up as a 2 × 2 partitioned matrix, where each diagonal block contains a traditional torsion-rotation Hamiltonian (as in the earlier program BELGI), and the two off-diagonal blocks contain tunneling terms. This hybrid formulation permits the use of the permutation-inversion group G6 (isomorphic to C(3v)) for terms in the two diagonal blocks but requires G12 for terms in the off-diagonal blocks. The first application of the new program is to 2-methylmalonaldehyde. Microwave data for this molecule were previously fit using an all-tunneling Hamiltonian formalism to treat both large-amplitude motions. For 2-methylmalonaldehyde, the hybrid program achieves the same quality of fit as was obtained with the all-tunneling program, but fits with the hybrid program eliminate a large discrepancy between internal rotation barriers in the OH and OD isotopologs of 2-methylmalonaldehyde that arose in fits with the all-tunneling program. This large isotopic shift in internal rotation barrier is thus almost certainly an artifact of the all-tunneling model. Other molecules for application of the hybrid program are mentioned.
| S-EPMC4758688 | biostudies-literature
Project description:PurposeIntensity-modulated proton therapy (IMPT) for lung tumors with a large tumor movement is challenging due to loss of robustness in the target coverage. Often an upper cut-off at 5-mm tumor movement is used for proton patient selection. In this study, we propose (1) a robust and easily implementable treatment planning strategy for lung tumors with a movement larger than 5 mm, and (2) a four-dimensional computed tomography (4DCT) robust evaluation strategy for evaluating the dose distribution on the breathing phases.Materials and methodsWe created a treatment planning strategy based on the internal target volume (ITV) concept (aim 1). The ITV was created as a union of the clinical target volumes (CTVs) on the eight 4DCT phases. The ITV expanded by 2 mm was the target during robust optimization on the average CT (avgCT). The clinical plan acceptability was judged based on a robust evaluation, computing the voxel-wise min and max (VWmin/max) doses over 28 error scenarios (range and setup errors) on the avgCT. The plans were created in RayStation (RaySearch Laboratories, Stockholm, Sweden) using a Monte Carlo dose engine, commissioned for our Mevion S250i Hyperscan system (Mevion Medical Systems, Littleton, MA, USA). We developed a new 4D robust evaluation approach (4DRobAvg; aim 2). The 28 scenario doses were computed on each individual 4DCT phase. For each scenario, the dose distributions on the individual phases were deformed to the reference phase and combined to a weighted sum, resulting in 28 weighted sum scenario dose distributions. From these 28 scenario doses, VWmin/max doses were computed. This new 4D robust evaluation was compared to two simpler 4D evaluation strategies: re-computing the nominal plan on each individual 4DCT phase (4DNom) and computing the robust VWmin/max doses on each individual phase (4DRobInd). The treatment planning and dose evaluation strategies were evaluated for 16 lung cancer patients with tumor movement of 4-26 mm.ResultsThe ratio of the ITV and CTV volumes increased linearly with the tumor amplitude, with an average ratio of 1.4. Despite large ITV volumes, a clinically acceptable plan fulfilling all target and organ at risk (OAR) constraints was feasible for all patients. The 4DNom and 4DRobInd evaluation strategies were found to under- or overestimate the dosimetric effect of the tumor movement, respectively. 4DRobInd showed target underdosage for five patients, not observed in the robust evaluation on the avgCT or in 4DRobAvg. The accuracy of dose deformation used in 4DRobAvg was quantified and found acceptable, with differences for the dose-volume parameters below 1 Gy in most cases.ConclusionThe proposed ITV-based planning strategy on the avgCT was found to be a clinically feasible approach with adequate tumor coverage and no OAR overdosage even for large tumor movement. The new proposed 4D robust evaluation, 4DRobAvg, was shown to give an easily interpretable understanding of the effect of respiratory motion dose distribution, and to give an accurate estimate of the dose delivered in the different breathing phases.
| S-EPMC8456954 | biostudies-literature
Project description:BackgroundHaplotypes extracted from human DNA can be used for gene mapping and other analysis of genetic patterns within and across populations. A fundamental problem is, however, that current practical laboratory methods do not give haplotype information. Estimation of phased haplotypes of unrelated individuals given their unphased genotypes is known as the haplotype reconstruction or phasing problem.ResultsWe define three novel statistical models and give an efficient algorithm for haplotype reconstruction, jointly called HaploRec. HaploRec is based on exploiting local regularities conserved in haplotypes: it reconstructs haplotypes so that they have maximal local coherence. This approach--not assuming statistical dependence for remotely located markers--has two useful properties: it is well-suited for sparse marker maps, such as those used in gene mapping, and it can actually take advantage of long maps.ConclusionOur experimental results with simulated and real data show that HaploRec is a powerful method for the large scale haplotyping needed in association studies. With sample sizes large enough for gene mapping it appeared to be the best compared to all other tested methods (Phase, fastPhase, PL-EM, Snphap, Gerbil; simulated data), with small samples it was competitive with the best available methods (real data). HaploRec is several orders of magnitude faster than Phase and comparable to the other methods; the running times are roughly linear in the number of subjects and the number of markers. HaploRec is publicly available at http://www.cs.helsinki.fi/group/genetics/haplotyping.html.
| S-EPMC1766938 | biostudies-literature
Project description:The exciting applications of molecular motion are still limited and are in urgent pursuit, although some fascinating concepts such as molecular motors and molecular machines have been proposed for years. Utilizing molecular motion in a nanoplatform for practical application has been scarcely explored due to some unconquered challenges such as how to achieve effective molecular motion in the aggregate state within nanoparticles. Here, we introduce a class of near infrared-absorbing organic molecules with intramolecular motion-induced photothermy inside nanoparticles, which enables most absorbed light energy to dissipate as heat. Such a property makes the nanoparticles a superior photoacoustic imaging agent compared to widely used methylene blue and semiconducting polymer nanoparticles and allow them for high-contrast photoacoustic imaging of tumours in live mice. This study not only provides a strategy for developing advanced photothermal/photoacoustic imaging nanoagents, but also enables molecular motion in a nanoplatform to find a way for practical application.
| S-EPMC6377612 | biostudies-literature
Project description:Establishing methods for controlling aspects of large amplitude submolecular movements is a prerequisite for the development of artificial devices that function through rotary motion at the molecular level. Here we demonstrate that the rate of rotation of the interlocked components of fumaramide-derived [2]rotaxanes can be accelerated, by >6 orders of magnitude, by isomerizing them to the corresponding maleamide [2]rotaxanes by using light.
| S-EPMC140866 | biostudies-literature
Project description:The relationship between the amplitude of motion and the accumulation of motion sickness in time is unclear. Here, we investigated this relationship at the individual and group level. Seventeen participants were exposed to four oscillatory motion stimuli, in four separate sessions, separated by at least 1 week to prevent habituation. Motion amplitude was varied between sessions at either 1, 1.5, 2, or 2.5 ms-2. Time evolution was evaluated within sessions applying: an initial motion phase for up to 60 min, a 10-min rest, a second motion phase up to 30 min to quantify hypersensitivity and lastly, a 5-min rest. At both the individual and the group level, motion sickness severity (MISC) increased linearly with respect to acceleration amplitude. To analyze the evolution of sickness over time, we evaluated three variations of the Oman model of nausea. We found that the slow (502 s) and fast (66.2 s) time constants of motion sickness were independent of motion amplitude, but varied considerably between individuals (slow STD = 838 s; fast STD = 79.4 s). We also found that the Oman model with output scaling following a power law with an exponent of 0.4 described our data much better as compared to the exponent of 2 proposed by Oman. Lastly, we showed that the sickness forecasting accuracy of the Oman model depended significantly on whether the participants had divergent or convergent sickness dynamics. These findings have methodological implications for pre-experiment participant screening, as well as online tuning of automated vehicle algorithms based on sickness susceptibility.
| S-EPMC9126086 | biostudies-literature
Project description:Accurate and efficient genotyping of simple sequence repeats (SSRs) constitutes the basis of SSRs as an effective genetic marker with various applications. However, the existing methods for SSR genotyping suffer from low sensitivity, low accuracy, low efficiency and high cost. In order to fully exploit the potential of SSRs as genetic marker, we developed a novel method for SSR genotyping, named as AmpSeq-SSR, which combines multiplexing polymerase chain reaction (PCR), targeted deep sequencing and comprehensive analysis. AmpSeq-SSR is able to genotype potentially more than a million SSRs at once using the current sequencing techniques. In the current study, we simultaneously genotyped 3105 SSRs in eight rice varieties, which were further validated experimentally. The results showed that the accuracies of AmpSeq-SSR were nearly 100 and 94% with a single base resolution for homozygous and heterozygous samples, respectively. To demonstrate the power of AmpSeq-SSR, we adopted it in two applications. The first was to construct discriminative fingerprints of the rice varieties using 3105 SSRs, which offer much greater discriminative power than the 48 SSRs commonly used for rice. The second was to map Xa21, a gene that confers persistent resistance to rice bacterial blight. We demonstrated that genome-scale fingerprints of an organism can be efficiently constructed and candidate genes, such as Xa21 in rice, can be accurately and efficiently mapped using an innovative strategy consisting of multiplexing PCR, targeted sequencing and computational analysis. While the work we present focused on rice, AmpSeq-SSR can be readily extended to animals and micro-organisms.
| S-EPMC5449614 | biostudies-literature
Project description:The global search stage of crystal structure prediction (CSP) methods requires a fine balance between accuracy and computational cost, particularly for the study of large flexible molecules. A major improvement in the accuracy and cost of the intramolecular energy function used in the CrystalPredictor II [Habgood et al. (2015). J. Chem. Theory Comput. 11, 1957-1969] program is presented, where the most efficient use of computational effort is ensured via the use of adaptive local approximate model (LAM) placement. The entire search space of the relevant molecule's conformations is initially evaluated using a coarse, low accuracy grid. Additional LAM points are then placed at appropriate points determined via an automated process, aiming to minimize the computational effort expended in high-energy regions whilst maximizing the accuracy in low-energy regions. As the size, complexity and flexibility of molecules increase, the reduction in computational cost becomes marked. This improvement is illustrated with energy calculations for benzoic acid and the ROY molecule, and a CSP study of molecule (XXVI) from the sixth blind test [Reilly et al. (2016). Acta Cryst. B72, 439-459], which is challenging due to its size and flexibility. Its known experimental form is successfully predicted as the global minimum. The computational cost of the study is tractable without the need to make unphysical simplifying assumptions.
| S-EPMC5134761 | biostudies-literature
Project description:We present an efficient method for calculating the electronic structure and total energy of strongly correlated electron systems. The method extends the traditional Gutzwiller approximation for one-particle operators to the evaluation of the expectation values of two particle operators in the many-electron Hamiltonian. The method is free of adjustable Coulomb parameters, and has no double counting issues in the calculation of total energy, and has the correct atomic limit. We demonstrate that the method describes well the bonding and dissociation behaviors of the hydrogen and nitrogen clusters, as well as the ammonia composed of hydrogen and nitrogen atoms. We also show that the method can satisfactorily tackle great challenging problems faced by the density functional theory recently discussed in the literature. The computational workload of our method is similar to the Hartree-Fock approach while the results are comparable to high-level quantum chemistry calculations.
| S-EPMC4551991 | biostudies-other
Project description:MotivationAs the number of publically available microarray experiments increases, the ability to analyze extremely large datasets across multiple experiments becomes critical. There is a requirement to develop algorithms which are fast and can cluster extremely large datasets without affecting the cluster quality. Clustering is an unsupervised exploratory technique applied to microarray data to find similar data structures or expression patterns. Because of the high input/output costs involved and large distance matrices calculated, most of the algomerative clustering algorithms fail on large datasets (30,000 + genes/200 + arrays). In this article, we propose a new two-stage algorithm which partitions the high-dimensional space associated with microarray data using hyperplanes. The first stage is based on the Balanced Iterative Reducing and Clustering using Hierarchies algorithm with the second stage being a conventional k-means clustering technique. This algorithm has been implemented in a software tool (HPCluster) designed to cluster gene expression data. We compared the clustering results using the two-stage hyperplane algorithm with the conventional k-means algorithm from other available programs. Because, the first stage traverses the data in a single scan, the performance and speed increases substantially. The data reduction accomplished in the first stage of the algorithm reduces the memory requirements allowing us to cluster 44,460 genes without failure and significantly decreases the time to complete when compared with popular k-means programs. The software was written in C# (.NET 1.1).AvailabilityThe program is freely available and can be downloaded from http://www.amdcc.org/bioinformatics/bioinformatics.aspx.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC2672630 | biostudies-literature