Exploration of Entropy Pair Functional Theory.

Evaluation of the entropy from molecular dynamics (MD) simulation remains an outstanding challenge. The standard approach requires thermodynamic integration across a series of simulations. Recent work Nicholson et al. demonstrated the ability to construct a functional that returns excess entropy, based on the pair correlation function (PCF); it was capable of providing, with acceptable accuracy, the absolute excess entropy of iron simulated with a pair potential in both fluid and crystalline states. In this work, the general applicability of the Entropy Pair Functional Theory (EPFT) approach is explored by applying it to three many-body interaction potentials. These potentials are state of the art for large scale models for the three materials in this study: Fe modelled with a modified embedded atom method (MEAM) potential, Cu modelled with an MEAM and Si modelled with a Tersoff potential. We demonstrate the robust nature of EPFT in determining excess entropy for diverse systems with many-body interactions. These are steps toward a universal Entropy Pair Functional, EPF, that can be applied with confidence to determine the entropy associated with sophisticated optimized potentials and first principles simulations of liquids, crystals, engineered structures, and defects.

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions.

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson's theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs-Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born-Oppenheimer molecular dynamics performed at a single temperature. This advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

Lung eQTLs to help reveal the molecular underpinnings of asthma.

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55 × 10(-151)). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases.

Molecular simulations of adsorption of RDX and TATP on IRMOF-1(Be).

The influence of different sorption sites of isoreticular metal-organic frameworks (IRMOFs) on interactions with explosive molecules is investigated. Different connector effects are taken into account by choosing IRMOF-1(Be) (IRMOF-1 with Zn replaced by Be), and two high explosive molecules: 1,3,5-trinitro-s-triazine (RDX) and triacetone triperoxide (TATP). The key interaction features (structural, electronic and energetic) of selected contaminants were analyzed by means of density functional calculations. The interaction of RDX and TATP with different IRMOF-1(Be) fragments is studied. The results show that physisorption is favored and occurs due to hydrogen bonding, which involves the C-H groups of both molecules and the carbonyl oxygen atoms of IRMOF-1(Be). Additional stabilization of RDX and TATP arises from weak electrostatic interactions. Interaction with IRMOF-1(Be) fragments leads to polarization of the target molecules. Of the molecular configurations we have studied, the Be-O-C cluster connected with six benzene linkers (1,4-benzenedicarboxylate, BDC), possesses the highest binding energy for the studied explosives (-16.4 kcal mol(-1) for RDX and -12.9 kcal mol(-1) for TATP). The main difference was discovered to be in the preferable adsorption site for adsorbates (RDX above the small and TATP placed above the big cage). Based on these results, IRMOF-1 can be suggested as an effective material for storage and also for separation of similar explosives. Hydration destabilizes most of the studied adsorption systems by 1-3 kcal mol(-1) but it leads to the same trend in the binding strength as found for the non-hydrated complexes.

A coarse-grained model for polyethylene glycol polymer.

A coarse-grained (CG) model of polyethylene glycol (PEG) was developed and implemented in CG molecular dynamics (MD) simulations of PEG chains with degree of polymerization (DP) 20 and 40. In the model, two repeat units of PEG are grouped as one CG bead. Atomistic MD simulation of PEG chains with DP = 20 was first conducted to obtain the bonded structural probability distribution functions (PDFs) and nonbonded pair correlation function (PCF) of the CG beads. The bonded CG potentials are obtained by simple inversion of the corresponding PDFs. The CG nonbonded potential is parameterized to the PCF using both an inversion procedure based on the Ornstein-Zernike equation with the Percus-Yevick approximation (OZPY(-1)) and a combination of OZPY(-1) with the iterative Boltzmann inversion (IBI) method (OZPY(-1)+IBI). As a simple one step method, the OZPY(-1) method possesses an advantage in computational efficiency. Using the potential from OZPY(-1) as an initial guess, the IBI method shows fast convergence. The coarse-grained molecular dynamics (CGMD) simulations of PEG chains with DP = 20 using potentials from both methods satisfactorily reproduce the structural properties from atomistic MD simulation of the same systems. The OZPY(-1)+IBI method yields better agreement than the OZPY(-1) method alone. The new CG model and CG potentials from OZPY(-1)+IBI method was further tested through CGMD simulation of PEG with DP = 40 system. No significant changes are observed in the comparison of PCFs from CGMD simulations of PEG with DP = 20 and 40 systems indicating that the potential is independent of chain length.

Applications of a general random-walk theory for confined diffusion.

