Use of a novel Hill-climbing genetic algorithm in protein folding simulations

We have developed a novel Hill-climbing genetic algorithm (GA) for simulation of protein folding. The program (written in C) builds a set of Cartesian points to represent an unfolded polypeptide's backbone. The dihedral angles determining the chain's configuration are stored in an array of chromosome structures that is copied and then mutated. The fitness of the mutated chain's configuration is determined by its radius of gyration. A four-helix bundle was used to optimise simulation conditions, and the program was compared with other, larger, genetic algorithms on a variety of structures. The program ran 50% faster than other GA programs. Overall, tests on 100 non-redundant structures gave comparable results to other genetic algorithms, with the Hill-climbing program running from between 20 and 50% faster. Examples including crambin, cytochrome c, cytochrome B and hemerythrin gave good secondary structure fits with overall alpha carbon atom rms deviations of between 5 and 5.6 Angstrom with an optimised hydrophobic term in the fitness function. (C) 2003 Elsevier Ltd. All rights reserved.

Use of essential and molecular dynamics to study gamma B-crystallin unfolding after non-enzymic post-translational modifications

Essential and Molecular Dynamics (ED/MD) have been used to model the conformational changes of a protein implicated in a conformational disease-cataract, the largest cause of blindness in the world-after non-enzymic post-translational modification. Cyanate modification did not significantly alter flexibility, while the Schiff's base adduct produced a more flexible N-terminal domain, and intra-secondary structure regions, than either the cyanate adduct or the native structure. Glycation also increased linker flexibility and disrupted the charge network. A number of post-translational adducts showed structural disruption around Cys15 and increased linker flexibility; this may be important in subsequent protein aggregation. Our modelling results are in accord with experimental evidence, and show that ED/MD is a useful tool in modelling conformational changes in proteins implicated in disease processes. (C) 2003 Published by Elsevier Ltd.

Computer modelling and estimation of recruitment patterns of non-branching coral colonies at three sites in Wakatobi Marine Park, S.E. Sulawesi, Indonesia; implications for coral reef conservation

We have studied growth and estimated recruitment of massive coral colonies at three sites, Kaledupa, Hoga and Sampela, separated by about 1.5 km in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. There was significantly higher species richness (P<0.05), coral cover (P<0.05) and rugosity (P<0.01) at Kaledupa than at Sampela. A model for coral reef growth has been developed based on a rational polynomial function, where dx/dt is an index of coral growth with time; W is the variable (for example, coral weight, coral length or coral area), up to the power of n in the numerator and m in the denominator; a1……an and b1…bm are constants. The values for n and m represent the degree of the polynomial, and can relate to the morphology of the coral. The model was used to simulate typical coral growth curves, and tested using published data obtained by weighing coral colonies underwater in reefs on the south-west coast of Curaçao [‘Neth. J. Sea Res. 10 (1976) 285’]. The model proved an accurate fit to the data, and parameters were obtained for a number of coral species. Surface area data was obtained on over 1200 massive corals at three different sites in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. The year of an individual's recruitment was calculated from knowledge of the growth rate modified by application of the rational polynomial model. The estimated pattern of recruitment was variable, with little numbers of massive corals settling and growing before 1950 at the heavily used site, Sampela, relative to the reef site with little or no human use, Kaledupa, and the intermediate site, Hoga. There was a significantly greater sedimentation rate at Sampela than at either Kaledupa (P<0.0001) or Hoga (P<0.0005). The relative mean abundance of fish families present at the reef crests at the three sites, determined using digital video photography, did not correlate with sedimentation rates, underwater visibility or lack of large non-branching coral colonies. Radial growth rates of three genera of non-branching corals were significantly lower at Sampela than at Kaledupa or at Hoga, and there was a high correlation (r=0.89) between radial growth rates and underwater visibility. Porites spp. was the most abundant coral over all the sites and at all depths followed by Favites (P<0.04) and Favia spp. (P<0.03). Colony ages of Porites corals were significantly lower at the 5 m reef flat on the Sampela reef than at the same depth on both other reefs (P<0.005). At Sampela, only 2.8% of corals on the 5 m reef crest are of a size to have survived from before 1950. The Scleractinian coral community of Sampela is severely impacted by depositing sediments which can lead to the suffocation of corals, whilst also decreasing light penetration resulting in decreased growth and calcification rates. The net loss of material from Sampela, if not checked, could result in the loss of this protective barrier which would be to the detriment of the sublittoral sand flats and hence the Sampela village.

