SEARCH · PubMed Health

Results for “Compaction simulator”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15Linked to original sources

PubMedDate not supplied

Compact and accurate linear and nonlinear autoregressive moving average model parameter estimation using laguerre functions.

A linear and nonlinear autoregressive moving average (ARMA) identification algorithm is developed for modeling time series data. The algorithm uses Laguerre expansion of kernals (LEK) to estimate Volterra-Wiener kernals. However, instead of estimating linear and nonlinear system dynamics via moving average models, as is the case for the Volterra-Wiener analysis, we propose an ARMA model-based approach. The proposed algorithm is essentially the same as LEK, but this algorithm is extended to include past values of the output as well. Thus, all of the advantages associated with using the Laguerre function remain with our algorithm; but, by extending the algorithm to the linear and nonlinear ARMA model, a significant reduction in the number of Laguerre functions can be made, compared with the Volterra-Wiener approach. This translates into a more compact system representation and makes the physiological interpretation of higher order kernels easier. Furthermore, simulation results show better performance of the proposed approach in estimating the system dynamics than LEK in certain cases, and it remains effective in the presence of significant additive measurement noise.

Algorithms↗

PubMed1997

Common structural features of UUCG and UACG tetraloops in very short hairpins determined by UV absorption, Raman, IR and NMR spectroscopies.

Thermodynamic and structural properties of two UNCG tetraloops in very short hairpin octamers, 5'-r(GCUUCGGC)-3' and 5'-r(GCUACGGC)-3', have been studied by means of various physical techniques. Melting profiles of both octamers, obtained from UV absorption spectra taken as a function of temperature, are consistent with a monophasic, progressive and completely reversible order-to-disorder transition and confirm their unusual structural stability (Tm > 51 degrees C). The 1H, 13C and 31P NMR chemical shifts and coupling constants of the UACG loop nucleotides are comparable with those reported previously for UUCG loops, i.e. 2'-endo/anti conformation of the second and third nucleotide of the loop as well as the syn orientation of the ultimate guanine base and the A-type double helical conformation of the hairpin stem. Simulation of quantitative NOESY volumes shows that the UACG octamer adopts a very rigid compact structure which is well represented by an average order parameter of 0.9. Three base-pairs and four additional strong hydrogen bonds are undoubtedly responsible for such limited flexibility. Raman and infrared spectra as a function of temperature reflect the order-to-disorder transition, as well. Vibrational conformational markers in low temperature spectra of both octamers indicate the hairpin structure as the major conformer in aqueous phase. These spectra further support the structural features of most of the nucleotides involved in the tetraloops and clearly demonstrate the structural similarities of the phosphodiester backbone in both hairpins. Consequently, on the basis of all present results, one can deduce that the conformational features of the UUCG and UACG tetraloops seem to be inherent to the UNCG type tetraloops, regardless of either the nature of the tetraloop second base or the stem length.

Base Sequence↗

PubMed1994

[Orthognathic surgery of dentofacial deformities post surgerical repairing of cleft palate and lip.].

10 cases with maxillary development limited post surgical repairing of cleft palate and lip were undergone modified Le Fort I segmental osteotomy,that is,put the maxilla forward without velopharyngeal function change and oro-nasal fistula.By Computer Aided Simulation System for orthognathic surgery(CASSOS),analysis of dental pulp,evaluation of velopharyngeal function,compacting profile pre-to-post operation, authors described that surgical procedure has the satisfactory effectiveness.

English Abstract↗

PubMed2002

Biophysical study of the globular organisation of interphase chromosomes.

The globular model of interphase chromosomes is studied using methods involving the statistical physics of polymers. An interphase chromatid is represented as a flexible chain of structural sub-units, or superdomains (SDs). Each SD is simulated as a number of chromatin fibre loops fixed at a nuclear matrix core. A chain of SDs is further folded in the nucleus in a compact conformation owing to volumetric interactions between SDs. The algorithm used is extended to incorporate the chain anchorage at different points. Excluded volume effects are taken into account in Monte Carlo simulation at both the SD and whole chromosome level. A variety of structures is observed in computer experiments. The simulation results correlate with the available experimental data.

Biophysics↗

PubMed2004

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein.

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.

Protein Folding↗

PubMed2006

Transformation pathways of silica under high pressure.

