PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

Microwave spectroscopic and theoretical studies on the phenylacetylene center dot center dot center dot H(2)O complex: C-H center dot center dot center dot O and O-H center dot center dot center dot pi hydrogen bonds as equal partners

Rotational spectra of five isotopologues of the title complex, C(6)H(5)CCH center dot center dot center dot H(2)O, C(6)H(5)CCH center dot center dot center dot HOD, C(6)H(5)CCH center dot center dot center dot D(2)O, C(6)H(5)CCH center dot center dot center dot H(2)(18)O and C(6)H(5)CCD center dot center dot center dot H(2)O, were measured and analyzed. The parent isotopologue is an asymmetric top with kappa = -0.73. The complex is effectively planar (ab inertial plane) and both `a' and `b' dipole transitions have been observed but no c dipole transition could be seen. All the transitions of the parent complex are split into two resulting from an internal motion interchanging the two H atoms in H(2)O. This is confirmed by the absence of such doubling for the C(6)H(5)CCH center dot center dot center dot HOD complex and a significant reduction in the splitting for the D(2)O analog. The rotational spectra, unambiguously, reveal a structure in which H(2)O has both O-H center dot center dot center dot pi (pi cloud of acetylene moiety) and C-H center dot center dot center dot O (ortho C-H group of phenylacetylene) interactions. This is in agreement with the structure deduced by IR-UV double resonance studies (Singh et al., J. Phys. Chem. A, 2008, 112, 3360) and also with the global minimum predicted by advanced electronic structure theory calculations (Sedlack et al., J. Phys. Chem. A, 2009, 113, 6620). Atoms in Molecule (AIM) theoretical analysis of the complex reveals the presence of both O-H center dot center dot center dot pi and C-H center dot center dot center dot O hydrogen bonds. More interestingly, based on the electron densities at the bond critical points, this analysis suggests that both these interactions are equally strong. Moreover, the presence of both these interactions leads to significant deviation from linearity of both hydrogen bonds.

Theoretical Studies on the Structure and Bonding of Metallacyclocumulenes, -cyclopentynes, and -cycloallenes

In this paper, we compare the electronic structure of the hafnacycloallene complex Cp(2)HfC(4)Rr'(2)R `' (5Hf), which was previously described by Erker et al., with those of the titanium, zirconium, and hafnium complexes Cp(2)M(eta(4)-RHC(4)HR) (3M; i.e. metallacyclopent-2,3,4-trienes, metallacyclocumulenes) and Cp(2)M(eta(2)-R(2)C(4)R(2)) (4M; i.e. 1-metallacyclopent-3-ynes) using density functional theory (BP86/LANL2DZ) calculations. Moreover, the eta(3)-phenylallenyl zirconocene complex 7Zr, which was synthesized by Wojcicki et al., is included for the comparison. These calculations and extended Huckel calculations show that the bonding in complex 5Hf is remarkably similar to that of complexes 4M and 7Zr. An analysis of the structural parameters and bonding reveals that the unique interaction of the internal carbon atoms along with the terminal carbon atoms with the bent-metallocene moiety is the reason for the unusual stability of these metallacycles. The molecular orbital analysis further suggests that complex 5Hf can react with another metal fragment to give the bimetallic complexes 9 and 10. The electronic structures of complexes 3M, 4M, 5Hf, and 7Zr have been comparatively studied to get a general understanding of the bonding in these metallacycles.

Novel ABO(3) oxides related to perovskite and YAlO3 structure types in the La-B-V-O (B = Ni, Cu) systems

