The molecular design of a new solvent for the absorption of carbon dioxide

Gas absorption, the removal of one or more constitutents from a gas mixture, is widely used in chemical processes. In many gas absorption processes, the gas mixture is already at high pressure and in recent years organic solvents have been developed for the process of physical absorption at high pressure followed by low pressure regeneration of the solvent and recovery of the absorbed gases. Until now the discovery of new solvents has usually been by expensive and time consuming trial and error laboratory tests. This work describes a new approach, whereby a solvent is selected from considerations of its molecular structure by applying recently published methods of predicting gas solubility from the molecular groups which make up the solvent molecule. The removal of the acid gases of carbon dioxide and hydrogen sulfide from methane or hydrogen was used as a commercially important example. After a preliminary assessment to identify promising moecular groups, more than eighty new solvent molecules were designed and evaluated by predicting gas solubility. The other important physical properties were also predicted by appropriate theoretical procedures, and a commercially promising new solvent was chosen to have a high solubility for acid gases, a low solubility for methane and hydrogen, a low vapour pressure, and a low viscosity. The solvent chosen, of molecular structure Ch3-COCH2-CH2-CO-CH3, was tested in the laboratory and shown to have physical properties, except for vapour pressures, close to those predicted. That is gas solubilities were within 10% but lower than predicted. Viscosity within 10% but higher than predicted and a vapour pressure significantly lower than predicted. A computer program was written to predict gas solubility in the new solvent at the high pressures (25 bar) used in practice. This is based on the group contribution method of Skold Jorgensen (1984). Before using this with the new solvent, Acetonyl acetone, the method was show to be sufficiently accurate by comparing predicted values of gas solubility with experimental solubilities from the literature for 14 systems up to 50 bar. A test of the commercial potential of the new solvent was made by means of two design studies which compared the size of plant and approximate relative costs of absorbing acid gases by means of the new solvent with other commonly used solvents. These were refrigerated methanol(Rectisol process) and Dimethyl Ether or Polyethylene Glycol(Selexol process). Both studies showed in terms of capital and operating cost some significant advantage for plant designed for the new solvent process.

The separation of hydrogen and carbon monoxide using polymer membranes

This work studies the development of polymer membranes for the separation of hydrogen and carbon monoxide from a syngas produced by the partial oxidation of natural gas. The CO product is then used for the large scale manufacture of acetic acid by reaction with methanol. A method of economic evaluation has been developed for the process as a whole and a comparison is made between separation of the H2/CO mixture by a membrane system and the conventional method of cryogenic distillation. Costs are based on bids obtained from suppliers for several different specifications for the purity of the CO fed to the acetic acid reactor. When the purity of the CO is set at that obtained by cryogenic distillation it is shown that the membrane separator offers only a marginal cost advantage. Cost parameters for the membrane separation systems have been defined in terms of effective selectivity and cost permeability. These new parameters, obtained from an analysis of the bids, are then used in a procedure which defines the optimum degree of separation and recovery of carbon monoxide for a minimum cost of manufacture of acetic acid. It is shown that a significant cost reduction is achieved with a membrane separator at the optimum process conditions. A method of "targeting" the properties of new membranes has been developed. This involves defining the properties for new (hypothetical -yet to be developed) membranes such that their use for the hydrogen/carbon monoxide separation will produce a reduced cost of acetic acid manufacture. The use of the targeting method is illustrated in the development of new membranes for the separation of hydrogen and carbon monoxide. The selection of polymeric materials for new membranes is based on molecular design methods which predict the polymer properties from the molecular groups making up the polymer molecule. Two approaches have been used. One method develops the analogy between gas solubility in liquids and that in polymers. The UNIFAC group contribution method is then used to predict gas solubility in liquids. In the second method the polymer Permachor number, developed by Salame, has been correlated with hydrogen and carbon monoxide permeabilities. These correlations are used to predict the permeabilities of gases through polymers. Materials have been tested for hydrogen and carbon monoxide permeabilities and improvements in expected economic performance have been achieved.

Molecular modelling assisted design and synthesis of cyclothialidine derivatives as DNA gyrase inhibitors