Multivariate relationships between groundwater chemistry and toxicity in an urban aquifer

Multivariate statistical methods were used to investigate file Causes of toxicity and controls on groundwater chemistry from 274 boreholes in an Urban area (London) of the United Kingdom. The groundwater was alkaline to neutral, and chemistry was dominated by calcium, sodium, and Sulfate. Contaminants included fuels, solvents, and organic compounds derived from landfill material. The presence of organic material in the aquifer caused decreases in dissolved oxygen, sulfate and nitrate concentrations. and increases in ferrous iron and ammoniacal nitrogen concentrations. Pearson correlations between toxicity results and the concentration of individual analytes indicated that concentrations of ammoinacal nitrogen, dissolved oxygen, ferrous iron, and hydrocarbons were important where present. However, principal component and regression analysis suggested no significant correlation between toxicity and chemistry over the whole area. Multidimensional Scaling was used to investigate differences in sites caused by historical use, landfill gas status, or position within the sample area. Significant differences were observed between sites with different historical land use and those with different gas status. Examination of the principal component matrix revealed that these differences are related to changes in the importance of reduced chemical species.

Some transition metal complexes derived from mono- and di-ethynyl perfluorobenzenes

Transition metal alkynyl complexes containing perfluoroaryl groups have been prepared directly from trimethylsilyl-protected mono- and di-ethynyl perfluoroarenes by simple desilylation/metallation reaction sequences. Reactions between Me3SiC CC6F5 and RuCl(dppe)Cp'[Cp' = Cp, Cp*] in the presence of KF in MeOH give the monoruthenium complexes Ru(C CC6F5)(dppe)Cp'[Cp' = Cp (2); Cp* (3)], which are related to the known compound Ru(C CC6F5)(PPh3)(2)Cp (1). Treatment of Me3SiC CC6F5 with Pt-2(mu-dppm)(2)Cl-2 in the presence of NaOMe in MeOH gave the bis(alkynyl) complex Pt-2(mu-dppm)(2)(C CC6F5)(2) (4). The Pd(0)/Cu(I)-catalysed reactions between Au(C CC6F5)(PPh3) and Mo( CBr)(CO)(2) Tp* [Tp* = hydridotris(3.5-dimethylpyrazoyl)borate], Co-3(mu(3)-CBr)(mu-dppm)(CO)(7) or IC CFc [Fc = (eta(5)-C5H4)FeCp] afford Mo( CC CC6F5)(CO)(2)Tp* (5), Co-3(mu 3-CC CC6F5)(mu-dppm)(CO)(7) (6) and FcC CC CC6F5 (7), respectively. The diruthenium complexes 1,4-{Cp'(PP)RuC C}(2)C6F4 [(PP)Cp'=(PPh3)(2)Cp (8); (dppe)Cp (9); (dppe)Cp* (10)] are prepared from 1,4-(Me3SiC C)(2)C6F4 in a manner similar to that described for the monoruthenium complexes 1-3. The non-fluorinated complexes 1,4-{Cp'(PP)RuC C}(2)C6H4 [(PP)Cp' = (PPh3)(2)Cp (11); ( dppe) Cp (12); ( dppe) Cp* (13)], prepared for comparison, are obtained from 1,4-(Me3SiC C)(2)C6H4. Spectro-electrochemical studies of the ruthenium aryl and arylene alkynyl complexes 2-3 and 8-13, together with DFT-based computational studies on suitable model systems, indicate that perfluorination of the aromatic ring has little effect on the electronic structures of these compounds, and that the frontier orbitals have appreciable diethynylphenylene character. Molecular structure determinations are reported for the fluoroaromatic complexes 1, 2, 3, 6 and 10.

