NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation ill each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA ar rd RNA duplex structures.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

MIPS: metal interactions in protein structures

MIPS (metal interactions in protein structures) is a database of metals in the three-dimensional acromolecular structures available in the Protein Data Bank. Bound metal ions in proteins have both catalytic and structural functions. The proposed database serves as an open resource for the analysis and visualization of all metals and their interactions with macromolecular (protein and nucleic acid) structures. MIPS can be searched via a user-friendly interface, and the interactions between metals and protein molecules, and the geometric parameters, can be viewed in both textual and graphical format using the freely available graphics plug-in Jmol. MIPS is updated regularly, by means of programmed scripts to find metal-containing proteins from newly released protein structures. The database is useful for studying the properties of coordination between metals and protein molecules. It also helps to improve understanding of the relationship between macromolecular structure and function. This database is intended to serve the scientific community working in the areas of chemical and structural biology, and is freely available to all users, around the clock, at http://dicsoft2.physics.iisc.ernet.in/mips/.

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I-turn forming tetrapeptide and an -helical tetradecapeptide

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt -helical conformations stabilized by 11 successive 5 1 hydrogen bonds. In addition, a single 4 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (, ) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle NCC() and the observed backbone , values. For > 106° , helices are observed, while fully extended structures are characterized by < 106° . The mean values for extended and folded conformations for the Dpg residue are 103.6° ± 1.7° and 109.9° ± 2.6° , respectively.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Parallel and antiparallel aggregation of alpha-helices. Crystal structures of two apolar decapeptides X-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe (X = Boc, Ac)

An apolar helical decapeptide with different end groups, Boc- or Ac-, crystallizes in a completely parallel fashion for the Boc-analog and in an antiparallel fashion for the Ac-analog. In both crystals, the packing motif consists of rows of parallel molecules. In the Boc-crystals, adjacent rows assemble with the helix axes pointed in the same direction. In the Ac-crystals, adjacent rows assemble with the helix axes pointed in opposite directions. The conformations of the molecules in both crystals are quite similar, predominantly alpha-helical, except for the tryptophanyl side chain where chi 1 congruent to 60 degrees in the Boc- analog and congruent to 180 degrees in the Ac-analog. As a result, there is one lateral hydrogen bond between helices, N(1 epsilon)...O(7), in the Ac-analog. The structures do not provide a ready rationalization of packing preference in terms of side-chain interactions and do not support a major role for helix dipole interactions in determining helix orientation in crystals. The crystal parameters are as follow. Boc-analog: C60H97N11O13.C3H7OH, space group Pl with a = 10.250(3) A, b = 12.451(4) A, c = 15.077(6) A, alpha = 96.55(3) degrees, beta = 92.31(3) degrees, gamma = 106.37(3) degrees, Z = 1, R = 5.5% for 5581 data ([F] greater than 3.0 sigma(F)), resolution 0.89 A. Ac-analog: C57H91N11O12, space group P2(1) with a = 9.965(1) A, b = 19.707(3) A, c = 16.648(3) A, beta = 94.08(1), Z = 2, R = 7.2% for 2530 data ([F] greater than 3.0 sigma(F)), resolution 1.00 A.

Modular design of synthetic protein mimics. Crystal structures, assembly, and hydration of two 15- and 16-residue apolar, leucyl-rich helical peptides

Two IS- and 16-residue peptides containing a-aminoisobutyric acid (Aib) have been synthesized, as part of a strategy to construct stereochemically rigid peptide helices, in a modular approach to design of protein mimics. The peptides Boc-(Val-Ala-Leu-Aib),-OMe ( I ) and Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-(Val-Ala-Leu-Aib()11z)- OhaMvee been crystallized.Both crystals are stable only in the presence of mother liquor or water. The crystal data are as follows. I: C78H140N16019~2H20,P2,, a = 16.391 (3) A, b = 16.860 (3) A, c = 18.428 (3) A, p = 103.02 (I)O, Z = 2, R = 9.6% for 3445 data with lFol >30(F), resolution 0.93 A. 11: C7,Hl,,N,S018.7.5H,0, C2221, a = 18.348 ( 5 ) A, b = 47.382 (1 1) A, c = 24.157 ( 5 ) A, Z =8, R = l0,6%, for 3147 data with lFol > 3a(F), resolution 1.00 A. The 15-residue peptide (11) is entirely a helical, while the 16-residue peptide ( I ) has a short segment of 310 helix at the N terminus. The packing of the helices in the crystals is rather incfficicnt with no particular attractions between Leu-Leu side chains, or any other pair. Both crystals have fairly large voids, which are filled with water molecules in a disordered fashion. Water molecule sites near the polar head-to-tail regions are well detcrmined, those closer to the hydrophobic side chains less so and a number of possible water sites in the remaining "empty" space are not determined. No interdigitation of Leu side chains is observed in the crystal as is hypothesized in the "leucine zipper" class of DNA binding proteins.

Crystal structures of a nonapeptide helix containing alpha,alpha-di-n-butylglycine (Dbg), Boc-Gly-Dbg-Ala-Val-Ala-Leu-Aib-Val-Leu-OMe

