Predicting metal-binding site residues in low-resolution structural models

The accurate prediction of the biochemical function of a protein is becoming increasingly important, given the unprecedented growth of both structural and sequence databanks. Consequently, computational methods are required to analyse such data in an automated manner to ensure genomes are annotated accurately. Protein structure prediction methods, for example, are capable of generating approximate structural models on a genome-wide scale. However, the detection of functionally important regions in such crude models, as well as structural genomics targets, remains an extremely important problem. The method described in the current study, MetSite, represents a fully automatic approach for the detection of metal-binding residue clusters applicable to protein models of moderate quality. The method involves using sequence profile information in combination with approximate structural data. Several neural network classifiers are shown to be able to distinguish metal sites from non-sites with a mean accuracy of 94.5%. The method was demonstrated to identify metal-binding sites correctly in LiveBench targets where no obvious metal-binding sequence motifs were detectable using InterPro. Accurate detection of metal sites was shown to be feasible for low-resolution predicted structures generated using mGenTHREADER where no side-chain information was available. High-scoring predictions were observed for a recently solved hypothetical protein from Haemophilus influenzae, indicating a putative metal-binding site.

Automatic prediction of functional site regions in low-resolution protein structures

World-wide structural genomics initiatives are rapidly accumulating structures for which limited functional information is available. Additionally, state-of-the art structural prediction programs are now capable of generating at least low resolution structural models of target proteins. Accurate detection and classification of functional sites within both solved and modelled protein structures therefore represents an important challenge. We present a fully automatic site detection method, FuncSite, that uses neural network classifiers to predict the location and type of functionally important sites in protein structures. The method is designed primarily to require only backbone residue positions without the need for specific side-chain atoms to be present. In order to highlight effective site detection in low resolution structural models FuncSite was used to screen model proteins generated using mGenTHREADER on a set of newly released structures. We found effective metal site detection even for moderate quality protein models illustrating the robustness of the method.

Characterization of the layered structure in main chain dibenzo-18-crown-6 ether polymers by simultaneous WAXS/MAXS-SAXS/DSC measurements

The structure and thermal properties of polymers containing dibenzo-18-crown-6 ether units in the main chain linked to an aliphatic spacer of different lengths (C10-C14) is reported. X-ray diffraction patterns of all the studied samples exhibit a peak in the medium angle region, revealing the existence of a lamellar structure. Simultaneous calorimetry and small, medium (SAXS-MAXS) and wide (WAXS) X-ray measurements during cooling and subsequent heating of the samples reveal that a layer phase is formed upon cooling. In the case of the homopolymers, this phase is almost simultaneously accompanied by the appearance of some reflections in the wide angle region as an indication of lateral crystallization. However, by copolymerization, the formation of the layer phase is decoupled from lateral crystallization, being stable in a wide temperature region.

Film-forming polymers containing in the main-chain dibenzo crown ethers with aliphatic (C-10-C-16), aliphatic-aromatic, or oxyindole spacers