Since cyclothialidine was discovered as the most active DNA gyrase inhibitor in 1994, enormous efforts have been devoted to make it into a commercial medicine by a number of pharmaceutical companies and research groups worldwide. However, no serious breakthrough has been made up to now. An essential problem involved with cyclothialidine is that though it demonstrated the potent inhibition of DNA gyrase, it showed little activity against bacteria. This probably is attributable to its inability to penetrate bacterial cell walls and membranes. We applied the TSAR programme to generate a QSAR equation to the gram-negative organisms. In that equation, LogP is profoundly indicated as the key factor influencing the cyclothialidine activity against bacteria. However, the synthesized new analogues have failed to prove that. In the structure based drug design stage, we designed a group of open chain cyclothialidine derivatives by applying the SPROUT programme and completed the syntheses. Improved activity is found in a few analogues and a 3D pharmacophore of the DNA gyrase B is proposed to lead to synthesis of the new derivatives for development of potent antibiotics.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

Design and synthesis of potential antimicrobial agents

Chorismate mutase is one of the essential enzymes in the shikimate pathway and is key to the survival of the organism Mycobacterium tuberculosis. The x-ray crystal structure of this enzyme from Mycobacterium tuberculosis was manipulated to prepare an initial set of in silico protein models of the active site. Known inhibitors of the enzyme were docked into the active site using the flexible ligand / flexible active site side chains approach implemented in CAChe Worksystem (Fujitsu Ltd). The resulting complexes were refined by molecular dynamics studies in explicit water using Amber 9. This yielded a further set of protein models that were used for additional rounds of ligand docking. A binding hypothesis was established for the enzyme and this was used to screen a database of commercially available drug-like compounds. From these results new potential ligands were designed that fitted appropriately into the active site and matched the functional groups and binding motifs founds therein. Some of these compounds and close analogues were then synthesized and submitted for biological evaluation. As a separate part of this thesis, analogues of very active anti-tuberculosis pyridylcarboxamidrazone were also prepared. This was carried out by the addition and the deletion of the substitutions from the lead compound thereby preparing heteroaryl carboxamidrazone derivatives and related compounds. All these compounds were initially evaluated for biological activity against various gram positive organisms and then sent to the TAACF (USA) for screening against Mycobacterium tuberculosis. Some of the new compounds proved to be at least as potent as the original lead compound but less toxic.

Structure and dynamics of multiple cationic vectors-siRNA complexation by all-atomic molecular dynamics simulations

Understanding the molecular mechanism of gene condensation is a key component to rationalizing gene delivery phenomena, including functional properties such as the stability of the gene-vector complex and the intracellular release of the gene. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with four cationic carrier systems of varying charge and surface topology at different charge ratios. At lower charge ratios, polymers bind quite effectively to siRNA, while at high charge ratios, the complexes are saturated and there are free polymers that are unable to associate with RNA. We also observed reduced fluctuations in RNA structures when complexed with multiple polymers in solution as compared to both free siRNA in water and the single polymer complexes. These novel simulations provide a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance progress toward rational design of nonviral gene delivery systems.

The effect of pH on PAMAM dendrimer-siRNA complexation:endosomal considerations as determined by molecular dynamics simulation

Intracellular degradation of genes, most notably within the endo-lysosomal compartment is considered a significant barrier to (non-viral) gene delivery in vivo. Previous reports based on in vitro studies claim that carriers possessing a mixture of primary, secondary and tertiary amines are able to buffer the acidic environment within the endosome, allowing for timely release of their contents, leading to higher transfection rates. In this report, we adopt an atomistic molecular dynamics (MD) simulation approach, comparing the complexation of 21-bp siRNA with low-generation polyamidoamine (PAMAM) dendrimers (G0 and G1) at both neutral and acidic pHs, the latter of which mimics the degradative environment within maturing 'late-endosomes'. Our simulations reveal that the time taken for the dendrimer-gene complex (dendriplex) to reach equilibrium is appreciably longer at low pH and this is accompanied by more compact packaging of the dendriplex, as compared to simulations performed at neutral pH. We also note larger absolute values of calculated binding free energies of the dendriplex at low pH, indicating a higher dendrimer-nucleic acid affinity in comparison with neutral pH. These novel simulations provide a more detailed understanding of low molecular-weight polymer-siRNA behavior, mimicking the endosomal environment and provide input of direct relevance to the "proton sponge theory", thereby advancing the rational design of non-viral gene delivery systems.

