Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Computer Graphics”

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 1,423 records · Page 79Linked to original sources

Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics.

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.

Aquaporins↗

Correlations between internal mobility and stability of globular proteins.

The recent work is surveyed which leads to the suggestions that the conformation of globular proteins in solution corresponds to a dynamic ensemble of rapidly interconverting spatial structures, that clusters of hydrophobic amino acid side chains have an important role in the architecture of protein molecules, and that mechanistic aspects of protein denaturation can be correlated with internal mobility seen in the native conformation. These conclusions resulted originally from high resolution 1H nuclear magnetic resonance (NMR) studies of aromatic ring mobility, exchange of interior amide protons and thermal denaturation of the basic pancreatic trypsin inhibitor and a group of related proteins. Various new approaches to further characterize proteins in solution have now been taken and preliminary data are presented. These include computer graphics to outline hydrophobic clusters in globular protein structures, high resolution 1H-NMR experiments at variable hydrostatic pressure and 13C-NMR relaxation measurements. At the present early stage of these new investigations it appears that the hydrophobic cluster model for globular proteins is compatible with the data obtained.

Amino Acids↗

Recognition and interactions controlling the assemblies of beta barrel domains.

We present a qualitative computer graphics approach to the characterization of forces important to the assembly of beta domains that should have general utility for examining protein interactions and assembly. In our approach, the nature of the molecular surface buried by the domain contacts, the specificity of the residue-to-residue interactions, and the identity of electrostatic, hydrophobic, and hydrophilic interactions are elucidated. These techniques are applied to the beta barrel domains of Cu, Zn superoxide dismutase (SOD), immunoglobulin Fab, and tomato bushy stunt virus coat protein (TBSV), a plant viral capsid protein. By looking at a set of proteins having different numbers of interacting beta domains, we have been able to see some of the variety and also some of the patterns common to these assembled domains. Strong beta domain interactions (identified by their biochemical integrity) are apparently due to chemical, electrostatic, and shape complementarity of the molecular surfaces buried from interaction with solvent molecules. Although the amount of hydrophobic buried surface area appears to correlate with the strength of the interaction, electrostatic forces appear to be important in both stabilizing and destabilizing specific contacts. In TBSV, analysis of electrostatic interactions may help explain mechanisms of subunit accommodation to different environments, particle expansion, and pathways of assembly. The possible molecular basis for observed differences in the stability and flexibility of the domain complexes is discussed.

Capsid↗

Orientation of the protonated retinal Schiff base group in bacteriorhodopsin from absorption linear dichroism.

Linear dichroism experiments are performed on light-adapted bacteriorhodopsin (BR568) films containing native retinal (A1) and its 3,4-dehydroretinal (A2) analogue to measure the angle between the chromophore transition dipole moment and the membrane normal. QCFF/pi calculations show that the angle between the transition moment and the long axis of the polyene is changed by 3.4 degrees when the C3-C4 bond is unsaturated. The difference vector between the two transition moments points in the same direction as the Schiff base (N----H) bond for the all-trans BR568 chromophore. Because the plane of the chromophore is perpendicular to the membrane plane, a comparison of the transition moment orientations in the A1- and A2-pigments enables us to determine the orientation of the N----H bond with respect to the absolute chromophore (N----C5 vector) orientation. The angles of the transition moments are 70.3 degrees +/- 0.4 degrees and 67.8 degrees +/- 0.4 degrees for the A1- and A2-pigments, respectively. The fact that the change in the transition moment angle (2.5 degrees) is close to the predicted 3.4 degrees supports the idea that the chromophore plane is nearly perpendicular to the membrane plane. The decreased transition moment angle in the A2-analogue requires that the N----H bond and the N----C5 vector point toward the same membrane surface. Available results indicate that the N----C5 vector points toward the exterior in BR568. With this assignment, we conclude that the N----H bond points toward the exterior surface and its most likely counterion Asp-212. This information makes possible the construction of a computer graphics model for the active site in BR568.

Bacteriorhodopsins↗

Polymorphism in the assembly of polyomavirus capsid protein VP1.

