Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MACROMOLECULAR SYSTEMS”

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 127 records · Page 7Linked to original sources

An efficient hybrid explicit/implicit solvent method for biomolecular simulations.

We present a new hybrid explicit/implicit solvent method for dynamics simulations of macromolecular systems. The method models explicitly the hydration of the solute by either a layer or sphere of water molecules, and the generalized Born (GB) theory is used to treat the bulk continuum solvent outside the explicit simulation volume. To reduce the computational cost, we implemented a multigrid method for evaluating the pairwise electrostatic and GB terms. It is shown that for typical ion and protein simulations our method achieves similar equilibrium and dynamical observables as the conventional particle mesh Ewald (PME) method. Simulation timings are reported, which indicate that the hybrid method is much faster than PME, primarily due to a significant reduction in the number of explicit water molecules required to model hydration effects.

Journal Article↗

A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions.

The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.

Algorithms↗

Automatic domain decomposition of proteins by a Gaussian Network Model.

Proteins are often comprised of domains of apparently independent folding units. These domains can be defined in various ways, but one useful definition divides the protein into substructures that seem to move more or less independently. The same methods that allow fairly accurate calculation of motion can be used to help classify these substructures. We show how the Gaussian Network Model (GNM), commonly used for determining motion, can also be adapted to automatically classify domains in proteins. Parallels between this physical network model and graph theory implementation are apparent. The method is applied to a nonredundant set of 55 proteins, and the results are compared to the visual assignments by crystallographers. Apart from decomposing proteins into structural domains, the algorithm can generally be applied to any large macromolecular system to decompose it into motionally decoupled sub-systems.

Algorithms↗

Proposed structural models of the prothrombinase (FXa-FVa) complex.

Activated coagulation factor V (FVa) functions as a cofactor to factor Xa (FXa) in the conversion of prothrombin (PT) to thrombin. This essential procoagulant reaction, despite being the subject of extensive investigation, is not fully understood structurally and functionally. To elucidate the structure of the FXa-FVa complex, we have performed protein:protein (Pr:Pr) docking simulation with the pseudo-Brownian Pr:Pr docking ICM package and with the shape-complementarity Pr:Pr docking program PPD. The docking runs were carried out using a new model of full-length human FVa and the X-ray structure of human FXa. Five representative models of the FXa-FVa complex were in overall agreement with some of the available experimental data, but only one model was found to be consistent with almost all of the reported experimental results. The use of hybrid docking approach (theoretical plus experimental) is definitively important to study such large macromolecular complexes. The FXa-FVa model we have created will be instrumental for further investigation of this macromolecular system and will guide future site directed mutagenesis experiments.

Factor Va↗

Treatment of electrostatic effects in proteins: multigrid-based Newton iterative method for solution of the full nonlinear Poisson-Boltzmann equation.

The nonlinear Poisson-Boltzmann equation (NPBE) provides a continuum description of the electrostatic field in an ionic medium around a macromolecule. Here, a novel approach to the solution of the full NPBE is developed. This robust and efficient algorithm combines multilevel techniques with a damped inexact Newton's method. The CPU time required for solution of the full NPBE, which is less than that for standard single-grid approaches in solving the corresponding linearized equation, is proportional to the number of unknowns enabling applications to very large macromolecular systems. Convergence of the method is demonstrated for a variety of protein systems. Comparison of the solutions to the linearized Poisson-Boltzmann equation shows that the damping of the electrostatic field around the charge is increased and that the potential scales logarithmically with charge. The inclusion of the full nonlinearity thus reduces the impact of highly charged residues on protein surfaces and provides a more realistic representation of electrostatic effects. This is demonstrated through calculation of potential around the active site regions of the 1,266-residue tryptophan synthase dimer and in the computation of rate constants from Brownian dynamics calculations in the superoxide dismutase-superoxide and antibody-antigen systems.

Algorithms↗

