Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Compaction simulator”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 145 records · Page 8Linked to original sources

Structural basis for the unusual properties of 2',5' nucleic acids and their complexes with RNA and DNA.

To provide insights into the unusual properties of 2',5' nucleic acids (iso nucleic acids), that includes their rejection by Nature as information molecules, modeling studies have been carried out to examine if they indeed possess the stereochemical ability to form helical duplexes and triplexes, just as their 3',5' linked constitutional isomers. The results show that the formation of helical duplexes with 2',5' linkages demands a mandatory displacement of the Watson and Crick base pairs from the helical axis, as a direct consequence of the lateral shift of the sugar-phosphate backbone from the periphery towards the interior of the helix. Thus, both duplexes and triplexes formed with a 2',5'-sugar-phosphate backbone possess this intrinsic trait, manifested normally only in A type duplexes of DNA and RNA. It was found that only a 10-fold symmetric parallel triplex with isomorphous T.AT triplets is stereochemically favorable for isoDNA with 'extended' nucleotide repeats, unlike the 12-fold symmetric triplex favored by DNA. The wider nature of a 12-fold triplex, concomitant with mandatory slide requirement for helix formation in isoDNA, demands even larger displacement, especially with 'extended' nucleotide structural repeats, thereby violating symmetry. However, a symmetric triplex possessing higher twist, can be naturally formed for isoDNA with a 'compact' nucleotide repeat. Two nanosecond molecular dynamics simulation of a 2',5'-B DNA duplex, formed with an intrinsic base pair displacement of -3.3 A, does not seem to favor a total transition to a typical A type duplex, although enhanced slide, X-displacement, decrease in helical rise and narrowing of the major groove during simulation seem to indicate a trend. Modeling of the interaction between the chimeric isoDNA.RNA duplex and E. coli RNase H has provided a structural basis for the inhibitory action of the enzyme. Interaction of residues Gln 80, Trp 81, Asn 16 and Lys 99, of E. coli RNase H with DNA of the DNA.RNA hybrid, are lost when the DNA backbone is replaced by isoDNA. Based on modeling and experimental observations, it is argued that 2',5' nucleic acids possess restricted conformational flexibility for helical polymorphism. The inability of isoDNA to favor the biologically relevant B form duplex and the associated topological inadequacies related to nucleic acid compaction and interactions with regulatory proteins may be some of the factors that might have led to the rejection of 2',5' links.

Base Pairing↗

Side-chain conformational disorder in a molten globule: molecular dynamics simulations of the A-state of human alpha-lactalbumin.

Molten globules are compact, partially folded forms of proteins consisting of an ensemble of interconverting conformers with disorder in the side-chain packing across the ensemble. Using insights from experimental data a strategy has been devised to describe this side-chain disorder in a molten globule. Structures of human alpha-lactalbumin have been generated with significantly different side-chain packings to those observed in the native protein by changing all the chi1 torsion angles. Molecular dynamics simulations have been carried out starting from these structures using conditions under which the molten globule is seen experimentally, i.e. at low pH and in the absence of calcium. In each simulation, low energy conformers were generated which are compact and retain substantial secondary structure and an overall native fold, but have very different side-chain contacts; there are all-atom root-mean-square differences of up to 0.5 nm between the different structures. The structures from the different simulations taken together provide an initial description of possible contributors to the ensemble of conformers in the molten globule state which is consistent, at least in general terms, with experimental data.

Calcium↗

Molecular dynamics simulation of protein denaturation: solvation of the hydrophobic cores and secondary structure of barnase.

The transition in barnase from the native state to a compact globule has been studied with high-temperature molecular dynamics simulations. A partial destruction of the alpha-helices and the outer strands of the beta-sheet is observed with water molecules replacing the hydrogen bonds of the secondary structural elements. Simultaneously, the main alpha-helix moves away from the beta-sheet and exposes the principal hydrophobic core, many of whose nonpolar side chains, beginning with the ones near the surface, become solvated by hydrogen-bonded water molecules. This step involves a significant increase in the solvent-exposed surface area; the resulting loss of stability due to the hydrophobic effect may be the major source of the activation barrier in the unfolding reaction. The detailed mechanism described here for the first stage of the denaturation of barnase, including the essential role of water molecules, is likely to be representative of protein denaturation, in general.