Design and construction of an ultrasonic probe for use in a Cryo-Magnet NMR Spectrometer

SINNMR (Sonically Induced Narrowing of the Nuclear Magnetic Resonance spectra of solids), is a novel technique that is being developed to enable the routine study of solids by nuclear magnetic resonance spectroscopy. SINNMR aims to narrow the broad resonances that are characteristic of solid state NMR by inducing rapid incoherent motion of solid particles suspended in a support medium, using high frequency ultrasound in the range 2-10 MHz. The width of the normal broad resonances from solids are due to incomplete averaging of several components of the total spin Hamiltonian caused by restrictions placed on molecular motion within a solid. At present Magic Angle Spinning (MAS) NMR is the classical solid state technique used to reduce line broadening, but: this has associated problems, not least of which is the appearance of many spinning side bands which confuse the spectra. It is hoped that SlNNMR will offer a simple alternative, particularly as it does not reveal spinning sidebands The fundamental question concerning whether the use of ultrasound within a cryo-magnet will cause quenching has been investigated with success, as even under the most extreme conditions of power, frequency and irradiator time, the magnet does not quench. The objective of this work is to design and construct a SINNMR probe for use in a super conducting cryo-magnet NMR spectrometer. A cell for such a probe has been constructed and incorporated into an adapted high resolution broadband probe. It has been proved that the cell is capable of causing cavitation, up to 10 MHz, by running a series of ultrasonic reactions within it and observing the reaction products. It was found that the ultrasound was causing the sample to be heated to unacceptable temperatures and this necessitated the incorporation of temperature stabilisation devices. Work has been performed on the investigation of the narrowing of the solid state 23Na spectrum of tri-sodium phosphate using high frequency ultrasound. Work has also been completed on the signal enhancement and T1 reduction of a liquid mixture and a pure compound using ultrasound. Some preliminary "bench" experiments have been completed on a novel ultrasonic device designed to help minimise sample heating. The concept involves passing the ultrasound through a temperature stabilised, liquid filled funnel that has a drum skin on the end that will enable the passage of ultrasound into the sample. Bench experiments have proved that acoustic attenuation is low and that cavitation in the liquid beyond the device is still possible.

Design and characterization of zeolite membranes for fluid separation

Experimental and theoretical methods have been used to study zeolite structures, properties and applications as membranes for separation purposes. Thin layers of silicalite-1 and Na-LTA zeolites have been synthesised onto carbon-graphite supports using a hydrothermal synthesis procedure. The separation behaviour of the composite membranes was characterized by gas permeation studies of pure, binary and ternary mixtures of methane, ethane and propane. The influence of temperature and feed gas mixture composition on the separation and selectivity performance of the membranes was also investigated. It was found that the silicalite-1 composite membranes synthesised onto the 4 hour oxidized carbon-graphite supports showed the most promising separation behaviour of all the composite membranes investigated. Molecular simulation methods were used to gain an understanding of how hydrocarbon molecules behave both within the pores and on the surfaces of silicalite-1, mordenite and LTA zeolites. Molecular dynamic simulations were used to investigate the influence of temperature and molecular loadings on the diffusional behaviour of hydrocarbons in zeolites. Both hydroxylated (surface termination with hydroxyl groups) and non-hydroxylated silicalite-1 and Na-mordenite surfaces were generated. For both zeolites the most stable surfaces correspond to the {010} surface. For the silicalite-1 {010} surface the adsorption of hydrocarbons and molecular water onto the hydroxylated surface showed a favourable exothermic adsorption process compared to adsorption on the non-hydroxylated surface. With the Na-mordenite {010} surface the adsorption of hydrocarbons onto both the hydroxylated and non-hydroxylated surfaces had a combination of favourable and non-favourable adsorption energies, while the adsorption of molecular water onto both types of surface was found to be a favourable adsorption process.

Liquid-phase adsorption of n-paraffins on type 5A molecular sieves

