Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, Engyodontium album BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca2? binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.

Molecular Characterisation and Phylogenetic Analysis of a Novel Isoform of Hepatic Antimicrobial Peptide, Hepcidin (Zc-hepc1), from the Coral Fish Moorish idol, Zanclus cornutus (Linnaeus,1758)

Hepcidin is a family of short cysteine-rich antimicrobial peptides (AMPs) participating in various physiological functions with inevitable role in host immune responses. Present study deals with identification and characterisation of a novel hepcidin isoform from coral fish Zanclus cornutus. The 81 amino acid (aa) preprohepcidin obtained from Z. cornutus consists of a hydrophobic aa rich 22 mer signal peptide, a highly variable proregion of 35 aa and a bioactive mature peptide with 8 conserved cysteine residues which contribute to the disulphide back bone. The mature hepcidin, Zc-hepc1 has a theoretical isoelectric point of 7.46, a predicted molecular weight of 2.43 kDa and a net positive charge of ?1. Phylogenetic analysis grouped Z. cornutus hepcidin with HAMP2 group hepcidins confirming the divergent evolution of hepcidin-like peptide in fishes. Zc-hepc1 can attain a b-hairpin-like structure with two antiparallel b-sheets. This is the first report of an AMP from the coral fish Z. cornutus.

Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, Engyodontium album BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca2? binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.

Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

A comparative analysis of two methods for the calculation of electric-field-induced perturbations to molecular vibration

Two common methods of accounting for electric-field-induced perturbations to molecular vibration are analyzed and compared. The first method is based on a perturbation-theoretic treatment and the second on a finite-field treatment. The relationship between the two, which is not immediately apparent, is made by developing an algebraic formalism for the latter. Some of the higher-order terms in this development are documented here for the first time. As well as considering vibrational dipole polarizabilities and hyperpolarizabilities, we also make mention of the vibrational Stark effec

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

High performance subgraph mining in molecular compounds

Structured data represented in the form of graphs arises in several fields of the science and the growing amount of available data makes distributed graph mining techniques particularly relevant. In this paper, we present a distributed approach to the frequent subgraph mining problem to discover interesting patterns in molecular compounds. The problem is characterized by a highly irregular search tree, whereby no reliable workload prediction is available. We describe the three main aspects of the proposed distributed algorithm, namely a dynamic partitioning of the search space, a distribution process based on a peer-to-peer communication framework, and a novel receiver-initiated, load balancing algorithm. The effectiveness of the distributed method has been evaluated on the well-known National Cancer Institute’s HIV-screening dataset, where the approach attains close-to linear speedup in a network of workstations.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

Analytical potentials for triatomic molecules IX. The prediction of anharmonic force constants from potential energy surfaces based on harmonic force fields and dissociation energies for SO2 and O3

Analytical potential energy functions which are valid at all dissociation limits have been derived for the ground states of SO2 and O3. The procedure involves minimizing the errors between the observed vibrational spectra and spectra calculated by a variational procedure. Good agreement is obtained between the observed and calculated spectra for both molecules. Comparisons are made between anharmonic force fields, previously determined from the spectral data, and the force fields obtained by differentiating the derived analytical functions at the equilibrium configurations.

General derivative relations for anharmonic force fields

The brace notation, introduced by Allen and Csaszar (1993, J. chem. Phys., 98, 2983), provides a simple and compact way to deal with derivatives of arbitrary non-tensorial quantities. One of its main advantages is that it builds the permutational symmetry of the derivatives directly into the formalism. The brace notation is applied to formulate the general nth-order Cartesian derivatives of internal coordinates, and to provide closed forms for general, nth-order transformation equations of anharmonic force fields, expressed as Taylor series, from internal to Cartesian or normal coordinate spaces.

Molecular guided catalytic hydrogenation by micelle-hosted Pd nanoparticle in supercritical CO2

Homogeneous dispersion of microemulsion containing palladium nanoparticles in scCO(2) is, for the first time, observed via sapphire window reactor and these particles show an unusual reluctance for double bond hydrogenation of citral aldehyde at hydrophobic end rather than hydrophilic end (high regioselectivity) owing to the unique micelle environment in supercritical carbon dioxide that guide a head-on attack of the molecule.

Chalcogenide-halides of niobium (V). 1. Gas-phase structures of NbOBr3, NbSBr3, and NbSCl3. 2. Matrix infrared spectra and vibrational force fields of NbOBr3, NbSBr3, NbSCl3, and NbOCl3

