Determining molecular structures and conformations directly from electron diffraction using a genetic algorithm.

A global optimization strategy, based upon application of a genetic algorithm (GA), is demonstrated as an approach for determining the structures of molecules possessing significant conformational flexibility directly from gas-phase electron diffraction data. In contrast to the common approach to molecular structure determination, based on trial-and-error assessment of structures available from quantum chemical calculations, the GA approach described here does not require expensive quantum mechanical calculations or manual searching of the potential energy surface of the sample molecule, relying instead upon simple comparison between the experimental and calculated diffraction pattern derived from a proposed trial molecular structure. Structures as complex as all-trans retinal and p-coumaric acid, both important chromophores in photosensing processes, may be determined by this approach. In the examples presented here, we find that the GA approach can determine the correct conformation of a flexible molecule described by 11 independent torsion angles. We also demonstrate applications to samples comprising a mixture of two distinct molecular conformations. With these results we conclude that applications of this approach are very promising in elucidating the structures of large molecules directly from electron diffraction data.     

Structural determinants of steroids for cytochrome P450 3A4-mediated metabolism

Cytochrome P450 3A4 (CYP3A4), a major member of cytochrome P450 isoenzymes, metabolizes the majority of steroids in 6β-position. For the purpose of determining requisite structural features of a series of structurally related steroids for CYP3A4-mediated metabolism, three-dimensional pharmacophore modeling as well as electrotopological state map were conducted for 15 steroids. Though prior studies speculated that the chemical reactivity of the allylic 6β-position might have a greater influence on CYP3A4 selective 6-hydroxylation than steric constraints in the enzyme, our results reveal that for CYP3A4 steroidal substrates, it is not the chemical reactivity of atoms at 6β-site, but the pharmacophoric features, i.e. the two hydrophobic rings together with two H-bond donors, that act as the key factors responsible for determining the CYP3A4 selective 6-hydroxylation of steroids.     

Nanocomposite materials from latex and cellulose whiskers

Composite materials were processed by casting a mixture of aqueous suspensions of latex and microfibrils. These microfibrils, or whiskers, are extracted from a sea animal and are monocrystals of cellulose, with an aspect ratio around 100 and an average diameter of 20 nm. It has been found that the mechanical properties (shear modulus) are increased by more than two orders of magnitude in the rubbery state of the polymeric matrix, when the whisker content was 6% (w/w). This very large effect is discussed on the basis of different types of mechanical models and it is concluded that these whiskers form a rigid network, probably linked by hydrogen bonds. The formation of this network is assumed to be governed by a percolation mechanism.     

An investigation into the construction of molecular models by the template joining method

Summary The results of a wide-ranging investigation into some of the different methods available for performing the joining of templates to build molecular models show that the choice of algorithm can significantly affect the quality of the results obtained, and different algorithms are most suited to particular categories of join.     

Cavity search: An algorithm for the isolation and display of cavity-like binding regions

Summary A set of algorithms designed to enhance the display of protein binding cavities is presented. These algorithms, collectively entitled CAVITY SEARCH, allow the user to isolate and fully define the extent of a particular cavity. Solid modeling techniques are employed to produce a detailed cast of the active site region, which can then be color-coded to show electrostatic and steric interactions between the protein cavity and a bound ligand.     

Engineering cytochrome P-450s

Cytochrome P-450 isozymes represent a critical component of nature’s spectrum of detoxification catalysts that could be exploited for bioremediation. The ethanol-inducible human cytochrome P-450 2E1 serves as a model eukaryotic P-450 that complements the bacterial P-450 cam in dehalogenation and detoxification of environmental pollutants. We explored the construction of novel chimeric P-450s using cytochrome P-450 camC and 2E1 genes. For construction of chimera 1 (478 amino acids, 55.14 kDa), 145 amino acids from the N-terminus of P-450 2E1 protein (493 amino acids, 56.84 kDa) were replaced with 130 amino acids from the N-terminus of P-450 camC protein (415 amino acids, 46.66 kDa). In chimera 2 (525 amino acids, 60.24 kDa) the strategy involves replacement of 28 amino acids in the C-terminus of chimera 1 with 75 amino acids from the C-terminus of P-450 camC gene. Homology models of both the chimeric proteins were developed using SWISS-MODEL based on the known crystal structure of cytochrome P-450 camC, BM-3, 1DT6A, and 2C17A. The models indicated that the proposed heme-binding site was intact, which is inevitable for catalytic activity of cytochrome P-450s. The expression of chimera 1 and 2 genes in Escherichia coli DH5α was evident from light-pink cell pellets, protein band in sodium dodecyl sulfate polyacrylamide gel electrophoresis, and diagnostic carbon monoxide-difference spectra. Our studies show that strategies can be developed to exploit the natural diversity of the P-450 superfamily to generate chimeric biocatalysts that would provide new templates amenable to directed evolution.     