Quantitative analysis of crystal growth. Tryptophanyl-tRNA synthetase crystal polymorphism and its relationship to catalysis.

We show that quantitative analysis of replicated, full-factorial crystal growth experiments and, by implication, similar studies of a wide variety of other phenomena, can be a powerful tool for analyzing macromolecular systems with complex, interacting dependencies on functionally significant factors. Bacillus stearothermophilus tryptophanyl-tRNA synthetase crystallizes in three different crystal forms depending on the ligands present under otherwise identical conditions. Comparison of crystallographic space groups for complexes with different ligands reveals that the three forms entail at least two very different families of packing arrangements that are correlated with specific changes in the enzyme ligation state. One is associated with the ligand-free enzyme, substrate ligands, and the binding of the activated amino acid; the other results from the presence of high ATP concentrations and/or the synthesis of the unusual acyl-transfer product, tryptophanyl-2'(3') ATP. Together with previous physico-chemical studies of aminoacyl-tRNA synthetases, these observations suggest that the two families are related, respectively, to the biochemical processes of amino acid activation and acyl transfer. Further evidence that the crystal polymorphism results from an underlying protein conformational polymorphism has now been obtained by quantitative analysis of how crystal growth depends on pH and the substrates tryptophan and ATP. The analysis consists first in showing that crystallization conditions for the unliganded protein are very favorable, suggesting that variation in crystal growth induced by pH and substrates under otherwise identical conditions is due to their effects on the protein conformation and not on incidental perturbations of crystal growth, per se. Next, crystal growth experiments are shown to be reproducible enough to support statistical analysis of quantitative scores assigned to the results. Finally, the observed variation in scores can be attributed at high confidence levels chiefly to three effects: that of pH alone, the synergistic effects of pH plus tryptophan, and of tryptophan plus ATP. These statistical inferences are consistent with other biochemical data, and support the conclusions based on crystal packing that representative stages of the enzyme mechanism have been trapped in the different crystal forms. The pH-tryptophan interaction implies that there is a pH-dependent conformational change favoring high affinity substrate binding at high pH. The pH-ATP interaction implies that a subsequent conformational change, not previously considered, occurs between tryptophan activation and acyl transfer.

Adenosine Triphosphate↗

Adsorption of biological molecules to a solid support for scanning probe microscopy.

Scanning probe microscopes are now established tools to study the surface structure of biological macromolecules under physiological conditions. Sample preparation methods for this microscopy all have the objective to attach the specimen firmly to a support. Here we analyse the commonly used method of adsorbing biological specimens to freshly cleaved mica. This is facilitated by adjusting the electrolyte concentration and the pH of the buffer solution. Native macromolecular systems absorbed to mica in this way can be reproducibly imaged at submolecular resolution.

Adsorption↗

HeteroTOCSY-based experiments for measuring heteronuclear relaxation in nucleic acids and proteins.

While both 31P and 113Cd are present at locations of interest in many different macromolecular systems, heteronuclear-detected relaxation measurements on these nuclei have been restrained by limitations in either resolution or signal-to-noise ratio. We have developed heteroTOCSY-based methods to overcome both of these problems. Two-dimensional versions of these experiments were utilized to measure 31P T1 and T2 values in DNA oligonucleotides; the additional resolution offered by a second dimension allowed determination of these values for most of the 31P resonances in a DNA dodecamer. The results from the experiments indicated that there was little significant variation in T1 values for the different phosphates in the DNA dodecamer; however, the T2 values showed a clear pattern, with lower values in the interior of the sequence than at the ends of the helix. Furthermore, a significant correlation between 31P chemical shifts and T2 values was observed. One-dimensional, frequency-selective versions of these experiments were also developed for use on systems containing a smaller number of heteronuclear spins. These methods were applied to investigate the heteronuclear relaxation properties of 113Cd in 113Cd2LAC9(61), a Cys6Zn2 DNA-binding domain. Data from the experiments confirm biochemical evidence that more significant differences occur in the metal-protein interactions between the two metal-binding sites than has been previously identified for proteins containing this motif.