Network-forming oxides with rigid polyhedral building blocks often possess significant capacity for densification under pressure owing to their open structures. The high-pressure behaviour of these oxides is key to the mechanical properties of engineering materials and geological processes in the Earth's interior. Concurrent molecular-dynamics simulations and first-principles calculations reveal that this densification follows a ubiquitous two-stage mechanism. First, a compact high-symmetry anion sublattice forms, as controlled by strong repulsion between the large oxygen anions, and second, cations redistribute onto the newly created interstices. The same mechanism is observed for two different polymorphs of silica, and in the particular case of cristobalite, is corroborated by the experimental finding of a previously unidentified metastable phase. Our simulations not only clarify the nature of this phase, but also identify its occurrence as key evidence in support of this densification mechanism.

Journal Article↗

PubMed2004

Pressure denaturation of staphylococcal nuclease studied by neutron small-angle scattering and molecular simulation.

We studied the pressure-induced folding/unfolding transition of staphylococcal nuclease (SN) over a pressure range of approximately 1-3 kilobars at 25 degrees C by small-angle neutron scattering and molecular dynamics simulations. We find that applying pressure leads to a twofold increase in the radius of gyration derived from the small-angle neutron scattering spectra, and P(r), the pair distance distribution function, broadens and shows a transition from a unimodal to a bimodal distribution as the protein unfolds. The results indicate that the globular structure of SN is retained across the folding/unfolding transition although this structure is less compact and elongated relative to the native structure. Pressure-induced unfolding is initiated in the molecular dynamics simulations by inserting water molecules into the protein interior and applying pressure. The P(r) calculated from these simulations likewise broadens and shows a similar unimodal-to-bimodal transition with increasing pressure. The simulations also reveal that the bimodal P(r) for the pressure-unfolded state arises as the protein expands and forms two subdomains that effectively diffuse apart during initial stages of unfolding. Hydrophobic contact maps derived from the simulations show that water insertions into the protein interior and the application of pressure together destabilize hydrophobic contacts between these two subdomains. The findings support a mechanism for the pressure-induced unfolding of SN in which water penetration into the hydrophobic core plays a central role.

Algorithms↗

PubMed2003

Structure of sphingomyelin bilayers: a simulation study.

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50 degrees C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed1992

Characterization of "native" apomyoglobin by molecular dynamics simulation.

We have used molecular dynamics simulation methods to study the structure and fluctuations of "native" apomyoglobin in aqueous solution for a period of greater than 0.5 nanosecond. This work was motivated by the recent attempts of Hughson et al. to characterize the structure and motion of both this molecule and the less compact, acid stabilized I stage, using methods of pulsed H/2H exchange. The study of these systems provides new insights into protein folding intermediates and our simulation has yielded a detailed model for structure and fluctuations in apomyoglobin which complements the experimental studies. We find that local (short-time) fluctuations agree well with fluctuations observed for the holoprotein in aqueous solution, as well as results from the crystallographic B-factors. In addition, the structural features we observe for native apomyoglobin are very similar to the holoprotein, in basic agreement with the findings of Hughson et al. By examining larger-scale motions, developing only over timescales in excess of a 100 picoseconds, we are able to identify conformationally "labile" and "non-labile" regions within native apomyoglobin. These regions correspond extremely well to those identified in the nuclear magnetic resonance experiments as unstable and stable "folding subdomains" in the I state of apomyoglobin. Overall we find that helices A, B, E, G and H show the least amount of motion and helices C, D and F move substantially over the timescales examined. The major motions, and the primary difference between the holo and apo structures as we have observed them, are due to the shifting motion of helices C, D and F into the vacant heme cavity. We also find that motions at the interface of helical segments can be large, with one important exception being the chain segment connecting helices G and H. This segment of chain interacts with the conformationally "non-labile" helix A to form a relatively rigid subdomain composed of helices A, G and H. We believe that these findings provide direct support for the suggestion of Hughson et al. that helices A, G and H constitute a compact subdomain that remains in a native-like conformation as the protein begins to unfold in environments of decreasing pH.

Apoproteins↗

PubMed1989

A computer model to dynamically simulate protein folding: studies with crambin.

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.

Computer Simulation↗

PubMed2005

A model of chromosome aberration induction: applications to space research.

