Predicting the behavior of heterogeneous nonequilibrium systems is currently analytically intractable. Consequently, complex biological systems have resisted unifying principles. Here, I introduce a ...mapping from dynamical systems to battery-resistor circuits. I show that in these transformed variables (i) arbitrary numbers of heterogeneous dynamical transitions can be reduced to a Thevenin equivalent resistor which is invariant to driving from equilibrium, (ii) resistors (together with the external driving sources) are sufficient to describe system behavior, and (iii) the resistor's directional symmetry leads to universal theorems of nonequilibrium behavior. This mapping is used to derive two general steady-state relations. First, for any cyclic process, the maximum amplification of any state is tightly bounded by the total dissipation of all states; experimental data are used to show that the master signal protein Ras achieves this bound. Second, for any process, the response of any reaction due to driving any other reaction is identical to the reciprocal response rescaled by the ratio of the corresponding Thevenin resistors. This result generalizes Onsager's reciprocal relation to the strongly driven regime and makes a testable prediction about how systems should be designed or evolved to maximize response. These analytic results represent a new perspective applicable to biological complexity and suggest that this mapping provides the natural variables to study heterogeneous nonequilibrium systems.
Tauopathies are neurodegenerative diseases characterized by intracellular amyloid deposits of tau protein. Missense mutations in the tau gene (MAPT) correlate with aggregation propensity and cause ...dominantly inherited tauopathies, but their biophysical mechanism driving amyloid formation is poorly understood. Many disease-associated mutations localize within tau's repeat domain at inter-repeat interfaces proximal to amyloidogenic sequences, such as
VQIVYK
. We use cross-linking mass spectrometry, recombinant protein and synthetic peptide systems, in silico modeling, and cell models to conclude that the aggregation-prone
VQIVYK
motif forms metastable compact structures with its upstream sequence that modulates aggregation propensity. We report that disease-associated mutations, isomerization of a critical proline, or alternative splicing are all sufficient to destabilize this local structure and trigger spontaneous aggregation. These findings provide a biophysical framework to explain the basis of early conformational changes that may underlie genetic and sporadic tau pathogenesis.
Protein fibril self-assembly is a universal transition implicated in neurodegenerative diseases. Although fibril structure/growth are well characterized, fibril nucleation is poorly understood. Here, ...we use a computational-experimental approach to resolve fibril nucleation. We show that monomer hairpin content quantified from molecular dynamics simulations is predictive of experimental fibril formation kinetics across a tau motif mutant library. Hairpin trimers are predicted to be fibril transition states; one hairpin spontaneously converts into the cross-beta conformation, templating subsequent fibril growth. We designed a disulfide-linked dimer mimicking the transition state that catalyzes fibril formation, measured by ThT fluorescence and TEM, of wild-type motif - which does not normally fibrillize. A dimer compatible with extended conformations but not the transition-state fails to nucleate fibril at any concentration. Tau repeat domain simulations show how long-range interactions sequester this motif in a mutation-dependent manner. This work implies that different fibril morphologies could arise from disease-dependent hairpin seeding from different loci.
The SNAP receptor (SNARE) proteins syntaxin-1, SNAP-25, and synaptobrevin mediate neurotransmitter release by forming tight SNARE complexes that fuse synaptic vesicles with the plasma membranes in ...microseconds. Membrane fusion is generally explained by the action of proteins on macroscopic membrane properties such as curvature, elastic modulus, and tension, and a widespread model envisions that the SNARE motifs, juxtamembrane linkers, and C-terminal transmembrane regions of synaptobrevin and syntaxin-1 form continuous helices that act mechanically as semirigid rods, squeezing the membranes together as they assemble ("zipper") from the N to the C termini. However, the mechanism underlying fast SNARE-induced membrane fusion remains unknown. We have used all-atom molecular dynamics simulations to investigate this mechanism. Our results need to be interpreted with caution because of the limited number and length of the simulations, but they suggest a model of membrane fusion that has a natural physicochemical basis, emphasizes local molecular events over general membrane properties, and explains extensive experimental data. In this model, the central event that initiates fast (microsecond scale) membrane fusion occurs when the SNARE helices zipper into the juxtamembrane linkers which, together with the adjacent transmembrane regions, promote encounters of acyl chains from both bilayers at the polar interface. The resulting hydrophobic nucleus rapidly expands into stalk-like structures that gradually progress to form a fusion pore, aided by the SNARE transmembrane regions and without clearly discernible intermediates. The propensity of polyunsaturated lipids to participate in encounters that initiate fusion suggests that these lipids may be important for the high speed of neurotransmitter release.
ABSTRACT
Fungi play many essential roles in ecosystems. They facilitate plant access to nutrients and water, serve as decay agents that cycle carbon and nutrients through the soil, water and ...atmosphere, and are major regulators of macro‐organismal populations. Although technological advances are improving the detection and identification of fungi, there still exist key gaps in our ecological knowledge of this kingdom, especially related to function. Trait‐based approaches have been instrumental in strengthening our understanding of plant functional ecology and, as such, provide excellent models for deepening our understanding of fungal functional ecology in ways that complement insights gained from traditional and ‐omics‐based techniques. In this review, we synthesize current knowledge of fungal functional ecology, taxonomy and systematics and introduce a novel database of fungal functional traits (FunFun). FunFun is built to interface with other databases to explore and predict how fungal functional diversity varies by taxonomy, guild, and other evolutionary or ecological grouping variables. To highlight how a quantitative trait‐based approach can provide new insights, we describe multiple targeted examples and end by suggesting next steps in the rapidly growing field of fungal functional ecology.
To investigate the concurrent validity and reliability of the Mini International Neuropsychiatric Interview for Children and Adolescents (MINI-KID), a short structured diagnostic interview for DSM-IV ...and ICD-10 psychiatric disorders in children and adolescents.
Participants were 226 children and adolescents (190 outpatients and 36 controls) aged 6 to 17 years. To assess the concurrent validity of the MINI-KID, participants were administered the MINI-KID and the Schedule for Affective Disorders and Schizophrenia for School Aged Children-Present and Lifetime Version (K-SADS-PL) by blinded interviewers in a counterbalanced order on the same day. Participants also completed a self-rated measure of disability. In addition, interrater (n = 57) and test-retest (n = 83) reliability data (retest interval, 1-5 days) were collected, and agreement between the parent version of the MINI-KID and the standard MINI-KID (n = 140) was assessed. Data were collected between March 2004 and January 2008.
Substantial to excellent MINI-KID to K-SADS-PL concordance was found for syndromal diagnoses of any mood disorder, any anxiety disorder, any substance use disorder, any ADHD or behavioral disorder, and any eating disorder (area under curve AUC = 0.81-0.96, kappa = 0.56-0.87). Results were more variable for psychotic disorder (AUC = 0.94, kappa = 0.41). Sensitivity was substantial (0.61-1.00) for 15/20 individual DSM-IV disorders. Specificity was excellent (0.81-1.00) for 18 disorders and substantial (> 0.73) for the remaining 2. The MINI-KID identified a median of 3 disorders per subject compared to 2 on the K-SADS-PL and took two-thirds less time to administer (34 vs 103 minutes). Interrater and test-retest kappas were substantial to almost perfect (0.64-1.00) for all individual MINI-KID disorders except dysthymia. Concordance of the parent version (MINI-KID-P) with the standard MINI-KID was good.
The MINI-KID generates reliable and valid psychiatric diagnoses for children and adolescents and does so in a third of the time as the K-SADS-PL.
How protein structure encodes functionality is not fully understood. For example, long-range intraprotein communication can occur without measurable conformational change and is often not captured by ...existing structural correlation functions. It is shown here that important functional information is encoded in the timing of protein motions, rather than motion itself. I introduce the conditional activity function to quantify such timing correlations among the degrees of freedom within proteins. For three proteins, the conditional activities between side-chain dihedral angles were computed using the output of microseconds-long atomistic simulations. The new approach demonstrates that a sparse fraction of side-chain pairs are dynamically correlated over long distances (spanning protein lengths up to 7 nm), in sharp contrast to structural correlations, which are short-ranged (<1 nm). Regions of high self- and inter-side-chain dynamical correlations are found, corresponding to experimentally determined functional modules and allosteric connections, respectively.
Allostery, the transfer of information between distant parts of a macromolecule, is a fundamental feature of protein function and regulation. However, allosteric mechanisms are usually not explained ...by protein structure, requiring information on correlated fluctuations uniquely accessible to molecular simulation. Existing work to extract allosteric pathways from molecular dynamics simulations has focused on thermodynamic correlations. Here, we show how kinetic correlations encode complementary information essential to explain observed variations in allosteric regulation. We applied kinetic and thermodynamic correlation analysis on atomistic simulations of H, K, and NRas isoforms in the apo, GTP, and GDP-bound states of Ras protein, with and without complexing to its downstream effector, Raf. We show that switch I and switch II are the primary components of thermodynamic and kinetic allosteric networks, consistent with the key roles of these two motifs. These networks connect the switches to an allosteric loop recently discovered from a crystal structure of HRas. This allosteric loop is inactive in KRas, but is coupled to the hydrolysis arm switch II in NRas and HRas. We find that the mechanism in the latter two isoforms are thermodynamic and kinetic, respectively. Binding of Raf-RBD further activates thermodynamic allostery in HRas and KRas but has limited effect on NRas. These results indicate that kinetic and thermodynamic correlations are both needed to explain protein function and allostery. These two distinct channels of allosteric regulation, and their combinatorial variability, may explain how subtle mutational differences can lead to diverse regulatory profiles among enzymatic proteins.
Synaptic vesicles are primed into a state that is ready for fast neurotransmitter release upon Ca
-binding to Synaptotagmin-1. This state likely includes trans-SNARE complexes between the vesicle and ...plasma membranes that are bound to Synaptotagmin-1 and complexins. However, the nature of this state and the steps leading to membrane fusion are unclear, in part because of the difficulty of studying this dynamic process experimentally. To shed light into these questions, we performed all-atom molecular dynamics simulations of systems containing trans-SNARE complexes between two flat bilayers or a vesicle and a flat bilayer with or without fragments of Synaptotagmin-1 and/or complexin-1. Our results need to be interpreted with caution because of the limited simulation times and the absence of key components, but suggest mechanistic features that may control release and help visualize potential states of the primed Synaptotagmin-1-SNARE-complexin-1 complex. The simulations suggest that SNAREs alone induce formation of extended membrane-membrane contact interfaces that may fuse slowly, and that the primed state contains macromolecular assemblies of trans-SNARE complexes bound to the Synaptotagmin-1 C
B domain and complexin-1 in a spring-loaded configuration that prevents premature membrane merger and formation of extended interfaces, but keeps the system ready for fast fusion upon Ca
influx.
Proteins can aggregate into disordered aggregates or ordered assemblies such as amyloid fibrils. These two distinct phases serve differing roles in function and disease. How protein sequence ...determines the preferred phase is unknown. Here we establish a statistical mechanical disorder‐to‐order transition condition for compact polymer aggregates, including proteins. The theory produces a simple universal equation determining the favored phase as a function of temperature, polymer length, and interaction energy variance. We show that the sequence‐dependent energy variance is efficiently calculated using atomistic simulations, so that the theory has no adjustable parameters. The equation accurately predicts experimental length‐dependent crystallization temperatures of synthetic polymers. The equation also predicts that all protein sequences that aggregate will also favor ordering. Consequently, energy must be expended to maintain the steady‐state disordered phase if it is not kinetically metastable on physiological timescales. More broadly, the theory suggests that aggregates of organic polymers will generally tend to order on habitable planets.