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

In Vivo CD8+ T-Cell Suppression of SIV Viremia Is Not Mediated by CTL Clearance of Productively Infected Cells

The CD8+ T-cell is a key mediator of antiviral immunity, potentially contributing to control of pathogenic lentiviral infection through both innate and adaptive mechanisms. We studied viral dynamics during antiretroviral treatment of simian immunodeficiency virus (SIV) infected rhesus macaques following CD8+ T-cell depletion to test the importance of adaptive cytotoxic effects in clearance of cells productively infected with SIV. As previously described, plasma viral load (VL) increased following CD8+ T-cell depletion and was proportional to the magnitude of CD8+ T-cell depletion in the GALT, confirming a direct relationship between CD8+ T-cell loss and viral replication. Surprisingly, first phase plasma virus decay following administration of antiretroviral drugs was not slower in CD8+ T-cell depleted animals compared with controls indicating that the short lifespan of the average productively infected cell is not a reflection of cytotoxic T-lymphocyte (CTL) killing. Our findings support a dominant role for non-cytotoxic effects of CD8+ T-cells on control of pathogenic lentiviral infection and suggest that cytotoxic effects, if present, are limited to early, pre-productive stages of the viral life cycle. These observations have important implications for future strategies to augment immune control of HIV

Persistence of viral reservoirs in multiple tissues after antiretroviral therapy suppression in a macaque RT-SHIV model

Although antiretroviral therapy (ART) can suppress HIV-1 replication sufficiently to eliminate measurable plasma viremia, infected cells remain and ensure viral recrudescence after discontinuation of ART. We used a macaque model of HIV-1/AIDS to evaluate the location of infected cells during ART. Twelve macaques were infected with RT-SHIVmne, a SIV containing HIV-1 reverse transcriptase, conferring sensitivity to non-nucleoside reverse transcriptase inhibitors (NNRTIs). Ten to fourteen weeks post-infection, 6 animals were treated with 3 or 4 antiretroviral drugs for 17-20 weeks; 6 control animals remained untreated. Viral DNA (vDNA) and RNA (vRNA) were measured in peripheral blood mononuclear cells (PBMC) and at necropsy in multiple tissues by quantitative PCR and RT-PCR. The majority of virally infected cells were located in lymphoid tissues with variable levels in the gastrointestinal tract of both treated and untreated animals. Tissue viral DNA levels correlated with week 1 plasma viremia, suggesting that tissues that harbor proviral DNA are established within the first week of infection. PBMC vDNA levels did not correlate with plasma viremia or tissue levels of vDNA. vRNA levels were high in lymphoid and gastrointestinal tissues of the untreated animals; animals on ART had little vRNA expressed in tissues and virus could not be cultured from lymph node resting CD4+ cells after 17-20 weeks on ART, indicating little or no ongoing viral replication. Strategies for eradication of HIV-1 will need to target residual virus in ART suppressed individuals, which may not be accurately reflected by frequencies of infected cells in blood. © 2013 Kline et al

Search for Gravitational Waves from Primordial Black Hole Binary Coalescences in the Galactic Halo

We use data from the second science run of the LIGO gravitational-wave detectors to search for the gravitational waves from primordial black hole (PBH) binary coalescence with component masses in the range 0.2--$1.0 M_\odot$. The analysis requires a signal to be found in the data from both LIGO observatories, according to a set of coincidence criteria. No inspiral signals were found. Assuming a spherical halo with core radius 5 kpc extending to 50 kpc containing non-spinning black holes with masses in the range 0.2--$1.0 M_\odot$, we place an observational upper limit on the rate of PBH coalescence of 63 per year per Milky Way halo (MWH) with 90% confidence.Comment: 7 pages, 4 figures, to be submitted to Phys. Rev.

Operationalizing frailty among older residents of assisted living facilities