Antimetabolites, Antineoplastic↗

Enhancement of etoposide-induced cytotoxicity by cyclosporin A.

Following the clinical observation of enhanced antineoplastic action of etoposide in the presence of cyclosporin A (CyA), we investigated this drug interaction in several in vitro and in vivo tumor systems. Macromolecular DNA damage induced by etoposide at drug levels comparable to plasma AUC values achieved in patients was increased not only in leukemic peripheral blood cells from patients but also in mononuclear peripheral blood cells from a healthy donor. Intracellular retention of radioactivity from 3H-etoposide was increased by a factor of 1.5 at the most in the presence of CyA. The cytotoxicity of etoposide and adriamycin to L 1210 leukemic cells was clearly enhanced, whereas CyA had no effect on the action of cisplatin or ionizing irradiation. At CyA blood levels not exceeding 1.44 microgram/ml, increased tumor inhibition of etoposide was observed in a human embryonal cancer xenograft, but there was also higher lethality in normal mice. We conclude from our own data and from other recent findings that with respect to chemosensitization the effects of CyA resemble those of calcium channel blockers or anticalmodulin agents. In contrast to calcium channel blockers, however, adequate plasma levels of CyA can well be achieved in patients.

Animals↗

On the mechanism of a radiation-induced change in enzymic differentiation during development.

In experiments on glucose-6-phosphatase and tyrosine aminotransferase it was shown that radiation induces changes in enzymic differentiation in perinatal rat liver. A study was made of the probable reasons for the observed changes. It was shown that the macromolecular system of the protein enzyme synthesis was not damaged by the radiation doses used. The observed decrease in glucose-6-phosphatase activity during late embryogenesis, after pre-irradiation at early organogenesis, is eliminated by administration of exogenous thyroxine. A radiation-induced rise in the tyrosine aminotransferase activity during the perinatal period correlated with the cyclic AMP system status. It is proposed that modification of enzymic differentiation after irradiation results from the change in the amount of inductors.

Animals↗

Molecular recognition in applied enzyme chemistry.

Molecular recognition impinges upon many fields of biological chemistry, especially those involving catalytic processes. This review gives examples from studies at Strathclyde of both small and macromolecular systems. Mechanism-based enzyme inhibitors are described with reference to dihydrofolate reductase, dihydroorotate dehydrogenase, and cholesterol metabolism. Applications of molecular recognition related to synthetic transformations are discussed in terms of aromatic substitution, chemically modified papain, and catalytic antibodies for Diels-Alder reactions.

Animals↗

Correlations between age and rat dermis modifications. Ultrastructural-morphometric evaluations and lysyl oxidase activity.

The extracellular matrix is a complex, integrated macromolecular system which plays a crucial role in the economy of each organ. In this study we focused our attention on the correlations between age and rat skin dermis. The latter was chosen as a model of the connective tissue, and was analyzed by means of electron microscopy and by measurement of the activity of lysyl oxidase, the enzyme involved in collagen and elastin crosslink formation. Ultrastructural and morphometric evaluations associated to body weight growth, showed a progressive increase in the amounts of extracellular components and a progressive reduction in the cell density. Skin from adult animals appeared characterized by a well organized matrix; by contrast, in old rats, we observed several degenerative features such as the disorganization of collagen bundles, the vacuolization of elastic fibers, and the atrophy of the mesenchimal cells. Morphometric evaluations in old animals showed a slight but significant reduction in the percentage of the total collagen measured, a fair stability in the area occupied by the elastin fibers, and an increase of the apparently non-structured matrix. The fact that lysyl oxidase activity was diminished in old rats does not corroborate the observation by several authors that increased collagen insolubility is a consequence of higher intra- and intermolecular crosslinking. This would suggest that other chemical modifications, such as crosslink oxidation or non enzymatic glycosylation, might be involved during the aging of connective tissue. The qualitative and quantitative modifications observed at all ages illustrate the correlation between connective tissue modifications and structural and/or functional properties of the skin.

