A gradient descent algorithm for minimizing amino acid coupling reactions when synthesizing cyclic-peptide libraries

Combinatorial chemistry has become an invaluable tool in medicinal chemistry for the identification of new drug leads. For example, libraries of predetermined sequences and head-to-tail cyclized peptides are routinely synthesized in our laboratory using the IRORI approach. Such libraries are used as molecular toolkits that enable the development of pharmacophores that define activity and specificity at receptor targets. These libraries can be quite large and difficult to handle, due to physical and chemical constraints imposed by their size. Therefore, smaller sub-libraries are often targeted for synthesis. The number of coupling reactions required can be greatly reduced if the peptides having common amino acids are grouped into the same sub-library (batching). This paper describes a schedule optimizer to minimize the number of coupling reactions by rotating and aligning sequences while simultaneously batching. The gradient descent method thereby reduces the number of coupling reactions required for synthesizing cyclic peptide libraries. We show that the algorithm results in a 75% reduction in the number of coupling reactions for a typical cyclic peptide library.

Molecular Sieve Silica (MSS) Membranes for Gas Separation and Reaction Processes

Weakly branched silica films formed by the two-step sol-gel process allow for the formation of high selectivity membranes for gas separation. 29Si NMR and gas permeation showed that reduced crosslinking leads to He/CH4 selectivity improvement from 300 to 1000. Applied in membrane reactor for cyclohexane conversion to benzene, conversions were achieved at 14 fold higher than a conventional reactor at 250°C. Hydrothermal stability studies showed that carbon templating of silica is required for hydrothermally stable membranes. From our work it was shown that with correct application of chemistry, practical membrane systems can be built to suit gas separation (e. g. hydrogen fuel) and reactor systems.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

On the validity of the shrinking core model in noncatalytic gas solid reaction

By using a matched asymptotic expansion technique, the shrinking core model (SCM) used in non-catalytic gas solid reactions with general kinetic expression is rigorously justified in this paper as a special case of the homogeneous model when the reaction rate is much faster than that of diffusion. The time-pendent velocity of the moving reacted-unreacted interface is found to be proportional to the gas flux at that interface for all geometries of solid particles, and the thickness order of the reaction zone and also the degree of chemical reaction at the interface is discussed in this paper.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Dissociation between skin conductance orienting and secondary task reaction time: Time course with a visual discrimination task

A dissociation between two putative measures of resource allocation skin conductance responding, and secondary task reaction time (RT), has been observed during auditory discrimination tasks. Four experiments investigated the time course of the dissociation effect with a visual discrimination task. participants were presented with circles and ellipses and instructed to count the number of longer-than-usual presentations of one shape (task-relevant) and to ignore presentations of the other shape (task-irrelevant). Concurrent with this task, participants made a speeded motor response to an auditory probe. Experiment 1 showed that skin conductance responses were larger during task-relevant stimuli than during task-irrelevant stimuli, whereas RT to probes presented at 150 ms following shape onset was slower during task-irrelevant stimuli. Experiments 2 to 4 found slower RT during task-irrelevant stimuli at probes presented at 300 ms before shape onset until 150 ms following shape onset. At probes presented 3,000 and 4,000 ms following shape onset probe RT was slower during task-relevant stimuli. The similarities between the observed time course and the so-called psychological refractory period (PRF) effect are discussed.

The influence of joint position on the dynamics of perception-action coupling

Six right-handed subjects performed rhythmic flexion and extension movements of the index finger in time with an auditory metronome. On each block of trials, the wrist of the response hand was placed in a extended, neutral or flexed position. In the flex-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger flexion with each beat of the metronome. In the extend-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger extension with each beat of the metronome. The frequency of the metronome was increased from 2.00 Hz to 3.75 Hz in 8 steps (8 s epochs) of 0.25 Hz. During trials prepared in the extend-on-the-beat pattern, all subjects exhibited transitions to either a flex-on-the-beat pattern or to phase wandering as the frequency of pacing was increased. The time at which these transitions occurred was reliably influenced by the position of the wrist. Four subjects exhibited qualitative departures from the flex-on-the-beat pattern at pacing frequencies that were greater than those at which the extend-on-the-beat pattern could be maintained. The lime at which these departures occurred was not influenced by the position of the wrist. These results are discussed with reference to the constraints imposed on the coordination dynamics by the intrinsic properties of the neuromuscular-skeletal system.

Sensitive and specific detection of Clavibacter xyli subsp. xyli, causal agent of ratoon stunting disease of sugarcane, with a polymerase chain reaction-based assay

A sensitive, specific polymerase chain reaction-based assay was developed for the detection of the causal agent of ratoon stunting disease of sugarcane, Clavibacter xyli subsp. xyli. This assay uses oligonucleotide primers derived from the internal transcribed spacer region between the 16S and 23S rRNA genes of the bacterial rRNA operon. The assay is specific for C. xyli subsp. xyli and does not produce an amplification product from the template of the closely related bacterium C. xyli subsp. cynodontis, nor from other bacterial species. The assay was successfully applied to the detection of C. xyli subsp. xyli in fibrovascular fluid extracted from sugarcane and was sensitive to approximately 22 cells per PCR assay. A multiplex PCR test was also developed which identified and differentiated C. xyli subsp. xyli and C. xyli subsp. cynodontis in a single PCR assay.

Three-dimensional quantum solitons with parametric coupling