The molecular structures of NbOBr3, NbSCl3, and NbSBr3 have been determined by gas-phase electron diffraction (GED) at nozzle-tip temperatures of 250 degreesC, taking into account the possible presence of NbOCl3 as a contaminant in the NbSCl3 sample and NbOBr3 in the NbSBr3 sample. The experimental data are consistent with trigonal-pyramidal molecules having C-3v symmetry. Infrared spectra of molecules trapped in argon or nitrogen matrices were recorded and exhibit the characteristic fundamental stretching modes for C-3v species. Well resolved isotopic fine structure (Cl-35 and Cl-37) was observed for NbSCl3, and for NbOCl3 which occurred as an impurity in the NbSCl3 spectra. Quantum mechanical calculations of the structures and vibrational frequencies of the four YNbX3 molecules (Y = O, S; X = Cl, Br) were carried out at several levels of theory, most importantly B3LYP DFT with either the Stuttgart RSC ECP or Hay-Wadt (n + 1) ECP VDZ basis set for Nb and the 6-311 G* basis set for the nonmetal atoms. Theoretical values for the bond lengths are 0.01-0.04 Angstrom longer than the experimental ones of type r(a), in accord with general experience, but the bond angles with theoretical minus experimental differences of only 1.0-1.5degrees are notably accurate. Symmetrized force fields were also calculated. The experimental bond lengths (r(g)/Angstrom) and angles (angle(alpha)/deg) with estimated 2sigma uncertainties from GED are as follows. NbOBr3: r(Nb=O) = 1.694(7), r(Nb-Br) = 2.429(2), angle(O=Nb-Br) = 107.3(5), angle(Br-Nb-Br) = 111.5(5). NbSBr3: r(Nb=S) = 2.134(10), r(Nb-Br) = 2.408(4), angle(S=Nb-Br) = 106.6(7), angle(Br-Nb-Br) = 112.2(6). NbSCl3: Nb=S) = 2.120(10), r(Nb-Cl) = 2.271(6), angle(S=Nb-Cl) = 107.8(12), angle(Cl-Nb-Cl) = 111.1(11).

Carbohydrate-based micelle clusters which enhance hydrophobic drug bioavailability by up to 1 order of magnitude

Amphiphilic chitosan-based polymers (M-w < 20 kDa) self-assemble in aqueous media at low micromolar concentrations to give previously unknown micellar clusters of 100-300 nm in size. Micellar clusters comprise smaller 10-30 nm aggregates, and the nanopolarity/drug incorporation efficiency of their hydrophobic domains can be tailored by varying the degree of lipidic derivatization and molecular weight of the carbohydrate. The extent of drug incorporation by these novel micellar clusters is 1 order of magnitude higher than is seen with triblock copolymers, with molar polymer/drug ratios of 1:48 to 1:67. On intravenous injection, the pharmacodynamic activity of a carbohydrate propofol formulation is increased by 1 order of magnitude when compared to a commercial emulsion formulation, and on topical ocular application of a carbohydrate prednisolone formulation, initial drug aqueous humor levels are similar to those found with a 10-fold dose of prednisolone suspension.

Molecular dynamics simulations of threadlike cetyltrimethylammonium chloride micelles: effects of sodium chloride and sodium salicylate salts

We use atomistic molecular dynamics simulations to probe the effects of added sodium chloride (NaCl) and sodium salicylate (NaSal) salts on the spherical-to-threadlike micelle shape transition in aqueous solutions of cetyltrimethylammonium chloride (CTAC) surfactants. Long threadlike micelles are found to be unstable and break into spherical micelles at low concentrations or NaCl, but remain stable for 20 ns above a threshold value of [NaCl] approximate to 3.0 M, which is about 2.5 times larger than the experimental salt concentration at which the transition between spherical and rodlike micelles occurs. The chloride counterions associate weakly oil the surface of the CTAC micelles with the degree of counterion dissociation decreasing slightly with increasing [NaCl] on spherical micelles, but dropping significantly on the threadlike micelles tit high [NaCl]. This effect indicates that the electrolyte ions drive the micellar shape transition by screening the electrostatic repulsions between the micellar headgroups, The aromatic salicylate counterions, on the other hand, penetrate inside the micelle with their hydrophilic groups staying in the surfactant headgroup region and the hydrophobic groups partially embedded into the hydrophobic core of the micelle. The strong association of the salicylate ions with the surfactant headgroups leads to dense packing of the surfactant molecules, which effectively reduces the surface area per surfactant, and increases intramicellar ordering of the surfactant headgroups, favoring the formation of long threadlike micelles. Simulation predictions of the geometric and electrostatic properties of the spherical and threadlike micelles are in good agreement with experiments.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.