Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

A Large Gene Network in Immature Erythroid Cells Is Controlled by the Myeloid and B Cell Transcriptional Regulator PU.1

PU.1 is a hematopoietic transcription factor that is required for the development of myeloid and B cells. PU.1 is also expressed in erythroid progenitors, where it blocks erythroid differentiation by binding to and inhibiting the main erythroid promoting factor, GATA-1. However, other mechanisms by which PU.1 affects the fate of erythroid progenitors have not been thoroughly explored. Here, we used ChIP-Seq analysis for PU.1 and gene expression profiling in erythroid cells to show that PU.1 regulates an extensive network of genes that constitute major pathways for controlling growth and survival of immature erythroid cells. By analyzing fetal liver erythroid progenitors from mice with low PU.1 expression, we also show that the earliest erythroid committed cells are dramatically reduced in vivo. Furthermore, we find that PU.1 also regulates many of the same genes and pathways in other blood cells, leading us to propose that PU.1 is a multifaceted factor with overlapping, as well as distinct, functions in several hematopoietic lineages

A comparison of probe-level and probeset models for small-sample gene expression data

<p>Abstract</p> <p>Background</p> <p>Statistical methods to tentatively identify differentially expressed genes in microarray studies typically assume larger sample sizes than are practical or even possible in some settings.</p> <p>Results</p> <p>The performance of several probe-level and probeset models was assessed graphically and numerically using three spike-in datasets. Based on the Affymetrix GeneChip, a novel nested factorial model was developed and found to perform competitively on small-sample spike-in experiments.</p> <p>Conclusions</p> <p>Statistical methods with test statistics related to the estimated log fold change tend to be more consistent in their performance on small-sample gene expression data. For such small-sample experiments, the nested factorial model can be a useful statistical tool. This method is implemented in freely-available R code (affyNFM), available with a tutorial document at <url>http://www.stat.usu.edu/~jrstevens</url>.</p

Diversity of HIV-1 Subtype B: Implications to the Origin of BF Recombinants

BACKGROUND: The HIV-1 subtype B epidemic in Brazil is peculiar because of the high frequency of isolates having the GWGR tetramer at V3 loop region. It has been suggested that GWGR is a distinct variant and less pathogenic than other subtype B isolates. METHODOLOGY/PRINCIPAL FINDINGS: Ninety-four percent of the HIV-1 subtype B worldwide sequences (7689/8131) obtained from the Los Alamos HIV database contain proline at the tetramer of the V3 loop of the env gene (GPGR) and only 0.74% (60/8131) have tryptophan (GWGR). By contrast, 48.4% (161/333) of subtype B isolates from Brazil have proline, 30.6% (102/333) contain tryptophan and 10.5% (35/333) have phenylalanine (F) at the second position of the V3 loop tip. The proportion of tryptophan and phenylalanine in Brazilian isolates is much higher than in worldwide subtype B sequences (chi-square test, p = 0.0001). The combined proportion of proline, tryptophan and phenylalanine (GPGR+GWGR+GFGR) of Brazilian isolates corresponds to 89% of all amino acids in the V3 loop. Phylogenetic analysis revealed that almost all subtype B isolates in Brazil have a common origin regardless of their motif (GWGR, GPGR, GGGR, etc.) at the V3 tetramer. This shared ancestral origin was also observed in CRF28_BF and CRF29_BF in a genome region (free of recombination) derived from parental subtype B. These results imply that tryptophan substitution (e.g., GWGR-to-GxGR), which was previously associated with the change in the coreceptor usage within the host, also occurs at the population level. CONCLUSIONS/SIGNIFICANCE: Based on the current findings and previous study showing that tryptophan and phenylalanine in the V3 loop are related with coreceptor usage, we propose that tryptophan and phenylalanine in subtype B isolates in Brazil are kept by selective mechanisms due to the distinct coreceptor preferences in target cells of GWGR, GFGR and GFGR viruses

Allyl Isothiocyanate that Induces GST and UGT Expression Confers Oxidative Stress Resistance on C. elegans, as Demonstrated by Nematode Biosensor

Electrophilic xenobiotics and endogenous products from oxidative stresses induce the glutathione S-transferases (GSTs), which form a large family within the phase II enzymes over both animal and plant kingdoms. The GSTs thus induced in turn detoxify these external as well as internal stresses. Because these stresses are often linked to ageing and damage to health, the induction of phase II enzymes without causing adverse effects would be beneficial in slowing down ageing and keeping healthy conditions. for use as a nematode biosensor. With the nematode biosensor, we found that AITC induced GST expression and conferred tolerance on the nematode against various oxidative stresses. We also present evidence that the transcription factor SKN-1 is involved in regulating the GST expression induced by AITC.We show the applicability of the nematode biosensor for discovering and evaluating functional food substances and chemicals that would provide anti-ageing or healthful benefits

Gammaherpesvirus-Driven Plasma Cell Differentiation Regulates Virus Reactivation from Latently Infected B Lymphocytes