A general random walk theory for diffusion in the presence of nanoscale confinement is developed and applied. The random-walk theory contains two parameters describing confinement: a cage size and a cage-to-cage hopping probability. The theory captures the correct nonlinear dependence of the mean square displacement (MSD) on observation time for intermediate times. Because of its simplicity, the theory also requires modest computational requirements and is thus able to simulate systems with very low diffusivities for sufficiently long time to reach the infinite-time-limit regime where the Einstein relation can be used to extract the self-diffusivity. The theory is applied to three practical cases in which the degree of order in confinement varies. The three systems include diffusion of (i) polyatomic molecules in metal organic frameworks, (ii) water in proton exchange membranes, and (iii) liquid and glassy iron. For all three cases, the comparison between theory and the results of molecular dynamics (MD) simulations indicates that the theory can describe the observed diffusion behavior with a small fraction of the computational expense. The confined-random-walk theory fit to the MSDs of very short MD simulations is capable of accurately reproducing the MSDs of much longer MD simulations. Furthermore, the values of the parameter for cage size correspond to the physical dimensions of the systems and the cage-to-cage hopping probability corresponds to the activation barrier for diffusion, indicating that the two parameters in the theory are not simply fitted values but correspond to real properties of the physical system.

Discovery of MK-0952, a selective PDE4 inhibitor for the treatment of long-term memory loss and mild cognitive impairment.

The structure-activity relationship of a novel series of 8-biarylnaphthyridinones acting as type 4 phosphodiesterase (PDE4) inhibitors for the treatment of long-term memory loss and mild cognitive impairment is described herein. The manuscript describes a new paradigm for the development of PDE4 inhibitor targeting CNS indications. This effort led to the discovery of the clinical candidate MK-0952, an intrinsically potent inhibitor (IC(50)=0.6 nM) displaying limited whole blood activity (IC(50)=555 nM). Supporting in vivo results in two preclinical efficacy tests and one test assessing adverse effects are also reported. The comparative profiles of MK-0952 and two other Merck compounds are described to validate the proposed hypothesis.

Use of the Ornstein-Zernike Percus-Yevick equation to extract interaction potentials from pair correlation functions.

In this work, we test the ability of the Ornstein-Zernike equation in the Percus-Yevick approximation (OZPY) to generate interaction potentials from pair correlation functions (PCFs) of monatomic and diatomic Lennard-Jones fluids. The PCFs are generated by solving OZPY equation (monatomic fluid) and molecular-dynamics (MD) simulations (diatomic fluid). Since the interaction potentials are inputs in the OZPY method and the MD simulation, the extraction of the potential from the PCFs using OZPY is a test of self-consistency. This test is necessary if the procedure is to be used to generate coarse-grained (CG) potentials from PCFs. We find that the procedure is completely self-consistent for the monatomic fluid in the whole range of densities studied (reduced density up to 0.55, under reduced temperature of 2.0). In the diatomic case, we find that the procedure is generally self-consistent under both low and high densities, although there is a systematic deviation at high densities. The method is able to reproduce the two parameters (ε and σ) of the input Lennard-Jones potential model to within about 1%. This CG potential generating procedure can be straightforwardly extended to more complicated molecules.

The discovery and synthesis of highly potent subtype selective phosphodiesterase 4D inhibitors.

The SAR study of a series of 6-aryloxymethyl-8-aryl substituted quinolines is described. Optimization of the series led to the discovery of compound 26b, a highly potent (IC50=0.6 nM) and selective PDE4D inhibitor with a 75-fold selectivity over the A, B, and C subtypes and over 18,000-fold selectivity against other PDE family members. Rat pharmacokinetics and tissue distribution are also summarized.

Effect of charge distribution on RDX adsorption in IRMOF-10.

Quantum mechanical (QM) calculations, classical grand canonical Monte Carlo (GCMC) simulations, and classical molecular dynamics (MD) simulations are performed to test the effect of charge distribution on hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) adsorption and diffusion in IRMOF-10. Several different methods for mapping QM electron distributions onto atomic point charges are explored, including the electrostatic potential (ESP) method, Mulliken population analysis, Lowdin population analysis, and natural bond orbital analysis. Classical GCMC and MD simulations of RDX in IRMOF-10 are performed using 15 combinations of charge sources of RDX and IRMOF-10. As the charge distributions vary, interaction potential energies, the adsorption loading, and the self-diffusivities are significantly different. None of the 15 combinations are able to quantitatively capture the dependence of the energy of adsorption on local configuration of RDX as observed in the QM calculations. We observe changes in the charge distributions of RDX and IRMOF-10 with the introduction of an RDX molecule into the cage. We also observe a large dispersion contribution to the interaction energy from QM calculations that is not reproduced in the classical simulations, indicating that the source of discrepancy may not lie exclusively with the assignment of charges.

Alkyl-bridged substituted 8-arylquinolines as highly potent PDE IV inhibitors.

