Global parameter search reveals design principles of the mammalian circadian clock

Background: Virtually all living organisms have evolved a circadian (~24 hour) clock that controls physiological and behavioural processes with exquisite precision throughout the day/night cycle. The suprachiasmatic nucleus (SCN), which generates these ~24 h rhythms in mammals, consists of several thousand neurons. Each neuron contains a gene-regulatory network generating molecular oscillations, and the individual neuron oscillations are synchronised by intercellular coupling, presumably via neurotransmitters. Although this basic mechanism is currently accepted and has been recapitulated in mathematical models, several fundamental questions about the design principles of the SCN remain little understood. For example, a remarkable property of the SCN is that the phase of the SCN rhythm resets rapidly after a 'jet lag' type experiment, i.e. when the light/ dark (LD) cycle is abruptly advanced or delayed by several hours. Results: Here, we describe an extensive parameter optimization of a previously constructed simplified model of the SCN in order to further understand its design principles. By examining the top 50 solutions from the parameter optimization, we show that the neurotransmitters' role in generating the molecular circadian rhythms is extremely important. In addition, we show that when a neurotransmitter drives the rhythm of a system of coupled damped oscillators, it exhibits very robust synchronization and is much more easily entrained to light/dark cycles. We were also able to recreate in our simulations the fast rhythm resetting seen after a 'jet lag' type experiment. Conclusion: Our work shows that a careful exploration of parameter space for even an extremely simplified model of the mammalian clock can reveal unexpected behaviours and non-trivial predictions. Our results suggest that the neurotransmitter feedback loop plays a crucial role in the robustness and phase resetting properties of the mammalian clock, even at the single neuron level

A Generic Program for Multistate Protein Design

Some protein design tasks cannot be modeled by the traditional single state design strategy of finding a sequence that is optimal for a single fixed backbone. Such cases require multistate design, where a single sequence is threaded onto multiple backbones (states) and evaluated for its strengths and weaknesses on each backbone. For example, to design a protein that can switch between two specific conformations, it is necessary to to find a sequence that is compatible with both backbone conformations. We present in this paper a generic implementation of multistate design that is suited for a wide range of protein design tasks and demonstrate in silico its capabilities at two design tasks: one of redesigning an obligate homodimer into an obligate heterodimer such that the new monomers would not homodimerize, and one of redesigning a promiscuous interface to bind to only a single partner and to no longer bind the rest of its partners. Both tasks contained negative design in that multistate design was asked to find sequences that would produce high energies for several of the states being modeled. Success at negative design was assessed by computationally redocking the undesired protein-pair interactions; we found that multistate design's accuracy improved as the diversity of conformations for the undesired protein-pair interactions increased. The paper concludes with a discussion of the pitfalls of negative design, which has proven considerably more challenging than positive design

Reversed halo sign in pneumocystis pneumonia: a case report

<p>Abstract</p> <p>Background</p> <p>The reversed halo sign may sometimes be seen in patients with cryptogenic organizing pneumonia, but is rarely associated with other diseases.</p> <p>Case presentation</p> <p>We present a case study of a 32-year-old male patient with acquired immunodeficiency syndrome, who had previously been treated with chemotherapy for non-Hodgkin's lymphoma. A chest X-ray showed bilateral patchy infiltrates. High-resolution computed tomography revealed the reversed halo sign in both upper lobes. The patient was diagnosed with pneumocystis pneumonia, which was successfully treated with sulfamethoxazole trimethoprim; the reversed halo sign disappeared, leaving cystic lesions. Cases such as this one are rare, but show that the reversed halo sign may occur in patients who do not have cryptogenic organizing pneumonia.</p> <p>Conclusion</p> <p>Physicians can avoid making an incorrect diagnosis and prescribing the wrong treatment by carefully evaluating all clinical criteria rather than assuming that the reversed halo sign only occurs with cryptogenic organizing pneumonia.</p

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

Anchored Design of Protein-Protein Interfaces

Few existing protein-protein interface design methods allow for extensive backbone rearrangements during the design process. There is also a dichotomy between redesign methods, which take advantage of the native interface, and de novo methods, which produce novel binders.Here, we propose a new method for designing novel protein reagents that combines advantages of redesign and de novo methods and allows for extensive backbone motion. This method requires a bound structure of a target and one of its natural binding partners. A key interaction in this interface, the anchor, is computationally grafted out of the partner and into a surface loop on the design scaffold. The design scaffold's surface is then redesigned with backbone flexibility to create a new binding partner for the target. Careful choice of a scaffold will bring experimentally desirable characteristics into the new complex. The use of an anchor both expedites the design process and ensures that binding proceeds against a known location on the target. The use of surface loops on the scaffold allows for flexible-backbone redesign to properly search conformational space.This protocol was implemented within the Rosetta3 software suite. To demonstrate and evaluate this protocol, we have developed a benchmarking set of structures from the PDB with loop-mediated interfaces. This protocol can recover the correct loop-mediated interface in 15 out of 16 tested structures, using only a single residue as an anchor

RosettaRemodel: A Generalized Framework for Flexible Backbone Protein Design