Novel, linear, soluble, high-molecular-weight, film-forming polymers and copolymers in which main-chain crown ether units alternate with aliphatic (C-10-C-16) units have been obtained for the first time from aromatic electrophilic substitution reactions of crown ethers by aliphatic dicarboxylic acids followed by reduction of the carbonyl groups. The crown ether unit is dibenzo-18-crown-6, dibenzo-21-crown-7, dibenzo-24-crown-8, or dibenzo-30-crown-10; the aliphatic spacer is derived from a dicarboxylic acid (sebacic, 1,12-dodecanedicarboxylic, hexadecanedioic or 1,4-phenylenediacetic acids). The reactions were performed at 35 degrees C in a mixture of methanesulfonic acid (MSA) with phosphorus pentoxide, 12:1 (w/w), (Eaton's reagent). The carbonyl groups in the polyketones obtained were completely reduced to methylene linkages by treatment at room temperature with triethylsilane in a mixture of trifluoroacetic acid and dichloromethane. Polymers containing in the main chain crown ethers alternating with oxyindole fragments were prepared by one-pot condensation of crown ethers with isatin in a medium of Eaton's reagent. A possible reaction mechanism is suggested. According to IR and NMR analyses, the polyacylation reactions lead to the formation of isomeric (syn/anti-substituted) crown ether units in the main chain. The polymers obtained were soluble in the common organic solvents, and flexible transparent films could be cast from the solutions. DSC and X-ray studies of the polymers with "symmetrical" crown ethers reveal the presence of the endotherms corresponding to the supramolecular assemblies.

Crystallization of random aromatic copolyesters containing flexible spacer chains and side-groups

The crystallization behaviour of a series of random copolymers of varying chemical composition is reported. For polymers containing a high proportion of alternating rigid aromatic units and flexible spacers, conventional liquid crystalline and crystalline phase behaviour is observed. The introduction of a substantial fraction of a second shorter rigid unit containing side-chains leads to a broad endotherm in the d.s.c. scan covering some 150°C. Subsequent isothermal crystallization at any point within the broad endotherm leads to the generation of sharp endotherms at temperatures just above the recrystallization temperature. We attribute this behaviour to the crystallization of clusters of molecules containing similar random sequences. Such crystals are non-periodic along the chain direction.

Molecular design of a discrete chain-folding polyimide for controlled inkjet deposition of supramolecular polymers

A supramolecular polymer based upon two complementary polymer components is formed by sequential deposition from solution in THF, using a piezoelectric drop-on-demand inkjet printer. Highly efficient cycloaddition or ‘click’ chemistry afforded a well-defined poly(ethylene glycol) featuring chain-folding diimide end groups, which possesses greatly enhanced solubility in THF relative to earlier materials featuring random diimide sequences. Blending the new polyimide with a complementary poly(ethylene glycol) system bearing pyrene end groups (which bind to the chain-folding diimide units) overcomes the limited solubility encountered previously with chain-folding polyimides in inkjet printing applications. The solution state properties of the resulting polymer blend were assessed via viscometry to confirm the presence of a supramolecular polymer before depositing the two electronically complementary polymers by inkjet printing techniques. The novel materials so produced offer an insight into ways of controlling the properties of printed materials through tuning the structure of the polymer at the (supra)molecular level.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

Effect of acyl chain length on transfection efficiency and toxicity of polyethylenimine

Polyethylenimine (PEI) is an efficient nonviral gene delivery vector because of its high buffering capacity and DNA condensation ability. In our study, the amino groups on the polymeric backbone were acylated using acetic or propionic anhydride to alter the protonation behaviour and the hydrophilic/hydrophobic balance of the polymer. The concentration of acylated primary amines was determined using trinitrobenzene sulphonic acid assay. Results showed that our modified polymers had lower buffering capacities in solutions compared to PEI. The polymers were complexed with plasmid encoding enhanced green fluorescent protein at three different ratios (1:1, 1:2 and 1:10 w/w DNA to polymer) to form polyplexes and their toxicities and transfection efficiencies were evaluated in HEK 293 cells. Acylation reduced the number of primary amines on the polymer and the surface charge, improving haemocompatibility and reducing cytotoxicity. The reduction in the concentration of amino groups helped to optimise DNA compaction and facilitated polyplex dissociation in the cell, which increased transfection efficiency of the modified polymers compared to the parent polymer. Polymers with buffering capacities greater than 50% and less than 80% relative to PEI, showed higher transfection efficiencies than PEI. The propionic anhydride modified polymers had appropriate interactions with DNA which provided both DNA compaction and polyplex dissociation. These systems interacted better with the cell membrane because of their slightly higher lipophilicity and formed polyplexes which were less cytotoxic than polyplexes of acetic anhydride modified polymers. Among the vectors tested, 1:0.3 mol/mol PEI:propionic anhydride in a 1:2 w/w DNA:polymer composition provided the best transfection system with improved transfection efficiency and reduced cytotoxicity.

Construction of coordination polymers of cadmium(II) with mixed hexamethylenetetramine and terephthalate or thiocyanate ligands

Three new polymeric complexes [Cd(hmt)(SCN)(2)(H2O)(2)](n) (1), [Cd-3(mu(2)-hmt)(2)(SCN)(6)(H2O)(2)](n) (2), and [Cd-2(hmt)(2)(tP)(2)(H2O)(6)](n) (3) [hmt = hexamethylenetetramine, tp = terephthalate] have been synthesized and characterized by single crystal X-ray diffraction. Both the compounds 1 and 2 are 1-D polymers where Cd units are linked by double end-to-end thiocyanate bridges but in 2 the chain is wider containing three cadmium atoms instead of one as found in 1. In both compounds the coordination environment around cadmium atom is distorted octahedral. Compound 3 is a three-dimensional polymer having water-filled microporous channels. Both tp and brut are mu(2)-bridged. One of the acid groups of tp is coordinated in chelating bidentate and the other in monodentate fashion. Half of its Cd atoms are hexa-coordinated and the rest are hepta-coordinated. Thermogravimetric analysis and X-ray diffraction study of 3 show that its framework remains intact upon removal of water molecules. The flexibility of coordination number around cadmium atoms (six or seven) probably plays an important role in establishing the rigidity of the framework. (C) 2003 Elsevier Ltd. All rights reserved.

Syntheses and crystal structures of two new coordination polymers of Cd(II) using organic spacer and inorganic-bridging ligands

Two new cadmium (II) complexes [Cd(hmt)(dca)(2)] (n) (1) and [Cd-3(hmt)(2)(SeCN)(6)(H2O)(2)] (n) (2) (hmt=hexamethylenetetramine, dca=dicyanamide) have been synthesized and characterized by X-ray single-crystal analysis. The complex 1 is a 2D rectangular grid of octahedral cadmium (II) with CdN6 chromophore where cadmium centers are doubly bridged by dicyanamide and hmt along a-axis, which are interlinked by dicyanamide running along c-axis. Whereas, complex 2 is a 1D chain of octahedral cadmium (II) with a three-leg ladder topology running along a-axis. The Cd(II) centers are doubly bridged through SeCN (infinite rail) along a-axis and singly bridged by hmt (two-step rung) along c-axis, having cadmium centers with CdSe2N3O and CdSe2N4 chromophores. The adjacent chains through H-bonding between coordinated water and hmt, and (SeSe)-Se-... interaction are extended to 2D supramolecular architecture.

Modeling entangled dynamics: comparison between stochastic single-chain and multichain models

To test the effectiveness of stochastic single-chain models in describing the dynamics of entangled polymers, we systematically compare one such model; the slip-spring model; to a multichain model solved using stochastic molecular dynamics(MD) simulations (the Kremer-Grest model). The comparison involves investigating if the single-chain model can adequately describe both a microscopic dynamical and a macroscopic rheological quantity for a range of chain lengths. Choosing a particular chain length in the slip-spring model, the parameter values that best reproduce the mean-square displacement of a group of monomers is determined by fitting toMDdata. Using the same set of parameters we then test if the predictions of the mean-square displacements for other chain lengths agree with the MD calculations. We followed this by a comparison of the time dependent stress relaxation moduli obtained from the two models for a range of chain lengths. After identifying a limitation of the original slip-spring model in describing the static structure of the polymer chain as seen in MD, we remedy this by introducing a pairwise repulsive potential between the monomers in the chains. Poor agreement of the mean-square monomer displacements at short times can be rectified by the use of generalized Langevin equations for the dynamics and resulted in significantly improved agreement.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Self-assembled molecular complexes and coordination polymers of CdII, hexamine, and monocarboxylates: structural analysis and theoretical studies of supramolecular interactions

Four new cadmium(II) complexes [Cd-2(bz)(4)(H2O)(4)(mu 2-hmt)]center dot Hbz center dot H2O (1), [Cd-3(bz)(6)(H2O)(6)(mu 2-hmt)(2)]center dot 6H(2)O (2), [Cd(pa)(2)(H2O)(mu(2)-hmt)](n) (3), and {[Cd-3(ac)(6)(H2O)(3)(mu(3)-hmt)(2)]center dot 6H(2)O}(n) (4) with hexamine (hmt) and monocarboxylate ions, benzoate (bz), phenylacetate (pa), or acetate (ac) have been synthesized and characterized structurally. Structure determinations reveal that 1 is dinuclear, 2 is trinuclear, 3 is a one-dimensional (1D) infinite chain, and 4 is a two-dimensional (2D) polymer with fused hexagonal rings consisting of Cd-II and hmt. All the Cd-II atoms in the four complexes (except one CdII in 2) possess seven-coordinate pentagonal bipyramidal geometry with the various chelating bidentate carboxylate groups in equatorial sites. One of the CdII ions in 2, a complex that contains two monodentate carboxylates is in a distorted octahedral environment. The bridging mode of hmt is mu 2- in complexes 1-3 but is mu 3- in complex 4. In all complexes, there are significant numbers of H-bonds, C-H/pi, and pi-pi interactions which play crucial roles in forming the supramolecular networks. The importance of the noncovalent interactions in terms of energies and geometries has been analyzed using high level ab initio calculations. The effect of the cadmium coordinated to hmt on the energetic features of the C-H/pi interaction is analyzed. Finally, the interplay between C-H/pi and pi-pi interactions observed in the crystal structure of 3 is also studied.