As a basis for the commercial separation of normal paraffins a detailed study has been made of factors affecting the adsorption of binary liquid mixtures of high molecular weight normal paraffins (C12, C16, and C20) from isooctane on type 5A molecular sieves. The literature relating to molecular sieve properties and applications, and to liquid-phase adsorption of high molecular weight normal paraffin compounds by zeolites, was reviewed. Equilibrium isotherms were determined experimentally for the normal paraffins under investigation at temperatures of 303oK, 323oK and 343oK and showed a non-linear, favourable- type of isotherm. A higher equilibrium amount was adsorbed with lower molecular weight normal paraffins. An increase in adsorption temperature resulted in a decrease in the adsorption value. Kinetics of adsorption were investigated for the three normal paraffins at different temperatures. The effective diffusivity and the rate of adsorption of each normal paraffin increased with an increase in temperature in the range 303 to 343oK. The value of activation energy was between 2 and 4 kcal/mole. The dynamic properties of the three systems were investigated over a range of operating conditions (i.e. temperature, flow rate, feed concentration, and molecular sieve size in the range 0.032 x 10-3 to 2 x 10-3m) with a packed column. The heights of adsorption zones calculated by two independent equations (one based on a constant width, constant velocity and adsorption zone and the second on a solute material balance within the adsorption zone) agreed within 3% which confirmed the validity of using the mass transfer zone concept to provide a simple design procedure for the systems under study. The dynamic capacity of type 5A sieves for n-eicosane was lower than for n-hexadecane and n-dodecane corresponding to a lower equilibrium loading capacity and lower overall mass transfer coefficient. The values of individual external, internal, theoretical and experimental overall mass transfer coefficient were determined. The internal resistance was in all cases rate-controlling. A mathematical model for the prediction of dynamic breakthrough curves was developed analytically and solved from the equilibrium isotherm and the mass transfer rate equation. The experimental breakthrough curves were tested against both the proposed model and a graphical method developed by Treybal. The model produced the best fit with mean relative percent deviations of 26, 22, and 13% for the n-dodecane, n-hexadecane, and n-eicosane systems respectively.

The design of novel antifolates against pneumocystis carinii

The pneumonia caused by Pneumocystis carinii is ultimately responsible for the death of many acquired immunodeficiency syndrome (AIDS) patients. Large doses of trimethoprim and pyrimethamine in combination with a sulphonamide and/or pentamidine suppress the infection but produce serious side-effects and seldom prevent recurrence after treatment withdrawal. However, the partial success of the aforementioned antifolates, and also trimetrexate used alone, does suggest dihydrofolate reductase (DHFR) as a target for the development of antipneumocystis agents. From the DHFR inhibitory activities of 3'-substituted pyrimethamine analogues it was suggested that the 3'-(3'',3''-dimethyltriazen-1''-yl) substituent may be responsible for the greater activity for the P.carinii over the mammalian enzyme. Crystallographic and molecular modeling studies revealed considerable geometrical and electronic differences between the triazene and the chemically related formamidine functions that may account for the differences in DHFR inhibitory profiles. Structural and electronic parameters calculated for a series of 3'-(3'',3''-disubstitutedtriazen-1''-yl) pyrimethamine analogues did not correlate with the DHFR inhibitory activities. However, the in vitro screening against P.carinii DHFR revealed that the 3''-hydroxyethyl-3''-benzyl analogue was the most active and selective. Models of the active sites of human and P.carinii DHFRs were constructed using DHFR sequence and structural homology data which had identified key residues involved in substrate and cofactor binding. Low energy conformations of the 3'',3''-dimethyl and 3''-hydroxyethyl-3''-benzyle analogues, determined from nuclear magnetic resonance studies and theoretical calculations, were docked by superimposing the diaminopyrimidine fragment onto a previously docked pyrimethamine analogue. Enzyme kinetic data supported the 3''-hydroxyethyl-3''-benzyl moiety being located in the NADPH binding groove. The 3''-benzyl substituent was able to locate to within 3 AA of a valine residue in the active site of P.carinii DHFR thereby producing a hydrophobic contact. The equivalent residue in human DHFR is threonine, more hydrophilic and less likely to be involved in such a contact. This difference may account for the greater inhibitory activity this analogue has for P.carinii DHFR and provide a basis for future drug design. From an in vivo model of PCP in immunosuppressed rats it was established that the 3"-hydroxyethyl-3"-benzyl analogue was able to reduce the.P.carinii burden more effectively with increasing doses, without causmg any visible signs of toxicity. However, equivalent doses were not as effective as pentamidine, a current treatment of choice for Pneumocystis carinii pneumonia.