Bacillus↗

Molecular-dynamics simulation of the beta domain of metallothionein with a semi-empirical treatment of the metal core.

The three-metal-containing beta domain of rat liver metallothionein-2 in aqueous solution was simulated with different metal contents. The Cd(3), the CdZn(2), and the Zn(3) variant were investigated using a conventional molecular dynamics simulation, as well as a simulation with a semi-empirical quantum-chemical description (MNDO and MNDO/d) of the metal core embedded in a classical environment. For the purely classical simulations, the standard GROMOS96 force-field parameters were used, and parameters were estimated for cadmium. The results of both kinds of simulations were compared to each other and to the corresponding experimental X-ray crystallographic and NMR solution data. The purely classical simulations were found to produce a too compact metal cluster with partially incorrect geometries, which affected the enfolding protein backbone structure. The inclusion of MNDO/d for the treatment of the metal cluster improved the results to give correct cluster geometries and an overall protein structure in agreement with the experiment. The metal cluster and the cysteine residues bound to it are structurally stable, while the irregular polypeptide backbone loops between the cysteines exhibit a considerable flexibility. MNDO without extension to d orbitals failed to maintain the structure of the metal core.

Animals↗

Unfolding and refolding of the native structure of bovine pancreatic trypsin inhibitor studied by computer simulations.

A new procedure for studying the folding and unfolding of proteins, with an application to bovine pancreatic trypsin inhibitor (BPTI), is reported. The unfolding and refolding of the native structure of the protein are characterized by the dimensions of the protein, expressed in terms of the three principal radii of the structure considered as an ellipsoid. A dynamic equation, describing the variations of the principal radii on the unfolding path, and a numerical procedure to solve this equation are proposed. Expanded and distorted conformations are refolded to the native structure by a dimensional-constraint energy minimization procedure. A unique and reproducible unfolding pathway for an intermediate of BPTI lacking the [30,51] disulfide bond is obtained. The resulting unfolded conformations are extended; they contain near-native local structure, but their longest principal radii are more than 2.5 times greater than that of the native structure. The most interesting finding is that the majority of expanded conformations, generated under various conditions, can be refolded closely to the native structure, as measured by the correct overall chain fold, by the rms deviations from the native structure of only 1.9-3.1 A, and by the energy differences of about 10 kcal/mol from the native structure. Introduction of the [30,51] disulfide bond at this stage, followed by minimization, improves the closeness of the refolded structures to the native structure, reducing the rms deviations to 0.9-2.0 A. The unique refolding of these expanded structures over such a large conformational space implies that the folding is strongly dictated by the interactions in the amino acid sequence of BPTI. The simulations indicate that, under conditions that favor a compact structure as mimicked by the volume constraints in our algorithm, the expanded conformations have a strong tendency to move toward the native structure; therefore, they probably would be favorable folding intermediates. The results presented here support a general model for protein folding, i.e., progressive formation of partially folded structural units, followed by collapse to the compact native structure. The general applicability of the procedure is also discussed.

Animals↗

Feature preserving refinement of surfaces for web-based surgical simulation.

In plastic surgery, 3D models of the affected part are often used for the purpose of visualizations and surgical simulations. The optimal models for web-based surgical simulations keep high accuracy in affected parts and keep low accuracy in other parts. Consequently, the data size becomes small. In this research, we propose a method to generate free-form surfaces based on Lattice Structure from polygonal meshes. The polygonal meshes are generated automatically from CT and MRI data using Marching Cubes. By changing the resolution of input images, the accuracy of output meshes is controlled. Free-form surfaces based on Lattice Structure are fitted to polygonal meshes. Lattice Structure is a method to manage a free-form surface with a simple base polygonal mesh. The data size is quite small because surface shape is converted and saved as a simple polygon. A free-form surface is quickly generated and high accuracy is maintained. Moreover, users can input character lines and they are reflected as boundaries of patches. The models generated with this method are partly accurate and compact. These data make it possible to simulate surgery on the WWW, because they can be quickly transferred

