Behavioral strategy shapes activation of the Vip-Sst disinhibitory circuit in visual cortex.

In complex environments, animals can adopt diverse strategies to find rewards. How distinct strategies differentially engage brain circuits is not well understood. Here, we investigate this question, focusing on the cortical Vip-Sst disinhibitory circuit between vasoactive intestinal peptide-postive (Vip) interneurons and somatostatin-positive (Sst) interneurons. We characterize the behavioral strategies used by mice during a visual change detection task. Using a dynamic logistic regression model, we find that individual mice use mixtures of a visual comparison strategy and a statistical timing strategy. Separately, mice also have periods of task engagement and disengagement. Two-photon calcium imaging shows large strategy-dependent differences in neural activity in excitatory, Sst inhibitory, and Vip inhibitory cells in response to both image changes and image omissions. In contrast, task engagement has limited effects on neural population activity. We find that the diversity of neural correlates of strategy can be understood parsimoniously as the increased activation of the Vip-Sst disinhibitory circuit during the visual comparison strategy, which facilitates task-appropriate responses.

Uncovering circuit mechanisms of current sinks and sources with biophysical simulations of primary visual cortex.

Local field potential (LFP) recordings reflect the dynamics of the current source density (CSD) in brain tissue. The synaptic, cellular, and circuit contributions to current sinks and sources are ill-understood. We investigated these in mouse primary visual cortex using public Neuropixels recordings and a detailed circuit model based on simulating the Hodgkin-Huxley dynamics of >50,000 neurons belonging to 17 cell types. The model simultaneously captured spiking and CSD responses and demonstrated a two-way dissociation: firing rates are altered with minor effects on the CSD pattern by adjusting synaptic weights, and CSD is altered with minor effects on firing rates by adjusting synaptic placement on the dendrites. We describe how thalamocortical inputs and recurrent connections sculpt specific sinks and sources early in the visual response, whereas cortical feedback crucially alters them in later stages. These results establish quantitative links between macroscopic brain measurements (LFP/CSD) and microscopic biophysics-based understanding of neuron dynamics and show that CSD analysis provides powerful constraints for modeling beyond those from considering spikes.

Simulations of cortical networks using spatially extended conductance-based neuronal models.

The Hodgkin-Huxley model of action potential generation and propagation, published in the Journal of Physiology in 1952, initiated the field of biophysically detailed computational modelling in neuroscience, which has expanded to encompass a variety of species and components of the nervous system. Here we review the developments in this area with a focus on efforts in the community towards modelling the mammalian neocortex using spatially extended conductance-based neuronal models. The Hodgkin-Huxley formalism and related foundational contributions, such as Rall's cable theory, remain widely used in these efforts to the current day. We argue that at present the field is undergoing a qualitative change due to new very rich datasets describing the composition, connectivity and functional activity of cortical circuits, which are being integrated systematically into large-scale network models. This trend, combined with the accelerating development of convenient software tools supporting such complex modelling projects, is giving rise to highly detailed models of the cortex that are extensively constrained by the data, enabling computational investigation of a multitude of questions about cortical structure and function.

EGFR oligomerization organizes kinase-active dimers into competent signalling platforms.

Epidermal growth factor receptor (EGFR) signalling is activated by ligand-induced receptor dimerization. Notably, ligand binding also induces EGFR oligomerization, but the structures and functions of the oligomers are poorly understood. Here, we use fluorophore localization imaging with photobleaching to probe the structure of EGFR oligomers. We find that at physiological epidermal growth factor (EGF) concentrations, EGFR assembles into oligomers, as indicated by pairwise distances of receptor-bound fluorophore-conjugated EGF ligands. The pairwise ligand distances correspond well with the predictions of our structural model of the oligomers constructed from molecular dynamics simulations. The model suggests that oligomerization is mediated extracellularly by unoccupied ligand-binding sites and that oligomerization organizes kinase-active dimers in ways optimal for auto-phosphorylation in trans between neighbouring dimers. We argue that ligand-induced oligomerization is essential to the regulation of EGFR signalling.

Membrane interaction of bound ligands contributes to the negative binding cooperativity of the EGF receptor.

The epidermal growth factor receptor (EGFR) plays a key role in regulating cell proliferation, migration, and differentiation, and aberrant EGFR signaling is implicated in a variety of cancers. EGFR signaling is triggered by extracellular ligand binding, which promotes EGFR dimerization and activation. Ligand-binding measurements are consistent with a negatively cooperative model in which the ligand-binding affinity at either binding site in an EGFR dimer is weaker when the other site is occupied by a ligand. This cooperativity is widely believed to be central to the effects of ligand concentration on EGFR-mediated intracellular signaling. Although the extracellular portion of the human EGFR dimer has been resolved crystallographically, the crystal structures do not reveal the structural origin of this negative cooperativity, which has remained unclear. Here we report the results of molecular dynamics simulations suggesting that asymmetrical interactions of the two binding sites with the membrane may be responsible (perhaps along with other factors) for this negative cooperativity. In particular, in our simulations the extracellular domains of an EGFR dimer spontaneously lay down on the membrane in an orientation in which favorable membrane contacts were made with one of the bound ligands, but could not be made with the other. Similar interactions were observed when EGFR was glycosylated, as it is in vivo.

Her2 activation mechanism reflects evolutionary preservation of asymmetric ectodomain dimers in the human EGFR family.