Two crystals structures of a nonapeptide (anhydrous and hydrated) containing the amino acid residue alpha, alpha-di-n-butylglycyl, reveal a mixed 3(10)/alpha-helical conformation. Residues 1-7 adopt phi, psi values in the helical region, with Val(8) being appreciably distorted. The Dbg residue has phi, psi values of -40, -37 degrees and -46, -40 degrees in two crystals with the two butyl side chains mostly extended in each. Peptide molecules in the crystals pack into helical columns. The crystal parameters are C50H91N9O12, space group P2(1), with a = 9.789(1) Angstrom, b = 20.240(2) Angstrom, c = 15.998(3) Angstrom, beta = 103.27(1); Z = 2, R = 10.3% for 1945 data observed >3 sigma(F) and C50H91N9O12. 3H(2)O, space group P2(1), with a = 9.747(3) Angstrom, b = 21.002(8) Angstrom, c = 15.885(6) Angstrom, beta = 102.22(3)degrees, Z = 2, R = 13.6% for 2535 data observed >3 sigma(F). The observation of a helical conformation at Dbg suggests that the higher homologs in the alpha, alpha-dialkylated glycine series also have a tendency to stabilize peptide helices. (C) Munksgaard 1996.

Disallowed Ramachandran Conformations of Amino Acid Residues in Protein Structures

An analysis of the nature and distribution of disallowed Ramachandran conformations of amino acid residues observed in high resolution protein crystal structures has been carried out. A data set consisting of 110 high resolution, non-homologous, protein crystal structures from the Brookhaven Protein Data Bank was examined. The data set consisted of a total of 18,708 non-Gly residues, which were characterized on the basis of their backbone dihedral angles (φ, ψ). Residues falling outside the defined “broad allowed limits” on the Ramachandran map were chosen and the reportedB-factor value of the α-carbon atom was used to further select well defined disallowed conformations. The conformations of the selected 66 disallowed residues clustered in distinct regions of the Ramachandran map indicating that specific φ, ψ angle distortions are preferred under compulsions imposed by local constraints. The distribution of various amino acid residues in the disallowed residue data set showed a predominance of small polar/charged residues, with bulky hydrophobic residues being infrequent. As a further check, for all the 66 cases non-hydrogen van der Waals short contacts in the protein structures were evaluated and compared with the ideal “Ala-dipeptide” constructed using disallowed dihedral angle (φ, ψ) values. The analysis reveals that short contacts are eliminated in most cases by local distortions of bond angles. An analysis of the conformation of the identified disallowed residues in related protein structures reveals instances of conservation of unusual stereochemistry.

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I 'beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I' beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 -> 1 hydrogen bonds. In addition, a single 4 -> 1 hydrogen bond is also observed at the N-terminus. All live Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C-alpha-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean r values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.

A numerical survey of relativistic rotating neutron star structures using the Hartle-Thorne formalism

A broad numerical survey of relativistic rotating neutron star structures was compiled using an exhaustive list of presently available equation of state models for neutron star matter. The structure parameters (spherical deformations in mass and radii, the moment of inertia and quadrupole moment, oblateness, and free precession) are calculated using the formalism proposed by Hartle and Thorne (1968). The results are discussed in relation to the relevant observational information. Binary pulsar data and X-ray burst sources provide information on the bulk properties of neutron stars, enabling the derivation of constraints that can be put on the structure of neutron stars and equation of state models.

Use of 2-pyrimidineamidoxime to generate polynuclear homo-/heterometallic assemblies: synthesis, crystal structures and magnetic study with theoretical investigations on the exchange mechanism

Three new transition metal complexes using 2-pyrimidineamidoxime (pmadH(2)) as multidentate chelating and/or bridging ligand have been synthesized and characterized. The ligand pmadH(2) has two potential bridging functional groups mu-O and mu-(N-O)] and consequently shows several coordination modes. While a polymeric 1D Cu-II complex Cu(pmadH(2))(2)(NO3)](NO3) (1) was obtained upon treatment of Cu(NO3)(2)center dot 3H(2)O with pmadH(2) at room temperature in the absence of base, a high temperature reaction in the presence of base yielded a tetranuclear Cu-II-complex Cu-4(pmad)(2)(pmadH)(2)(NO3)](NO3)(H2O) (2). One of the Cu-II centers is in a square pyramidal environment while the other three are in a square planar geometry. Reaction of the same ligand with an equimolar mixture of both Cu(NO3)(2)center dot 3H(2)O and NiCl2 center dot 6H(2)O yielded a tetranuclear heterometallic (Cu2Ni2II)-Ni-II complex Cu2Ni2(pmad)(2)(pmadH)(2)Cl-2]center dot H2O (3) containing both square planar (Ni-II) and square pyramidal (Cu-II) metal centers. Complexes 1-3 represent the first examples of polynuclear metal complexes of 2-pyrimidineamidoxime. The analysis of variable temperature magnetic susceptibility data of 2 reveals that both ferromagnetic and antiferromagnetic interactions exist in this complex (J(1) = +10.7 cm(-1) and J(2) = -2.7 cm(-1) with g = 2.1) leading to a resultant ferromagnetic behavior. Complex 3 shows expected antiferromagnetic interaction between two Cu-II centers through -N-O- bridging pathway with J(1) = -3.4 cm(-1) and g = 2.08. DFT calculations have been used to corroborate the magnetic results.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Reactivity Of Pph3 Toward Ruiiruiiicl(O2car)4 - Syntheses, Molecular-Structures, And Spectroscopic And Electrochemical Properties Of Ruiiruiii(Oh2)Cl(Mecn)(O2car)4(Pph3)2 And Ruii2(Oh2)(Mecn)2(O2car)4(Pph3)2

