Marriage Matching with Correlated Preferences

Authors of experimental, empirical, theoretical and computational studies of two-sided matching markets have recognized the importance of correlated preferences. We develop a general method for the study of the effect of correlation of preferences on the outcomes generated by two-sided matching mechanisms. We then illustrate our method by using it to quantify the effect of correlation of preferences on satisfaction with the men-propose Gale-Shapley matching for a simple one-to-one matching problem.

NMR determination of the major solution conformation of a peptoid pentamer with chiral side chains

Polymers of N-substituted glycines (“peptoids”) containing chiral centers at the α position of their side chains can form stable structures in solution. We studied a prototypical peptoid, consisting of five para-substituted (S)-N-(1-phenylethyl)glycine residues, by NMR spectroscopy. Multiple configurational isomers were observed, but because of extensive signal overlap, only the major isomer containing all cis-amide bonds was examined in detail. The NMR data for this molecule, in conjunction with previous CD spectroscopic results, indicate that the major species in methanol is a right-handed helix with cis-amide bonds. The periodicity of the helix is three residues per turn, with a pitch of ≈6 Å. This conformation is similar to that anticipated by computational studies of a chiral peptoid octamer. The helical repeat orients the amide bond chromophores in a manner consistent with the intensity of the CD signal exhibited by this molecule. Many other chiral polypeptoids have similar CD spectra, suggesting that a whole family of peptoids containing chiral side chains is capable of adopting this secondary structure motif. Taken together, our experimental and theoretical studies of the structural properties of chiral peptoids lay the groundwork for the rational design of more complex polypeptoid molecules, with a variety of applications, ranging from nanostructures to nonviral gene delivery systems.

Effects of unique ion chemistry on thin-film growth by plasma–surface interactions

Plasma processing is a standard industrial method for the modification of material surfaces and the deposition of thin films. Polyatomic ions and neutrals larger than a triatomic play a critical role in plasma-induced surface chemistry, especially in the deposition of polymeric films from fluorocarbon plasmas. In this paper, low energy CF3+ and C3F5+ ions are used to modify a polystyrene surface. Experimental and computational studies are combined to quantify the effect of the unique chemistry and structure of the incident ions on the result of ion-polymer collisions. C3F5+ ions are more effective at growing films than CF3+, both at similar energy/atom of ≈6 eV/atom and similar total kinetic energies of 25 and 50 eV. The composition of the films grown experimentally also varies with both the structure and kinetic energy of the incident ion. Both C3F5+ and CF3+ should be thought of as covalently bound polyatomic precursors or fragments that can react and become incorporated within the polystyrene surface, rather than merely donating F atoms. The size and structure of the ions affect polymer film formation via differing chemical structure, reactivity, sticking probabilities, and energy transfer to the surface. The different reactivity of these two ions with the polymer surface supports the argument that larger species contribute to the deposition of polymeric films from fluorocarbon plasmas. These results indicate that complete understanding and accurate computer modeling of plasma–surface modification requires accurate measurement of the identities, number densities, and kinetic energies of higher mass ions and energetic neutrals.

Surface geochemistry of the clay minerals

Clay minerals are layer type aluminosilicates that figure in terrestrial biogeochemical cycles, in the buffering capacity of the oceans, and in the containment of toxic waste materials. They are also used as lubricants in petroleum extraction and as industrial catalysts for the synthesis of many organic compounds. These applications derive fundamentally from the colloidal size and permanent structural charge of clay mineral particles, which endow them with significant surface reactivity. Unraveling the surface geochemistry of hydrated clay minerals is an abiding, if difficult, topic in earth sciences research. Recent experimental and computational studies that take advantage of new methodologies and basic insights derived from the study of concentrated ionic solutions have begun to clarify the structure of electrical double layers formed on hydrated clay mineral surfaces, particularly those in the interlayer region of swelling 2:1 layer type clay minerals. One emerging trend is that the coordination of interlayer cations with water molecules and clay mineral surface oxygens is governed largely by cation size and charge, similarly to a concentrated ionic solution, but the location of structural charge within a clay layer and the existence of hydrophobic patches on its surface provide important modulations. The larger the interlayer cation, the greater the influence of clay mineral structure and hydrophobicity on the configurations of adsorbed water molecules. This picture extends readily to hydrophobic molecules adsorbed within an interlayer region, with important implications for clay–hydrocarbon interactions and the design of catalysts for organic synthesis.