The receptor tyrosine kinase Her2, an intensely pursued drug target, differs from other members of the EGFR family in that it does not bind EGF-like ligands, relying instead on heterodimerization with other (ligand-bound) EGFR-family receptors for activation. The structural basis for Her2 heterodimerization, however, remains poorly understood. The unexpected recent finding of asymmetric ectodomain dimer structures of Drosophila EGFR (dEGFR) suggests a possible structural basis for Her2 heterodimerization, but all available structures for dimers of human EGFR family ectodomains are symmetric. Here, we report results from long-timescale molecular dynamics simulations indicating that a single ligand is necessary and sufficient to stabilize the ectodomain interface of Her2 heterodimers, which assume an asymmetric conformation similar to that of dEGFR dimers. This structural parallelism suggests a dimerization mechanism that has been conserved in the evolution of the EGFR family from Drosophila to human. DOI:http://dx.doi.org/10.7554/eLife.00708.001.

Transitions to catalytically inactive conformations in EGFR kinase.

The epidermal growth factor receptor (EGFR) is a key protein in cellular signaling, and its kinase domain (EGFR kinase) is an intensely pursued target of small-molecule drugs. Although both catalytically active and inactive conformations of EGFR kinase have been resolved crystallographically, experimental characterization of the transitions between these conformations remains difficult. Using unbiased, all-atom molecular dynamics simulations, we observed EGFR kinase spontaneously transition from the active to the so-called "Src-like inactive" conformation by way of two sets of intermediate conformations: One corresponds to a previously identified locally disordered state and the other to previously undescribed "extended" conformations, marked by the opening of the ATP-binding site between the two lobes of the kinase domain. We also simulated the protonation-dependent transition of EGFR kinase to another ["Asp-Phe-Gly-out" ("DFG-out")] inactive conformation and observed similar intermediate conformations. A key element observed in the simulated transitions is local unfolding, or "cracking," which supports a prediction of energy landscape theory. We used hydrogen-deuterium (H/D) exchange measurements to corroborate our simulations and found that the simulated intermediate conformations correlate better with the H/D exchange data than existing active or inactive EGFR kinase crystal structures. The intermediate conformations revealed by our simulations of the transition process differ significantly from the existing crystal structures and may provide unique possibilities for structure-based drug discovery.

Conformational coupling across the plasma membrane in activation of the EGF receptor.

How the epidermal growth factor receptor (EGFR) activates is incompletely understood. The intracellular portion of the receptor is intrinsically active in solution, and to study its regulation, we measured autophosphorylation as a function of EGFR surface density in cells. Without EGF, intact EGFR escapes inhibition only at high surface densities. Although the transmembrane helix and the intracellular module together suffice for constitutive activity even at low densities, the intracellular module is inactivated when tethered on its own to the plasma membrane, and fluorescence cross-correlation shows that it fails to dimerize. NMR and functional data indicate that activation requires an N-terminal interaction between the transmembrane helices, which promotes an antiparallel interaction between juxtamembrane segments and release of inhibition by the membrane. We conclude that EGF binding removes steric constraints in the extracellular module, promoting activation through N-terminal association of the transmembrane helices.

Architecture and membrane interactions of the EGF receptor.

Dimerization-driven activation of the intracellular kinase domains of the epidermal growth factor receptor (EGFR) upon extracellular ligand binding is crucial to cellular pathways regulating proliferation, migration, and differentiation. Inactive EGFR can exist as both monomers and dimers, suggesting that the mechanism regulating EGFR activity may be subtle. The membrane itself may play a role but creates substantial difficulties for structural studies. Our molecular dynamics simulations of membrane-embedded EGFR suggest that, in ligand-bound dimers, the extracellular domains assume conformations favoring dimerization of the transmembrane helices near their N termini, dimerization of the juxtamembrane segments, and formation of asymmetric (active) kinase dimers. In ligand-free dimers, by holding apart the N termini of the transmembrane helices, the extracellular domains instead favor C-terminal dimerization of the transmembrane helices, juxtamembrane segment dissociation and membrane burial, and formation of symmetric (inactive) kinase dimers. Electrostatic interactions of EGFR's intracellular module with the membrane are critical in maintaining this coupling.

Oncogenic mutations counteract intrinsic disorder in the EGFR kinase and promote receptor dimerization.

The mutation and overexpression of the epidermal growth factor receptor (EGFR) are associated with the development of a variety of cancers, making this prototypical dimerization-activated receptor tyrosine kinase a prominent target of cancer drugs. Using long-timescale molecular dynamics simulations, we find that the N lobe dimerization interface of the wild-type EGFR kinase domain is intrinsically disordered and that it becomes ordered only upon dimerization. Our simulations suggest, moreover, that some cancer-linked mutations distal to the dimerization interface, particularly the widespread L834R mutation (also referred to as L858R), facilitate EGFR dimerization by suppressing this local disorder. Corroborating these findings, our biophysical experiments and kinase enzymatic assays indicate that the L834R mutation causes abnormally high activity primarily by promoting EGFR dimerization rather than by allowing activation without dimerization. We also find that phosphorylation of EGFR kinase domain at Tyr845 may suppress the intrinsic disorder, suggesting a molecular mechanism for autonomous EGFR signaling.

Molecular dynamics simulations of the complete satellite tobacco mosaic virus.

This work presents an all-atom molecular dynamics simulation of a complete virus, the satellite tobacco mosaic virus. Simulations with up to 1 million atoms for over 50 ns demonstrate the stability of the entire virion and of the RNA core alone, while the capsid without RNA exhibits a pronounced instability. Physical properties of the simulated virus particle including electrostatic potential, radial distribution of viral components, and patterns of correlated motion are analyzed, and the implications for the assembly and infection mechanism of the virus are discussed.