Inter- and intramolecular potential for the N-formylglycinamide-water system. A comparison between theoretical modeling and empirical force fields

An intramolecular NEMO potential is presented for the N-formylglycinamide molecule together with an intermolecular potential for the N-formylglycinamide-water system. The intramolecular N-formylglycinamide potential can be used as a building block for the backbone of polypeptides and proteins. Two intramolecular minima have been obtained. One, denoted as C5, is stabilized by a hydrogen bonded five member ring, and the other, denoted as C7, corresponds to a seven membered ring. The interaction between one water molecule and the N-formylglycinamide system is also studied and compared with Hartree-Fock SCF calculations and with the results obtained for some of the more commonly used force fields. The agreement between the NEMO and SCF energies for the complexes is in general superior to that of the other force fields. In the C7 region the surfaces obtained from the intramolecular part of the commonly used force fields are too flat compared to the NEMO potential and the ab initio calculations. We further analyze the possibility of using a charge distribution obtained from one conformation to describe the charge distribution of other conformations. We have found that the use of polarizabilities and generic dipoles can model most of the changes in charge density due to the different geometry of the new conformations, but that one can expect additional errors in the interaction energies that are of the order of 1 kcal/mol.     

Plasma Thermal Conversion of Methane to Acetylene

This paper describes a re-examination of a known process for the direct plasma thermal conversion of methane to acetylene. Conversion efficiencies (% methane converted) approached 100% and acetylene yields in the 90–95% range with 2–4% solid carbon production were demonstrated. Specificity for acetylene was higher than in prior work. Improvements in conversion efficiency, yield, and specificity were due primarily to improved injector design and reactant mixing, and minimization of temperature gradients and cold boundary layers. At the 60-kilowatt scale cooling by wall heat transfer appears to be sufficient to quench the product stream and prevent further reaction of acetylene resulting in the formation of heavier hydrocarbon products or solid carbon. Significantly increasing the quenching rate by aerodynamic expansion of the products through a converging–diverging nozzle led to a reduction in the yield of ethylene but had little effect on the yield of other hydrocarbon products. While greater product selectivity for acetylene has been demonstrated, the specific energy consumption per unit mass of acetylene produced was not improved upon. A kinetic model that includes the reaction mechanisms resulting in the formation of acetylene and heavier hydrocarbons, through benzene, is described.     

Computer simulation of conformational movement based on interconversion phenomena

A computer method has been developed which is an alternative to molecular dynamics in the sense that it pictures conformational motion. It simulates propagation in conformational movement based on conformational interconversion phenomena. The method starts with the knowledge of the topology of the conformational potential energy hyper-surface, which is described by the minima and the transition states. The topology is obtained by the recently described software . The simulation of conformational motion is based on the Boltzmann statistics of movement between the minima and the transition states at a given temperature. The method is illustrated for methyl--glucopyranoside and - -galNAc(1-3)[- -Fuc(1-2)]Gal-O-Me molecules. Conformational transitions of hydroxyl groups as well as glycosidic linkages are discussed.     

Performance of SuSi: a method for generating atomistic models of amorphous polymers based on a random search of energy minima

The performance of a recently developed method to generate representative atomistic models of amorphous polymers has been investigated. This method, which is denoted SuSi, can be defined as a random generator of energy minima. The effects produced by different parameters used to define the size of the system and the characteristics of the generation algorithm have been examined. Calculations have been performed on poly(L,D-lactic) acid (rho = 1.25 g/cm3) and nylon 6 (rho = 1.084 g/cm(3)), which are important commercial polymers.