We consider the quantum field theory of two bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium, this corresponds to the process of second-harmonic generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. The quantum solitons or energy eigenstates (bound-state solutions) are obtained exactly in the simplest case of two-particle binding, in one, two, and three space dimensions. We also investigate three-particle binding in one space dimension. The results indicate that the exact quantum solitons of this field theory have a singular, pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. To estimate the physically accessible radii and binding energies of the bound states, we impose a momentum cutoff on the nonlinear couplings. In the case of nonlinear optical interactions, the resulting radii and binding energies of these photonic particlelike excitations in highly nonlinear parametric media appear to be close to physically observable values.

Coupling of a Jahn-Teller pseudorotation with a hindered internal rotation in an isolated molecule: 9-hydroxytriptycene

The irregular vibronic structure in the S-1<--S-0 resonant two-photon ionization (R2PI) spectrum of supersonically cooled triptycene is a result of a classic Exe Jahn-Teller effect [A. Furlan et al., J. Chem. Phys. 96, 7306 (1992)]. This is well characterized and can be used as an effective probe of intramolecular perturbations. Here we examine the S-1<--S-0 R2PI spectrum of 9-hydroxytriptycene and the fluorescence from various excited state vibronic levels. In this system the pseudorotation of the Jahn-Teller vibration is strongly coupled to the torsional motion of the bridgehead hydroxy group. This torsional motion results in a tunneling splitting in both the ground and excited states. The population of the upper level in the ground electronic state results in additional vibronic transitions becoming symmetry allowed in the R2PI spectrum that are forbidden in the bare triptycene molecule. The assignment of the R2PI and fluorescence spectra allows the potential energy surfaces of these vibrational modes to be accurately quantified. The full C-3v vibronic point group must be used to interpret the spectra. The time scale of the internal rotation of the-OH group and the butterfly flapping of the Jahn-Teller pseudorotation are of similar magnitude. The tunneling between the nine minima on the three dimensional potential energy surface is such that the Jahn-Teller pseudorotation occurs in concert with the-OH internal rotation. The Berry phase that is acquired during this motion is discussed. The simple physical picture emerges of the angle between two of the three benzene moieties opening in three equivalent ways in the S-1 electronic state. This geometry follows the position of the hydroxy group, which preferentially orients itself to point between these two rings. (C) 1998 American Institute of Physics. [S0021-9606(98)02348-4].

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Nickel complexes with tris(2-aminoethyl)amine (tren): [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2)center dot 6H(2)O (ox = oxalate), {[Ni-2(tren)(3)][ClO4](4)center dot H2O}(n), and [Ni-2(tren)(2)(aepd)][ClO4](4)center dot 2H(2)O (aepd = N-(2-aminoethyl)pyrrolidi

Reaction of K-3[Cr(ox)(3)] (ox = oxalate) with nickel(II) and tris(2-aminoethyl)amine (tren) in aqueous solution resulted in isolation of the bimetallic assembly [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2). 6H(2)O. The polymeric complex {[Ni-2(tren)(3)][ClO4](4). H2O}(n) has been prepared by reaction of nickel(II) perchlorate and tren in aqueous solution. From the same reaction mixture the complex [Ni-2(tren)(2)(aepd)][ClO4](4). 2H(2)O (aepd = N-(2-aminoethyl)pyrrolidine-3,4-diamine), in which a bridging tren ligand contains a carbon-carbon bond between two arms forming a substituted pyrrolidine, has been isolated. The complexes have been characterized by X-ray crystallography. The magnetic susceptibility (300-4.2 K) and magnetization data (2, 4 K, H = 0-5 T) for {[Ni-2(tren)(3)][ClO4](4). H2O}(n) (300 K , 4.23 mu(B)) exhibit evidence of weak antiferromagnetic coupling and zero field splitting (2J = -1.8 cm(-1); \ D\ = 2 cm(-1)) at low temperature. For [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2). 6H(2)O the susceptibility data at 300 K are indicative of uncoupled nickel(II) and chromium(III) sites with zero-field splitting and intramolecular antiferromagnetic coupling predicted at low temperature.

Synthesis of [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)] (PF6)(2), a hydroxo-bridged copper(II) complex of N,N '-dimethyl-1-thia-4,7-diazacyclononane (Me-2[9]aneN(2)S). X-ray structural analysis, magnetic susceptibility, EPR and visible spectroscopy

The bis(mu-hydroxo) complex [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)](PF6)(2) (Me-2[9]aneN(2)S = N,N'-dimethyl-1-thia-4,7-diazacyclononane) results after reaction of [Cu(Me-2[9]aneN(2)S)(MeCN)] (PF6) with dioxygen at -78 degrees C in acetonitrile. The complex has been characterized by X-ray crystallography: orthorhombic, space group Pnma, with a 18.710(3), b 16.758(2), c 9.593(2) Angstrom, and Z = 4. The structure refined to a final R value of 0.051. The complex contains two copper(II) ions bridged by two hydroxo groups with Cu ... Cu 2.866(1) Angstrom. The solid-state magnetic susceptibility study reveals ferromagnetic coupling, the fitting parameters being J = +46+/-5 cm(-1), g = 2.01+/-0.01 and theta = -0.58+/-0.03 K. The frozen-solution e.p.r. spectrum in dimethyl sulfoxide is characteristic of a monomeric copper(II) ion (g(parallel to) 2.300, g(perpendicular to) 2.063; A(parallel to) 156.2 x 10(-4) cm(-1), A(perpendicular to) 9.0 x 10(-4) cm(-1)) with an N2O2 donor set. Thioether coordination to the copper(II) in solution is supported by the presence of an intense absorption assigned to a sigma(S)-->Cu-II LMCT transition at c. 34000 cm(-1). The single-crystal spectrum of [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)] (PF6)(2) (273 K) reveals d-->d transitions at 14500 and 18300 cm(-1) and a weak pi(S)-->Cu-II charge-transfer band at approximately 25000 cm(-1).