A mechanistic model and Monte Carlo code simulating chromosome aberration induction in human lymphocytes is presented. The model is based on the assumption that aberrations arise from clustered DNA lesions and that only the free ends of clustered lesions created in neighboring chromosome territories or in the same territory can join and produce exchanges. The lesions are distributed in the cell nucleus according to the radiation track structure. Interphase chromosome territories are modeled as compact intranuclear regions with volumes proportional to the chromosome DNA contents. Both Giemsa staining and FISH painting can be simulated, and background aberrations can be taken into account. The good agreement with in vitro data provides validation of the model in terms of both the assumptions adopted and the simulation techniques. As an application in the field of space research, the model predictions were compared with aberration yields measured among crew members of long-term missions on board Mir and ISS, assuming an average radiation quality factor of 2.4. The agreement obtained also validated the model for in vivo exposure scenarios and suggested possible applications to the prediction of other relevant aberrations, typically translocations.

Chromosome Aberrations↗

PubMed2004

Practicability of protontherapy using compact laser systems.

Protontherapy is a well-established approach to treat cancer due to the favorable ballistic properties of proton beams. Nevertheless, this treatment is today only possible with large scale accelerator facilities which are very difficult to install at existing hospitals. In this article we report on a new approach for proton acceleration up to energies within the therapeutic window between 60 and 200 MeV by using modern, high intensity and compact laser systems. By focusing such laser beams onto thin foils we obtained on target intensities of 6 x 10(19) W/cm2, which is sufficient to produce a well-collimated proton beam with an energy of up to 10 MeV. These results are in agreement with numerical simulations and indicate that proton energies within the therapeutic window should be obtained in the very near future using such economical and very compact laser systems. Hence, this approach could revolutionize cancer treatment by bringing the "lab to the hospital-rather than the hospital to the lab".

Biophysical Phenomena↗

PubMed1998

Calculations on folding of segment B1 of streptococcal protein G.

We present an investigation of the folding thermodynamics and mechanism of segment B1 of streptococcal protein G. Molecular dynamics simulations of the fully solvated protein are used to probe thermodynamically significant states at different stages of folding. We performed several unfolding simulations to generate a database of initial conditions. The database is analyzed and clustered. The cluster centers extracted from this database were then used as starting points for umbrella sampling of the folding free energy landscape under folding conditions. The resulting sampling was combined with the weighted histogram analysis method. One and two-dimensional free energy surfaces were constructed along several order parameters and used to analyze the folding process. Our findings indicate that an initial collapse precedes the formation of significant native structure. Elements of local structure originate in the regions of the protein shown to have higher H/2H exchange protection factors in early stages of folding. A non-native contact, observed experimentally at the N terminus of the alpha-helix in a peptide excised from the protein, is seen to pre-organize the chain in early stages of folding. Collapse and early structure formation yields a compact globule with a significant number of water molecules present. Desolvation of the protein core is coincident with the final stages of folding from the compact state.

Bacterial Proteins↗

PubMed2006

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Algorithms↗

PubMed1999

Additional studies on the distribution of stresses during vertical compaction of gutta-percha in the root canal.

OBJECTIVE: This study was designed to investigate the effect of certain pathological alterations of the dental structures (diminishing bone support, internal resorption, root perforation, periapical lesion) on stress distribution during root canal filling procedures by the warm vertical compaction technique. DESIGN: The computer stress analyses were done for a maxillary canine tooth model which was based on dimensions recovered from a human cadaveric maxilla scanned by CT. METHODS: The finite element method was used to calculate the stresses generated during root canal filling procedures by warm vertical compaction technique. Patterns of stress distribution associated with various alterations in dental structures were investigated. For this purpose 60 cases were simulated. The hypothetical force of 10 N is taken as a unit representation. For other magnitudes of applied force, the corresponding stresses would be scaled directly because the calculations were made for linear materials. RESULTS AND CONCLUSION: It is found that, when diminishing bone support and internal resorption are concurrently simulated, a marked increase in stress magnitudes occur (maximum von Mises stress 5.37 N/mm2). However, these values still remain much below the most frequently reported tensile strength of dentine (50-100 N/mm2). If dentist's handwork is transformed into equivalent edge tractions on gutta-percha, then stresses in dentine, even when they are corrected for 3-kg applied force, appear to remain below fracture strengths of this material. This result leads us to conclude that when warm vertical compaction technique is skillfully performed and inadvertent undue force is not applied, a premature root fracture in a large rooted maxillary anterior tooth with straight root canal anatomy is not likely to occur, even for the unfavourable conditions simulated in our model. This result, like all results derived from modelling applications, is of course contingent upon agreement between the way in which the clinical operations are performed and the way in which they are mirrored for computer representation. We believe that the approach described here avoids the spurious stresses that have been reported in similar investigations.