<p>Abstract</p> <p>Background</p> <p>Frailty in later life is viewed as a state of heightened vulnerability to poor outcomes. The utility of frailty as a measure of vulnerability in the assisted living (AL) population remains unexplored. We examined the feasibility and predictive accuracy of two different interpretations of the Cardiovascular Health Study (CHS) frailty criteria in a population-based sample of AL residents.</p> <p>Methods</p> <p>CHS frailty criteria were operationalized using two different approaches in 928 AL residents from the Alberta Continuing Care Epidemiological Studies (ACCES). Risks of one-year mortality and hospitalization were estimated for those categorized as frail or pre-frail (compared with non-frail). The prognostic significance of individual criteria was explored, and the area under the ROC curve (AUC) was calculated for select models to assess the utility of frailty in predicting one-year outcomes.</p> <p>Results</p> <p>Regarding feasibility, complete CHS criteria could not be assessed for 40% of the initial 1,067 residents. Consideration of supplementary items for select criteria reduced this to 12%. Using absolute (CHS-specified) cut-points, 48% of residents were categorized as frail and were at greater risk for death (adjusted risk ratio [RR] 1.75, 95% CI 1.08-2.83) and hospitalization (adjusted RR 1.54, 95% CI 1.20-1.96). Pre-frail residents defined by absolute cut-points (48.6%) showed no increased risk for mortality or hospitalization compared with non-frail residents. Using relative cut-points (derived from AL sample), 19% were defined as frail and 55% as pre-frail and the associated risks for mortality and hospitalization varied by sex. Frail (but not pre-frail) women were more likely to die (RR 1.58 95% CI 1.02-2.44) and be hospitalized (RR 1.53 95% CI 1.25-1.87). Frail and pre-frail men showed an increased mortality risk (RR 3.21 95% CI 1.71-6.00 and RR 2.61 95% CI 1.40-4.85, respectively) while only pre-frail men had an increased risk of hospitalization (RR 1.58 95% CI 1.15-2.17). Although incorporating either frailty measure improved the performance of predictive models, the best AUCs were 0.702 for mortality and 0.633 for hospitalization.</p> <p>Conclusions</p> <p>Application of the CHS criteria for frailty was problematic and only marginally improved the prediction of select adverse outcomes in AL residents. Development and validation of alternative approaches for detecting frailty in this population, including consideration of female/male differences, is warranted.</p

Does weight loss improve semen quality and reproductive hormones? results from a cohort of severely obese men

<p>Abstract</p> <p>Background</p> <p>A high body mass index (BMI) has been associated with reduced semen quality and male subfecundity, but no studies following obese men losing weight have yet been published. We examined semen quality and reproductive hormones among morbidly obese men and studied if weight loss improved the reproductive indicators.</p> <p>Methods</p> <p>In this pilot cohort study, 43 men with BMI > 33 kg/m²were followed through a 14 week residential weight loss program. The participants provided semen samples and had blood samples drawn, filled in questionnaires, and had clinical examinations before and after the intervention. Conventional semen characteristics as well as sperm DNA integrity, analysed by the sperm chromatin structure assay (SCSA) were obtained. Serum levels of testosterone, estradiol, sex hormone-binding globulin (SHBG), luteinizing hormone (LH), follicle-stimulating hormone (FSH), anti-Müllerian hormone (AMH) and inhibin B (Inh-B) were measured.</p> <p>Results</p> <p>Participants were from 20 to 59 years of age (median = 32) with BMI ranging from 33 to 61 kg/m². At baseline, after adjustment for potential confounders, BMI was inversely associated with sperm concentration (p = 0.02), total sperm count (p = 0.02), sperm morphology (p = 0.04), and motile sperm (p = 0.005) as well as testosterone (p = 0.04) and Inh-B (p = 0.04) and positively associated to estradiol (p < 0.005). The median (range) percentage weight loss after the intervention was 15% (3.5 - 25.4). Weight loss was associated with an increase in total sperm count (p = 0.02), semen volume (p = 0.04), testosterone (p = 0.02), SHBG (p = 0.03) and AMH (p = 0.02). The group with the largest weight loss had a statistically significant increase in total sperm count [193 millions (95% CI: 45; 341)] and normal sperm morphology [4% (95% CI: 1; 7)].</p> <p>Conclusion</p> <p>This study found obesity to be associated with poor semen quality and altered reproductive hormonal profile. Weight loss may potentially lead to improvement in semen quality. Whether the improvement is a result of the reduction in body weight per se or improved lifestyles remains unknown.</p