Computer Graphics↗

[Impact of localized compaction and ridge fertilization on field nitrate transport and nitrate use efficiency].

The primary season of nitrate leaching in the Loess Plateau region is the monsoon, which is caused by heavy rainfall during the growth season of corn (Zea mays L.). Nitrate leaching to groundwater is an increasing concern in agriculture, and is one of the major nitrogen losing ways in dryland farming system. Localized compaction and ridge fertilization is a method for nitrogen fertilizer application, by which, less fertilizer leaching would occur. The NO3(-)-N transport in soil profile, corn yields and nitrogen use efficiency under localized compaction and ridge fertilization were investigated through two years field study. The factors that affect NO3(-)-N transport under localized compaction and ridge fertilization were studied, combined with simulated experiment. The results showed that NO3(-)-N was leached to below 90 cm in plat fertilization in the year of about 370 mm rainfall, a mean precipitation during the season, while the NO3(-)-N leakage of the fertilizer zone was reduced by localized compaction and ridge fertilization, as a result that the NO3(-)-N concentration below 60 cm was less than 10 mg.kg-1, and NO3(-)-N accumulated in 20-40 cm with a concentration 80-90 mg.kg-1. There was no significant difference in yield between application methods with 240.0 kg N.hm-2. However, the absorbed amount of nitrogen was improved significantly by localized compaction and ridge fertilization, and the nitrogen use efficiency was increased by 9%. The bulk density of the barriers had an evident effect on NO3(-)-N transport under localized compaction and ridge fertilization, but the effect of ridge slope was insignificant.

Biological Transport↗

RNA tertiary structure of the HIV RRE domain II containing non-Watson-Crick base pairs GG and GA: molecular modeling studies.

We have used molecular modeling techniques to model the RNA tertiary structure of the viral RNA element (referred to as domain II of Rev responsive element, RRE) bound by the Rev protein of HIV. In this study, the initial three-dimensional model was built from its established RNA secondary structure, including three non-Watson-Crick G:G, G:A and G:U base pairs. Molecular dynamics (MD) simulations were performed with hydrated or unhydrated sodium ions. Our results indicate that the non-Watson-Crick base pairs in the simulation with unhydrated sodium ions and water are more stable than those with hydrated sodium ions only. The RNA can maintain its compact double helical structure throughout the course of the MD simulations with water and unhydrated sodium ions, although the non-Watson-Crick base pairs and two bulge loops show much more flexibility and conformational distortion than the classical RNA helical region. The distinct distortion of the sugar-phosphate backbone significantly widens the RNA major groove so that the major groove is readily accessible for hydrogen bonding by specific Rev binding. This model emphasizes the importance of specific hydrogen bonding in the stabilization of the three-dimensional structure of the HIV Rev core binding element, not only between the nucleotide bases, but also among the ribose hydroxyls, phosphate anionic oxygens, base oxygens and nitrogens, and bridging water molecules. Moreover, our results suggest that sodium ions play an important role in the formation of base pairs G:G and G:A of the RRE by a manner similar to the arginine of the Rev-RRE complex.

Adenine↗

Morphology and mechanical properties of surfactant aggregates at water-silica interfaces: molecular dynamics simulations.

