Models of gas-liquid solubilities

A recent method for phase equilibria, the AGAPE method, has been used to predict activity coefficients and excess Gibbs energy for binary mixtures with good accuracy. The theory, based on a generalised London potential (GLP), accounts for intermolecular attractive forces. Unlike existing prediction methods, for example UNIFAC, the AGAPE method uses only information derived from accessible experimental data and molecular information for pure components. Presently, the AGAPE method has some limitations, namely that the mixtures must consist of small, non-polar compounds with no hydrogen bonding, at low moderate pressures and at conditions below the critical conditions of the components. Distinction between vapour-liquid equilibria and gas-liquid solubility is rather arbitrary and it seems reasonable to extend these ideas to solubility. The AGAPE model uses a molecular lattice-based mixing rule. By judicious use of computer programs a methodology was created to examine a body of experimental gas-liquid solubility data for gases such as carbon dioxide, propane, n-butane or sulphur hexafluoride which all have critical temperatures a little above 298 K dissolved in benzene, cyclo-hexane and methanol. Within this methodology the value of the GLP as an ab initio combining rule for such solutes in very dilute solutions in a variety of liquids has been tested. Using the GLP as a mixing rule involves the computation of rotationally averaged interactions between the constituent atoms, and new calculations have had to be made to discover the magnitude of the unlike pair interactions. These numbers have been seen as significant in their own right in the context of the behaviour of infinitely-dilute solutions. A method for extending this treatment to "permanent" gases has also been developed. The findings from the GLP method and from the more general AGAPE approach have been examined in the context of other models for gas-liquid solubility, both "classical" and contemporary, in particular those derived from equations-of-state methods and from reference solvent methods.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Hydrotalcite intercalated siRNA:computational characterization of the interlayer environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

Quinols as novel therapeutic agents. 2.1 4-(1-arylsulfonylindol- 2-yl)-4-hydroxycyclohexa-2,5-dien-1-ones and related agents as potent and selective antitumor agents

A series of substituted 4-(1-arylsulfonylindol-2-yl)-4-hydroxycyclohexa-2, 5-dien-1-ones (indolylquinols) has been synthesized on the basis of the discovery of lead compound 1a and screened for antitumor activity. Synthesis of this novel series was accomplished via the "one-pot" addition of lithiated (arylsulfonyl)indoles to 4,4-dimethoxycyclohexa-2,5-dienone followed by deprotection under acidic conditions. Similar methodology gave rise to the related naphtho-, 1H-indole-, and benzimidazole-substituted quinols. A number of compounds in this new series were found to possess in vitro human tumor cell line activity substantially more potent than the recently reported antitumor 4-substituted 4-hydroxycyclohexa-2,5-dien-1-ones1 with similar patterns of selectivity against colon, renal, and breast cell lines. The most potent compound in the series in vitro, 4-(1-benzenesulfonyl-6-fluoro-1H-indol- 2-yl)-4-hydroxycyclohexa-2,5-dienone (1h), exhibits a mean GI50 value of 16 nM and a mean LC50 value of 2.24 μM in the NCI 60-cell-line screen, with LC50 activity in the HCT 116 human colon cancer cell line below 10 nM. The crystal structure of the unsubstituted indolylquinol 1a exhibits two independent molecules, both participating in intermolecular hydrogen bonds from quinol OH to carbonyl O, but one OH group also interacts intramolecularly with a sulfonyl O atom. This interaction, which strengthens upon ab initio optimization, may influence the chemical environment of the bioactive quinol moiety. In vivo, significant antitumor activity was recorded (day 28) in mice bearing subcutaneously implanted MDA-MB-435 xenografts, following intraperitoneal treatment of mice with compound 1a at 50 mg/kg.

X-ray crystallographic and theoretical studies of an anticonvulsant enaminone:methyl 4-(4'-bromophenyl)amino-6-methyl-2-oxocyclohex-3-en-1-oate

Objective: The aims of this study were to establish the structure of the potent anticonvulsant enaminone methyl 4-(4′-bromophenyl)amino-6-methyl-2- oxocyclohex-3-en-1-oate (E139), and to determine the energetically preferred conformation of the molecule, which is responsible for the biological activity. Materials and Methods: The structure of the molecule was determined by X-ray crystallography. Theoretical ab initio calculations with different basis sets were used to compare the energies of the different enantiomers and to other structurally related compounds. Results: The X-ray crystal structure revealed two independent molecules of E139, both with absolute configuration C11(S), C12(R), and their inverse. Ab initio calculations with the 6-31G, 3-21G and STO-3G basis sets confirmed that the C11(S), C12(R) enantiomer with both substituents equatorial had the lowest energy. Compared to relevant crystal structures, the geometry of the theoretical structures shows a longer C-N and shorter C=O distance with more cyclohexene ring puckering in the isolated molecule. Conclusion: Based on a pharmacophoric model it is suggested that the enaminone system HN-C=C-C=O and the 4-bromophenyl group in E139 are necessary to confer anticonvulsant property that could lead to the design of new and improved anticonvulsant agents. Copyright © 2003 S. Karger AG, Basel.

Simple one-electron invariants of molecular chirality

Pseudoscalar measures of electronic chirality for molecular systems are derived using the spectral moment theory applied to the frequency-dependent rotational susceptibility. In this scheme a one-electron chirality operator κ^ naturally emerges as a quantum counterpart of the triple scalar product, involving velocity, acceleration and second acceleration. Averaging κ^ over an electronic state vector gives rise to an additive chirality invariant (κ-index), considered as a quantitative measure of chirality. A simple computational technique for quick calculation of the κ-index is developed and various structural classes (cyclic hydrocarbons, cage-shaped systems, etc.) are studied. Reasonable behaviour of the chirality index is demonstrated. The chirality changes during the β-turn formation in Leu-Enkephalin is presented as a useful example of the chirality analysis for conformational transitions.