9 resultados para crystal structure and molecular dynamics
em Brock University, Canada
Resumo:
The crystal structure of Cu(PM)2(N03hoH20 (where PM is pyridoxamine, CSHI2N202) has been determined from three dimensional x-ray diffraction data. The crystals are triclinic, space group pI, a = 14.248 (2), b = 8.568 (1), c = 9.319 (1) 1, a = 94.08 (1), e = 89.73 (1), y~~ 99.18 (1)°, z = 2, jl(MoK) = 10.90 em-I, Po = 1.61 g/cm3 and Pc = 1.61 g/em3• The structure a was solved by Patterson techniques from data collected on a Picker 4-circle diffractometer to 26max = 45°. All atoms, including hydrogens, have been located. Anisotropic thermal parameters have been refined for all nonhydrogen atoms. For the 2390 independent reflections with F ? 3cr(F) , R = 0.0408. The results presented here provide the first detailed structural information of a metal complex with PM itself. The copper atoms are located on centres of symmetry and each is chela ted by two PM zwitterions through the amino groups and phenolate oxygen atoms. The zwitterionic form found in this structure involves the loss of a proton from the phenolate group and protonation of the pyridine ring nitrogen atoms. The two independent Cu(PM)2 moieties are symmetrically bridged by a single oxygen atom from one of the nitrate groups. The second nitrate group is not coordinated to the copper atoms but is central to an extensive hydrogen bonding network involving the water molecule and uncoordinated functional groups of PM.
Resumo:
The x-ray crystal structure of thiamine hydroiodide,C1ZH18N40S12' has been determined. The unit cell parameters are a = 13.84 ± 0.03, o b = 7.44 ± 0.01, c = 20.24 ± 0.02 A, 8 = 120.52 ± 0.07°, space group P2/c, z = 4. A total of 1445 reflections having ,2 > 2o(F2), 26 < 40° were collected on a Picker four-circle diffractometer with MoKa radiation by the 26 scan technique. The structure was solved by the heavy atom method. The iodine and sulphur atoms were refined anisotropically; only the positional parameters were refined for the hydrogen atoms. Successive least squares cycles yielded an unweighted R factor of 0.054. The site of protonation of the pyrimidine ring is the nitrogen opposite the amino group. The overall structure conforms very closely to the structures of other related thiamine compounds. The bonding surrounding the iodine atoms is distorted tetrahedral. The iodine atoms make several contacts with surrounding atoms most of them at or near the van der Waal's distances A thiaminium tetrachlorocobaltate salt was produced whose molecular and crystal structure was j~dged to be isomorphous to thiaminium tetrachlorocadmate.
Resumo:
Hg(18-Crown-6)C12 and Cd(18-Crown-6)C12 are isostructura1, space group Cl~ Z = 2. For the mercury compound, a = 10.444(2) A° , b = 11. 468(1) A° , c = 7.754(1) A° , a = 90.06(1)°, B = 82.20(1)°, Y = 90.07(1)°, Dobs = 1.87, Dca1c = 1.93, V = 920.05 13, R = 4.66%. For the cadmium compound, 000 a = 10.374(1) A, b = 11.419(2) A, c = 7.729(1) A, a = 89.95(1)°, B = 81.86(2)°, Y = 89.99(1)°, Dobs = 1.61, Dcalc = 1.64, V = 906.4613, R = 3.95%. The mercury and cadmium ions exhibit hexagonal bipyramidal coordination, with the metal ion located on a centre of symmetry in the plane of the oxygen atoms. The main differences between the two structures are an increase in the metal-oxygen distance and a reduction in the metalchloride distance when the central ion changes from Cd2+ to Hg2+. These differences may be explained in terms of the differences in hardness or softness of the metal ions and the donor atoms.
Resumo:
Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.
Resumo:
Polyglutamine is a naturally occurring peptide found within several proteins in neuronal cells of the brain, and its aggregation has been implicated in several neurodegenerative diseases, including Huntington's disease. The resulting aggregates have been demonstrated to possess ~-sheet structure, and aggregation has been shown to start with a single misfolded peptide. The current project sought to computationally examine the structural tendencies of three mutant poly glutamine peptides that were studied experimentally, and found to aggregate with varying efficiencies. Low-energy structures were generated for each peptide by simulated annealing, and were analyzed quantitatively by various geometry- and energy-based methods. According to the results, the experimentally-observed inhibition of aggregation appears to be due to localized conformational restraint placed on the peptide backbone by inserted prolines, which in tum confines the peptide to native coil structure, discouraging transition towards the ~sheet structure required for aggregation. Such knowledge could prove quite useful to the design of future treatments for Huntington's and other related diseases.
Resumo:
We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.
Resumo:
The nucleotide sequence of a genomic DNA fragment thought previously to contain the dihydrofolate reductase gene (DFR1) of Saccharomyces cerevisiae by genetic criteria was determined. This DNA fragment of 1784' basepairs contains a large open reading frame from position 800 to 1432, which encodes a enzyme with a predicted molecular weight of 24,229.8 Daltons. Analysis of the amino acid sequence of this protein revealed that the yeast polypep·tide contained 211 amino acids, compared to the 186 residues commonly found in the polypeptides of other eukaryotes. The difference in size of the gene product can be attributed mainly to an insert in the yeast gene. Within this region, several consensus sequences required for processing of yeast nuclear and class II mitochondrial introns were identified, but appear not sufficient for the RNA splicing. The primary structure of the yeast DHFR protein has considerable sequence homology with analogous polypeptides from other organisms, especially in the consensus residues involved in cofactor and/or inhibitor binding. Analysis of the nucleotide sequence also revealed the presence of a number of canonical sequences identified in yeast as having some function in the regulation of gene expression. These include UAS elements (TGACTC) required for tIle amino acid general control response, and "TATA H boxes as well as several consensus sequences thought to be required for transcriptional termination and polyadenylation. Analysis of the codon usage of the yeast DFRl coding region revealed a codon bias index of 0.0083. this valve very close to zero suggestes 3 that the gene is expressed at a relatively low level under normal physiological conditions. The information concerning the organization of the DFRl were used to construct a variety of fusions of its 5' regulatory region with the coding region of the lacZ gene of E. coli. Some of such fused genes encoded a fusion product that expressed in E.coli and/or in yeast under the control of the 5' regulatory elements of the DFR1. Further studies with these fusion constructions revealed that the beta-galactosidase activity encoded on multicopy plasmids was stimulated transiently by prior exposure of yeast host cells to UV light. This suggests that the yeast PFRl gene is indu.ced by UV light and nlay in1ply a novel function of DHFR protein in the cellular responses to DNA damage. Another novel f~ature of yeast DHFR was revealed during preliminary studies of a diploid strain containing a heterozygous DFRl null allele. The strain was constructed by insertion of a URA3 gene within the coding region of DFR1. Sporulation of this diploid revealed that meiotic products segregated 2:0 for uracil prototrophy when spore clones were germinated on medium supplemented with 5-formyltetrahydrofolate (folinic acid). This finding suggests that, in addition to its catalytic activity, the DFRl gene product nlay play some role in the anabolisln of folinic acid. Alternatively, this result may indicate that Ura+ haploid segregants were inviable and suggest that the enzyme has an essential cellular function in this species.
Resumo:
The interaction of biological molecules with water is an important determinant of structural properties both in molecular assemblies, and in conformation of individual macromolecules. By observing the effects of manipulating the activity of water (which can be accomplished by limiting its concentration or by adding additional solutes, "osmotic stress"), one can learn something about intrinsic physical properties of biological molecules as well as measure an energetic contribution of closely associated water molecules to overall equilibria in biological reactions. Here two such studies are reported. The first of these examines several species of lysolipid which, while present in relatively low concentrations in biomembranes, have been shown to affect many cellular processes involving membrane-protein or membrane-membrane interactions. Monolayer elastic constants were determined by combining X-ray diffraction and the osmotic stress technique. Spontaneous radii of curvature of lysophosphatidylcholines were determined to be positive and in the range +30A to +70A, while lysophosphatidylethanolamines proved to be essentially flat. Neither lysolipid significantly affected the bending modulus of the monolayer in which it was incorporated. The second study examines the role of water in theprocess of polymerization of actin into filaments. Water activity was manipulated by adding osmolytes and the effect on the equilibrium dissociation constant (measured as the criticalmonomer concentration) was determined. As water activity was decreased, the critical concentration was reduced for Ca-actin but not for Mg-actin, suggesting that 10-12 fewer water molecules are associated with Ca-actin in the polymerized state. Thisunexpectedly small amount of water is discussed in the context of the common structural motif of a nucleotide binding cleft.
Resumo:
The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.