The synthesis, structure and magnetic properties of mixed-metal oxides of ABO(3) composition in the La-B-V-O (B = Ni, Cu) systems are described in the present paper. While the B = Ni oxides adopt GdFeO3-like perovskite structure containing disordered nickel and vanadium at the octahedral B site, La3Cu2VO9 crystallizes in a YAlO3-type structure. A detailed investigation of the superstructure of nominal La3Cu2VO9 by WDS analysis and Rietveld refinement of powder XRD data reveal that the likely composition of the phase is La13Cu9V4O38.5, where the Cu and V atoms are ordered in a root13a(h) (a(h) = hexagonal a parameter of YAlO3-like subcell) superstructure. Magnetic susceptibility data support the proposed superstructure consisting of triangular Cu-3 clusters. At low temperatures, the magnetic moment corresponds to S = 1/2 per Cu-3 cluster, while at high temperatures the behavior is Curie-Weiss like, showing S = 1/2 per copper. The present work reveals the contrasting behavior of La-Cu-V-O and La-Ni-V-O systems: while a unique line-phase related to YAlO3 structure is formed around La3Cu2VO9 Composition in the copper system, a continuous series of perovskite-GdFeO3 solid solutions, LaNi1-xVxO3 for 0 less than or equal to x less than or equal to 1/3 seems to be obtained in the nickel system, where the oxidation state of nickel varies from 3+ to 2+.

The central region of meso-3,4-diphenylhexane-2,5-dione

As with 1,2-diphenylethane (dpe), X-ray crystallographic methods measure the central bond in meso-3,4-diphenylhexane-2,5-done (dphd) as significantly shorter than normal for an sp(3)-sp(3) bond. The same methods measure the benzylic (ethane C-Ph) bonds in dphd as unusually long for sp(3)-sp(2) liaisons. Torsional motions of the phenyl rings about the C-Ph bonds have been proposed as the artifacts behind the result of a 'short' central bond in dpe. While a similar explanation can, presumably, hold for the even 'shorter' central bond in dphd, it cannot account for the 'long' C-Ph bonds. The phenyl groups, departing much from regular hexagonal shape, adopt highly skewed conformations with respect to the plane constituted by the four central atoms. It is thought that-the thermal motions of the phenyl rings, conditioned by the potential wells in which they are ensconced in the unit cell, are largely libratory around their normal axes. In what appears to be a straightforward explanation under the 'rigid-body' concept, it appears that these libratory motions of the phenyl rings, that account, at the same time, for the 'short' central bond, are the artifacts behind the 'long' measurement of the C-Ph bonds. These motions could be superimposed on torsional motions analogous to those proposed in the case of dpe. An inspection of the ORTEP diagram from the 298 K data on dphd clearly suggests these possibilities. Supportive evidence for these qualitative explanations from an analysis of the differences between the mean square displacements of C(1) and C(7)/C(1a) and C(7a) based on the 'rigid-body model' is discussed. (C) 2002 Elsevier Science B.V. All rights reserved.

Visualizing lone pairs in compounds containing heavier congeners of the carbon and nitrogen group elements

In this mini-review, I discuss some recent work on the stereochemistry and bonding of lone pairs of electrons in divalent compounds of the heavier carbon group elements (SnII, PbII) and in trivalent compounds of the heavier nitrogen group elements (BiIII). Recently developed methods that permit the real-space visualization of bonding patterns on the basis of density functional calculations of electronic structure, reveal details of the nature of s electron lone pairs in compounds of the heavier main group elements – their stereochemistry and their inertness (or lack thereof). An examination of tetragonal P4/nmm SnO, a-PbO and BiOF, and cubic Fm3m PbS provides a segue into perovskite phases of technological significance, including ferroelectric PbTiO3 and antiferroelectric/piezoelectric PbZrO3, in both of which the lone pairs on Pb atoms play a pivotal rôle.

Cyclometalated complexes derived from calix[4] arene bisphosphites and their catalytic applications in cross-coupling reactions