Gammaherpesviruses chronically infect their host and are tightly associated with the development of lymphoproliferative diseases and lymphomas, as well as several other types of cancer. Mechanisms involved in maintaining chronic gammaherpesvirus infections are poorly understood and, in particular, little is known about the mechanisms involved in controlling gammaherpesvirus reactivation from latently infected B cells in vivo. Recent evidence has linked plasma cell differentiation with reactivation of the human gammaherpesviruses EBV and KSHV through induction of the immediate-early viral transcriptional activators by the plasma cell-specific transcription factor XBP-1s. We now extend those findings to document a role for a gammaherpesvirus gene product in regulating plasma cell differentiation and thus virus reactivation. We have previously shown that the murine gammaherpesvirus 68 (MHV68) gene product M2 is dispensable for virus replication in permissive cells, but plays a critical role in virus reactivation from latently infected B cells. Here we show that in mice infected with wild type MHV68, virus infected plasma cells (ca. 8% of virus infected splenocytes at the peak of viral latency) account for the majority of reactivation observed upon explant of splenocytes. In contrast, there is an absence of virus infected plasma cells at the peak of latency in mice infected with a M2 null MHV68. Furthermore, we show that the M2 protein can drive plasma cell differentiation in a B lymphoma cell line in the absence of any other MHV68 gene products. Thus, the role of M2 in MHV68 reactivation can be attributed to its ability to manipulate plasma cell differentiation, providing a novel viral strategy to regulate gammaherpesvirus reactivation from latently infected B cells. We postulate that M2 represents a new class of herpesvirus gene products (reactivation conditioners) that do not directly participate in virus replication, but rather facilitate virus reactivation by manipulating the cellular milieu to provide a reactivation competent environment

The tuberculosis necrotizing toxin kills macrophages by hydrolyzing NAD.

Mycobacterium tuberculosis (Mtb) induces necrosis of infected cells to evade immune responses. Recently, we found that Mtb uses the protein CpnT to kill human macrophages by secreting its C-terminal domain, named tuberculosis necrotizing toxin (TNT), which induces necrosis by an unknown mechanism. Here we show that TNT gains access to the cytosol of Mtb-infected macrophages, where it hydrolyzes the essential coenzyme NAD(+). Expression or injection of a noncatalytic TNT mutant showed no cytotoxicity in macrophages or in zebrafish zygotes, respectively, thus demonstrating that the NAD(+) glycohydrolase activity is required for TNT-induced cell death. To prevent self-poisoning, Mtb produces an immunity factor for TNT (IFT) that binds TNT and inhibits its activity. The crystal structure of the TNT-IFT complex revealed a new NAD(+) glycohydrolase fold of TNT, the founding member of a toxin family widespread in pathogenic microorganisms

Proctitis following stereotactic body radiation therapy for prostate cancer

Background Proctitis after radiation therapy for prostate cancer remains an ongoing clinical challenge and critical quality of life issue. SBRT could minimize rectal toxicity by reducing the volume of rectum receiving high radiation doses and offers the potential radiobiologic benefits of hypofractionation. This study sought to evaluate the incidence and severity of proctitis following SBRT for prostate cancer. Methods Between February 2008 and July 2011, 269 men with clinically localized prostate cancer were treated definitively with SBRT monotherapy at Georgetown University Hospital. All patients were treated to 35-36.25Gy in 5 fractions delivered with the CyberKnife Radiosurgical System (Accuray). Rectal bleeding was recorded and scored using the CTCAE v.4. Telangiectasias were graded using the Vienna Rectoscopy Score (VRS). Proctitis was assessed via the Bowel domain of the Expanded Prostate Index Composite (EPIC)-26 at baseline and at 1, 3, 6, 9, 12, 18 and 24 months post-SBRT. Results The median age was 69 years with a median prostate volume of 39 cc. The median follow-up was 3.9 years with a minimum follow-up of two years. The 2-year actuarial incidence of late rectal bleeding ≥ grade 2 was 1.5%. Endoscopy revealed VRS Grade 2 rectal telangiectasias in 11% of patients. All proctitis symptoms increased at one month post-SBRT but returned to near-baseline with longer follow-up. The most bothersome symptoms were bowel urgency and frequency. At one month post-SBRT, 11.2% and 8.5% of patients reported a moderate to big problem with bowel urgency and frequency, respectively. The EPIC bowel summary scores declined transiently at 1 month and experienced a second, more protracted decline between 6 months and 18 months before returning to near-baseline at two years post-SBRT. Prior to treatment, 4.1% of men felt their bowel function was a moderate to big problem which increased to 11.5% one month post-SBRT but returned to near-baseline at two years post-SBRT. Conclusions In this single institution cohort, the rate and severity of proctitis observed following SBRT is low. QOL decreased on follow-up; however, our results compare favorably to those reported for patients treated with alternative radiation modalities. Future prospective randomized studies are needed to confirm these observations