Accelerating self-consistent field convergence with the augmented Roothaan-Hall energy function.

Based on Pulay's direct inversion iterative subspace (DIIS) approach, we present a method to accelerate self-consistent field (SCF) convergence. In this method, the quadratic augmented Roothaan-Hall (ARH) energy function, proposed recently by Høst and co-workers [J. Chem. Phys. 129, 124106 (2008)], is used as the object of minimization for obtaining the linear coefficients of Fock matrices within DIIS. This differs from the traditional DIIS of Pulay, which uses an object function derived from the commutator of the density and Fock matrices. Our results show that the present algorithm, abbreviated ADIIS, is more robust and efficient than the energy-DIIS (EDIIS) approach. In particular, several examples demonstrate that the combination of ADIIS and DIIS ("ADIIS+DIIS") is highly reliable and efficient in accelerating SCF convergence.

Convergence of numerical time-averaging and stationary measures via Poisson equations

Numerical approximation of the long time behavior of a stochastic di.erential equation (SDE) is considered. Error estimates for time-averaging estimators are obtained and then used to show that the stationary behavior of the numerical method converges to that of the SDE. The error analysis is based on using an associated Poisson equation for the underlying SDE. The main advantages of this approach are its simplicity and universality. It works equally well for a range of explicit and implicit schemes, including those with simple simulation of random variables, and for hypoelliptic SDEs. To simplify the exposition, we consider only the case where the state space of the SDE is a torus, and we study only smooth test functions. However, we anticipate that the approach can be applied more widely. An analogy between our approach and Stein's method is indicated. Some practical implications of the results are discussed. Copyright © by SIAM. Unauthorized reproduction of this article is prohibited.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

Impact of economic, regulatory, and patent policies on innovation in cancer chemoprevention.

Chemoprevention agents are an emerging new scientific area that holds out the promise of delaying or avoiding a number of common cancers. These new agents face significant scientific, regulatory, and economic barriers, however, which have limited investment in their research and development (R&D). These barriers include above-average clinical trial scales, lengthy time frames between discovery and Food and Drug Administration approval, liability risks (because they are given to healthy individuals), and a growing funding gap for early-stage candidates. The longer time frames and risks associated with chemoprevention also cause exclusivity time on core patents to be limited or subject to significant uncertainties. We conclude that chemoprevention uniquely challenges the structure of incentives embodied in the economic, regulatory, and patent policies for the biopharmaceutical industry. Many of these policy issues are illustrated by the recently Food and Drug Administration-approved preventive agents Gardasil and raloxifene. Our recommendations to increase R&D investment in chemoprevention agents include (a) increased data exclusivity times on new biological and chemical drugs to compensate for longer gestation periods and increasing R&D costs; chemoprevention is at the far end of the distribution in this regard; (b) policies such as early-stage research grants and clinical development tax credits targeted specifically to chemoprevention agents (these are policies that have been very successful in increasing R&D investment for orphan drugs); and (c) a no-fault liability insurance program like that currently in place for children's vaccines.

Using the institutional grammar tool to understand regulatory compliance: The case of Colorado aquaculture

What is the relationship between the design of regulations and levels of individual compliance? To answer this question, Crawford and Ostrom's institutional grammar tool is used to deconstruct regulations governing the aquaculture industry in Colorado, USA. Compliance with the deconstructed regulatory components is then assessed based on the perceptions of the appropriateness of the regulations, involvement in designing the regulations, and intrinsic and extrinsic motivations. The findings suggest that levels of compliance with regulations vary across and within individuals regarding various aspects of the regulatory components. As expected, the level of compliance is affected by the perceived appropriateness of regulations, participation in designing the regulations, and feelings of guilt and fear of social disapproval. Furthermore, there is a strong degree of interdependence among the written components, as identified by the institutional grammar tool, in affecting compliance levels. The paper contributes to the regulation and compliance literature by illustrating the utility of the institutional grammar tool in understanding regulatory content, applying a new Q-Sort technique for measuring individual levels of compliance, and providing a rare exploration into feelings of guilt and fear outside of the laboratory setting. © 2012 Blackwell Publishing Asia Pty Ltd.

A C-terminal motif found in the beta2-adrenergic receptor, P2Y1 receptor and cystic fibrosis transmembrane conductance regulator determines binding to the Na+/H+ exchanger regulatory factor family of PDZ proteins.

The Na+/H+ exchanger regulatory factor (NHERF) binds to the tail of the beta2-adrenergic receptor and plays a role in adrenergic regulation of Na+/H+ exchange. NHERF contains two PDZ domains, the first of which is required for its interaction with the beta2 receptor. Mutagenesis studies of the beta2 receptor tail revealed that the optimal C-terminal motif for binding to the first PDZ domain of NHERF is D-S/T-x-L, a motif distinct from those recognized by other PDZ domains. The first PDZ domain of NHERF-2, a protein that is 52% identical to NHERF and also known as E3KARP, SIP-1, and TKA-1, exhibits binding preferences very similar to those of the first PDZ domain of NHERF. The delineation of the preferred binding motif for the first PDZ domain of the NHERF family of proteins allows for predictions for other proteins that may interact with NHERF or NHERF-2. For example, as would be predicted from the beta2 receptor tail mutagenesis studies, NHERF binds to the tail of the purinergic P2Y1 receptor, a seven-transmembrane receptor with an intracellular C-terminal tail ending in D-T-S-L. NHERF also binds to the tail of the cystic fibrosis transmembrane conductance regulator, which ends in D-T-R-L. Because the preferred binding motif of the first PDZ domain of the NHERF family of proteins is found at the C termini of a variety of intracellular proteins, NHERF and NHERF-2 may be multifunctional adaptor proteins involved in many previously unsuspected aspects of intracellular signaling.