Polyomavirus major capsid protein VP1, purified after expression of the recombinant gene in Escherichia coli, forms stable pentamers in low-ionic strength, neutral, or alkaline solutions. Electron microscopy showed that the pentamers, which correspond to viral capsomeres, can be self-assembled into a variety of polymorphic aggregates by lowering the pH, adding calcium, or raising the ionic strength. Some of the aggregates resembled the 500-A-diameter virus capsid, whereas other considerably larger or smaller capsids were also produced. The particular structures formed on transition to an environment favoring assembly depended on the pathway of the solvent changes as well as on the final conditions. Mass measurements from cryoelectron micrographs and image analysis of negatively stained specimens established that a distinctive 320-A-diameter particle consists of 24 close-packed pentamers arranged with octahedral symmetry. Comparison of this unexpected octahedral assembly with a 12-capsomere icosahedral aggregate and the 72-capsomere icosahedral virus capsid by computer graphics methods indicates that similar connections are made among trimers of pentamers in these shells of different size. The polymorphism in the assembly of VP1 pentamers can be related to the switching in bonding specificity required to build the virus capsid.

Capsid↗

Characterization of precrystallization aggregation of canavalin by dynamic light scattering.

The aggregation processes leading to crystallization and precipitation of canavalin have been investigated by dynamic light scattering (DLS) in photon correlation spectroscopy (PCS) mode. The sizes of aggregates formed under various conditions of pH, salt concentration, and protein concentrations were deduced from the correlation functions generated by the fluctuating intensity of light scattered by the solutions of the protein. Results obtained indicate that the barrier to crystallization of canavalin is the formation of the trimer, a species that has been characterized by x-ray crystallographic studies (McPherson, A. 1980. J. Biol. Chem. 255:10472-10480). The dimensions of the trimer in solution are in good agreement with those obtained both from the crystal (McPherson, A. 1980. J. Biol. Chem. 255:10472-10480) and from a low angle x-ray scattering study in solution (Plietz, P., P. Damaschun, J. J. Müller, and B. Schlener. 1983. FEBS [Fed. Eur. Biochem. Soc.] Lett. 162:43-46). Furthermore, under conditions known to lead to the formation of rhombohedral crystals of canavalin, a limiting size is reached at high concentrations of canavalin. The size measured corresponds to an aggregate of trimers making a unit rhombohedral cell consistent with x-ray crystallographic data (McPherson, A. 1980. J. Biol. Chem. 255:10472-10480). Presumably, such aggregates are the nuclei from which crystal growth proceeds. The present study was undertaken primarily to test the potential of DLS (PCS) as a tool for rapid, routine screening to determine the ultimate fate of protein solutions (i.e., crystallization or amorphous precipitation) at an early stage, therefore eliminating the need for long-term visual observation. Achieving this goal would constitute amajor advance in the practive of protein crystallization. Delays imposed by visual observation would be considerably reduced, and a more systematic approach could be adopted to select experimental conditions.Our findings with canavalin demonstrate that DLS(PCS) is, indeed, a selective and sensitive probe of precrystallization conditions. Other advantages of this technique include the facts that it is noninvasive, nondestructive,universal, and does not require calibration.

Computer Graphics↗

Effect of geometrical irregularities on propagation delay in axonal trees.

Multiple successive geometrical inhomogeneities, such as extensive arborization and terminal varicosities, are usual characteristics of axons. Near such regions the velocity of the action potential (AP) changes. This study uses AXONTREE, a modeling tool developed in the companion paper for two purposes: (a) to gain insights into the consequence of these irregularities for the propagation delay along axons, and (b) to simulate the propagation of APs along a reconstructed axon from a cortical cell, taking into account information concerning the distribution of boutons (release sites) along such axons to estimate the distribution of arrival times of APs to the axons release sites. We used Hodgkin and Huxley (1952) like membrane properties at 20 degrees C. Focusing on the propagation delay which results from geometrical changes along the axon (and not from the actual diameters or length of the axon), the main results are: (a) the propagation delay at a region of a single geometrical change (a step change in axon diameter or a branch point) is in the order of a few tenths of a millisecond. This delay critically depends on the kinetics and the density of the excitable channels; (b) as a general rule, the lag imposed on the AP propagation at a region with a geometrical ratio GR greater than 1 is larger than the lead obtained at a region with a reciprocal of that GR value; (c) when the electronic distance between two successive geometrical changes (Xdis) is small, the delay is not the sum of the individual delays at each geometrical change, when isolated. When both geometrical changes are with GR greater than 1 or both with GR less than 1, this delay is supralinear (larger than the sum of individual delays). The two other combinations yield a sublinear delay; and (d) in a varicose axon, where the diameter changes frequently from thin to thick and back to thin, the propagation velocity may be slower than the velocity along a uniform axon with the thin diameter. Finally, we computed propagation delays along a morphologically characterized axon from layer V of the somatosensory cortex of the cat. This axon projects mainly to area 4 but also sends collaterals to areas 3b and 3a. The model predicts that, for this axon, areas 3a, 3b, and the proximal part of area 4 are activated approximately 2 ms before the activation of the distal part of area 4.

