Successful molecular dynamics simulation of two zinc complexes bridged by a hydroxide in phosphotriesterase using the cationic dummy atom method

I report herein two 2.0 ns (1.0 fs time step) MD simulations of two zinc complexes bridged by a hydroxide in phosphotriesterase (PTE) employing the nonbonded method and the cationic dummy atom method that uses virtual atoms to impose orientational requirement for zinc ligands. The cationic dummy atom method was able to simulate the four-ligand coordination of the two zinc complexes in PTE. The distance (3.39 +/- 0.07A) between two nearby zinc ions in the time-average structure of PTE derived from the MD simulation using the cationic dummy atoms matched that in the X-ray structure (3.31 +/- 0.001A). Unequivocally, the time-average structure of PTE was able to fit into the experimentally determined difference electron density map of the corresponding X-ray structure. The results demonstrate the practicality of the cationic dummy atom method for MD simulations of zinc proteins bound with multiple zinc ions. In contrast, a 2.0 ns (1.0 fs time step) MD simulation using the nonbonded method revealed a striking difference in the active site between the X-ray structure and the time-average structure that was unable to fit into the density map of PTE. The results suggest that caution should be used in the MD simulations using the nonbonded method.     

Computer modeling and molecular dynamics simulations of ligand bound complexes of bovine angiogenin: dinucleotide topology at the active site of RNase a family proteins

We have undertaken the modeling of substrate-bound structures of angiogenin. In our recent study, we modeled the dinucleotide ligand binding to human angiogenin. In the present study, the substrates CpG, UpG, and CpA were docked onto bovine angiogenin. This was achieved by overcoming the problem of an obstruction to the B1 site by the C-terminus and identifying residues that bind to the second base. The modeled complexes retain biochemically important interactions. The docked models were subjected to 1 ns of molecular dynamics, and structures from the simulation were refined by using simulated annealing. Our models explained the enzyme's specificity for both B1 and B2 bases as observed experimentally. The nature of binding of the dinucleotide substrate was compared with that of the mononucleotide product. The models of these complexes were also compared with those obtained earlier with human angiogenin. On the basis of the simulations and annealed structures, we came up with a consensus topology of dinucleotide ligands that binds to human and bovine angiogenins. This dinucleotide conformation can serve as a starting model for ligand-bound complex structures for RNase A family of proteins. We demonstrated this capability by generating the complex structure of CpA bound to eosinophil-derived neurotoxin (EDN) by fitting the consensus topology of CpA to the crystal structure of native EDN.     

Kinetic nonoptimality and vibrational stability of proteins

Scaling of folding times in Go models of proteins and of decoy structures with the Lennard-Jones potentials in the native contacts reveal power law trends when studied under optimal folding conditions. The power law exponent depends on the type of native geometry. Its value indicates lack of kinetic optimality in the model proteins. In proteins, mechanical and thermodynamic stabilities are correlated.     

Bioremediation of DDT in soil by genetically improved strains of soil fungus Fusarium solani

Bioremediation of DDT in soil by genetically improved recombinants of the soil fungus Fusarium solani was studied. The parent strains were isolated from soil enriched with DDD or DDE (immediate anaerobic and aerobic degradation products of DDT), as further degradation of these products are slow processes compared to the parent compound. These naturally occurring strains isolated from soil, however, are poor degraders of DDT and differed in their capability to degrade its metabolites such as DDD, DDE, DDOH and DBP and other organochlorine pesticides viz. kelthane and lindane. Synergistic effect was shown by some of these strains, when grown together in the medium containing DDD and kelthane under mixed culture condition. No synergism in DDE degradation was observed with the strains isolated from enriched soil. DDD-induced proteins extracted from individual culture filtrate (exo-enzyme) when subjected to SDS-Polyacrylamide Gel Electrophoresis (SDS-PAGE) showed complementary polypeptide bands in these strains i.e., each strain produced distinct DDD degrading polypeptide bands and the recombinant or hybrid strains produced all of the bands of the two parents and degraded DDD better than the parental strains. Recombinant hybrid strains with improved dehalogenase activity were raised by parasexual hybridisation of two such complementary isolates viz. isolate 1(P-1) and 4(P-2) showing highest complementation and are compatible for hyphal fusion inducing heterokaryosis. These strains are genetically characterised as Kel⁺Ben^RDBP^-Lin^- and Kel^-Ben^rDBP⁺Lin⁺ respectively.Recombinants with mixed genotype, i.e., Kel⁺Ben^RDBP⁺Lin⁺ showing superior degradation quality for DDT were selected for bioremediation study. Recombination was confirmed by polypeptide band analysis of DDD induced exo-proteins from culture filtrate usingSDS-Polyacrylamide Gel Electrophoresis (PAGE) and RAPD (Random Amplified Polymorphic DNA) of genomic DNA using PCR (Polymerase Chain Reaction) technique. SDS-PAGE showed combination of DDD induced polypeptide bands characteristic of both the parents in the recombinants or the hybrids. PCR study showed the parent specific bands in the recombinant strains confirming gene transformation.     

Investigation of disaccharide recognition by molecularly imprinted polymers

The selectivity of carbohydrate-imprinted polymers for several disaccharides, namely cellobiose, maltose, lactose and gentiobiose, is investigated. An ternary ligand–Cu(II)–carbohydrate complex was formed in alkaline solution and captured afterwards in the polymer. The accessibility of the polymer matrix for disaccharides was investigated by HPLC analysis, refractometry and ¹H NMR spectroscopy applying excess of the original template during rebinding experiments under saturation conditions in unbuffered, aqueous solution at neutral pH and 20°C. The selective discrimination of the - and -glycosidic linkage of cellobiose and maltose is demonstrated. It is further shown, that the disaccharide-imprinted polymers slightly distinguish between the 1,4-- and the 1,6--glycosidic linkage of cellobiose and gentiobiose, while cellobiose and lactose are not selectively recognized. Due to the weak apparent binding constant of the functional Cu(II) monomers with the targeted disaccharides at physiological pH, the recognition process is dominated by the shape of the created imprinted cavity under the applied conditions.     

