Nanoactuator for Neuronal Optoporation

Light-driven modulation of neuronal activity at high spatial-temporal resolution is becoming of high interest in neuroscience. In addition to optogenetics, nongenetic membrane-targeted nanomachines that alter the electrical state of the neuronal membranes are in demand. Here, we engineered and characterized a photoswitchable conjugated compound (BV-1) that spontaneously partitions into the neuronal membrane and undergoes a charge transfer upon light stimulation. The activity of primary neurons is not affected in the dark, whereas millisecond light pulses of cyan light induce a progressive decrease in membrane resistance and an increase in inward current matched to a progressive depolarization and action potential firing. We found that illumination of BV-1 induces oxidation of membrane phospholipids, which is necessary for the electrophysiological effects and is associated with decreased membrane tension and increased membrane fluidity. Time-resolved atomic force microscopy and molecular dynamics simulations performed on planar lipid bilayers revealed that the underlying mechanism is a light-driven formation of pore-like structures across the plasma membrane. Such a phenomenon decreases membrane resistance and increases permeability to monovalent cations, namely, Na+, mimicking the effects of antifungal polyenes. The same effect on membrane resistance was also observed in nonexcitable cells. When sustained light stimulations are applied, neuronal swelling and death occur. The light-controlled pore-forming properties of BV-1 allow performing "on-demand" light-induced membrane poration to rapidly shift from cell-attached to perforated whole-cell patch-clamp configuration. Administration of BV-1 to ex vivo retinal explants or in vivo primary visual cortex elicited neuronal firing in response to short trains of light stimuli, followed by activity silencing upon prolonged light stimulations. BV-1 represents a versatile molecular nanomachine whose properties can be exploited to induce either photostimulation or space-specific cell death, depending on the pattern and duration of light stimulation.

Case Report: Novel biallelic moderately damaging variants in RTTN in a patient with cerebellar dysplasia