Palladium and platinum dichloride complexes of a series of symmetrically and unsymmetrically substituted 25,26;27,28-dibridged p-tert-butyl-calix[4]arene bisphosphites in which two proximal phenolic oxygen atoms of p-tert-butyl-or p-H-calix[4]arene are connected to a P(OR) ( R = substituted phenyl) moiety have been synthesized. The palladium dichloride complexes of calix[4]arene bisphosphites bearing sterically bulky aryl substituents undergo cyclometalation by C-C or C-H bond scission. An example of cycloplatinated complex is also reported. The complexes have been characterized by NMR spectroscopic and single crystal X-ray diffraction studies. During crystallization of the palladium dichloride complex of a symmetrically substituted calix[4]arene bisphosphite in dichloromethane, insertion of oxygen occurs into the Pd-P bond to give a P,O-coordinated palladium dichloride complex. The calix[4]arene framework in these bisphosphites and their metal complexes adopt distorted cone conformation; the cone conformation is more flattened in the metal complexes than in the free calix[4]arene bisphosphites. Some of these cyclometalated complexes proved to be active catalysts for Heck and Suzuki C-C cross-coupling reactions but, on an average, the yields are only modest. (C) 2011 Elsevier B.V. All rights reserved.

Dicopper(II) complexes showing DNA hydrolase activity and monomeric adduct formation with bis(4-nitrophenyl)phosphate

Ferromagnetic dicopper(II) complexes [Cu(2)(mu-O(2)CCH(3))(mu-OH)(L)(2)(mu-L(1))](PF(6))(2), where L = 1,10-phenanthroline (phen), L(1) = H(2)O in 1 and L = dipyrido[3,2-d:2',3'-f]quinoxaline (dpq), L(1) = CH(3)CN in 2, are prepared and structurally characterized. Crystals of 1 and 2 belong to the monoclinic space group of P2(1)/n and P2(1)/m, respectively. The copper(II) centers display distorted square-pyramidal geometry having a phenanthroline base and two oxygen atoms of the bridging hydroxo and acetate group in the basal plane. The fifth coordination site has weak axially bound bridging solvent molecule H(2)O in 1 and CH(3)CN in 2. The Cu center dot center dot center dot Cu distances are 3.034 and 3.046 angstrom in 1 and 2, respectively. The complexes show efficient hydrolytic cleavage of supercoiled pUC19 DNA as evidenced from the mechanistic studies that include T4 DNA ligase experiments. The binuclear complexes form monomeric copper(II) adducts [Cu(L)(2)(BNPP)](PF(6)) (L = phen, 3; dpq, 4) with bis(4-nitrophenyl)phosphate (BNPP) as a model phosphodiester. The crystal structures of 3 and 4 reveal distorted trigonal bipyramidal geometry in which BNPP binds through the oxygen atom of the phosphate. The kinetic data of the DNA cleavage reactions of the binuclear complexes under pseudo- and true-Michaelis-Menten conditions indicate remarkable enhancement in the DNA hydrolysis rate in comparison to the control data. (C) 2011 Elsevier B.V. All rights reserved.

Substitutional Self-Assembly of Alkanethiol and Selenol SAMs from a Lying-Down Doubly Tethered Butanedithiol SAM on Gold

Substitutional self-assembly of thiol and selenol SAMs from a lying-down phase of butanedithiol (C4DT) (SAM) were investigated using thiols, disulfide, and diselenide molecules. The intent was to address the question if formation of a lying-down dithiol phase is an impediment to formation of standing-up dithiol phases as it has been assumed. It is demonstrated that this is not the case, and the C4DT SAM, where both the sulfur atoms are chemisorbed on gold, is removed and replaced in all cases. Differences in substitution kinetics are observed.

3-(2-Bromoacetyl)-6-fluoro-2H-chromen-2-one

The non-H atoms of the title compound, C(11)H(6)BrFO(3), are essentially coplanar (r.m.s. deviation for all non-H atoms = 0.074 angstrom). In the crystal, the molecules are linked by C-H center dot center dot center dot O and C-H center dot center dot center dot Br interactions.