Can Non-lytic CD8+T Cells Drive HIV-1 Escape?

The CD8+ T cell effector mechanisms that mediate control of HIV-1 and SIV infections remain poorly understood. Recent work suggests that the mechanism may be primarily non-lytic. This is in apparent conflict with the observation that SIV and HIV-1 variants that escape CD8+ T cell surveillance are frequently selected. Whilst it is clear that a variant that has escaped a lytic response can have a fitness advantage compared to the wild-type, it is less obvious that this holds in the face of non-lytic control where both wild-type and variant infected cells would be affected by soluble factors. In particular, the high motility of T cells in lymphoid tissue would be expected to rapidly destroy local effects making selection of escape variants by non-lytic responses unlikely. The observation of frequent HIV-1 and SIV escape poses a number of questions. Most importantly, is the consistent observation of viral escape proof that HIV-1- and SIV-specific CD8+ T cells lyse infected cells or can this also be the result of non-lytic control? Additionally, the rate at which a variant strain escapes a lytic CD8+ T cell response is related to the strength of the response. Is the same relationship true for a non-lytic response? Finally, the potential anti-viral control mediated by non-lytic mechanisms compared to lytic mechanisms is unknown. These questions cannot be addressed with current experimental techniques nor with the standard mathematical models. Instead we have developed a 3D cellular automaton model of HIV-1 which captures spatial and temporal dynamics. The model reproduces in vivo HIV-1 dynamics at the cellular and population level. Using this model we demonstrate that non-lytic effector mechanisms can select for escape variants but that outgrowth of the variant is slower and less frequent than from a lytic response so that non-lytic responses can potentially offer more durable control

The Timing of the Cognitive Cycle

We propose that human cognition consists of cascading cycles of recurring brain events. Each cognitive cycle senses the current situation, interprets it with reference to ongoing goals, and then selects an internal or external action in response. While most aspects of the cognitive cycle are unconscious, each cycle also yields a momentary “ignition” of conscious broadcasting. Neuroscientists have independently proposed ideas similar to the cognitive cycle, the fundamental hypothesis of the LIDA model of cognition. High-level cognition, such as deliberation, planning, etc., is typically enabled by multiple cognitive cycles. In this paper we describe a timing model LIDA's cognitive cycle. Based on empirical and simulation data we propose that an initial phase of perception (stimulus recognition) occurs 80–100 ms from stimulus onset under optimal conditions. It is followed by a conscious episode (broadcast) 200–280 ms after stimulus onset, and an action selection phase 60–110 ms from the start of the conscious phase. One cognitive cycle would therefore take 260–390 ms. The LIDA timing model is consistent with brain evidence indicating a fundamental role for a theta-gamma wave, spreading forward from sensory cortices to rostral corticothalamic regions. This posteriofrontal theta-gamma wave may be experienced as a conscious perceptual event starting at 200–280 ms post stimulus. The action selection component of the cycle is proposed to involve frontal, striatal and cerebellar regions. Thus the cycle is inherently recurrent, as the anatomy of the thalamocortical system suggests. The LIDA model fits a large body of cognitive and neuroscientific evidence. Finally, we describe two LIDA-based software agents: the LIDA Reaction Time agent that simulates human performance in a simple reaction time task, and the LIDA Allport agent which models phenomenal simultaneity within timeframes comparable to human subjects. While there are many models of reaction time performance, these results fall naturally out of a biologically and computationally plausible cognitive architecture