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.

Apoptosis and caspases regulate death and inflammation in sepsis.

Although the prevailing concept has been that mortality in sepsis results from an unbridled hyper-inflammatory cytokine-mediated response, the failure of more than 30 clinical trials to treat sepsis by controlling this cytokine response requires a 'rethink' of the molecular mechanism underpinning the development of sepsis. As we discuss here, remarkable new studies indicate that most deaths from sepsis are actually the result of a substantially impaired immune response that is due to extensive death of immune effector cells. Rectification of this apoptotic-inflammatory imbalance using modulators of caspases and other components of the cell-death pathway have shown striking efficacy in stringent animal models of sepsis, indicating an entirely novel path forward for the clinical treatment of human sepsis.

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.

Enhanced bacterial clearance and sepsis resistance in caspase-12-deficient mice.

Caspases function in both apoptosis and inflammatory cytokine processing and thereby have a role in resistance to sepsis. Here we describe a novel role for a caspase in dampening responses to bacterial infection. We show that in mice, gene-targeted deletion of caspase-12 renders animals resistant to peritonitis and septic shock. The resulting survival advantage was conferred by the ability of the caspase-12-deficient mice to clear bacterial infection more efficiently than wild-type littermates. Caspase-12 dampened the production of the pro-inflammatory cytokines interleukin (IL)-1beta, IL-18 (interferon (IFN)-gamma inducing factor) and IFN-gamma, but not tumour-necrosis factor-alpha and IL-6, in response to various bacterial components that stimulate Toll-like receptor and NOD pathways. The IFN-gamma pathway was crucial in mediating survival of septic caspase-12-deficient mice, because administration of neutralizing antibodies to IFN-gamma receptors ablated the survival advantage that otherwise occurred in these animals. Mechanistically, caspase-12 associated with caspase-1 and inhibited its activity. Notably, the protease function of caspase-12 was not necessary for this effect, as the catalytically inactive caspase-12 mutant Cys299Ala also inhibited caspase-1 and IL-1beta production to the same extent as wild-type caspase-12. In this regard, caspase-12 seems to be the cFLIP counterpart for regulating the inflammatory branch of the caspase cascade. In mice, caspase-12 deficiency confers resistance to sepsis and its presence exerts a dominant-negative suppressive effect on caspase-1, resulting in enhanced vulnerability to bacterial infection and septic mortality.

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.

Recruitment and activation of caspase-8 by the Huntingtin-interacting protein Hip-1 and a novel partner Hippi.

In Huntington disease, polyglutamine expansion of the protein huntingtin (Htt) leads to selective neurodegenerative loss of medium spiny neurons throughout the striatum by an unknown apoptotic mechanism. Binding of Hip-1, a protein normally associated with Htt, is reduced by polyglutamine expansion. Free Hip-1 binds to a hitherto unknown polypeptide, Hippi (Hip-1 protein interactor), which has partial sequence homology to Hip-1 and similar tissue and subcellular distribution. The availability of free Hip-1 is modulated by polyglutamine length within Htt, with disease-associated polyglutamine expansion favouring the formation of pro-apoptotic Hippi-Hip-1 heterodimers. This heterodimer can recruit procaspase-8 into a complex of Hippi, Hip-1 and procaspase-8, and launch apoptosis through components of the 'extrinsic' cell-death pathway. We propose that Htt polyglutamine expansion liberates Hip-1 so that it can form a caspase-8 recruitment complex with Hippi. This novel non-receptor-mediated pathway for activating caspase-8 might contribute to neuronal death in Huntington disease.

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.

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.

Differential efficacy of caspase inhibitors on apoptosis markers during sepsis in rats and implication for fractional inhibition requirements for therapeutics.