By the reaction of Ru2Cl(O2CAr)4 (1) and PPh3 in MeCN-H2O the diruthenium(II,III) and diruthenium(II) compounds of the type Ru2(OH2)Cl(MeCN)(O2CAr)4(PPh3)2 (2) and Ru2(OH2)(MeCN)2(O2CAr)4(PPh3)2 (3) were prepared and characterized by analytical, spectral, and electrochemical data (Ar is an aryl group, C6H4-p-X; X = H, OMe, Me, Cl, NO2). The molecular structure of Ru2(OH2)Cl(MeCN)(O2CC6H4-p-OMe)4(PPh3)2 was determined by X-ray crystallography. Crystal data are as follows: triclinic, P1BAR, a = 13.538 (5) angstrom, b = 15.650 (4) angstrom, c = 18.287 (7) angstrom, alpha = 101.39 (3)-degrees, beta = 105.99 (4)-degrees, gamma = 97.94 (3)-degrees, V = 3574 angstrom 3, Z = 2. The molecule is asymmetric, and the two ruthenium centers are clearly distinguishable. The Ru(III)-Ru(II), Ru(III)-(mu-OH2), and Ru(II)-(mu-OH2) distances and the Ru-(mu-OH2)-Ru angle in [{Ru(III)Cl(eta-1-O2CC6H4-p-OMe)(PPh3)}(mu-OH2)(mu-O2CC6H4-p-OMe)2{Ru(II)(MeCN)(eta-1-O2CC6H4-p-OMe)(PPh3)}] are 3.604 (1), 2.127 (8), and 2.141 (10) angstrom and 115.2 (5)-degrees, respectively. The compounds are paramagnetic and exhibit axial EPR spectra in the polycrystalline form. An intervalence transfer (IT) transition is observed in the range 900-960 nm in chloroform in these class II type trapped mixed-valence species 2. Compound 2 displays metal-centered one-electron reduction and oxidation processes near -0.4 and +0.6 V (vs SCE), respectively in CH2Cl2-TBAP. Compound 2 is unstable in solution phase and disproportionates to (mu-aquo)diruthenium(II) and (mu-oxo)diruthenium(III) complexes. The mechanistic aspects of the core conversion are discussed. The molecular structure of a diruthenium(II) compound, Ru2(OH2)(MeCN)2(O2CC6H4-p-NO2)4(PPh3)2.1.5CH2Cl2, was obtained by X-ray crystallography. The compound crystallizes in the space group P2(1)/c with a = 23.472 (6) angstrom, b = 14.303 (3) angstrom, c = 23.256 (7) angstrom, beta = 101.69 (2)-degrees, V = 7645 angstrom 3, and Z = 4. The Ru(II)-Ru(II) and two Ru(II)-(mu-OH2) distances and the Ru(II)-(mu-OH2)-Ru(II) angle in [{(PPh3)-(MeCN)(eta-1-O2CC6H4-p-NO2)Ru}2(mu-OH2)(mu-O2CC6H4-p-NO2)2] are 3.712 (1), 2.173 (9), and 2.162 (9) angstrom and 117.8 (4)-degrees, respectively. In both diruthenium(II,III) and diruthenium(II) compounds, each metal center has three facial ligands of varying pi-acidity and the aquo bridges are strongly hydrogen bonded with the eta-1-carboxylato facial ligands. The diruthenium(II) compounds are diamagnetic and exhibit characteristic H-1 NMR spectra in CDCl3. These compounds display two metal-centered one-electron oxidations near +0.3 and +1.0 V (vs SCE) in CH2Cl2-TBAP. The overall reaction between 1 and PPh3 in MeCN-H2O through the intermediacy of 2 is of the disproportionation type. The significant role of facial as well as bridging ligands in stabilizing the core structures is observed from electrochemical studies.

Empirical torsional potential functions from protiesn structure DATA Phi- and Psi-Potentials for Non-glycyl Amino Acid Residues

The torsional potential functions Vt(φ) and Vt(ψ) around single bonds N–Cα and Cα-C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (φ, ψ)-plane with the value of Vtot(φ, ψ), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in l-configuration, are Vt(φ) = – 1.0 cos (φ + 60°); Vt(ψ) = – 0.5 cos (ψ + 60°) – 1.0 cos (2ψ + 30°) – 0.5 cos (3ψ + 30°). The dipeptide energy maps Vtot(φ, ψ) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line ψ = 0°. These functions, derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.