Computer Simulation↗

PubMed2004

Uniform sampling of steady-state flux spaces: means to design experiments and to interpret enzymopathies.

Reconstruction of genome-scale metabolic networks is now possible using multiple different data types. Constraint-based modeling is an approach to interrogate capabilities of reconstructed networks by constraining possible cellular behavior through the imposition of physicochemical laws. As a result, a steady-state flux space is defined that contains all possible functional states of the network. Uniform random sampling of the steady-state flux space allows for the unbiased appraisal of its contents. Monte Carlo sampling of the steady-state flux space of the reconstructed human red blood cell metabolic network under simulated physiologic conditions yielded the following key results: 1), probability distributions for the values of individual metabolic fluxes showed a wide variety of shapes that could not have been inferred without computation; 2), pairwise correlation coefficients were calculated between all fluxes, determining the level of independence between the measurement of any two fluxes, and identifying highly correlated reaction sets; and 3), the network-wide effects of the change in one (or a few) variables (i.e., a simulated enzymopathy or fixing a flux range based on measurements) were computed. Mathematical models provide the most compact and informative representation of a hypothesis of how a cell works. Thus, understanding model predictions clearly is vital to driving forward the iterative model-building procedure that is at the heart of systems biology. Taken together, the Monte Carlo sampling procedure provides a broadening of the constraint-based approach by allowing for the unbiased and detailed assessment of the impact of the applied physicochemical constraints on a reconstructed network.

Animals↗

PubMed2000

Estimation of synteny conservation and genome compaction between pufferfish (Fugu) and human.

BACKGROUND: Knowledge of the amount of gene order and synteny conservation between two species gives insights to the extent and mechanisms of divergence. The vertebrate Fugu rubripes (pufferfish) has a small genome with little repetitive sequence which makes it attractive as a model genome. Genome compaction and synteny conservation between human and Fugu were studied using data from public databases. METHODS: Intron length and map positions of human and Fugu orthologues were compared to analyse relative genome compaction and synteny conservation respectively. The divergence of these two genomes by genome rearrangement was simulated and the results were compared to the real data. RESULTS: Analysis of 199 introns in 22 orthologous genes showed an eight-fold average size reduction in Fugu, consistent with the ratio of total genome sizes. There was no consistent pattern relating the size reduction in individual introns or genes to gene base composition in either species. For genes that are neighbours in Fugu (genes from the same cosmid or GenBank entry), 40-50% have conserved synteny with a human chromosome. This figure may be underestimated by as much as two-fold, due to problems caused by incomplete human genome sequence data and the existence of dispersed gene families. Some genes that are neighbours in Fugu have human orthologues that are several megabases and tens of genes apart. This is probably caused by small inversions or other intrachromosomal rearrangements. CONCLUSIONS: Comparison of observed data to computer simulations suggests that 4000-16 000 chromosomal rearrangements have occurred since Fugu and human shared a common ancestor, implying a faster rate of rearrangement than seen in human/mouse comparisons.

Animals↗

PubMed2026

Marine air promotes structural compaction and coating growth of soot aerosols after long-range transport from East Asia.

Soot aerosol, a key global warming contributor, undergoes morphological and chemical transformations during atmospheric transport, particularly in humidified marine environments. This study investigates morphology, mixing state, and aging mechanisms of soot particles collected in the Bohai Sea and Yellow Sea. Transmission electron microscopy analyses reveal that coated soot particles dominate the marine atmosphere, accounting for over 98 % of soot-containing particles, with a mean mixing state index (χ) of 0.83. The fractal dimension (Df) of soot particles is 1.84 ± 0.05 in the Northern Yellow Sea, 1.90 ± 0.08 in the Bohai Sea, and 1.96 ± 0.07 in the Southern Yellow Sea, indicating structural compaction during long-range transport. Correspondingly, the average Dp/Dcore ratios (particle to core size ratio) are 5.3 in the Bohai Sea, 4.2 in the Northern Yellow Sea, and 3.9 in the Southern Yellow Sea. Notably, those ratios are higher in marine environments compared to those observed during continental regional transport from northern to southern China (3.54), suggesting enhanced coating growth in humid marine air. The results highlight the important role of marine atmospheres in accelerating soot aging, which in turn leads to significantly stronger light absorption compared to soot in continental air. Our results highlight the necessity of incorporating compact morphologies, uniform mixing states, and thick coatings into optical models for accurate radiative forcing simulations.

Aerosols↗