RosettaScripts: A Scripting Language Interface to the Rosetta Macromolecular Modeling Suite

Macromolecular modeling and design are increasingly useful in basic research, biotechnology, and teaching. However, the absence of a user-friendly modeling framework that provides access to a wide range of modeling capabilities is hampering the wider adoption of computational methods by non-experts. RosettaScripts is an XML-like language for specifying modeling tasks in the Rosetta framework. RosettaScripts provides access to protocol-level functionalities, such as rigid-body docking and sequence redesign, and allows fast testing and deployment of complex protocols without need for modifying or recompiling the underlying C++ code. We illustrate these capabilities with RosettaScripts protocols for the stabilization of proteins, the generation of computationally constrained libraries for experimental selection of higher-affinity binding proteins, loop remodeling, small-molecule ligand docking, design of ligand-binding proteins, and specificity redesign in DNA-binding proteins

Characterization of a Novel Interaction between Bcl-2 Members Diva and Harakiri

Interactions within proteins of the Bcl-2 family are key in the regulation of apoptosis. The death-inducing members control apoptotic mechanisms partly by antagonizing the prosurvival proteins through heterodimer formation. Structural and biophysical studies on these complexes are providing important clues to understand their function. To help improve our knowledge on protein-protein interactions within the Bcl-2 family we have studied the binding between two of its members: mouse Diva and human Harakiri. Diva has been shown to perform both prosurvival and killing activity. In contrast, Harakiri induces cell death by interacting with antiapoptotic Bcl-2 members. Here we show using ELISA and NMR that Diva and Harakiri can interact in vitro. Combining the NMR data with the previously reported three-dimensional structure of Diva we find that Harakiri binds to a specific region in Diva. This interacting surface is equivalent to the known binding area of prosurvival Bcl-2 members from the reported structures of the complexes, suggesting that Diva could function at the structural level similarly to the antiapoptotic proteins of the Bcl-2 family. We illustrate this result by building a structural model of the heterodimer using molecular docking and the NMR data as restraints. Moreover, combining circular dichroism and NMR we also show that Harakiri is largely unstructured with residual (13%) α-helical conformation. This result agrees with intrinsic disorder previously observed in other Bcl-2 members. In addition, Harakiri constructs of different length were studied to identify the region critical for the interaction. Differential affinity for Diva of these constructs suggests that the amino acid sequence flanking the interacting region could play an important role in binding

Intrinsic Order and Disorder in the Bcl-2 Member Harakiri: Insights into Its Proapoptotic Activity

Harakiri is a BH3-only member of the Bcl-2 family that localizes in membranes and induces cell death by binding to prosurvival Bcl-xL and Bcl-2. The cytosolic domain of Harakiri is largely disorder with residual α-helical conformation according to previous structural studies. As these helical structures could play an important role in Harakiri's function, we have used NMR and circular dichroism to fully characterize them at the residue-atomic level. In addition, we report structural studies on a peptide fragment spanning Harakiri's C-terminal hydrophobic sequence, which potentially operates as a transmembrane domain. We initially checked by enzyme immunoassays and NMR that peptides encompassing different lengths of the cytosolic domain are functional as they bind Bcl-xL and Bcl-2. The structural data in water indicate that the α-helical conformation is restricted to a 25-residue segment comprising the BH3 domain. However, structure calculation was precluded because of insufficient NMR restraints. To bypass this problem we used alcohol-water mixture to increase structure population and confirmed by NMR that the conformation in both milieus is equivalent. The resulting three-dimensional structure closely resembles that of peptides encompassing the BH3 domain of BH3-only members in complex with their prosurvival partners, suggesting that preformed structural elements in the disordered protein are central to binding. In contrast, the transmembrane domain forms in micelles a monomeric α-helix with a population close to 100%. Its three-dimensional structure here reported reveals features that explain its function as membrane anchor. Altogether these results are used to propose a tentative structural model of how Harakiri works

Common Genetic Polymorphisms Influence Blood Biomarker Measurements in COPD