Molecular aluminum hydrides identified by inelastic neutron scattering during H-2 regeneration of catalyst-doped NaAlH4

Catalyst-doped sodium aluminum hydrides have been intensively studied as solid hydrogen carriers for onboard proton-exchange membrane (PEM) fuel cells. Although the importance of catalyst choice in enhancing kinetics for both hydrogen uptake and release of this hydride material has long been recognized, the nature of the active species and the mechanism of catalytic action are unclear. We have shown by inelastic neutron scattering (INS) spectroscopy that a volatile molecular aluminum hydride is formed during the early stage of H-2 re-eneration of a depleted, catalyst-doped sodium aluminum hydride. Computational modeling of the INS spectra suggested the formation of AlH3 and oligomers (AlH3)(n) (Al2H6, Al3H9, and Al4H12 clusters), which are pertinent to the mechanism of hydrogen storage. This paper demonstrates, for the first time, the existence of these volatile species.

Synthesis, structure and solution chemistry of mixed-ligand oxovanadium(IV) and oxovanadium(V) complexes incorporating tridentate ONO donor hydrazone ligands

In the title family, the ONO donor ligands are the acetylhydrazones of salicylaidehyde (H2L1) and 2-hydroxyacetophenone (H2L2) (general abbreviation, H2L). The reaction of bis(acetylacetonato)oxovanadium(IV) with a mixture of tridentate H2L and a bidentate NN donor [e.g., 2,2'-bipyridine(bpy) or 1,10-phenanthroline(phen), hereafter B] ligands in equimolar ratio afforded the tetravalent complexes of the type [(VO)-O-IV(L)(B)]; complexes (1)-(4) whereas, if B is replaced by 8-hydroxyquinoline(Hhq) (which is a bidentate ON donor ligand), the above reaction mixture yielded the pentavalent complexes of the type [(VO)-O-V(L)(hq)]; complexes (5) and (6). Aerial oxygen is most likely the oxidant (for the oxidation of V-IV -> V-V) in the synthesis of pentavalent complexes (5) and (6). [(VO)-O-IV(L)(B)] complexes are one electron paramagnetic and display axial EPR spectra, while the [(VO)-O-V(L)(hq)] complexes are diamagnetic. The X-ray structure of [(VO)-O-V(L-2)(hq)] (6) indicates that H2L2 ligand is bonded with the vanadium meridionally in a tridentate dinegative fashion through its phenolic-O, enolic-O and imine-N atoms. The general bond length order is: oxo < phenolato < enolato. The V-O (enolato) bond is longer than V-O (phenolato) bond by similar to 0.07 angstrom and is identical with V-O (carboxylate) bond. H-1 NMR spectrum of (6) in CDCl3 solution indicates that the binding nature in the solid state is also retained in solution. Complexes (1)(4) display two ligand-field transitions in the visible region near 820 and 480 nm in DMF solution and exhibit irreversible oxidation peak near +0.60 V versus SCE in DMSO solution, while complexes (5) and (6) exhibit only LMCT band near 535 nm and display quasi-reversible one electron reduction peak near -0.10 V versus SCE in CH2Cl2 solution. The VO3+-VO2+ E-1/2 values shift considerably to more negative values when neutral NN donor is replaced by anionic ON donor species and it also provides better VO3+ binding via phenolato oxygen. For a given bidentate ligand, E-1/2 increases in the order: (L-2)(2-) < (L-1)(2-). (c) 2004 Elsevier B.V. All rights reserved.