A rodent model of sepsis was used to establish the relationship between caspase inhibition and inhibition of apoptotic cell death in vivo. In this model, thymocyte cell death was blocked by Bcl-2 transgene, indicating that apoptosis was predominantly dependent on the mitochondrial pathway that culminates in caspase-3 activation. Caspase inhibitors, including the selective caspase-3 inhibitor M867, were able to block apoptotic manifestations both in vitro and in vivo but with strikingly different efficacy for different cell death markers. Inhibition of DNA fragmentation required substantially higher levels of caspase-3 attenuation than that required for blockade of other apoptotic events such as spectrin proteolysis and phosphatidylserine externalization. These data indicate a direct relationship between caspase inhibition and some apoptotic manifestations but that small quantities of uninhibited caspase-3 suffice to initiate genomic DNA breakdown, presumably through the escape of catalytic quantities of caspase-activated DNase. These findings suggest that putative caspase-independent apoptosis may be overestimated in some systems since blockade of spectrin proteolysis and other cell death markers does not accurately reflect the high degrees of caspase-3 inhibition needed to prevent DNA fragmentation. Furthermore, this requirement presents substantial therapeutic challenges owing to the need for persistent and complete caspase blockade.

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.

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.

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.

Differential modulation of endotoxin responsiveness by human caspase-12 polymorphisms.

Caspases mediate essential key proteolytic events in inflammatory cascades and the apoptotic cell death pathway. Human caspases functionally segregate into two distinct subfamilies: those involved in cytokine maturation (caspase-1, -4 and -5) and those involved in cellular apoptosis (caspase-2, -3, -6, -7, -8, -9 and -10). Although caspase-12 is phylogenetically related to the cytokine maturation caspases, in mice it has been proposed as a mediator of apoptosis induced by endoplasmic reticulum stress including amyloid-beta cytotoxicity, suggesting that it might contribute to the pathogenesis of Alzheimer's disease. Here we show that a single nucleotide polymorphism in caspase-12 in humans results in the synthesis of either a truncated protein (Csp12-S) or a full-length caspase proenzyme (Csp12-L). The read-through single nucleotide polymorphism encoding Csp12-L is confined to populations of African descent and confers hypo-responsiveness to lipopolysaccharide-stimulated cytokine production in ex vivo whole blood, but has no significant effect on apoptotic sensitivity. In a preliminary study, we find that the frequency of the Csp12-L allele is increased in African American individuals with severe sepsis. Thus, Csp12-L attenuates the inflammatory and innate immune response to endotoxins and in doing so may constitute a risk factor for developing sepsis.

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.

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.

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.

Axonal dynactin p150Glued transports caspase-8 to drive retrograde olfactory receptor neuron apoptosis.

Olfactory receptor neurons (ORNs) undergo caspase-mediated retrograde apoptosis after target removal (bulbectomy), in which axonal caspase-9 and caspase-3 activation leads to terminal apoptosis in ORN soma of the olfactory epithelium. Here, we show that caspase-8 can act as an initiator of ORN apoptosis after bulbectomy and also after synaptic instability is induced by NMDA-mediated excitotoxic death of ORN target neurons in the olfactory bulb. Caspase-8 and caspase-3 are sequentially activated within ORN presynaptic terminals, and caspase-8 complexes with dynactin p150Glued, (a retrograde motor protein) and is transported retrogradely, preceding axonal caspase-3 activation and apoptosis of ORN cell bodies. Focal in vivo inhibition of initiator caspase activation or microtubule-dependent transport (with Taxol) at the lesioned axon terminus results in a significant reduction in retrograde axonal caspase-8 and caspase-3 activation and inhibition of retrograde ORN death. Caspase-8 activation and retrograde transport after NMDA lesion is similarly reduced in mice null for p75, the low-affinity nerve growth factor receptor. The retrograde apoptosis of ORNs thus involves a novel mechanism that used p75 in the local activation of caspase-8. Once caspase-8 is maximally activated in the presynaptic terminal, it is transported retrogradely by the motor complex dynactin/dynein, a process that can be inhibited focally to inhibit ORN apoptosis after acute axonal lesion. These data have revealed a novel mechanism of retrograde apoptosis, in which caspase-8 complexes directly with axonal dynactin p150Glued to reveal a differential vulnerability of subpopulations of ORNs to undergo apoptosis after axonal damage and the loss of olfactory bulb target neurons.