Action Potentials↗

Proposed mechanisms for binding of apo[a] kringle type 9 to apo B-100 in human lipoprotein[a].

The protein component of human lipoprotein[a] consists primarily of two apolipoproteins, apo[a] and apo B-100, linked through a cystine disulfide(s). In the amino acid sequence of apo bd, Cys4057 located within a plasminogen kringle 4-like repeat sequence (3991-4068) is believed to form a disulfide bond with a specific cysteine residue in apo B-100. Our fluorescence-labeling experiments and molecular modeling studies have provided evidence for possible interactions between this apo[a] kringle type and apo B-100. The fluorescent probe, fluorescein-5-maleimide, was used in parallel experiments to label free sulfhydryl moieties in lipoprotein[a] and low-density lipoprotein (LDL). In apo B-100 of LDL, Cys3734 was labeled with the probe, but this site was not labeled in autologous lipoprotein[a]. The result strongly implicates Cys3734 of apo B-100 as the residue forming the disulfide linkage with Cys4057 of apo[a]. To explore possible noncovalent interactions between apo B-100 and apo[a], the crystallographic coordinates for plasminogen kringle 4 were used to generate molecular models of the apo[a] kringle-repeat sequence (3991-4068, LPaK9), the only plasminogen kringle 4 type repeat in apo[a] having an extra cysteine residue not involved in an intramolecular disulfide bond. The Cys4057 residue (henceforth designated as Cys67 in the LPaK9 sequence) is believed to form an intermolecular disulfide bond with a cysteine of apo B-100. In computer graphics molecular models of LPaK9, Cys67 is located on the surface of the kringle near the lysine ligand binding site. Selected segments of the LDL apo B-100 sequence that contain free sulfhydryl cysteines were subjected to energy minimization and docking with the ligand binding site and adjacent regions of the LPaK9 model. In the docking experiments, apo B-100 segment 3732-3745 (PSCKLDFREIQIYK) displayed the best fit and the largest number of van der Waals contacts with models of LPaK9. Other apo B-100 peptides with sulfhydryl cysteine were found to be less compatible when minimized with this kringle. These results support and extend previously suggested mechanisms for a complex interaction between apo[a] and apo B-100 that involve more than a simple covalent disulfide bond.

Amino Acid Sequence↗

XRayView: a teaching aid for X-ray crystallography.

A software package, XRayView, has been developed that uses interactive computer graphics to introduce basic concepts of x-ray diffraction by crystals, including the reciprocal lattice, the Ewald sphere construction, Laue cones, the wavelength dependence of the reciprocal lattice, primitive and centered lattices and systematic extinctions, rotation photography. Laue photography, space group determination and Laue group symmetry, and the alignment of crystals by examination of reciprocal space. XRayView is designed with "user-friendliness" in mind, using pull-down menus to control the program. Many of the experiences of using real x-ray diffraction equipment to examine crystalline diffraction can be simulated. Exercises are available on-line to guide the users through many typical x-ray diffraction experiments.

Biophysics↗

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes.

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.

Base Sequence↗

Influence of cis double bonds in the sn-2 acyl chain of phosphatidylethanolamine on the gel-to-liquid crystalline phase transition.