Curtin-Hammett versus non-Curtin-Hammett frameworks is optimizing the enantioselective binolam/titanium(IV)-catalyzed cyanobenzoylation of aldehydes: Part 2

Extensive experimental and computational studies have been carried out on the enantioselective titanium(IV)-catalyzed cyanobenzoylation of aldehydes using 1:n Binolam:Ti(OiPr)4 mixtures as precatalysts, with the purpose of identifying the key mechanistic aspects governing enantioselectivity. HCN and isopropyl benzoate were detected in the reacting mixtures. This, as well as the reaction’s response to the presence of an exogenous base, and the failure to react in the presence of Binol:Ti(OiPr)4 mixtures, led us to propose not a direct cyanobenzoylation but an indirect process involving enantioselective hydrocyanation followed by O-benzoylation. Computational work provided positive evidence for the intervention of both indirect and direct cyanobenzoylation routes, the former being faster. However, the standard Curtin–Hammett-based optimization search ended with unsatisfactory results. Experimental and computational DFT studies (B3LYP/6-31G*) led us to conclude that: (1) the overall cyanobenzoylation of aldehydes catalyzed by 1:n Binolam:Ti(OiPr)4 mixtures involves an enantioselective hydrocyanation followed by an stereochemically inert O-benzoylation; (2) the initial complexes prevailing in a 1:1 Binolam:Ti(OiPr)4 mixture are the solvated mononuclear monomer 5·2(iPrOH) and solvated dinuclear dimer 9·2(iPrOH), whereas 9·2(iPrOH) is the major component in a 1:2 or higher 1:n mixture; (3) since the slowest step is that of benzoylation of ligated iPrOH which yields the actual catalysts 5–9, the catalytic system fits into a non-Curtin–Hammett framework, the final products deriving from a kinetic quench of the competing routes; and (4) accordingly, catalysis by 1:1 Binolam:Ti(OiPr)4 mixtures should involve cyanobenzoylations promoted by mononuclear 5, contaminated with those promoted by some dinuclear open dimer 9, whereas cyanobenzoylations catalyzed by a 1:2 and higher 1:n mixtures should be the result of catalysis promoted by the large amounts of dinuclear open dimer 9.

Binap-gold (I) versus Binap-silver trifluoroacetate complexes as catalysts in 1,3-dipolar cycloadditions of axomethine ylides

The 1,3-dipolar cycloaddition between azomethine ylides and alkenes is efficiently catalysed by [{(Sa)-Binap-Au(tfa)}2] (Binap=2,2′-bis(diphenylphosphino)-1,1′-binaphthyl; tfa=trifluoroacetyl). Maleimides, 1,2-bis(phenylsulfonyl)ethylene, chalcone and nitrostyrene were suitable dipolarophiles even when using sterically hindered 1,3-dipole precursors. The results obtained in these transformations improve the analogous ones obtained in the same reactions catalysed by [Binap–Ag(tfa)]. In addition, computational studies have also been carried out to demonstrate both the high enantioselectivity exhibited by the chiral gold(I) complex, and the non-linear effect observed in this transformation.

Regio and diastereoselective multicomponent 1,3-dipolar cycloadditions between prolinate hydrochlorides, aldehydes and dipolarophiles for the direct synthesis of pyrrolizidines