Impaired formation of beta-adrenergic receptor-nucleotide regulatory protein complexes in pseudohypoparathyroidism.

Decreased activity of the guanine nucleotide regulatory protein (N) of the adenylate cyclase system is present in cell membranes of some patients with pseudohypoparathyrodism (PHP-Ia) whereas others have normal activity of N (PHP-Ib). Low N activity in PHP-Ia results in a decrease in hormone (H)-stimulatable adenylate cyclase in various tissues, which might be due to decreased ability to form an agonist-specific high affinity complex composed of H, receptor (R), and N. To test this hypothesis, we compared beta-adrenergic agonist-specific binding properties in erythrocyte membranes from five patients with PHP-Ia (N = 45% of control), five patients with PHP-Ib (N = 97%), and five control subjects. Competition curves that were generated by increasing concentrations of the beta-agonist isoproterenol competing with [125I]pindolol were shallow (slope factors less than 1) and were computer fit to a two-state model with corresponding high and low affinity for the agonist. The agonist competition curves from the PHP-Ia patients were shifted significantly (P less than 0.02) to the right as a result of a significant (P less than 0.01) decrease in the percent of beta-adrenergic receptors in the high affinity state from 64 +/- 22% in PHP-Ib and 56 +/- 5% in controls to 10 +/- 8% in PHP-Ia. The agonist competition curves were computer fit to a "ternary complex" model for the two-step reaction: H + R + N in equilibrium HR + N in equilibrium HRN. The modeling was consistent with a 60% decrease in the functional concentration of N, and was in good agreement with the biochemically determined decrease in erythrocyte N protein activity. These in vitro findings in erythrocytes taken together with the recent observations that in vivo isoproterenol-stimulated adenylate cyclase activity is decreased in patients with PHP (Carlson, H. E., and A. S. Brickman, 1983, J. Clin. Endocrinol. Metab. 56:1323-1326) are consistent with the notion that N is a bifunctional protein interacting with both R and the adenylate cyclase. It may be that in patients with PHP-Ia a single molecular and genetic defect accounts for both decreased HRN formation and decreased adenylate cyclase activity, whereas in PHP-Ib the biochemical lesion(s) appear not to affect HRN complex formation.

Cutting edge: distinct glycolytic and lipid oxidative metabolic programs are essential for effector and regulatory CD4+ T cell subsets.

Stimulated CD4(+) T lymphocytes can differentiate into effector T cell (Teff) or inducible regulatory T cell (Treg) subsets with specific immunological roles. We show that Teff and Treg require distinct metabolic programs to support these functions. Th1, Th2, and Th17 cells expressed high surface levels of the glucose transporter Glut1 and were highly glycolytic. Treg, in contrast, expressed low levels of Glut1 and had high lipid oxidation rates. Consistent with glycolysis and lipid oxidation promoting Teff and Treg, respectively, Teff were selectively increased in Glut1 transgenic mice and reliant on glucose metabolism, whereas Treg had activated AMP-activated protein kinase and were dependent on lipid oxidation. Importantly, AMP-activated protein kinase stimulation was sufficient to decrease Glut1 and increase Treg generation in an asthma model. These data demonstrate that CD4(+) T cell subsets require distinct metabolic programs that can be manipulated in vivo to control Treg and Teff development in inflammatory diseases.

Rapamycin Interferes With Postdepletion Regulatory T Cell Homeostasis and Enhances DSA Formation Corrected by CTLA4-Ig.

Previously, we demonstrated that alemtuzumab induction with rapamycin as sole maintenance therapy is associated with an increased incidence of humoral rejection in human kidney transplant patients. To investigate the role of rapamycin in posttransplant humoral responses after T cell depletion, fully MHC mismatched hearts were transplanted into hCD52Tg mice, followed by alemtuzumab treatment with or without a short course of rapamycin. While untreated hCD52Tg recipients acutely rejected B6 hearts (n = 12), hCD52Tg recipients treated with alemtuzumab alone or in conjunction with rapamycin showed a lack of acute rejection (MST > 100). However, additional rapamycin showed a reduced beating quality over time and increased incidence of vasculopathy. Furthermore, rapamycin supplementation showed an increased serum donor-specific antibodies (DSA) level compared to alemtuzumab alone at postoperation days 50 and 100. Surprisingly, additional rapamycin treatment significantly reduced CD4(+) CD25(+) FoxP3(+) T reg cell numbers during treatment. On the contrary, ICOS(+) PD-1(+) CD4 follicular helper T cells in the lymph nodes were significantly increased. Interestingly, CTLA4-Ig supplementation in conjunction with rapamycin corrected rapamycin-induced accelerated posttransplant humoral response by directly modulating Tfh cells but not Treg cells. This suggests that rapamycin after T cell depletion could affect Treg cells leading to an increase of Tfh cells and DSA production that can be reversed by CTLA4-Ig.