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.

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.

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.

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

The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav 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 Nav1.2 model and the cryo-EM Nav1.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 Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav 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 NH2-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 NH2-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 versus Nav1.1. channels. The results identify PRRT2 as a dual-domain protein in which the NH2-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.

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.

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 Ca2+ binding site, whose physiological role is unexplored. Molecular dynamics simulations revealed that the E373K point mutation irreversibly deletes Ca2+ 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 Ca2+ chelators and associated with a constitutive increase in resting terminal Ca2+ 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 Ca2+ buffer in excitatory terminals. Blocking Ca2+ binding to SynI results in higher constitutive Ca2+ 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 Ca2+ binding to SynI in the regulation of SV clustering and trafficking in nerve terminals.

Multiscale modelling of claudin-based assemblies: A magnifying glass for novel structures of biological interfaces.

Claudins (Cldns) define a family of transmembrane proteins that are the major determinants of the tight junction integrity and tissue selectivity. They promote the formation of either barriers or ion-selective channels at the interface between two facing cells, across the paracellular space. Multiple Cldn subunits form complexes that include cis- (intracellular) interactions along the membrane of a single cell and trans- (intercellular) interactions across adjacent cells. The first description of Cldn assemblies was provided by electron microscopy, while electrophysiology, mutagenesis and cell biology experiments addressed the functional role of different Cldn homologs. However, the investigation of the molecular details of Cldn subunits and complexes are hampered by the lack of experimental native structures, currently limited to Cldn15. The recent implementation of computer-based techniques greatly contributed to the elucidation of Cldn properties. Molecular dynamics simulations and docking calculations were extensively used to refine the first Cldn multimeric model postulated from the crystal structure of Cldn15, and contributed to the introduction of a novel, alternative, arrangement. While both these multimeric assemblies were found to account for the physiological properties of some family members, they gave conflicting results for others. In this review, we illustrate the major findings on Cldn-based systems that were achieved by using state-of-the-art computational methodologies. The information provided by these results could be useful to improve the characterization of the Cldn properties and help the design of new efficient strategies to control the paracellular transport of drugs or other molecules.

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.

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.

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.

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.

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 KD and dissociation rates koff 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.

Expanding the Nude SCID/CID Phenotype Associated with FOXN1 Homozygous, Compound Heterozygous, or Heterozygous Mutations.

Human nude SCID is a rare autosomal recessive inborn error of immunity (IEI) characterized by congenital athymia, alopecia, and nail dystrophy. Few cases have been reported to date. However, the recent introduction of newborn screening for IEIs and high-throughput sequencing has led to the identification of novel and atypical cases. Moreover, immunological alterations have been recently described in patients carrying heterozygous mutations. The aim of this paper is to describe the extended phenotype associated with FOXN1 homozygous, compound heterozygous, or heterozygous mutations. We collected clinical and laboratory information of a cohort of 11 homozygous, 2 compound heterozygous, and 5 heterozygous patients with recurrent severe infections. All, except one heterozygous patient, had signs of CID or SCID. Nail dystrophy and alopecia, that represent the hallmarks of the syndrome, were not always present, while almost 50% of the patients developed Omenn syndrome. One patient with hypomorphic compound heterozygous mutations had a late-onset atypical phenotype. A SCID-like phenotype was observed in 4 heterozygous patients coming from the same family. A spectrum of clinical manifestations may be associated with different mutations. The severity of the clinical phenotype likely depends on the amount of residual activity of the gene product, as previously observed for other SCID-related genes. The severity of the manifestations in this heterozygous family may suggest a mechanism of negative dominance of the specific mutation or the presence of additional mutations in noncoding regions.

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.

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.