Substituted 8-arylquinoline analogs bearing alkyl-linked side chain were identified as potent inhibitors of type 4 phophodiesterase. These compounds address the potential liabilities of the clinical candidate L-454560. The pharmacokinetic profile of the best analogs and the in vivo efficacy in an ovalbumin-induced bronchoconstriction assay in conscious guinea pigs are reported.

Gender differences in expression of the human caspase-12 long variant determines susceptibility to Listeria monocytogenes infection.

Inflammatory caspases are important effectors of innate immunity. Caspase-12, of the inflammatory caspase subfamily, is expressed in all mammals tested to date, but has acquired deleterious mutation in humans. A single-nucleotide polymorphism introduces a premature stop codon in caspase-12 in the majority of the population. However, in 20% of African descendants, caspase-12 is expressed and sensitizes to infections and sepsis. Here, we examined the modalities by which human caspase-12 confers susceptibility to infection. We have generated a fully humanized mouse that expresses the human caspase-12 rare variant (Csp-12L) in a mouse casp-12(-/-) background. Characterization of the humanized mouse uncovered sex differences in Csp-12L expression and gender disparity in innate immunity to Listeria monocytogenes infection. The Csp-12L transgene completely reversed the knockout resistance-to-infection phenotype in casp-12(-/-) males. In contrast, it had a marginal effect on the response of female mice. We found that estrogen levels modulated the expression of caspase-12. Csp-12L was expressed in male mice but its expression was repressed in female mice. Administration of 17-beta-estradiol (E2) to humanized male mice had a direct suppressive effect on Csp-12L expression and conferred relative resistance to infection. Chromatin immunoprecipitation experiments revealed that caspase-12 is a direct transcriptional target of the estrogen receptor alpha (ERalpha) and mapped the estrogen response element (ERE) to intron 7 of the gene. We propose that estrogen-mediated inhibition of Csp-12L expression is a built-in mechanism that has evolved to protect females from infection.

Characterization of alpha(4)beta(1) (CD49d/CD29) on equine leukocytes: potential utility of a potent alpha(4)beta(1) (CD49d/CD29) receptor antagonist in the treatment of equine heaves (recurrent airway obstruction).

The purpose of this study was to characterize the alpha(4)beta(1) receptor (CD49d/CD29, very late antigen-4, VLA-4) on circulating equine leukocytes and to evaluate the intrinsic potency of an alpha(4)beta(1) receptor antagonist (Compound B) in the horse. Ultimately, these studies would allow us to determine the suitability of treating recurrent airway obstruction (RAO; heaves) affected horses by blocking the cellular recruitment of lymphocytes and neutrophils into the lung. The data demonstrates the alpha(4)beta(1) integrin is present on horse lymphocytes and neutrophils (fluorescence-assisted cell sorter, FACS) and can bind low molecular weight alpha(4)beta(1) antagonists (Compounds A and B) with high affinity. K(D) values for the binding of Compound A to non-activated alpha(4)beta(1) on isolated horse PBMCs (peripheral blood mononuclear cells) and activated neutrophils were 17 pM and 27 pM, respectively. Compound B was identified as a suitable antagonist for performing a series of in vivo experiments. Compound B was found to possess excellent potency in horse whole blood, possessing IC(50) and IC(90) values of 39 pM and 172 pM, respectively. This represents a 3.9-fold molar excess of drug over the alpha(4)beta(1) concentration in blood. Following oral administration of Compound B (5 mg/kg) to beagle dogs and rhesus monkeys, rapid and sustained alpha(4)beta(1) receptor occupancy (>80%) was achieved and maintained for a period of 24 h. When Compound B was administered intravenously to the horse, by either a slow or rapid infusion at a dose of 0.3 mg/kg, receptor blockade of >80% was observed out to 24 h with a concomitant leukocytosis. We believe that Compound B possesses suitable intrinsic and pharmacological properties to be evaluated clinically in horses affected by RAO.

Huntingtin interacting protein 1 can regulate neurogenesis in Drosophila.

Huntington's disease (HD) is associated with a range of cellular consequences including selective neuronal death and decreased levels of neurogenesis. Ultimately, these altered processes are dependent upon proteins that interact with Huntingtin (Htt) such as the Huntingtin-interacting protein 1 (Hip1) which has a reduced binding preference to expanded Htt. These effects are similar to those observed with modified Notch signal transduction. As Hip1 plays a key role in endocytosis and intracellular transport, and activation of the Notch signal requires both, we investigated putative links between Hip1 and Notch signaling in flies. We have identified two forms of Hip1 that may be produced through the use of alternative first exons: a version of Hip1 with a lipid-binding ANTH domain and Hip1DeltaANTH lacking this domain. The directed expression of Hip1 decreases, while expression of Hip1DeltaANTH increases, the density of sensory microchaetae on the dorsal notum, a classical model of neurogenesis. A reduction in microchaetae density associated with Notch(Microchaetae Deficient (MCD)) (N(MCD) ) alleles is sensitive to both Hip1 and Hip1DeltaANTH levels, as are the bristle phenotypes generated by misexpression of deltex, a key mediator of Notch signaling. Genetic studies further demonstrate that the observed effects of Hip1 and of Hip1DeltaANTH are sensitive to achaete gene dosage while insensitive to the levels of E(Spl), suggesting a non-canonical Notch neurogenic signal through a deltex-dependent pathway. The novel role we describe for Hip1 in Notch-mediated neurogenesis provides a functional link between Notch signaling and proteins related to HD.