Dilute and concentrated surfactant systems at the solid-liquid interface are examined using classical molecular dynamics simulations. Particular emphasis is placed on understanding how surfactants aggregate and form the micellar structure, how micelles change shape at high concentrations in aqueous media and in the presence of hydrophilic surfaces, and at what force this micellar structure breaks apart during indentation of micelle-covered surfaces with a proximal probe microscope tip. The specific system of interest is C12TAB (n-dodecyltrimethylammonium bromide) surfactant in an aqueous medium that is modeled with empirical potentials. The simulations predict that the micelle structure in water is compact and either spherical or elliptical in shape. In the presence of a hydrophilic surface of silica, the structure evolves into a flat elliptical shape, in agreement with experimental findings. The simulated indentation of the micelle/silica system causes the micelle to break apart at an indentation force of about 1 nN and form a surfactant monolayer. The predicted force curve is in excellent agreement with experimental measurements.

Journal Article↗

Two-state expansion and collapse of a polypeptide.

The initial phase of folding for many proteins is presumed to be the collapse of the polypeptide chain from expanded to compact, but still denatured, conformations. Theory and simulations suggest that this collapse may be a two-state transition, characterized by barrier-crossing kinetics, while the collapse of homopolymers and random heteropolymers is continuous and multi-phasic. A new rapid-mixing flow technique has been used to resolve the late stages of polypeptide collapse, at time scales >/=45 microseconds. We have used a laser temperature-jump with fluorescence spectroscopy to resolve the complete time-course of the collapse of denatured cytochrome c with nanosecond time resolution. We find the process to be exponential in time and thermally activated, with an apparent activation energy approximately 9 k(B)T (after correction for solvent viscosity). These results indicate that polypeptide collapse is kinetically a two-state transition. Because of the observed free energy barrier, the time scale of polypeptide collapse is dramatically slower than is predicted by Langevin models for homopolymer collapse.

Animals↗

Two-state expansion and collapse of a polypeptide.

The initial phase of folding for many proteins is presumed to be the collapse of the polypeptide chain from expanded to compact, but still denatured, conformations. Theory and simulations suggest that this collapse may be a two-state transition, characterized by barrier-crossing kinetics, while the collapse of homopolymers is continuous and multi-phasic. We have used a laser temperature-jump with fluorescence spectroscopy to measure the complete time-course of the collapse of denatured cytochrome c with nanosecond time resolution. We find the process to be exponential in time and thermally activated, with an apparent activation energy approximately 9 k(B)T (after correction for solvent viscosity). These results indicate that polypeptide collapse is kinetically a two-state transition. Because of the observed free energy barrier, the time scale of polypeptide collapse is dramatically slower than is predicted by Langevin models for homopolymer collapse.

Animals↗

Exact discrete compactlike traveling kinks and pulses in phi(4) nonlinear lattices.

We show that by properly choosing the analytical form of a solitary wave solution of discrete phi(4) models we can calculate the parameters of the potential which allow the propagation of compact (kink and pulses) solutions. Our numerical simulations show that narrow kinks and pulses with finite extent can propagate freely, and that discrete breathers with finite but long lifetime, can emerge from their collisions. Moreover, our numerical simulations reveal that the propagation of two successive pulses at a relative distance of two lattice spacings propagate freely, i.e., without interaction.

Journal Article↗

Renormalization group study of a kinetically constrained model for strong glasses.

We derive a dynamic field theory for a kinetically constrained model, based on the Fredrickson-Andersen model, which we expect to describe the properties of an Arrhenius (strong) supercooled liquid at the coarse-grained level. We study this field theory using the renormalization group. For mesoscopic length and time scales, and for space dimension d>/=2 , the behavior of the model is governed by a zero-temperature dynamical critical point in the directed percolation universality class. We argue that in d=1 its behavior is that of compact directed percolation. We perform detailed numerical simulations of the corresponding Fredrickson-Andersen model on the lattice in various dimensions, and find reasonable quantitative agreement with the field theory predictions.

Journal Article↗

N-body simulations of galaxies and groups of galaxies with the Marseille GRAPE systems.

I review the Marseille GRAPE systems and the N-body simulations done with them. First I briefly describe the available hardware and software, their possibilities and their limitations. I then describe work done on interacting galaxies and groups of galaxies. This includes simulations of the formation of ring galaxies, simulations of bar destruction by a massive compact satellite, of merging in compact groups and of the formation of brightest members in clusters of galaxies.