Structure and Bonding in Stannadiphospholes and their Dianions SnC(2)P(2)R(2)(m) (R=H, tBu m=0,-2): A Comparative Study with C(5)H(5)(+) and C(5)H(5)(-) Analogues

The potential energy surfaces of both neutral and dianionic SnC(2)P(2)R(2) (R=H, tBu) ring systems have been explored at the B3PW91/LANL2DZ (Sn) and 6-311 + G* (other atoms) level. In the neutral isomers the global minimum is a nido structure in which a 1,2-diphosphocyclobutadiene ring (1,2-DPCB) is capped by the Sn. Interestingly, the structure established by Xray diffraction analysis, for R=tBu, is a 1,3-DPCB ring capped by Sn and it is 2.4 kcal mol(-1) higher in energy than the 1,2-DPCB ring isomer. This is possibly related to the kinetic stability of the 1,3-DPCB ring, which might originate from the synthetic precursor ZrCp(2)tBu(2)C(2)P(2). In the case of the dianionic isomers we observe only a 6 pi-electron aromatic structure as the global minimum, similarly to the cases of our previously reported results with other types of heterodiphospholes.([1,4,19]) The existence of large numbers of cluster-type isomers in neutral and 6 pi-planar structures in the dianions SnC(2)P(2)R(2)(2-) (R=H, tBu) is due to 3D aromaticity in neutral clusters and to 2D pi aromaticity of the dianionic rings. Relative energies of positional isomers mainly depend on: 1) the valency and coordination number of the Sn centre, 2) individual bond strengths, and 3) the steric effect of tBu groups. A comparison of neutral stannadiphospholes with other structurally related C(5)H(5)(+) analogues indicates that Sn might be a better isolobal analogue to P(+) than to BH or CH(+). The variation in global minima in these C(5)H(5)(+) analogues is due to characteristic features such as 1) the different valencies of C, B, P and Sn, 2) the electron deficiency of B, 3) weaker p pi-p pi bonding by P and Sn atoms, and 4) the tendency of electropositive elements to donate electrons to nido clusters. Unlike the C5H5+ systems, all C(5)H(5)(-) analogues have 6 pi-planar aromatic structures as global minima. The differences in the relative ordering of the positional isomers and ligating properties are significant and depend on 1) the nature of the pi orbitals involved, and 2) effective overlap of orbitals.

Synthesis, Structure, Negative Thermal Expansion, and Photocatalytic Property of Mo Doped ZrV(2)O(7)

A new series of compounds identified in the phase diagram of ZrO(2)-V(2)O(8)-MoO(3) have been synthesized via the solution combustion method. Single crystals of one of the compounds in the series, ZrV(1.50)Mo(0.50)O(7.25), were grown by the melt-cool technique from the starting materials with double the MoO(3) quantity. The room temperature average crystal structure of the grown crystals was solved using the single crystal X-ray diffraction technique. The crystals belong to the cubic crystal system, space group Pa (3) over bar (No. 205) with a = 8.8969 (4) angstrom, V = 704.24 (6) angstrom(3), and Z = 4. The final R(1) value of 0.0213 was achieved for 288 independent reflections during the structure refinement. The Zr(4+) occupies the special position (4a) whereas V(5+) and Mo(6+) occupy two unique (8c) Wyckoff positions. Two fully occupied O atoms, (24d) and (4b), one partially occupied 0 atom (8c) have been identified for this molybdovanadate, which is a unique feature for these crystals. The structure is related to both ZrV(2)O(7) and cubic ZrMo(2)O(8). The temperature dependent single crystal studies show negative thermal expansion above 370 K. The compounds have been characterized by powder X-ray diffraction, solid-state UV-vis diffuse reflectance spectra, scanning electron microscopy (SEM), and transmission electron microscopy (TEM). The photocatalytic activity of these compounds has been investigated for the degradation of various dyes, and these compounds show specificity toward the degradation of non-azoic dyes.