A general synthesis of highly substituted pyrrolizidines can be performed by a multicomponent 1,3-dipolar cycloaddition using proline ester hydrochlorides, aldehydes and dipolarophiles, at room temperature without catalysts or in the presence of AgOAc (5 mol %). In the case of (2S,4R)-4-hydroxyproline derivatives it is possible to obtain enantioenriched pyrrolizidines with high control of the regio- and diastereoselectivity affording the adducts 2,4-trans-2,5-trans according to an endo-approach and a S-dipole geometry of the in situ generated azomethine ylide. For proline esters a similar regioselectivity and endo-diastereoselectivity are observed when the dipole promotes an α-attack. However, when ethyl glyoxylate is used as aldehyde component the γ-attack of the S-ylide takes place preferentially giving rise the opposite regioselectivity for acrylic dipolarophiles, being crucial the role of silver acetate. In this case, the exo-adducts with a 2,3-cis-2,5-trans relative configuration are diastereoselectively obtained. In addition, computational studies have also been carried out to shed light on the origins of the diastereo- and regioselectivity observed for the described 1,3-dipolar cycloadditions.

Thermodynamics of binding of regulatory ligands to tissue transglutaminase

The transamidating activity of tissue transglutaminase is regulated by the ligands calcium and GTP, via conformational changes which facilitate or interfere with interaction with the peptidyl-glutamine substrate. We have analysed binding of these ligands by calorimetric and computational approaches. In the case of GTP we have detected a single high affinity site (K (D) approximately 1 muM), with moderate thermal effects suggestive that binding GTP involves replacement of GDP, normally bound to the protein. On line with this possibility no significant binding was observed during titration with GDP and computational studies support this view. Titration with calcium at a high cation molar excess yielded a complex binding isotherm with a number of "apparent binding sites" in large excess over those detectable by equilibrium dialysis (6 sites). This binding pattern is ascribed to occurrence of additional thermal contributions, beyond those of binding, due to the occurrence of conformational changes and to catalysis itself (with protein self-crosslinking). In contrast only one site for binding calcium with high affinity (K (D) approximately 0.15 muM) is observed with samples of enzyme inactivated by alkylation at the active site (to prevent enzyme crosslinkage and thermal effects of catalysis). These results indicate an intrinsic ability of tissue transglutaminase to bind calcium with high affinity and the necessity of careful reassessment of the enzyme regulatory pattern in relation to the concentrations of ligands in living cells, taking also in account effects of ligands on protein subcellular compartimentation.

The dawn of the RNA World:toward functional complexity through ligation of random RNA oligomers

A main unsolved problem in the RNA World scenario for the origin of life is how a template-dependent RNA polymerase ribozyme emerged from short RNA oligomers obtained by random polymerization on mineral surfaces. A number of computational studies have shown that the structural repertoire yielded by that process is dominated by topologically simple structures, notably hairpin-like ones. A fraction of these could display RNA ligase activity and catalyze the assembly of larger, eventually functional RNA molecules retaining their previous modular structure: molecular complexity increases but template replication is absent. This allows us to build up a stepwise model of ligation- based, modular evolution that could pave the way to the emergence of a ribozyme with RNA replicase activity, step at which information-driven Darwinian evolution would be triggered. Copyright © 2009 RNA Society.

On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10 molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies-characterized by a fixed number of paired nucleotides-and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited. © 2008 Elsevier Ltd. All rights reserved.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

Towards understanding promiscuity in multidrug efflux pumps

Drug export from cells is a major factor in the acquisition of cellular resistance to antimicrobial and cancer chemotherapy, and poses a significant threat to future clinical management of disease. Many of the proteins that catalyse drug efflux do so with remarkably low substrate specificity, a phenomenon known as multidrug transport. For these reasons we need a greater understanding of drug recognition and transport in multidrug pumps to inform research that attempts to circumvent their action. Structural and computational studies have been heralded as being great strides towards a full elucidation of multidrug recognition and transport. In this review we summarise these advances and ask how close we are to a molecular understanding of this remarkable phenomenon. © 2013 Elsevier Ltd.

Realization of an Optimal Schedule for the Two-machine Flow-Shop with Interval Job Processing Times