We describe RosettaRemodel, a generalized framework for flexible protein design that provides a versatile and convenient interface to the Rosetta modeling suite. RosettaRemodel employs a unified interface, called a blueprint, which allows detailed control over many aspects of flexible backbone protein design calculations. RosettaRemodel allows the construction and elaboration of customized protocols for a wide range of design problems ranging from loop insertion and deletion, disulfide engineering, domain assembly, loop remodeling, motif grafting, symmetrical units, to de novo structure modeling

Rotational strength, range of motion, and function in people with unaffected shoulders from various stages of life

Background: Different measurements are used to assess shoulder function, including range of motion, strength, functional performance and self-report function. To understand disablement, it is necessary to understand the relationship between impairments and function in persons without shoulder problems. This study was conducted to enhance existing comparative data in subjects without upper extremity pathology, and to assess the relationships between impairments (range of motion, strength) and selfreported or measured function/disability. The impact of age, gender and dominance was determined. Methods: Two-hundred ninety-four subjects with unaffected shoulders were recruited. The subjects (mean age: 37 years old) were divided into three subgroups, 18–39, 40–59, and over 60 years of age. During a single session, at least two of the following variables were measured: self-reported function (shoulder disability scales), range of motion, isometric rotational strength, or upper limb functional performance (FIT-HaNSA). Two-way analysis of variance was used to determine, for each variable, the effects of age and gender. The relationship between the outcomes was established using Pearson product correlations. Results: Men were significantly stronger than women for all age categories. There was an age-related decline in strength in men in the over-60 age category. Significant negative correlations between strength and range of motion were demonstrated (-0.22 \u3c -0.32). Women had a significantly higher range of motion than men for external rotation in the 40–59 age category. Furthermore, the subjects in the over- 60 age category experienced a decrease of range of motion. There was minimal disability reported in all age groups on self-report scales. Only the Simple Shoulder Test demonstrated significant decreases in the over-60 age category and correlated with age (r = -0.202). Conclusion: Self-reported disability was low in individuals without upper extremity problems, although recruitment of such individuals was difficult in the older age groups due to the high prevalence of shoulder pathology. A low correlation between self-report disability and strength/range of motion in these unaffected subjects reflects the lack of disability reported by all subjects without pathology despite normal variations in strength and motion

Atomic-accuracy prediction of protein loop structures through an RNA-inspired ansatz

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. This article introduces a modeling strategy based on a 'stepwise ansatz', recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth RNA-puzzle competition. These results establish all-atom enumeration as a systematic approach to protein structure that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic resolution.Comment: Identity of four-loop blind test protein and parts of figures 5 have been omitted in this preprint to ensure confidentiality of the protein structure prior to its public releas

Constant light enhances synchrony among circadian clock cells and promotes behavioral rhythms in VPAC(2)-signaling deficient mice

Individual neurons in the suprachiasmatic nuclei (SCN) contain an intracellular molecular clock and use intercellular signaling to synchronize their timekeeping activities so that the SCN can coordinate brain physiology and behavior. The neuropeptide vasoactive intestinal polypeptide (VIP) and its VPAC2 receptor form a key component of intercellular signaling systems in the SCN and critically control cellular coupling. Targeted mutations in either the intracellular clock or intercellular neuropeptide signaling mechanisms, such as VIP-VPAC2 signaling, can lead to desynchronization of SCN neuronal clocks and loss of behavioral rhythms. An important goal in chronobiology is to develop interventions to correct deficiencies in circadian timekeeping. Here we show that extended exposure to constant light promotes synchrony among SCN clock cells and the expression of ~24 h rhythms in behavior in mice in which intercellular signaling is disrupted through loss of VIP-VPAC2 signaling. This study highlights the importance of SCN synchrony for the expression of rhythms in behavior and reveals how non-invasive manipulations in the external environment can be used to overcome neurochemical communication deficits in this important brain system

Maturation of GABAergic Inhibition Promotes Strengthening of Temporally Coherent Inputs among Convergent Pathways

Spike-timing-dependent plasticity (STDP), a form of Hebbian plasticity, is inherently stabilizing. Whether and how GABAergic inhibition influences STDP is not well understood. Using a model neuron driven by converging inputs modifiable by STDP, we determined that a sufficient level of inhibition was critical to ensure that temporal coherence (correlation among presynaptic spike times) of synaptic inputs, rather than initial strength or number of inputs within a pathway, controlled postsynaptic spike timing. Inhibition exerted this effect by preferentially reducing synaptic efficacy, the ability of inputs to evoke postsynaptic action potentials, of the less coherent inputs. In visual cortical slices, inhibition potently reduced synaptic efficacy at ages during but not before the critical period of ocular dominance (OD) plasticity. Whole-cell recordings revealed that the amplitude of unitary IPSCs from parvalbumin positive (Pv+) interneurons to pyramidal neurons increased during the critical period, while the synaptic decay time-constant decreased. In addition, intrinsic properties of Pv+ interneurons matured, resulting in an increase in instantaneous firing rate. Our results suggest that maturation of inhibition in visual cortex ensures that the temporally coherent inputs (e.g. those from the open eye during monocular deprivation) control postsynaptic spike times of binocular neurons, a prerequisite for Hebbian mechanisms to induce OD plasticity