Size dependent microstructure for Ag-Ni nanoparticles

The Ag-Ni system is characterized by large differences in atomic sizes (14%) and a positive heat of mixing (+23 kJ mol(-1)). The binary equilibrium diagram for this system therefore exhibits a large miscibility gap in both solid and liquid state. This paper explores the size-dependent changes in microstructure and the suppression of the miscibility gap which occurs when free alloy particles of nanometer size are synthesized by co-reduction of Ag and Ni metal precursors. The paper reports that complete mixing between Ag and Ni atoms could be achieved for smaller nanoparticles (<7 nm). These particles exhibit a single-phase solid solution with face-centered cubic (fcc) structure. With increase in size, the nanoparticles revealed two distinct regions. One of the regions is composed of pure Ag. This region partially surrounds a region of fcc solid solution at an early stage of decomposition. Experimental observations were compared with the results obtained from the thermodynamic calculations, which compared the free energies corresponding to a physical mixture of pure Ag and Ni phases and a fcc Ag-Ni solid solution for different particle sizes. Results from the theoretical calculations revealed that, for the Ag-Ni system, solid solution was energetically preferred over the physical mixture configuration for particle sizes of 7 nm and below. The experimentally observed two-phase microstructure for larger particles was thus primarily due to the growth of Ag-rich regions epitaxially on initially formed small fcc Ag-Ni nanoparticles. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Ferromagnetism Exhibited by Nanoparticles of Noble Metals

Gold nanoparticles with average diameters in the range 2.515 nm, prepared at the organic/aqueous interface by using tetrakis( hydroxymethyl) phosphonium chloride (THPC) as reducing agent, exhibit ferromagnetism whereby the saturation magnetization M(S) increases with decreasing diameter and varies linearly with the fraction of surface atoms. The value of M(S) is higher when the particles are present as a film instead of as a sol. Capping with strongly interacting ligands such as alkane thiols results in a higher M(S) value, which varies with the strength of the metal-sulfur bond. Ferromagnetism is also found in Pt and Ag nanoparticles prepared as sols, and the M(S) values vary as Pt > Au > Ag. A careful study of the temperature variation of the magnetization of Au nanoparticles, along with certain other observations, suggests that small bare nanoparticles of noble metals could indeed possess ferromagnetism, albeit weak, which is accentuated in the presence of capping agents, specially alkane thiols which form strong metal-sulfur bonds.

Microwave spectrum and structure of C(6)H(5)CCH center dot center dot center dot H(2)S complex

Rotational spectra of C(6)H(5)CCH center dot center dot center dot H(2)S, C(6)H(5)CCH center dot center dot center dot H(2)(34)S, C(6)H(5)CCH center dot center dot center dot HDS, C(6)H(5)CCH center dot center dot center dot D(2)S and C(6) H(5)CCD center dot center dot center dot H(2)S complexes have been observed using a pulsed nozzle Fourier transform microwave spectrometer. The observed spectrum is consistent with a structure in which hydrogen sulfide is located over the phenyl ring pi cloud and the distance between the centers of masses of the two monomers is 3.74 +/- 0.01 angstrom. In the complex, the H(2)S unit is shifted from the phenyl ring center towards the acetylene group. The vibrationally averaged structure has an effective Cs symmetry. Ab initio calculations were performed at MP2/aug-cc-pVDZ level of theory to locate the possible geometries of the complex. The calculations reveal the experimentally observed structure to be more stable than a coplanar arrangement of the monomers, which was observed for the C(6)H(5)CCH center dot center dot center dot H(2)O complex. Atoms in molecule theoretical analysis shows the presence of S-H center dot center dot center dot pi hydrogen bond. For the parent isotopologue, each transition frequency was found to split into two resulting from an interchange of the equivalent hydrogens of H(2)S unit in the complex. (C) 2011 Elsevier Inc. All rights reserved.