Age Factors↗

Pressure effects on water-swollen elastin. A model for hydrophobic interactions in proteins.

The effect of pressure on the swelling of elastin in pure water was investigated. Because elastin is a very non-polar protein, and because the swelling of elastin can be directly related to changes in the strength of hydrophobic interactions, we have used elastin as a model to study the effect of pressure on hydrophobic interactions in proteins. The elastin swelling model is particularly useful because it is based on a macromolecular system very similar to a globular protein and not on dilute aqueous solutions of small non-polar compounds. Increased pressure causes elastin to increase its swollen volume, and the observed swelling changes were analyzed in terms of the Flory-Rehner theory (Flory, P.J. and Rehner, Jr., J (1943) J. Chem. Phys. 11, 521--526) for the swelling of kinetically free, random polymer networks. Our calculations provide a measure of the volume change for the process of transferring 1 mol of an average non-polar amino acid side chain from a region where the side chains are surrounded by other non-polar groups and have no contact with water (i.e. a hydrophobic region) into contact with water. The results indicate that there is a small, negative volume change associated with this process, and quantitative estimates indicate that the volume change is of the order of --6 ml/mol side chain. The results support the hypothesis that the free energy required to transfer a non-polar side chain from a hydrophobic region into water becomes less positive (i.e. the hydrophobic interaction becomes weaker) as hydrostatic pressure is increased.

Elastin↗

Protein docking and complementarity.

Predicting the structures of protein-protein complexes is a difficult problem owing to the topographical and thermodynamic complexity of these structures. Past efforts in this area have focussed on fitting the interacting proteins together using rigid body searches, usually with the conformations of the proteins as they occur in crystal structure complexes. Here we present work which uses a rigid body docking method to generate the structures of three known protein complexes, using both the bound and unbound conformations of the interacting molecules. In all cases we can regenerate the geometry of the crystal complexes to high accuracy. We also are able to find geometries that do not resemble the crystal structure but nevertheless are surprisingly reasonable both mechanistically and by some simple physical criteria. In contrast to previous work in this area, we find that simple methods for evaluating the complementarity at the protein-protein interface cannot distinguish between the configurations that resemble the crystal structure complex and those that do not. Methods that could not distinguish between such similar and dissimilar configurations include surface area burial, solvation free energy, packing and mechanism-based filtering. Evaluations of the total interaction energy and the electrostatic interaction energy of the complexes were somewhat better. Of the techniques that we tried, energy minimization distinguished most clearly between the "true" and "false" positives, though even here the energy differences were surprisingly small. We found the lowest total interaction energy from amongst all of the putative complexes generated by docking was always within 5 A root-mean-square of the crystallographic structure. There were, however, several putative complexes that were very dissimilar to the crystallographic structure but had energies that were close to that of the low energy structure. The magnitude of the error in energy calculations has not been established in macromolecular systems, and thus the reliability of the small differences in energy remains to be determined. The ability of this docking method to regenerate the crystallographic configurations of the interacting proteins using their unbound conformations suggests that it will be a useful tool in predicting the structures of unsolved complexes.

Algorithms↗

Perturbation of molecular species distribution in steady states supported by a flow of energy. Models analogous to Ca2(+)-dependent ATPase and phosphorylase b.

When an allosteric macromolecular system is capable of existing in two conformations, both of which may be converted into energy-storing forms by the binding of a substrate or by the absorption of radiant energy, then a kinetic process may occur, such as an enzymic conversion of the substrate into products, which liberates energy and selectively depletes one or more of the forms of the macromolecule. Upon a continuous supply of energy, a steady state, or pseudoequilibrium, is reached during which the selective depletion of molecular species results effectively in a directional flow of energy through the system. This perturbs the distribution of the various molecular species. This effect may simulate both positive and negative binding cooperativity, and mimic the presence of multiple binding sites with different affinities even in monomeric, monovalent systems. Specific model systems are presented analogous to the transport of Ca2+ by sarcoplasmic reticulum and the allosteric behavior of phosphorylase b.