We have semisynthesized 19 species of mixed-chain phosphatidylethanolamines (PEs) in which the sn-1 acyl chain is derived from saturated fatty acids with varying chain lengths and the sn-2 acyl chain has different chain lengths but contains 0, 1, and 2 cis double bond(s). The gel-to-liquid crystalline phase transition temperatures (Tm) of lipid bilayers prepared from these 19 mixed-chain PEs were determined calorimetrically. When the Tm values are compared with those of saturated and monounsaturated counterparts, a common Tm profile is observed in the plot of Tm versus the number of cis double bonds. Specifically, a marked stepwise decrease in Tm is detected as the number of cis double bonds in the sn-2 acyl chain of the mixed-chain PE is successively increased from 0 to 1 and then to 2. The large Tm-lowering effect of the acyl chain unsaturation can be attributed to the increase in Gibbs free energy of the gel-state bilayer as a result of weaker lateral chain-chain interactions. In addition, we have applied molecular mechanics calculations to simulate the molecular structure of dienoic mixed-chain C(X):C(Y:2 delta n,n+3)PE in the gel-state bilayer, thus enabling the three independent structural parameters (N, delta C, and LS) to be calculated in terms of X, Y, and n, which are intrinsic quantities of C(X):C(Y:2 delta n,n+3)PE. When the Tm values and the corresponding N and delta C values of all dienoic mixed-chain PEs under study are first codified and then analyzed statistically by multiple regressions, the dependence of Tm on the structural parameters can be described quantitatively by a simple and general equation. The physical meaning and the usefulness of this simple and general equation are explained.

Calorimetry, Differential Scanning↗

Homology modeling of cephalopod lens S-crystallin: a natural mutant of sigma-class glutathione transferase with diminished endogenous activity.

The soluble S-crystallin constitutes the major lens protein in cephalopods. The primary amino acid sequence of S-crystallin shows an overall 41% identity with the digestive gland sigma-class glutathione transferase (GST) of cephalopod. However, the lens S-crystallin fails to bind to the S-hexylglutathione affinity column and shows very little GST activity in the nucleophilic aromatic substitution reaction between GSH and 1-chloro-2,4-dinitrobenzene. When compared with other classes of GST, the S-crystallin has an 11-amino acid residues insertion between the conserved alpha4 and alpha5 helices. Based on the crystal structure of squid sigma-class GST, a tertiary structure model for the octopus lens S-crystallin is constructed. The modeled S-crystallin structure has an overall topology similar to the squid sigma-class GST, albeit with longer alpha4 and alpha5 helical chains, corresponding to the long insertion. This insertion, however, makes the active center region of S-crystallin to be in a more closed conformation than the sigma-class GST. The active center region of S-crystallin is even more shielded and buried after dimerization, which may explain for the failure of S-crystallin to bind to the immobilized-glutathione in affinity chromatography. In the active site region, the electrostatic potential surface calculated from the modeled structure is quite different from that of squid GST. The positively charged environment, which contributes to stabilize the negatively charged Meisenheimer complex, is altered in S-crystallin probably because of mutation of Asn99 in GST to Asp101 in S-crystallin. Furthermore, the important Phe106 in authentic GST is changed to His108 in S-crystallin. Combining the topological differences as revealed by computer graphics and sequence variation at these structurally relevant residues provide strong structural evidences to account for the much decreased GST activity of S-crystallin as compared with the authentic GST of the digestive gland.

Algorithms↗

Representation and searching of carbohydrate structures using graph-theoretic techniques.

This paper describes how the carbohydrate structures in the Complex Carbohydrate Structure Database (CCSD) can be represented by labelled graphs, in which the nodes and edges of a graph are used to denote the residues and the inter-residue linkages, respectively, of a carbohydrate. These graph representations are then searched with a subgraph-isomorphism algorithm. We describe the use of one such algorithm, that due to Ullmann, and demonstrate that it provides a very precise way of searching the structures in CCSD. We also describe the use of screening techniques that can eliminate many of the CCSD structures from the subgraph-isomorphism search, with a consequent increase in the speed of the search.

Algorithms↗

The binding of CC-1065 to thymidine and deoxyadenosine oligonucleotides and to poly(dA).poly(dT).