Implementing precision medicine for complex diseases such as chronic obstructive lung disease (COPD) will require extensive use of biomarkers and an in-depth understanding of how genetic, epigenetic, and environmental variations contribute to phenotypic diversity and disease progression. A meta-analysis from two large cohorts of current and former smokers with and without COPD [SPIROMICS (N = 750); COPDGene (N = 590)] was used to identify single nucleotide polymorphisms (SNPs) associated with measurement of 88 blood proteins (protein quantitative trait loci; pQTLs). PQTLs consistently replicated between the two cohorts. Features of pQTLs were compared to previously reported expression QTLs (eQTLs). Inference of causal relations of pQTL genotypes, biomarker measurements, and four clinical COPD phenotypes (airflow obstruction, emphysema, exacerbation history, and chronic bronchitis) were explored using conditional independence tests. We identified 527 highly significant (p 10% of measured variation in 13 protein biomarkers, with a single SNP (rs7041; p = 10−392) explaining 71%-75% of the measured variation in vitamin D binding protein (gene = GC). Some of these pQTLs [e.g., pQTLs for VDBP, sRAGE (gene = AGER), surfactant protein D (gene = SFTPD), and TNFRSF10C] have been previously associated with COPD phenotypes. Most pQTLs were local (cis), but distant (trans) pQTL SNPs in the ABO blood group locus were the top pQTL SNPs for five proteins. The inclusion of pQTL SNPs improved the clinical predictive value for the established association of sRAGE and emphysema, and the explanation of variance (R2) for emphysema improved from 0.3 to 0.4 when the pQTL SNP was included in the model along with clinical covariates. Causal modeling provided insight into specific pQTL-disease relationships for airflow obstruction and emphysema. In conclusion, given the frequency of highly significant local pQTLs, the large amount of variance potentially explained by pQTL, and the differences observed between pQTLs and eQTLs SNPs, we recommend that protein biomarker-disease association studies take into account the potential effect of common local SNPs and that pQTLs be integrated along with eQTLs to uncover disease mechanisms. Large-scale blood biomarker studies would also benefit from close attention to the ABO blood group

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk.

Blood pressure is a heritable trait influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (≥140 mm Hg systolic blood pressure or  ≥90 mm Hg diastolic blood pressure). Even small increments in blood pressure are associated with an increased risk of cardiovascular events. This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.Peer reviewe

31st Annual Meeting and Associated Programs of the Society for Immunotherapy of Cancer (SITC 2016) : part two

Background The immunological escape of tumors represents one of the main ob- stacles to the treatment of malignancies. The blockade of PD-1 or CTLA-4 receptors represented a milestone in the history of immunotherapy. However, immune checkpoint inhibitors seem to be effective in specific cohorts of patients. It has been proposed that their efficacy relies on the presence of an immunological response. Thus, we hypothesized that disruption of the PD-L1/PD-1 axis would synergize with our oncolytic vaccine platform PeptiCRAd. Methods We used murine B16OVA in vivo tumor models and flow cytometry analysis to investigate the immunological background. Results First, we found that high-burden B16OVA tumors were refractory to combination immunotherapy. However, with a more aggressive schedule, tumors with a lower burden were more susceptible to the combination of PeptiCRAd and PD-L1 blockade. The therapy signifi- cantly increased the median survival of mice (Fig. 7). Interestingly, the reduced growth of contralaterally injected B16F10 cells sug- gested the presence of a long lasting immunological memory also against non-targeted antigens. Concerning the functional state of tumor infiltrating lymphocytes (TILs), we found that all the immune therapies would enhance the percentage of activated (PD-1pos TIM- 3neg) T lymphocytes and reduce the amount of exhausted (PD-1pos TIM-3pos) cells compared to placebo. As expected, we found that PeptiCRAd monotherapy could increase the number of antigen spe- cific CD8+ T cells compared to other treatments. However, only the combination with PD-L1 blockade could significantly increase the ra- tio between activated and exhausted pentamer positive cells (p= 0.0058), suggesting that by disrupting the PD-1/PD-L1 axis we could decrease the amount of dysfunctional antigen specific T cells. We ob- served that the anatomical location deeply influenced the state of CD4+ and CD8+ T lymphocytes. In fact, TIM-3 expression was in- creased by 2 fold on TILs compared to splenic and lymphoid T cells. In the CD8+ compartment, the expression of PD-1 on the surface seemed to be restricted to the tumor micro-environment, while CD4 + T cells had a high expression of PD-1 also in lymphoid organs. Interestingly, we found that the levels of PD-1 were significantly higher on CD8+ T cells than on CD4+ T cells into the tumor micro- environment (p < 0.0001). Conclusions In conclusion, we demonstrated that the efficacy of immune check- point inhibitors might be strongly enhanced by their combination with cancer vaccines. PeptiCRAd was able to increase the number of antigen-specific T cells and PD-L1 blockade prevented their exhaus- tion, resulting in long-lasting immunological memory and increased median survival