Allosteric Regulation↗

Microcalorimetry of interaction of dihydro-imidazo-phenanthridinium (DIP)-based compounds with duplex DNA.

Isothermal titration (ITC) and differential scanning calorimetry (DSC) have been used to screen the binding thermodynamics of a family of DNA intercalators based on the dihydro-imidazo-phenanthridinium (DIP) framework. All members of this DIP-based ligand family bind to both genomic (calf thymus and/or salmon testes) and a synthetic dodecamer d(CGCGAATTCGCG) duplex DNA with broadly similar affinities regardless of side chain size or functionality. Viscosity measurements confirm that binding satisfies standard criteria for intercalation. Binding is exothermic but with an additional favourable positive entropy contribution in most cases at 25 degrees C, although a significant negative heat capacity effect (DeltaC(p)) means that both DeltaH(0) and DeltaS(0) decrease with increasing temperature. DIP-ligand binding to DNA also shows significant entropy-enthalpy compensation effects that are now almost standard in such situations, probably reflecting the conformational flexibility of macromolecular systems involving a multiplicity of weak non-covalent interactions. This ability to vary side chain functionality without compromising DNA binding suggests that the DIP framework should be a promising basis for more adventurous chemistry at the DNA level.

Calorimetry↗

Carbon-13 labeling for improved tracer depth profiling of organic materials using secondary ion mass spectrometry.

13C labeling is introduced as an alternative to deuterium labeling for analysis of organic materials using secondary ion mass spectrometry (SIMS). A model macromolecular system composed of polystyrene (PS) and poly(methyl methacrylate) (PMMA) was used to compare the effects of isotopic labeling using both deuterium substitution (dPS) and 13C labeling (13C-PS). Clear evidence is shown that deuterium labeling does introduce changes in the thermodynamic properties of the system, with the observation of segregation of dPS to an hPS:dPS/hPMMA interface. This type of behavior could significantly impact many types of investigations due to the potential for improper interpretation of experimental results as a consequence of labeling-induced artifacts. 13C labeling is shown to provide a true tracer for analysis using SIMS.

Carbon Isotopes↗

Biosorption of cadmium by the protonated macroalga Sargassum muticum: binding analysis with a nonideal, competitive, and thermodynamically consistent adsorption (NICCA) model.

Protonated biomass of the seaweed Sargassum muticum was investigated for its ability to remove cadmium(II) from aqueous solutions. In this work, a nonideal, semiempirical, thermodynamically consistent (NICCA) isotherm was proposed to fit the experimental ion binding data obtained in NaNO3 0.05 mol L(-1). This model describes the competition between protons and metal ions satisfactorily. Moreover, it reflects the complexity of the macromolecular systems that take part in biosorption considering the heterogeneity of the sorbent. It was demonstrated in this work that the NICCA isotherm constitutes a great improvement with respect to a simpler Langmuir competitive equation, which was not able to describe all the experimental data satisfactorily. Potentiometric acid-base titrations in the absence of cadmium were made to estimate the maximum amount of acid functional groups (2.61 mmol g(-1)) and the conditional proton binding parameters, logK (H) (3.8) and m(H) (0.54). The values of the binding parameters for the cadmium ion were chosen to provide the best simultaneous description of the isotherm at pH 4.5, as well as the dependence of cadmium adsorption on pH. Values of logK (Cd) (3.1), n(Cd) (1.8), and p (0.19) in the case of the NICCA isotherm or logK(Cd) (2.94-3.4) for Langmuir competitive models were obtained. Kinetic experiments were performed at two different pH values (3.0 and 4.5), establishing the time dependence that represents the sorption of cadmium with a pseudo-second-order kinetic model. It was observed that 4 h is enough to ensure that the equilibrium uptake was reached.

Adsorption↗