Non-preemptive two-machine flow-shop scheduling problem with uncertain processing times of n jobs is studied. In an uncertain version of a scheduling problem, there may not exist a unique schedule that remains optimal for all possible realizations of the job processing times. We find necessary and sufficient conditions (Theorem 1) when there exists a dominant permutation that is optimal for all possible realizations of the job processing times. Our computational studies show the percentage of the problems solvable under these conditions for the cases of randomly generated instances with n ≤ 100 . We also show how to use additional information about the processing times of the completed jobs during optimal realization of a schedule (Theorems 2 – 4). Computational studies for randomly generated instances with n ≤ 50 show the percentage of the two- machine flow-shop scheduling problems solvable under the sufficient conditions given in Theorems 2 – 4.

Visualizing Rare Watson-Crick-Like Tautomeric and Anionic Mismatches in DNA and RNA

The central dogma of molecular biology relies on the correct Watson-Crick (WC) geometry of canonical deoxyribonucleic acid (DNA) dG•dC and dA•dT base pairs to replicate and transcribe genetic information with speed and an astonishing level of fidelity. In addition, the Watson-Crick geometry of canonical ribonucleic acid (RNA) rG•rC and rA•rU base pairs is highly conserved to ensure that proteins are translated with high fidelity. However, numerous other potential nucleobase tautomeric and ionic configurations are possible that can give rise to entirely new pairing modes between the nucleotide bases. Very early on, James Watson and Francis Crick recognized their importance and in 1953 postulated that if bases adopted one of their less energetically disfavored tautomeric forms (and later ionic forms) during replication it could lead to the formation of a mismatch with a Watson-Crick-like geometry and could give rise to “natural mutations.”

Since this time numerous studies have provided evidence in support of this hypothesis and have expanded upon it; computational studies have addressed the energetic feasibilities of different nucleobases’ tautomeric and ionic forms in siico; crystallographic studies have trapped different mismatches with WC-like geometries in polymerase or ribosome active sites. However, no direct evidence has been given for (i) the direct existence of these WC-like mismatches in canonical DNA duplex, RNA duplexes, or non-coding RNAs; (ii) which, if any, tautomeric or ionic form stabilizes the WC-like geometry. This thesis utilizes nuclear magnetic resonance (NMR) spectroscopy and rotating frame relaxation dispersion (R1ρ RD) in combination with density functional theory (DFT), biochemical assays, and targeted chemical perturbations to show that (i) dG•dT mismatches in DNA duplexes, as well as rG•rU mismatches RNA duplexes and non-coding RNAs, transiently adopt a WC-like geometry that is stabilized by (ii) an interconnected network of rapidly interconverting rare tautomers and anionic bases. These results support Watson and Crick’s tautomer hypothesis, but additionally support subsequent hypotheses invoking anionic mismatches and ultimately tie them together. This dissertation shows that a common mismatch can adopt a Watson-Crick-like geometry globally, in both DNA and RNA, and whose geometry is stabilized by a kinetically linked network of rare tautomeric and anionic bases. The studies herein also provide compelling evidence for their involvement in spontaneous replication and translation errors.

Insight into the mechanism of galactokinase: role of a critical glutamate residue and helix/coil transitions

Galactokinase, the enzyme which catalyses the first committed step in the Leloir pathway, has attracted interest due to its potential as a biocatalyst and as a possible drug target in the treatment of type I galactosemia. The mechanism of the enzyme is not fully elucidated. Molecular dynamics (MD) simulations of galactokinase with the active site residues Arg-37 and Asp-186 altered predicted that two regions (residues 174-179 and 231-240) had different dynamics as a consequence. Interestingly, the same two regions were also affected by alterations in Arg-105, Glu-174 and Arg- 228. These three residues were identified as important in catalysis in previous computational studies on human galactokinase. Alteration of Arg-105 to methionine resulted in a modest reduction in activity with little change in stability. When Arg-228 was changed to methionine, the enzyme’s interaction with both ATP and galactose was affected. This variant was significantly less stable than the wild-type protein. Changing Glu-174 to glutamine (but not to aspartate) resulted in no detectable activity and a less stable enzyme. Overall, these combined in silico and in vitro studies demonstrate the importance of a negative charge at position 174 and highlight the critical role of the dynamics in to key regions of the protein. We postulate that these regions may be critical for mediating the enzyme’s structure and function.