Rotatin, encoded by the RTTN gene, is a centrosomal protein with multiple, emerging functions, including left-right specification, ciliogenesis, and neuronal migration. Recessive variants in RTTN are associated with a neurodevelopmental disorder with microcephaly and malformations of cortical development known as "Microcephaly, short stature, and polymicrogyria with seizures" (MSSP, MIM #614833). Affected individuals show a wide spectrum of clinical manifestations like intellectual disability, poor/absent speech, short stature, microcephaly, and congenital malformations. Here, we report a subject showing a distinctive neuroradiological phenotype and harboring novel biallelic variants in RTTN: the c.5500A>G, p.(Asn1834Asp), (dbSNP: rs200169343, ClinVar ID:1438510) and c.19A>G, p.(Ile7Val), (dbSNP: rs201165599, ClinVar ID:1905275) variants. In particular brain magnetic resonance imaging (MRI) showed a peculiar pattern, with cerebellar hypo-dysplasia, and multiple arachnoid cysts in the lateral cerebello-medullary cisterns, in addition to left Meckel cave. Thus, we compare his phenotypic features with current literature, speculating a possible role of newly identified RTTN variants in his clinical picture, and supporting a relevant variability in this emerging condition.

Charting Nanocluster Structures via Convolutional Neural Networks

A general method to obtain a representation of the structural landscape of nanoparticles in terms of a limited number of variables is proposed. The method is applied to a large data set of parallel tempering molecular dynamics simulations of gold clusters of 90 and 147 atoms, silver clusters of 147 atoms, and copper clusters of 147 atoms, covering a plethora of structures and temperatures. The method leverages convolutional neural networks to learn the radial distribution functions of the nanoclusters and distills a low-dimensional chart of the structural landscape. This strategy is found to give rise to a physically meaningful and differentiable mapping of the atom positions to a low-dimensional manifold in which the main structural motifs are clearly discriminated and meaningfully ordered. Furthermore, unsupervised clustering on the low-dimensional data proved effective at further splitting the motifs into structural subfamilies characterized by very fine and physically relevant differences such as the presence of specific punctual or planar defects or of atoms with particular coordination features. Owing to these peculiarities, the chart also enabled tracking of the complex structural evolution in a reactive trajectory. In addition to visualization and analysis of complex structural landscapes, the presented approach offers a general, low-dimensional set of differentiable variables that has the potential to be used for exploration and enhanced sampling purposes.

Missense mutations in the membrane domain of PRRT2 affect its interaction with Nav1.2 voltage-gated sodium channels

PRRT2 is a neuronal protein that controls neuronal excitability and network stability by modulating voltage-gated Na⁺ channel (Nav). PRRT2 pathogenic variants cause pleiotropic syndromes including epilepsy, paroxysmal kinesigenic dyskinesia and episodic ataxia attributable to loss-of-function pathogenetic mechanism. Based on the evidence that the transmembrane domain of PRRT2 interacts with Nav1.2/1.6, we focused on eight missense mutations located within the domain that show expression and membrane localization similar to the wild-type protein. Molecular dynamics simulations showed that the mutants do not alter the structural stability of the PRRT2 membrane domain and preserve its conformation. Using affinity assays, we found that the A320V and V286M mutants displayed respectively decreased and increased binding to Nav1.2. Accordingly, surface biotinylation showed an increased Nav1.2 surface exposure induced by the A320V mutant. Electrophysiological analysis confirmed the lack of modulation of Nav1.2 biophysical properties by the A320V mutant with a loss-of-function phenotype, while the V286M mutant displayed a gain-of-function with respect to wild-type PRRT2 with a more pronounced left-shift of the inactivation kinetics and delayed recovery from inactivation. The data confirm the key role played by the PRRT2-Nav interaction in the pathogenesis of the PRRT2-linked disorders and suggest an involvement of the A320 and V286 residues in the interaction site. Given the similar clinical phenotype caused by the two mutations, we speculate that circuit instability and paroxysmal manifestations may arise when PRRT2 function is outside the physiological range.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

The intramembrane COOH-terminal domain of PRRT2 regulates voltage-dependent Na+ channels

Proline-rich transmembrane protein 2 (PRRT2) is the single causative gene for pleiotropic paroxysmal syndromes including epilepsy, kinesigenic dyskinesia, episodic ataxia and migraine. PRRT2 is a neuron-specific type-2 membrane protein with a COOH-terminal intramembrane domain and a long proline-rich NH₂-terminal cytoplasmic region. A large array of experimental data indicates that PRRT2 is a neuron stability gene that negatively controls intrinsic excitability by regulating surface membrane localization and biophysical properties of voltage-dependent Na⁺ channels Nav1.2 and Nav1.6, but not Nav1.1. To further investigate the regulatory role of PRRT2, we studied the structural features of this membrane protein with molecular dynamics simulations, and its structure-function relationships with Nav1.2 channels by biochemical and electrophysiological techniques. We found that the intramembrane COOH-terminal region maintains a stable conformation over time, with the first transmembrane domain forming a helix-loop-helix motif within the bilayer. The unstructured NH₂-terminal cytoplasmic region bound to the Nav1.2 better than the isolated COOH-terminal intramembrane domain, mimicking full-length PRRT2, while the COOH-terminal intramembrane domain was able to modulate Na⁺ current and channel biophysical properties, still maintaining the striking specificity for Nav1.2 vs Nav1.1. channels. The results identify PRRT2 as a dual-domain protein in which the NH₂-terminal cytoplasmic region acts as a binding antenna for Na⁺ channels, while the COOH-terminal membrane domain regulates channel exposure on the membrane and its biophysical properties.

The impact of pathogenic and artificial mutations on Claudin-5 selectivity from molecular dynamics simulations

Tight-junctions (TJs) are multi-protein complexes between adjacent endothelial or epithelial cells. In the blood-brain-barrier (BBB), they seal the paracellular space and the Claudin-5 (Cldn5) protein forms their backbone. Despite the fundamental role in brain homeostasis, little is known on Cldn5-based TJ assemblies. Different structural models were suggested, with Cldn5 protomers generating paracellular pores that restrict the passage of ions and small molecules. Recently, the first Cldn5 pathogenic mutation, G60R, was identified and shown to induce Cl^--selective channels and Na⁺ barriers in BBB TJs, providing an excellent opportunity to validate the structural models. Here, we used molecular dynamics to study the permeation of ions and water through two distinct G60R-Cldn5 paracellular architectures. Only the so-called Pore I reproduces the functional modification observed in experiments, displaying a free energy (FE) minimum for Cl^- and a barrier for Na⁺ consistent with anionic selectivity. We also studied the artificial Q57D and Q63D mutations in the constriction region, Q57 being conserved in Cldns except for cation permeable homologs. In both cases, we obtain FE profiles consistent with facilitated passage of cations. Our calculations provide the first in-silico description of a Cldn5 pathogenic mutation, further assessing the TJ Pore I model and yielding new insight on BBB's paracellular selectivity.

Computational study of ion permeation through claudin-4 paracellular channels

Claudins (Cldns) form a large family of protein homologs that are essential for the assembly of paracellular tight junctions (TJs), where they form channels or barriers with tissue-specific selectivity for permeants. In contrast to several family members whose physiological role has been identified, the function of claudin 4 (Cldn4) remains elusive, despite experimental evidence suggesting that it can form anion-selective TJ channels in the renal epithelium. Computational approaches have recently been employed to elucidate the molecular basis of Cldns' function, and hence could help in clarifying the role of Cldn4. In this work, we use structural modeling and all-atom molecular dynamics simulations to transfer two previously introduced structural models of Cldn-based paracellular complexes to Cldn4 to reproduce a paracellular anion channel. Free energy calculations for ionic transport through the pores allow us to establish the thermodynamic properties driving the ion-selectivity of the structures. While one model shows a cavity permeable to chloride and repulsive to cations, the other forms barrier to the passage of all the major physiological ions. Furthermore, our results confirm the charge selectivity role of the residue Lys65 in the first extracellular loop of the protein, rationalizing Cldn4 control of paracellular permeability.

In Silico Conformational Features of Botulinum Toxins A1 and E1 According to Intraluminal Acidification

Although botulinum neurotoxins (BoNTs) are among the most toxic compounds found in nature, their molecular mechanism of action is far from being elucidated. A key event is the conformational transition due to acidification of the interior of synaptic vesicles, leading to translocation of the BoNT catalytic domain into the neuronal cytosol. To investigate these conformational variations, homology modeling and atomistic simulations are combined to explore the internal dynamics of the sub-types BoNT/A1 (the most-used sub-type in medical applications) and BoNT/E1 (the most kinetically efficient sub-type). This first simulation study of di-chain BoNTs in closed and open states considers the effects of both neutral and acidic pH. The conformational mobility is driven by domain displacements of the ganglioside-binding site in the receptor binding domain, the translocation domain (HCNT) switch, and the belt α-helix, which present multiple conformations, depending on the primary sequence and the pH. Fluctuations of the belt α-helix are observed for closed conformations of the toxins and at acidic pH, while patches of more solvent-accessible residues appear under the same conditions in the core translocation domain HCNT. These findings suggest that, during translocation, the higher mobility of the belt could be transmitted to HCNT, leading to the favorable interaction of HCNT residues with the non-polar membrane environment.

Phenotypic and genetic spectrum of ATP6V1A encephalopathy: a disorder of lysosomal homeostasis

Vacuolar-type H+-ATPase (V-ATPase) is a multimeric complex present in a variety of cellular membranes that acts as an ATP-dependent proton pump and plays a key role in pH homeostasis and intracellular signalling pathways. In humans, 22 autosomal genes encode for a redundant set of subunits allowing the composition of diverse V-ATPase complexes with specific properties and expression. Sixteen subunits have been linked to human disease. Here we describe 26 patients harbouring 20 distinct pathogenic de novo missense ATP6V1A variants, mainly clustering within the ATP synthase α/β family-nucleotide-binding domain. At a mean age of 7 years (extremes: 6 weeks, youngest deceased patient to 22 years, oldest patient) clinical pictures included early lethal encephalopathies with rapidly progressive massive brain atrophy, severe developmental epileptic encephalopathies and static intellectual disability with epilepsy. The first clinical manifestation was early hypotonia, in 70%; 81% developed epilepsy, manifested as developmental epileptic encephalopathies in 58% of the cohort and with infantile spasms in 62%; 63% of developmental epileptic encephalopathies failed to achieve any developmental, communicative or motor skills. Less severe outcomes were observed in 23% of patients who, at a mean age of 10 years and 6 months, exhibited moderate intellectual disability, with independent walking and variable epilepsy. None of the patients developed communicative language. Microcephaly (38%) and amelogenesis imperfecta/enamel dysplasia (42%) were additional clinical features. Brain MRI demonstrated hypomyelination and generalized atrophy in 68%. Atrophy was progressive in all eight individuals undergoing repeated MRIs. Fibroblasts of two patients with developmental epileptic encephalopathies showed decreased LAMP1 expression, Lysotracker staining and increased organelle pH, consistent with lysosomal impairment and loss of V-ATPase function. Fibroblasts of two patients with milder disease, exhibited a different phenotype with increased Lysotracker staining, decreased organelle pH and no significant modification in LAMP1 expression. Quantification of substrates for lysosomal enzymes in cellular extracts from four patients revealed discrete accumulation. Transmission electron microscopy of fibroblasts of four patients with variable severity and of induced pluripotent stem cell-derived neurons from two patients with developmental epileptic encephalopathies showed electron-dense inclusions, lipid droplets, osmiophilic material and lamellated membrane structures resembling phospholipids. Quantitative assessment in induced pluripotent stem cell-derived neurons identified significantly smaller lysosomes. ATP6V1A-related encephalopathy represents a new paradigm among lysosomal disorders. It results from a dysfunctional endo-lysosomal membrane protein causing altered pH homeostasis. Its pathophysiology implies intracellular accumulation of substrates whose composition remains unclear, and a combination of developmental brain abnormalities and neurodegenerative changes established during prenatal and early postanal development, whose severity is variably determined by specific pathogenic variants.

Non-Markov-Type Analysis and Diffusion Map Analysis for Molecular Dynamics Trajectory of Chignolin at a High Temperature

We apply the non-Markov-type analysis of state-to-state transitions to nearly microsecond molecular dynamics (MD) simulation data at a folding temperature of a small artificial protein, chignolin, and we found that the time scales obtained are consistent with our previous result using the weighted ensemble simulations, which is a general path-sampling method to extract the kinetic properties of molecules. Previously, we also applied diffusion map (DM) analysis, which is one of a manifold of learning techniques, to the same trajectory of chignolin in order to cluster the conformational states and found that DM and relaxation mode analysis give similar results for the eigenvectors. In this paper, we divide the same trajectory into shorter pieces and further apply DM to such short-length trajectories to investigate how the obtained eigenvectors are useful to characterize the conformational change of chignolin.

Computational Assessment of Different Structural Models for Claudin-5 Complexes in Blood-Brain Barrier Tight Junctions

The blood-brain barrier (BBB) strictly regulates the exchange of ions and molecules between the blood and the central nervous system. Tight junctions (TJs) are multimeric structures that control the transport through the paracellular spaces between the adjacent brain endothelial cells of the BBB. Claudin-5 (Cldn5) proteins are essential for TJ formation and assemble into multiprotein complexes via cis-interactions within the same cell membrane and trans-interactions across two contiguous cells. Despite the relevant biological function of Cldn5 proteins and their role as targets of brain drug delivery strategies, the molecular details of their assembly within TJs are still unclear. Two different structural models have been recently introduced, in which Cldn5 dimers belonging to opposite cells join to generate paracellular pores. However, a comparison of these models in terms of ionic transport features is still lacking. In this work, we used molecular dynamics simulations and free energy (FE) calculations to assess the two Cldn5 pore models and investigate the thermodynamic properties of water and physiological ions permeating through them. Despite different FE profiles, both structures present single/multiple FE barriers to ionic permeation, while being permissive to water flux. These results reveal that both models are compatible with the physiological role of Cldn5 TJ strands. By identifying the protein-protein surface at the core of TJ Cldn5 assemblies, our computational investigation provides a basis for the rational design of synthetic peptides and other molecules capable of opening paracellular pores in the BBB.

Ca2+ binding to synapsin I regulates resting Ca2+ and recovery from synaptic depression in nerve terminals

Synapsin I (SynI) is a synaptic vesicle (SV)-associated phosphoprotein that modulates neurotransmission by controlling SV trafficking. The SynI C-domain contains a highly conserved ATP binding site mediating SynI oligomerization and SV clustering and an adjacent main Ca²⁺ binding site, whose physiological role is unexplored. Molecular dynamics simulations revealed that the E373K point mutation irreversibly deletes Ca²⁺ binding to SynI, still allowing ATP binding, but inducing a destabilization of the SynI oligomerization interface. Here, we analyzed the effects of this mutation on neurotransmitter release and short-term plasticity in excitatory and inhibitory synapses from primary hippocampal neurons. Patch-clamp recordings showed an increase in the frequency of miniature excitatory postsynaptic currents (EPSCs) that was totally occluded by exogenous Ca²⁺ chelators and associated with a constitutive increase in resting terminal Ca²⁺ concentrations. Evoked EPSC amplitude was also reduced, due to a decreased readily releasable pool (RRP) size. Moreover, in both excitatory and inhibitory synapses, we observed a marked impaired recovery from synaptic depression, associated with impaired RRP refilling and depletion of the recycling pool of SVs. Our study identifies SynI as a novel Ca²⁺ buffer in excitatory terminals. Blocking Ca²⁺ binding to SynI results in higher constitutive Ca²⁺ levels that increase the probability of spontaneous release and disperse SVs. This causes a decreased size of the RRP and an impaired recovery from depression due to the failure of SV reclustering after sustained high-frequency stimulation. The results indicate a physiological role of Ca²⁺ binding to SynI in the regulation of SV clustering and trafficking in nerve terminals.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K⁺ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants K_D and dissociation rates k_off were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.

Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations

Activation of voltage-gated ion channels is regulated by conformational changes of the voltage sensor domains (VSDs), four water- and ion-impermeable modules peripheral to the central, permeable pore domain. Anomalous currents, defined as ω-currents, have been recorded in response to mutations of residues on the VSD S4 helix and associated with ion fluxes through the VSDs. In humans, gene defects in the potassium channel Kv7.2 result in a broad range of epileptic disorders, from benign neonatal seizures to severe epileptic encephalopathies. Experimental evidence suggests that the R207Q mutation in S4, associated with peripheral nerve hyperexcitability, induces ω-currents at depolarized potentials, but the fine structural details are still elusive. In this work, we use atom-detailed molecular dynamics simulations and a refined model structure of the Kv7.2 VSD in the active conformation in a membrane/water environment to study the effect of R207Q and four additional mutations of proven clinical importance. Our results demonstrate that the R207Q mutant shows the most pronounced increase of hydration in the internal VSD cavity, a feature favoring the occurrence of ω-currents. Free energy and kinetics calculations of sodium permeation through the native and mutated VSD indicate as more favorable the formation of a cationic current in the latter. Overall, our simulations establish a mechanistic linkage between genetic variations and their physiological outcome, by providing a computational description that includes both thermodynamic and kinetic features of ion permeation associated with ω-currents.

Isobaric Labeling Proteomics Allows a High-Throughput Investigation of Protein Corona Orientation

The formation of the biomolecular corona represents a crucial factor in controlling the biological interactions and trafficking of nanomaterials. In this context, the availability of key epitopes exposed on the surface of the corona, and able to engage the biological machinery, is important to define the biological fate of the material. While the full biomolecular corona composition can be investigated by conventional bottom-up proteomics, the assessment of the spatial orientation of proteins in the corona in a high-throughput fashion is still challenging. In this work, we show that labeling corona proteins with isobaric tags in their native conditions and analyzing the MS/MS spectra of tryptic peptides allow an easy and high-throughput assessment of the inner/outer orientation of the corresponding proteins in the original corona. We put our results in the context of what is currently known of the protein corona of graphene-based nanomaterials. Our conclusions are in line with previous data and were confirmed by in silico calculations.

Genotype-phenotype correlations in patients with de novo KCNQ2 pathogenic variants

Objective

Early identification of de novo KCNQ2 variants in patients with epilepsy raises prognostic issues toward optimal management. We analyzed the clinical and genetic information from a cohort of patients with de novo KCNQ2 pathogenic variants to dissect genotype-phenotype correlations.

Methods

Patients with de novo KCNQ2 pathogenic variants were identified from Italy, Denmark, and Belgium. Atomic resolution Kv7.2 structures were also generated using homology modeling to map the variants.

Results

We included 34 patients with a mean age of 4.7 years. Median seizure onset was 2 days, mainly with focal seizures with autonomic signs. Twenty-two patients (65%) were seizure free at the mean age of 1.2 years. More than half of the patients (17/32) displayed severe/profound intellectual disability; however, 4 (13%) of them had a normal cognitive outcome.

A total of 28 de novo pathogenic variants were identified, most missense (25/28), and clustered in conserved regions of the protein; 6 variants recurred, and 7 were novel. We did not identify a relationship between variant position and seizure offset or cognitive outcome in patients harboring missense variants. Besides, recurrent variants were associated with overlapping epilepsy features but also variable evolution regarding the intellectual outcome.

Conclusions

We highlight the complexity of variant interpretation to assess the impact of a class of de novo KCNQ2 mutations. Genetic modifiers could be implicated, but the study paradigms to successfully address the impact of each single mutation need to be developed.

Thermodynamics and Kinetics of Ion Permeation in Wild-Type and Mutated Open Active Conformation of the Human α7 Nicotinic Receptor

Molecular studies of human pentameric ligand-gated ion channels (LGICs) expressed in neurons and at neuromuscular junctions are of utmost importance in the development of therapeutic strategies for neurological disorders. We focus here on the nicotinic acetylcholine receptor nAChR-α7, a homopentameric channel widely expressed in the human brain, with a proven role in a wide spectrum of disorders including schizophrenia and Alzheimer’s disease. By exploiting an all-atom structural model of the full (transmembrane and extracellular) protein in the open, agonist-bound conformation we recently developed, we evaluate the free energy and the mean first passage time of single-ion permeation using molecular dynamics simulations and the milestoning method with Voronoi tessellation. The results for the wild-type channel provide the first available mapping of the potential of mean force in the full-length α7 nAChR, reveal its expected cationic nature, and are in good agreement with simulation data for other channels of the LGIC family and with experimental data on nAChRs. We then investigate the role of a specific mutation directly related to ion selectivity in LGICs, the E-1′ → A-1′ substitution at the cytoplasmatic selectivity filter. We find that the mutation strongly affects sodium and chloride permeation in opposite directions, leading to a complete inversion of selectivity, at variance with the limited experimental results available that classify this mutant as cationic. We thus provide structural determinants for the observed cationic-to-anionic inversion, revealing a key role of the protonation state of residue rings far from the mutation, in the proximity of the hydrophobic channel gate.

Engineering REST-Specific Synthetic PUF Proteins to Control Neuronal Gene Expression: A Combined Experimental and Computational Study

Regulation of gene transcription is an essential mechanism for differentiation and adaptation of organisms. A key actor in this regulation process is the repressor element 1 (RE1)-silencing transcription factor (REST), a transcriptional repressor that controls more than 2000 putative target genes, most of which are neuron-specific. With the purpose of modulating REST expression, we exploited synthetic, ad hoc designed, RNA binding proteins (RBPs) able to specifically target and dock to REST mRNA. Among the various families of RBPs, we focused on the Pumilio and FBF (PUF) proteins, present in all eukaryotic organisms and controlling a variety of cellular functions. Here, a combined experimental and computational approach was used to design and test 8- and 16-repeat PUF proteins specific for REST mRNA. We explored the conformational properties and atomic features of the PUF-RNA recognition code by Molecular Dynamics simulations. Biochemical assays revealed that the 8- and 16-repeat PUF-based variants specifically bind the endogenous REST mRNA without affecting its translational regulation. The data also indicate a key role of stacking residues in determining the binding specificity. The newly characterized REST-specific PUF-based constructs act as excellent RNA-binding modules and represent a versatile and functional platform to specifically target REST mRNA and modulate its endogenous expression.

Neuronal firing modulation by a membrane-targeted photoswitch

Optical technologies allowing modulation of neuronal activity at high spatio-temporal resolution are becoming paramount in neuroscience. In this respect, azobenzene-based photoswitches are promising nanoscale tools for neuronal photostimulation. Here we engineered a light-sensitive azobenzene compound (Ziapin2) that stably partitions into the plasma membrane and causes its thinning through trans-dimerization in the dark, resulting in an increased membrane capacitance at steady state. We demonstrated that in neurons loaded with the compound, millisecond pulses of visible light induce a transient hyperpolarization followed by a delayed depolarization that triggers action potential firing. These effects are persistent and can be evoked in vivo up to 7 days, proving the potential of Ziapin2 for the modulation of membrane capacitance in the millisecond timescale, without directly affecting ion channels or local temperature.

Synapsins are expressed at neuronal and non-neuronal locations in Octopus vulgaris

Synapsins are a family of phosphoproteins fundamental to the regulation of neurotransmitter release. They are typically neuron-specific, although recent evidence pointed to their expression in non-neuronal cells where they play a role in exocytosis and vesicle trafficking. In this work, we characterized synapsin transcripts in the invertebrate mollusk Octopus vulgaris and present evidence of their expression not only in the brain but also in male and female reproductive organs. We identified three synapsin isoforms phylogenetically correlated to that of other invertebrates and with a modular structure characteristic of mammalian synapsins with a central, highly conserved C domain, important for the protein functions, and less conserved A, B and E domains. Our molecular modeling analysis further provided a solid background for predicting synapsin functional binding to ATP, actin filaments and secretory vesicles. Interestingly, we found that synapsin expression in ovary and testis increased during sexual maturation in cells with a known secretory role, potentially matching the occurrence of a secretion process. This might indicate that its secretory role has evolved across animals according to cell activity in spite of cell identity. We believe that this study may yield insights into the convergent evolution of ubiquitously expressed proteins between vertebrates and invertebrates.

De novo mutations of the ATP6V1A gene cause developmental encephalopathy with epilepsy

V-type proton (H⁺) ATPase (v-ATPase) is a multi-subunit proton pump that regulates pH homeostasis in all eukaryotic cells; in neurons, v-ATPase plays additional and unique roles in synapse function. Through whole exome sequencing, we identified de novo heterozygous mutations (p.Pro27Arg, p.Asp100Tyr, p.Asp349Asn, p.Asp371Gly) in ATP6V1A, encoding the A subunit of v-ATPase, in four patients with developmental encephalopathy with epilepsy. Early manifestations, observed in all patients, were developmental delay and febrile seizures, evolving to encephalopathy with profound delay, hypotonic/dyskinetic quadriparesis and intractable multiple seizure types in two patients (p.Pro27Arg, p.Asp100Tyr), and to moderate delay with milder epilepsy in the other two (p.Asp349Asn, p.Asp371Gly). Modelling performed on the available prokaryotic and eukaryotic structures of v-ATPase predicted p.Pro27Arg to perturb subunit interaction, p.Asp100Tyr to cause steric hindrance and destabilize protein folding, p.Asp349Asn to affect the catalytic function and p.Asp371Gly to impair the rotation process, necessary for proton transport. We addressed the impact of p.Asp349Asn and p.Asp100Tyr mutations on ATP6V1A expression and function by analysing ATP6V1A-overexpressing HEK293T cells and patients’ lymphoblasts. The p.Asp100Tyr mutant was characterized by reduced expression due to increased degradation. Conversely, no decrease in expression and clearance was observed for p.Asp349Asn. In HEK293T cells overexpressing either pathogenic or control variants, p.Asp349Asn significantly increased LysoTracker® fluorescence with no effects on EEA1 and LAMP1 expression. Conversely, p.Asp100Tyr decreased both LysoTracker® fluorescence and LAMP1 levels, leaving EEA1 expression unaffected. Both mutations decreased v-ATPase recruitment to autophagosomes, with no major impact on autophagy. Experiments performed on patients’ lymphoblasts using the LysoSensor™ probe revealed lower pH of endocytic organelles for p.Asp349Asn and a reduced expression of LAMP1 with no effect on the pH for p.Asp100Tyr. These data demonstrate gain of function for p.Asp349Asn characterized by an increased proton pumping in intracellular organelles, and loss of function for p.Asp100Tyr with decreased expression of ATP6V1A and reduced levels of lysosomal markers. We expressed p.Asp349Asn and p.Asp100Tyr in rat hippocampal neurons and confirmed significant and opposite effects in lysosomal labelling. However, both mutations caused a similar defect in neurite elongation accompanied by loss of excitatory inputs, revealing that altered lysosomal homeostasis markedly affects neurite development and synaptic connectivity. This study provides evidence that de novo heterozygous ATP6V1A mutations cause a developmental encephalopathy with a pathomechanism that involves perturbations of lysosomal homeostasis and neuronal connectivity, uncovering a novel role for v-ATPase in neuronal development.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Interfacing Graphene-Based Materials With Neural Cells

The scientific community has witnessed an exponential increase in the applications of graphene and graphene-based materials in a wide range of fields, from engineering to electronics to biotechnologies and biomedical applications. For what concerns neuroscience, the interest raised by these materials is two-fold. On one side, nanosheets made of graphene or graphene derivatives (graphene oxide, or its reduced form) can be used as carriers for drug delivery. Here, an important aspect is to evaluate their toxicity, which strongly depends on flake composition, chemical functionalization and dimensions. On the other side, graphene can be exploited as a substrate for tissue engineering. In this case, conductivity is probably the most relevant amongst the various properties of the different graphene materials, as it may allow to instruct and interrogate neural networks, as well as to drive neural growth and differentiation, which holds a great potential in regenerative medicine. In this review, we try to give a comprehensive view of the accomplishments and new challenges of the field, as well as which in our view are the most exciting directions to take in the immediate future. These include the need to engineer multifunctional nanoparticles (NPs) able to cross the blood-brain-barrier to reach neural cells, and to achieve on-demand delivery of specific drugs. We describe the state-of-the-art in the use of graphene materials to engineer three-dimensional scaffolds to drive neuronal growth and regeneration in vivo, and the possibility of using graphene as a component of hybrid composites/multi-layer organic electronics devices. Last but not least, we address the need of an accurate theoretical modeling of the interface between graphene and biological material, by modeling the interaction of graphene with proteins and cell membranes at the nanoscale, and describing the physical mechanism(s) of charge transfer by which the various graphene materials can influence the excitability and physiology of neural cells.

A refined model of claudin-15 tight junction paracellular architecture by molecular dynamics simulations

Tight-junctions between epithelial cells of biological barriers are specialized molecular structures that regulate the flux of solutes across the barrier, parallel to cell walls. The tight-junction backbone is made of strands of transmembrane proteins from the claudin family, but the molecular mechanism of its function is still not completely understood. Recently, the crystal structure of a mammalian claudin-15 was reported, displaying for the first time the detailed features of transmembrane and extracellular domains. Successively, a structural model of claudin-15-based paracellular channels has been proposed, suggesting a putative assembly that illustrates how claudins associate in the same cell (via cis interactions) and across adjacent cells (via trans interactions). Although very promising, the model offers only a static conformation, with residues missing in the most important extracellular regions and potential steric clashes. Here we present detailed atomic models of paracellular single and double pore architectures, obtained from the putative assembly and refined via structural modeling and all-atom molecular dynamics simulations in double membrane bilayer and water environment. Our results show an overall stable configuration of the complex with a fluctuating pore size. Extracellular residue loops in trans interaction are able to form stable contacts and regulate the size of the pore, which displays a stationary radius of 2.5-3.0 Å at the narrowest region. The side-by-side interactions of the cis configuration are preserved via stable hydrogen bonds, already predicted by cysteine crosslinking experiments. Overall, this work introduces an improved version of the claudin-15-based paracellular channel model that strengthens its validity and that can be used in further computational studies to understand the structural features of tight-junctions regulation.

Optogenetic Modulation of Intracellular Signalling and Transcription: Focus on Neuronal Plasticity

Several fields in neuroscience have been revolutionized by the advent of optogenetics, a technique that offers the possibility to modulate neuronal physiology in response to light stimulation. This innovative and far-reaching tool provided unprecedented spatial and temporal resolution to explore the activity of neural circuits underlying cognition and behaviour. With an exponential growth in the discovery and synthesis of new photosensitive actuators capable of modulating neuronal networks function, other fields in biology are experiencing a similar re-evolution. Here, we review the various optogenetic toolboxes developed to influence cellular physiology as well as the diverse ways in which these can be engineered to precisely modulate intracellular signalling and transcription. We also explore the processes required to successfully express and stimulate these photo-actuators in vivo before discussing how such tools can enlighten our understanding of neuronal plasticity at the systems level.

A Novel Topology of Proline-rich Transmembrane Protein 2 (PRRT2): HINTS FOR AN INTRACELLULAR FUNCTION AT THE SYNAPSE

Proline-rich transmembrane protein 2 (PRRT2) has been identified as the single causative gene for a group of paroxysmal syndromes of infancy, including epilepsy, paroxysmal movement disorders, and migraine. On the basis of topology predictions, PRRT2 has been assigned to the recently characterized family of Dispanins, whose members share the two-transmembrane domain topology with a large N terminus and short C terminus oriented toward the outside of the cell. Because PRRT2 plays a role at the synapse, it is important to confirm the exact orientation of its N and C termini with respect to the plasma membrane to get clues regarding its possible function. Using a combination of different experimental approaches, including live immunolabeling, immunogold electron microscopy, surface biotinylation and computational modeling, we demonstrate a novel topology for this protein. PRRT2 is a type II transmembrane protein in which only the second hydrophobic segment spans the plasma membrane, whereas the first one is associated with the internal surface of the membrane and forms a helix-loop-helix structure without crossing it. Most importantly, the large proline-rich N-terminal domain is not exposed to the extracellular space but is localized intracellularly, and only the short C terminus is extracellular (N_cyt/C_exo topology). Accordingly, we show that PRRT2 interacts with the Src homology 3 domain-bearing protein Intersectin 1, an intracellular protein involved in synaptic vesicle cycling. These findings will contribute to the clarification of the role of PRRT2 at the synapse and the understanding of pathogenic mechanisms on the basis of PRRT2-related neurological disorders.

Extended Phase-Space Methods for Enhanced Sampling in Molecular Simulations: A Review

Molecular Dynamics simulations are a powerful approach to study biomolecular conformational changes or protein-ligand, protein-protein, and protein-DNA/RNA interactions. Straightforward applications, however, are often hampered by incomplete sampling, since in a typical simulated trajectory the system will spend most of its time trapped by high energy barriers in restricted regions of the configuration space. Over the years, several techniques have been designed to overcome this problem and enhance space sampling. Here, we review a class of methods that rely on the idea of extending the set of dynamical variables of the system by adding extra ones associated to functions describing the process under study. In particular, we illustrate the Temperature Accelerated Molecular Dynamics (TAMD), Logarithmic Mean Force Dynamics (LogMFD), and Multiscale Enhanced Sampling (MSES) algorithms. We also discuss combinations with techniques for searching reaction paths. We show the advantages presented by this approach and how it allows to quickly sample important regions of the free-energy landscape via automatic exploration.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Temperature Accelerated Molecular Dynamics with Soft-Ratcheting Criterion Orients Enhanced Sampling by Low-Resolution Information

Many proteins exhibit an equilibrium between multiple conformations, some of them being characterized only by low-resolution information. Visiting all conformations is a demanding task for computational techniques performing enhanced but unfocused exploration of collective variable (CV) space. Otherwise, pulling a structure toward a target condition biases the exploration in a way difficult to assess. To address this problem, we introduce here the soft-ratcheting temperature-accelerated molecular dynamics (sr-TAMD), where the exploration of CV space by TAMD is coupled to a soft-ratcheting algorithm that filters the evolving CV values according to a predefined criterion. Any low resolution or even qualitative information can be used to orient the exploration. We validate this technique by exploring the conformational space of the inactive state of the catalytic domain of the adenyl cyclase AC from Bordetella pertussis. The domain AC gets activated by association with calmodulin (CaM), and the available crystal structure shows that in the complex the protein has an elongated shape. High-resolution data are not available for the inactive, CaM-free protein state, but hydrodynamic measurements have shown that the inactive AC displays a more globular conformation. Here, using as CVs several geometric centers, we use sr-TAMD to enhance CV space sampling while filtering for CV values that correspond to centers moving close to each other, and we thus rapidly visit regions of conformational space that correspond to globular structures. The set of conformations sampled using sr-TAMD provides the most extensive description of the inactive state of AC up to now, consistent with available experimental information.

Temperature-accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X-ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature-accelerated molecular dynamics (TAMD) simulations to generate putative all-atom models of globular conformations sampled by CaM-free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter-domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations.

Behavioural and physiological indicators of shelter dogs' welfare: reflections on the no-kill policy on free-ranging dogs in Italy revisited on the basis of 15 years of implementation

The Italian National Law 281 of 1991 forbids the euthanatization of free-ranging dogs, unless they have an incurable illness or are proved to be dangerous. Without neglecting the undeniable benefits of the "no-kill" policy, nevertheless it has brought about a chronic overpopulation in shelters and, as a result, higher costs of management and welfare problems since some dogs remain in the shelter for life. In 2004-2008, the Istituto Zooprofilattico Sperimentale of the Lazio and Tuscany regions carried out a survey in the Lazio Region to verify the effects of the Italian National Law 281/91 on free-ranging dog management following 15 years from its implementation. One of the aims of the study was an assessment of the welfare of dogs in a shelter sample (8 shelters out of 47 censused in the Lazio Region). 97 mixed-breed dogs were selected, their behaviour was studied and a blood sample was taken for each dog in order to determine the individual blood concentration of cortisol and the amount of oxidative damage (level of dRoms), as well as the amount of antioxidants to cope with it. Moreover, the total leukocyte count (leukogram) was accomplished. We ran general backward stepwise regression models using "level of antioxidant", "level of dRoms" and "level of serum cortisol" as dependent variables respectively. The results showed that the most important variable that improved the level of welfare of dogs consisted in having the opportunity to regularly go out of the cage for a walk, whereas other variables like gender, size of the cage (small, medium, large), being alone in the cage, and being neutered/entire, had no significant effect on the physiological indicators of welfare. Dogs that enjoyed the regular walk had a higher total antioxidant capacity, and performed a lower frequency of displacing activities and stereotyped behaviour. Moreover, oxidative stress parameters seem to be indicators well matched with behavioural indicators of stress. Thus, for the first time, markers of oxidative status are utilised for the welfare evaluation in the domestic dog. Furthermore, the results of this paper give some suggestion about how small steps can help to improve shelters and, furthermore, this paper intends to solicit the debate on the no-kill policy.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Comparison between Mean Forces and Swarms-of-Trajectories String Methods

The original formulation of the string method in collective variable space is compared with a recent variant called string method with swarms-of-trajectories. The assumptions made in the original method are revisited and the significance of the minimum free energy path (MFEP) is discussed in the context of reactive events. These assumptions are compared to those made in the string method with swarms-of-trajectories, and shown to be equivalent in a certain regime: in particular an expression for the path identified by the swarms-of-trajectories method is given and shown to be closely related to the MFEP. Finally, the algorithmic aspects of both methods are compared.

Intermediate state trapping of a voltage sensor

Voltage sensor domains (VSDs) regulate ion channels and enzymes by undergoing conformational changes depending on membrane electrical signals. The molecular mechanisms underlying the VSD transitions are not fully understood. Here, we show that some mutations of I241 in the S1 segment of the Shaker Kv channel positively shift the voltage dependence of the VSD movement and alter the functional coupling between VSD and pore domains. Among the I241 mutants, I241W immobilized the VSD movement during activation and deactivation, approximately halfway between the resting and active states, and drastically shifted the voltage activation of the ionic conductance. This phenotype, which is consistent with a stabilization of an intermediate VSD conformation by the I241W mutation, was diminished by the charge-conserving R2K mutation but not by the charge-neutralizing R2Q mutation. Interestingly, most of these effects were reproduced by the F244W mutation located one helical turn above I241. Electrophysiology recordings using nonnatural indole derivatives ruled out the involvement of cation-Π interactions for the effects of the Trp inserted at positions I241 and F244 on the channel’s conductance, but showed that the indole nitrogen was important for the I241W phenotype. Insight into the molecular mechanisms responsible for the stabilization of the intermediate state were investigated by creating in silico the mutations I241W, I241W/R2K, and F244W in intermediate conformations obtained from a computational VSD transition pathway determined using the string method. The experimental results and computational analysis suggest that the phenotype of I241W may originate in the formation of a hydrogen bond between the indole nitrogen atom and the backbone carbonyl of R2. This work provides new information on intermediate states in voltage-gated ion channels with an approach that produces minimum chemical perturbation.

Direct imaging of DNA fibers: the visage of double helix

Direct imaging becomes important when the knowledge at few/single molecule level is requested and where the diffraction does not allow to get structural and functional information. Here we report on the direct imaging of double stranded (ds) λ-DNA in the A conformation, obtained by combining a novel sample preparation method based on super hydrophobic DNA molecules self-aggregation process with transmission electron microscopy (TEM). The experimental breakthrough is the production of robust and highly ordered paired DNA nanofibers that allowed its direct TEM imaging and the double helix structure revealing.

Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer

An extremely scalable computational strategy is described for calculations of the potential of mean force (PMF) in multidimensions on massively distributed supercomputers. The approach involves coupling thousands of umbrella sampling (US) simulation windows distributed to cover the space of order parameters with a Hamiltonian molecular dynamics replica-exchange (H-REMD) algorithm to enhance the sampling of each simulation. In the present application, US/H-REMD is carried out in a two-dimensional (2D) space and exchanges are attempted alternatively along the two axes corresponding to the two order parameters. The US/H-REMD strategy is implemented on the basis of parallel/parallel multiple copy protocol at the MPI level, and therefore can fully exploit computing power of large-scale supercomputers. Here the novel technique is illustrated using the leadership supercomputer IBM Blue Gene/P with an application to a typical biomolecular calculation of general interest, namely the binding of calcium ions to the small protein Calbindin D9k. The free energy landscape associated with two order parameters, the distance between the ion and its binding pocket and the root-mean-square deviation (rmsd) of the binding pocket relative the crystal structure, was calculated using the US/H-REMD method. The results are then used to estimate the absolute binding free energy of calcium ion to Calbindin D9k. The tests demonstrate that the 2D US/H-REMD scheme greatly accelerates the configurational sampling of the binding pocket, thereby improving the convergence of the potential of mean force calculation.

"DFG-flip" in the insulin receptor kinase is facilitated by a helical intermediate state of the activation loop

We have characterized a large-scale inactive-to-active conformational change in the activation-loop of the insulin receptor kinase domain at the atomistic level via untargeted temperature-accelerated molecular dynamics (TAMD) and free-energy calculations using the string method. TAMD simulations consistently show folding of the A-loop into a helical conformation followed by unfolding to an active conformation, causing the highly conserved DFG-motif (Asp(1150), Phe(1151), and Gly(1152)) to switch from the inactive "D-out/F-in" to the nucleotide-binding-competent "D-in/F-out" conformation. The minimum free-energy path computed from the string method preserves these helical intermediates along the inactive-to-active path, and the thermodynamic free-energy differences are consistent with previous work on various other kinases. The mechanisms revealed by TAMD also suggest that the regulatory spine can be dynamically assembled/disassembled either by DFG-flip or by movement of the αC-helix. Together, these findings both broaden our understanding of kinase activation and point to intermediates as specific therapeutic targets.