Journal Article↗

Noninvasive localization of electromagnetic epileptic activity. I. Method descriptions and simulations.

This paper considers the solution of the bioelectromagnetic inverse problem with particular emphasis on focal compact sources that are likely to arise in epileptic data. Two linear inverse methods are proposed and evaluated in simulations. The first method belongs to the class of distributed inverse solutions, capable of dealing with multiple simultaneously active sources. This solution is based on a Local Auto Regressive Average (LAURA) model. Since no assumption is made about the number of activated sources, this approach can be applied to data with multiple sources. The second method, EPIFOCUS, assumes that there is only a single focal source. However, in contrast to the single dipole model, it allows the source to have a spatial extent beyond a single point and avoids the non-linear optimization process required by dipole fitting. The performance of both methods is evaluated with synthetic data in noisy and noise free conditions. The simulation results demonstrate that LAURA and EPIFOCUS increase the number of sources retrieved with zero dipole localization error and produce lower maximum error and lower average error compared to Minimum Norm, Weighted Minimum Norm and Minimum Laplacian (LORETA). The results show that EPIFOCUS is a robust and powerful tool to localize focal sources. Alternatives to localize data generated by multiple sources are discussed. A companion paper (Lantz et al. 2001, this issue) illustrates the application of LAURA and EPIFOCUS to the analysis of interictal data in epileptic patients.

Artifacts↗

Bioadhesive dosage form for peroral administration of timolol base.

Timolol base was prepared from its maleate salt and checked for purity, pKa (9.03) and n-octanol/phosphate buffer pH 6.6 partition coefficient (1.72). The rate of swelling of sodium carboxymethylcellulose, Carbopol 934 (CP) and hydroxypropylcellulose (HPC) was examined prior to use in bioadhesive compacts with the model drug methylene blue to study their influence on device integrity and drug diffusion. A final compact containing a core of timolol base and Precirol, a bioadhesive layer of CP and HPC, and a cap of magnesium stearate gave sustained release of the drug in simulated saliva pH 6.6. After preliminary evaluation in dogs, the compact was evaluated in a panel of humans in whom it was shown that the flux of drug could be increased from 70 to 127 micrograms mm-2 h-1 by inclusion of the penetration enhancer sodium lauryl-sulphate into the core formulation.

Adhesiveness↗

Protein self-association in the cell: a mechanism for fine tuning the level of macromolecular crowding?

A new role for protein self-association in the cell is discussed. An argument is advanced that when cellular protein is in its associated state the excluded volume component of the solution is minimized. Conversely, when cellular protein is in its dissociated state the excluded volume component of the solution is maximized. For proteins that make up a substantial fraction of the intracellular protein concentration, control of the self-association event thus presents itself as a means of regulating cellular processes that are influenced by different levels of volume exclusion. In this communication we examine how the control of protein association/dissociation might influence one such important process, namely the folding of a protein to a compact state.

Computer Simulation↗

Collisional cooling of large ions in electrospray mass spectrometry.

Collisional cooling of ions in the rf-only multipole guides has become a method of choice for coupling electrospray sources to various mass analyzers. Normally parameters of such ion guides (length, pressure) provide enough thermalization and focusing for ions in a wide mass range. Noncovalent complexes, however, have more compact conformations than denatured biomolecules of similar mass and, therefore may not be transmitted efficiently through standard ion guides, as demonstrated by theoretical analysis, simulations, and experiments. Several methods of improving collisional cooling for large compact ions have been developed on a quadrupole time-of-flight instrument, which include operating the ion guides at higher pressure and trapping ions to increase the cooling time. Improved transmission of heavy ions obtained with those methods is studied in experiments with proteasome 20S, an oligomeric protein noncovalent complex with molecular weight around 692,000, and a few other compounds.

Journal Article↗