In this work, we report on the binding of the novel antitumor agent CC-1065 to poly(dA).poly(dT) and to mixtures of dA and dT oligomers as determined by electronic absorption and circular dichroism (CD) methods. In addition, the DNA binding properties of CC-1065 and its binding mechanism are compared to those of netropsin. CC-1065 binds to the polymer by at least three mechanisms to produce one irreversibly and two reversibly bound species. One reversibly bound species is moderately stable, but in time (days), it converts to the irreversibly bound species. Both of these species bind within the minor groove of the polymer and exhibit intense CC-1065 induced CD spectra. The other reversibly bound species does not acquire an induced CD. CC-1065 forces B-form duplex formation between mixtures of single strand dA and dT oligomers and binds irreversibly to the duplexes without showing the presence of an intermediate, reversibly bound species. The induced CD increases with increasing length of the oligomer, from the 5-mer (barely detectable CD) to the 14-mer (intense CD). The 7-, 10- and 14-mer mixtures bind about 1, between 1 and 2, and between 2 and 3 CC-1065 molecules, respectively. Computer graphic models of the CC-1065-DNA complex show that the covalent adduct of CC-1065 and unreacted CC-1065 can attain the same close van der Waals contacts between adenine C2 hydrogens and antibiotic CH groups that were observed in the crystal structure of the netropsin-DNA complex. These contacts may account for the dA-dT base pair binding specificity of CC-1065 and for the stability of the reversibly bound CC-1065 species.

Binding Sites↗

Numerical approach to plug-flow activated sludge reactor kinetics.

In this study, general relationships applicable to a wide variety of microbial-mediated treatment processes are developed using a number of continuous-flow stirred tank reactors-in-series. Simultaneous equations resulting from mass balance on substrate and biomass are solved numerically taking the longitudinal biomass gradient into account. A relationship between substrate and biomass concentrations is developed as a function of hydraulic residence time. Monod kinetics and mathematical models of reactors-in-series are used to represent the actual conditions resulting from varying degrees of axial dispersion and wastewater quality. Dimensionless quantities are used to reduce the number of parameters to be taken into account. Computer techniques are applied to express the results generally. By means of the numerical approach, the variation of the ratio between the hydraulic residence time of a given reactor and that of an equivalent plug-flow tank with the same inlet and outlet conditions is investigated. Results are expressed in the form of graphs to characterize the plug-flow activated sludge systems and to provide a basis of design. Experimental data reported in the literature are also evaluated to demonstrate the cases where the existing classical solution to the problem differs from the actual results obtained from the computer program and the design graphs given.

Algorithms↗

Topical cardiac cooling--computer simulation of myocardial temperature changes.

Topical cardiac cooling (TC) is often used in cardiac surgery. We used a computer simulation to study temperature changes in the heart, especially in the right ventricular wall. A three-dimensional computer heart model, derived from Visible Human Data set, National Library of Medicine was used. The model is made from cubes, with spatial resolution of 1mm. Explicit Finite Different method and temperature diffusion equation were used to calculate new temperatures. Three different simulations were performed and simulated temperatures were drawn on a cross-section of heart model in different colors. The results show that areas not immersed into TC solution are less protected against ischemia. It is important not to rely solely on topical cooling but use appropriate method of myocardial protection technique.

Anatomy, Cross-Sectional↗

Three-dimensional visualization of protein interaction networks.

Protein interaction networks provide us with contextual information within which protein function can be interpreted and will assist many biomedical studies. We have developed a new force-directed layout algorithm for visualizing protein interactions in three-dimensional space. Our algorithm divides nodes into three groups based on their interacting properties: bi-connected sub-graph in the center, terminal nodes at the outermost region, and the rest in between them. Experimental results show that our algorithm efficiently generates a clear and aesthetically pleasing drawing of large-scale protein interaction networks and that it is an order of magnitude faster than other force-directed layouts.

Algorithms↗

Coarse-grain parallel computing for very large scale neural simulations in the NEXUS simulation environment.

We describe a neural simulator designed for simulating very large scale models of cortical architectures. This simulator, NEXUS, uses coarse-grain parallel computing by distributing computation and data onto multiple conventional workstations connected via a local area network. Coarse-grain parallel computing offers natural advantages in simulating functionally segregated neural processes. We partition a complete model into modules with locally dense connections--a module may represent a cortical area, column, layer, or functional entity. Asynchronous data communications among workstations are established through the Network File System, which, together with the implicit modularity, decreases communications overhead, and increases overall performance. Coarse-grain parallelism also benefits from the standardization of conventional workstations and LAN, including portability between generations and vendors.

Computer Communication Networks↗