Confinement of caspase-12 proteolytic activity to autoprocessing.

Caspase-12 is a dominant-negative regulator of caspase-1 (IL-1beta-converting enzyme) and an attenuator of cytokine responsiveness to septic infections. This molecular role for caspase-12 appears to be akin to the role of cFLIP in regulating caspase-8 in the extrinsic cell death pathway; however, unlike cFLIP/Usurpin, we demonstrate here that caspase-12 is catalytically competent. To examine these catalytic properties, rat caspase-12 was cloned, and the recombinant enzyme was used to examine the cleavage of macromolecular and synthetic fluorogenic substrates. Although caspase-12 could mediate autoproteolytic maturation of its own proenzyme, in both cis and trans, it was not able to cleave any other polypeptide substrate, including other caspase proenzymes, apoptotic substrates, cytokine precursors, or proteins in the endoplasmic reticulum that normally undergo caspase-mediated proteolysis. The dearth of potential substrates for caspase-12 also was confirmed by whole-cell diagonal-gel analysis. Autolytic cleavage within the caspase-12 proenzyme was mapped to a single site at the large-small subunit junction, ATAD(319), and this motif was recognized by caspase-12 when incorporated into synthetic fluorogenic substrates. The specific activity of caspase-12 with these substates was several orders of magnitude lower than caspases-1 and -3, highlighting its relative catalytic paucity. In intact cells, caspase-12 autoproteolysis occurred in the inhibitory complex containing caspase-1. We propose that the proteolytic activity of caspase-12 is confined to its own proenzyme and that autocleavage within the caspase-1 complex may be a means for temporal limitation of the inhibitory effects of caspase-12 on proinflammatory cytokine maturation.

Design, synthesis, and biological evaluation of 8-biarylquinolines: a novel class of PDE4 inhibitors.

The structure-activity relationship of a novel series of 8-biarylquinolines acting as type 4 phosphodiesterase (PDE4) inhibitors is described herein. Prototypical compounds from this series are potent and non-selective inhibitors of the four distinct PDE4 (IC(50)<10 nM) isozymes (A-D). In a human whole blood in vitro assay, they inhibit (IC(50)<0.5 microM) the LPS-induced release of the cytokine TNF-alpha. Optimized inhibitors were evaluated in vivo for efficacy in an ovalbumin-induced bronchoconstriction model in conscious guinea pigs. Their propensity to produce an emetic response was evaluated by performing pharmacokinetic studies in squirrel monkeys. This work has led to the identification of several compounds with excellent in vitro and in vivo profiles, including a good therapeutic window of efficacy over emesis.

Solid-phase analogue synthesis of caspase-3 inhibitors via palladium-catalyzed amination of 3-bromopyrazinones.

Caspase-3 is a cysteinyl protease that mediates apoptotic cell death. Its inhibition may have an important impact on the treatment of several degenerative diseases. Here we report the synthesis of reversible inhibitors via a solid-support palladium-catalyzed amination of 3-bromopyrazinones and the discovery of a pan-caspase reversible inhibitor.

Hippi is essential for node cilia assembly and Sonic hedgehog signaling.

Hippi functions as an adapter protein that mediates pro-apoptotic signaling from poly-glutamine-expanded huntingtin, an established cause of Huntington disease, to the extrinsic cell death pathway. To explore other functions of Hippi we generated Hippi knock-out mice. This deletion causes randomization of the embryo turning process and heart looping, which are hallmarks of defective left-right (LR) axis patterning. We report that motile monocilia normally present at the surface of the embryonic node, and proposed to initiate the break in LR symmetry, are absent on Hippi-/- embryos. Furthermore, defects in central nervous system development are observed. The Sonic hedgehog (Shh) pathway is downregulated in the neural tube in the absence of Hippi, which results in failure to establish ventral neural cell fate. Together, these findings demonstrate a dual role for Hippi in cilia assembly and Shh signaling during development, in addition to its proposed role in apoptosis signal transduction in the adult brain under pathogenically stressful conditions.