Relationship between enantioselectivity of alternative molecularly imprinted polymeric membranes and species of amino acid residues composing chiral recognition sites

Molecularly imprinted polymeric membranes with tetrapeptide residue H-Asp(OcHex)-Asp(OcHex)-Asp(OcHex)-Asp(OcHex)-CH₂- (DDDD) or H-Glu(OBzl)-Glu(OBzl)-Glu(OBzl)-Glu(OBzl)-C H₂- (EEEE) were prepared during membrane preparation (casting) processing in the presence of print molecules. The Boc-L-Trp imprinted polymeric membranes thus obtained showed adsorption selectivity toward Ac-L-Trp from its racemic mixtures. From adsorption isotherms of Ac-Trp, the chiral recognition site, that had been formed by the presence of print molecules in the membrane preparation process, exclusively recognized Ac-L-Trp that possessed the same configuration of the print molecule. The affinity constants between chiral recognition sites in the membrane and Ac-L-Trp was determined to be 1.00 × 10⁴ mol^–1 dm³ and 1.08 × 10⁴ mol^–1 dm³ for the DDDD and EEEE membranes, respectively. Enantioselective electrodialysis could be attained by applying an optimum potential difference to give permselectivity, with a value close to its adsorption selectivity.     

Genetic analysis of rabies virus isolates in the Philippines

To determine the genetic characteristics of the rabies virus in the Philippines, 59 rabies virus isolates were obtained from domestic rabid dogs and their partial nucleotide sequences of nucleoprotein (N) gene were compared. Based on comparison with reported sequences, phylogenetic analysis revealed that all isolates from the Philippines had close genetic relations and formed two subgroups. The Philippines isolates belonged to a different lineage from other Asian isolates but were closer to them than to terrestrial isolates and laboratory strains. Several specific nucleotide and amino acid substitutions were observed among the Philippines isolates. Our results suggest that rabies viruses in the Philippines might have a characteristic evolution.     

Interaction of the heart-specific LIM domain protein, FHL2, with DNA-binding nuclear protein, hNP220

Using a yeast two-hybrid library screen, we have identified that the heart specific FHL2 protein, four-and-a-half LIM protein 2, interacted with human DNA-binding nuclear protein, hNP220. Domain studies by the yeast two-hybrid interaction assay revealed that the second LIM domain together with the third and the fourth LIM domains of FHL2 were responsible to the binding with hNP220. Using green fluorescent protein (GFP)-FHL2 and blue fluorescent protein (BFP)-hNP220 fusion proteins co-expressed in the same cell, we demonstrated a direct interaction between FHL2 and hNP220 in individual nucleus by two-fusion Fluorescence Resonance Energy Transfer (FRET) assay. Besides, Western blot analysis using affinity-purified anti-FHL2 antipeptide antibodies confirmed a 32-kDa protein of FHL2 in heart only. Virtually no expression of FHL2 protein was detected in brain, liver, lung, kidney, testis, skeletal muscle, and spleen. Moreover, the expression of FHL2 protein was also detectable in the human diseased heart tissues. Our results imply that FHL2 protein can shuttle between cytoplasm and nucleus and may act as a molecular adapter to form a multicomplex with hNP220 in the nucleus, thus we speculate that FHL2 may be particularly important for heart muscle differentiation and the maintenance of the heart phenotype.     

Development of a novel diffusion-based method to estimate the size of the aggregated Abeta species responsible for neurotoxicity

beta-Amyloid peptide (Abeta) is the primary protein component of senile plaques in Alzheimer's disease and is believed to be responsible for the neurodegeneration associated with the disease. Abeta is toxic only when aggregated, however, the size and structure of the aggregated species associated with toxicity is unknown. In the present study, we developed a diffusion-based method to simultaneously separate and detect the biological activity of toxic Abeta oligomers and used the method to examine the relationship between size of aggregated protein and toxicity to SH-SY5Y cells. From these measurements, the effective diffusivity and hydrodynamic radius of the toxic oligomeric species of Abeta could be determined. A sensitivity analysis was performed to examine the effects of model assumptions used in data analysis on the effective diffusivity calculated. The method provides a new estimate of the size of small toxic Abeta species associated with fibril formation. This work contributes to our understanding of the relationship between Abeta structure and toxicity and with further refinements may aid in our ability to design agents which alter the Abeta aggregation/dissociation processes associated with neurotoxicity.     

Testing the relationship between morphological and molecular rates of change along phylogenies

Abstract.— Molecular evolution has been considered to be essentially a stochastic process, little influenced by the pace of phenotypic change. This assumption was challenged by a study that demonstrated an association between rates of morphological and molecular change estimated for "total-evidence" phylogenies, a finding that led some researchers to challenge molecular date estimates of major evolutionary radiations. Here we show that Omland's (1997) result is probably due to methodological bias, particularly phylogenetic nonindependence, rather than being indicative of an underlying evolutionary phenomenon. We apply three new methods specifically designed to overcome phylogenetic bias to 13 published phylogenetic datasets for vertebrate taxa, each of which includes both morphological characters and DNA sequence data. We find no evidence of an association between rates of molecular and morphological rates of change.