The implementation of a design of experiments strategy to increase recombinant protein yields in yeast (review)

Biological processes are subject to the influence of numerous factors and their interactions, which may be non-linear in nature. In a recombinant protein production experiment, understanding the relative importance of these factors, and their influence on the yield and quality of the recombinant protein being produced, is an essential part of its optimisation. In many cases, implementing a design of experiments (DoE) approach has delivered this understanding. This chapter aims to provide the reader with useful pointers in applying a DoE strategy to improve the yields of recombinant yeast cultures.

Proteins accessible to immune surveillance show significant T-cell epitope depletion:implications for vaccine design

T cell activation is the final step in a complex pathway through which pathogen-derived peptide fragments can elicit an immune response. For it to occur, peptides must form stable complexes with Major Histocompatibility Complex (MHC) molecules and be presented on the cell surface. Computational predictors of MHC binding are often used within in silico vaccine design pathways. We have previously shown that, paradoxically, most bacterial proteins known experimentally to elicit an immune response in disease models are depleted in peptides predicted to bind to human MHC alleles. The results presented here, derived using software proven through benchmarking to be the most accurate currently available, show that vaccine antigens contain fewer predicted MHC-binding peptides than control bacterial proteins from almost all subcellular locations with the exception of cell wall and some cytoplasmic proteins. This effect is too large to be explained from the undoubted lack of precision of the software or from the amino acid composition of the antigens. Instead, we propose that pathogens have evolved under the influence of the host immune system so that surface proteins are depleted in potential MHC-binding peptides, and suggest that identification of a protein likely to contain a single immuno-dominant epitope is likely to be a productive strategy for vaccine design.

Molecular and genetic evidence for abnormalities in the nodes of Ranvier in schizophrenia

Context: Genetic, neuroimaging, and molecular neurobiological evidence support the hypothesis that the disconnectivity syndrome in schizophrenia (SZ) could arise from failures of saltatory conduction and abnormalities at the nodes of Ranvier (NOR) interface where myelin and axons interact. Objective: To identify abnormalities in the expression of oligodendroglial genes and proteins that participate in the formation, maintenance, and integrity of the NOR in SZ. Design: The messenger RNA (mRNA) expression levels of multiple NOR genes were quantified in 2 independent postmortem brain cohorts of individuals with SZ, and generalizability to protein expression was confirmed. The effect of the ANK3 genotype on the mRNA expression level was tested in postmortem human brain. Case-control analysis tested the association of the ANK3 genotype with SZ. The ANK3 genotype's influence on cognitive task performance and functional magnetic resonance imaging activation was tested in 2 independent cohorts of healthy individuals. Setting: Research hospital. Patients: Postmortem samples from patients with SZ and healthy controls were used for the brain expression study (n=46) and the case-control analysis (n=272). Healthy white men and women participated in the cognitive (n=513) and neuroimaging (n=52) studies. Main Outcome Measures: The mRNA and protein levels in postmortem brain samples, genetic association with schizophrenia, cognitive performance, and blood oxygenation level-dependent functional magnetic resonance imaging. Results: The mRNA expression of multiple NOR genes was decreased in schizophrenia. The ANK3 rs9804190 C allele was associated with lower ANK3 mRNA expression levels, higher risk for SZ in the case-control cohort, and poorer working memory and executive function performance and increased prefrontal activation during a working memory task in healthy individuals. Conclusions: These results point to abnormalities in the expression of genes and protein associated with the integrity of the NOR and suggest them as substrates for the disconnectivity syndrome in SZ. The association of ANK3 with lower brain mRNA expression levels implicates a molecular mechanism for its genetic, clinical, and cognitive associations with SZ. ©2012 American Medical Association. All rights reserved.

Molecular dynamics simulations:bring biomolecular structures alive on a computer

The molecular dynamics (MD) simulations play a very important role in science today. They have been used successfully in binding free-energy calculations and rational design of drugs and vaccines. MD simulations can help visualize and understand structures and dynamics at an atomistic level when combined with molecular graphics programs. The molecular and atomistic properties can be displayed on a computer in a time-dependent way, which opens a road toward a better understanding of the relationship of structure, dynamics, and function. In this chapter, the basics of MD are explained, together with a step-by-step description of setup and running an MD simulation.