Interplay between VSD, pore, and membrane lipids in electromechanical coupling in HCN channels.

Hyperpolarized-activated and cyclic nucleotide-gated (HCN) channels are the only members of the voltage-gated ion channel superfamily in mammals that open upon hyperpolarization, conferring them pacemaker properties that are instrumental for rhythmic firing of cardiac and neuronal cells. Activation of their voltage-sensor domains (VSD) upon hyperpolarization occurs through a downward movement of the S4 helix bearing the gating charges, which triggers a break in the alpha-helical hydrogen bonding pattern at the level of a conserved Serine residue. Previous structural and molecular simulation studies had however failed to capture pore opening that should be triggered by VSD activation, presumably because of a low VSD/pore electromechanical coupling efficiency and the limited timescales accessible to such techniques. Here, we have used advanced modeling strategies, including enhanced sampling molecular dynamics simulations exploiting comparisons between non-domain swapped voltage-gated ion channel structures trapped in closed and open states to trigger pore gating and characterize electromechanical coupling in HCN1. We propose that the coupling mechanism involves the reorganization of the interfaces between the VSD helices, in particular S4, and the pore-forming helices S5 and S6, subtly shifting the balance between hydrophobic and hydrophilic interactions in a 'domino effect' during activation and gating in this region. Remarkably, our simulations reveal state-dependent occupancy of lipid molecules at this emergent coupling interface, suggesting a key role of lipids in hyperpolarization-dependent gating. Our model provides a rationale for previous observations and a possible mechanism for regulation of HCN channels by the lipidic components of the membrane.

Ivacaftor-Mediated Potentiation of ABCB4 Missense Mutations Affecting Critical Motifs of the NBDs: Repositioning Perspectives for Hepatobiliary Diseases.

ABCB4 (ATP-binding cassette subfamily B member 4) is a hepatocanalicular floppase involved in biliary phosphatidylcholine (PC) secretion. Variations in the ABCB4 gene give rise to several biliary diseases, including progressive familial intrahepatic cholestasis type 3 (PFIC3), an autosomal recessive disease that can be lethal in the absence of liver transplantation. In this study, we investigated the effect and potential rescue of ten ABCB4 missense variations in NBD1:NBD2 homologous positions (Y403H/Y1043H, K435M/K1075M, E558K/E1200A, D564G/D1206G and H589Y/H1231Y) all localized at the conserved and functionally critical motifs of ABC transporters, six of which are mutated in patients. By combining structure analysis and in vitro studies, we found that all ten mutants were normally processed and localized at the canalicular membrane of HepG2 cells, but showed dramatically impaired PC transport activity that was significantly rescued by treatment with the clinically approved CFTR potentiator ivacaftor. Our results provide evidence that functional ABCB4 mutations are rescued by ivacaftor, paving the way for the repositioning of this potentiator for the treatment of selected patients with PFIC3 caused by mutations in the ATP-binding sites of ABCB4.

Acting on the CFTR Membrane-Spanning Domains Interface Rescues Some Misfolded Mutants.

ABC transporters are large membrane proteins sharing a complex architecture, which comprises two nucleotide-binding domains (NBDs) and two membrane-spanning domains (MSDs). These domains are susceptible to mutations affecting their folding and assembly. In the CFTR (ABCC7) protein, a groove has been highlighted in the MSD1 at the level of the membrane inner leaflet, containing both multiple mutations affecting folding and a binding site for pharmaco-chaperones that stabilize this region. This groove is also present in ABCB proteins, however it is covered by a short elbow helix, while in ABCC proteins it remains unprotected, due to a lower position of the elbow helix in the presence of the ABCC-specific lasso motif. Here, we identified a MSD1 second-site mutation located in the vicinity of the CFTR MSD1 groove that partially rescued the folding defect of cystic fibrosis causing mutations located within MSD1, while having no effect on the most frequent mutation, F508del, located within NBD1. A model of the mutated protein 3D structure suggests additional interaction between MSD1 and MSD2, strengthening the assembly at the level of the MSD intracellular loops. Altogether, these results provide insightful information in understanding key features of the folding and function of the CFTR protein in particular, and more generally, of type IV ABC transporters.

Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium.

SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins - the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening.

Pharmacological chaperones improve intra-domain stability and inter-domain assembly via distinct binding sites to rescue misfolded CFTR.

Protein misfolding is involved in a large number of diseases, among which cystic fibrosis. Complex intra- and inter-domain folding defects associated with mutations in the cystic fibrosis transmembrane regulator (CFTR) gene, among which p.Phe508del (F508del), have recently become a therapeutical target. Clinically approved correctors such as VX-809, VX-661, and VX-445, rescue mutant protein. However, their binding sites and mechanisms of action are still incompletely understood. Blind docking onto the 3D structures of both the first membrane-spanning domain (MSD1) and the first nucleotide-binding domain (NBD1), followed by molecular dynamics simulations, revealed the presence of two potential VX-809 corrector binding sites which, when mutated, abrogated rescue. Network of amino acids in the lasso helix 2 and the intracellular loops ICL1 and ICL4 allosterically coupled MSD1 and NBD1. Corrector VX-445 also occupied two potential binding sites on MSD1 and NBD1, the latter being shared with VX-809. Binding of both correctors on MSD1 enhanced the allostery between MSD1 and NBD1, hence the increased efficacy of the corrector combination. These correctors improve both intra-domain folding by stabilizing fragile protein-lipid interfaces and inter-domain assembly via distant allosteric couplings. These results provide novel mechanistic insights into the rescue of misfolded proteins by small molecules.

Insights into the Role of the Discontinuous TM7 Helix of Human Ferroportin through the Prism of the Asp325 Residue.

The negatively charged Asp325 residue has proved to be essential for iron export by human (HsFPN1) and primate Philippine tarsier (TsFpn) ferroportin, but its exact role during the iron transport cycle is still to be elucidated. It has been posited as being functionally equivalent to the metal ion-coordinating residue His261 in the C-lobe of the bacterial homolog BbFpn, but the two residues arise in different sequence motifs of the discontinuous TM7 transmembrane helix. Furthermore, BbFpn is not subject to extracellular regulation, contrary to its mammalian orthologues which are downregulated by hepcidin. To get further insight into the molecular mechanisms related to iron export in mammals in which Asp325 is involved, we investigated the behavior of the Asp325Ala, Asp325His, and Asp325Asn mutants in transiently transfected HEK293T cells, and performed a comparative structural analysis. Our biochemical studies clearly distinguished between the Asp325Ala and Asp325His mutants, which result in a dramatic decrease in plasma membrane expression of FPN1, and the Asp325Asn mutant, which alters iron egress without affecting protein localization. Analysis of the 3D structures of HsFPN1 and TsFpn in the outward-facing (OF) state indicated that Asp325 does not interact directly with metal ions but is involved in the modulation of Cys326 metal-binding capacity. Moreover, models of the architecture of mammalian proteins in the inward-facing (IF) state suggested that Asp325 may form an inter-lobe salt-bridge with Arg40 (TM1) when not interacting with Cys326. These findings allow to suggest that Asp325 may be important for fine-tuning iron recognition in the C-lobe, as well as for local structural changes during the IF-to-OF transition at the extracellular gate level. Inability to form a salt-bridge between TM1 and TM7b during iron translocation could lead to protein instability, as shown by the Asp325Ala and Asp325His mutants.

Structure and Sequence-based Computational Approaches to Allosteric Signal Transduction: Application to Electromechanical Coupling in Voltage-gated Ion Channels.

Allosteric signaling underlies the function of many biomolecules, including membrane proteins such as ion channels. Experimental methods have enabled specific quantitative insights into the coupling between the voltage sensing domain (VSD) and the pore gate of voltage-gated ion channels, located tens of Ångström apart from one another, as well as pinpointed specific residues and domains that participate in electromechanical signal transmission. Nevertheless, an overall atomic-level resolution picture is difficult to obtain from these methods alone. Today, thanks to the cryo-EM resolution revolution, we have access to high resolution structures of many different voltage-gated ion channels in various conformational states, putting a quantitative description of the processes at the basis of these changes within our close reach. Here, we review computational methods that build on structures to detect and characterize allosteric signaling and pathways. We then examine what has been learned so far about electromechanical coupling between VSD and pore using such methods. While no general theory of electromechanical coupling in voltage-gated ion channels integrating results from all these methods is available yet, we outline the types of insights that could be achieved in the near future using the methods that have not yet been put to use in this field of application.

New insights into structure and function of bis-phosphinic acid derivatives and implications for CFTR modulation.

C407 is a compound that corrects the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) protein carrying the p.Phe508del (F508del) mutation. We investigated the corrector effect of c407 and its derivatives on F508del-CFTR protein. Molecular docking and dynamics simulations combined with site-directed mutagenesis suggested that c407 stabilizes the F508del-Nucleotide Binding Domain 1 (NBD1) during the co-translational folding process by occupying the position of the p.Phe1068 side chain located at the fourth intracellular loop (ICL4). After CFTR domains assembly, c407 occupies the position of the missing p.Phe508 side chain. C407 alone or in combination with the F508del-CFTR corrector VX-809, increased CFTR activity in cell lines but not in primary respiratory cells carrying the F508del mutation. A structure-based approach resulted in the synthesis of an extended c407 analog G1, designed to improve the interaction with ICL4. G1 significantly increased CFTR activity and response to VX-809 in primary nasal cells of F508del homozygous patients. Our data demonstrate that in-silico optimized c407 derivative G1 acts by a mechanism different from the reference VX-809 corrector and provide insights into its possible molecular mode of action. These results pave the way for novel strategies aiming to optimize the flawed ICL4-NBD1 interface.

A dynamic view of DNA structure within the nucleosome: Biological implications.

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNANCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNANCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNANCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNANCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNANCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

Functional and Pharmacological Characterization of the Rare CFTR Mutation W361R.

Understanding the functional consequence of rare cystic fibrosis (CF) mutations is mandatory for the adoption of precision therapeutic approaches for CF. Here we studied the effect of the very rare CF mutation, W361R, on CFTR processing and function. We applied western blot, patch clamp and pharmacological modulators of CFTR to study the maturation and ion transport properties of pEGFP-WT and mutant CFTR constructs, W361R, F508del and L69H-CFTR, expressed in HEK293 cells. Structural analyses were also performed to study the molecular environment of the W361 residue. Western blot showed that W361R-CFTR was not efficiently processed to a mature band C, similar to F508del CFTR, but unlike F508del CFTR, it did exhibit significant transport activity at the cell surface in response to cAMP agonists. Importantly, W361R-CFTR also responded well to CFTR modulators: its maturation defect was efficiently corrected by VX-809 treatment and its channel activity further potentiated by VX-770. Based on these results, we postulate that W361R is a novel class-2 CF mutation that causes abnormal protein maturation which can be corrected by VX-809, and additionally potentiated by VX-770, two FDA-approved small molecules. At the structural level, W361 is located within a class-2 CF mutation hotspot that includes other mutations that induce variable disease severity. Analysis of the 3D structure of CFTR within a lipid environment indicated that W361, together with other mutations located in this hotspot, is at the edge of a groove which stably accommodates lipid acyl chains. We suggest this lipid environment impacts CFTR folding, maturation and response to CFTR modulators.

Targeting different binding sites in the CFTR structures allows to synergistically potentiate channel activity.

Recent evidence shows that combination of correctors and potentiators, such as the drug ivacaftor (VX-770), can significantly restore the functional expression of mutated Cystic Fibrosis Transmembrane conductance Regulator (CFTR), an anion channel which is mutated in cystic fibrosis (CF). The success of these combinatorial therapies highlights the necessity of identifying a broad panel of specific binding mode modulators, occupying several distinct binding sites at structural level. Here, we identified two small molecules, SBC040 and SBC219, which are two efficient cAMP-independent potentiators, acting at low concentration of forskolin with EC50 close to 1 µM and in a synergic way with the drug VX-770 on several CFTR mutants of classes II and III. Molecular dynamics simulations suggested potential SBC binding sites at the vicinity of ATP-binding sites, distinct from those currently proposed for VX-770, outlining SBC molecules as members of a new family of potentiators.

Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

Ferroportin 1 (FPN1) is a major facilitator superfamily transporter that is essential for proper maintenance of human iron homeostasis at the systemic and cellular level. FPN1 dysfunction leads to the progressive accumulation of iron in reticuloendothelial cells, causing hemochromatosis type 4A (or ferroportin disease), an autosomal dominant disorder that displays large phenotypic heterogeneity. Although crystal structures have unveiled the outward- and inward-facing conformations of the bacterial homolog Bdellovibrio bacteriovorus Fpn (or Bd2019) and calcium has recently been identified as an essential cofactor, our molecular understanding of the iron transport mechanism remains incomplete. Here, we used a combination of molecular modeling, molecular dynamics simulations, and Ala site-directed mutagenesis, followed by complementary in vitro functional analyses, to explore the structural architecture of the human FPN1 intracellular gate. We reveal an interdomain network that involves 5 key amino acids and is likely very important for stability of the iron exporter facing the extracellular milieu. We also identify inter- and intradomain interactions that rely on the 2 Asp84 and Asn174 critical residues and do not exist in the bacterial homolog. These interactions are thought to play an important role in the modulation of conformational changes during the transport cycle. We interpret these results in the context of hemochromatosis type 4A, reinforcing the idea that different categories of loss-of-function mutations exist. Our findings provide an unprecedented view of the human FPN1 outward-facing structure and the particular function of the so-called "gating residues" in the mechanism of iron export.-Guellec, J., Elbahnsi, A., Le Tertre, M., Uguen, K., Gourlaouen, I., Férec, C., Ka, C., Callebaut, I., Le Gac, G. Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

Combining theoretical and experimental data to decipher CFTR 3D structures and functions.

Cryo-electron microscopy (cryo-EM) has recently provided invaluable experimental data about the full-length cystic fibrosis transmembrane conductance regulator (CFTR) 3D structure. However, this experimental information deals with inactive states of the channel, either in an apo, quiescent conformation, in which nucleotide-binding domains (NBDs) are widely separated or in an ATP-bound, yet closed conformation. Here, we show that 3D structure models of the open and closed forms of the channel, now further supported by metadynamics simulations and by comparison with the cryo-EM data, could be used to gain some insights into critical features of the conformational transition toward active CFTR forms. These critical elements lie within membrane-spanning domains but also within NBD1 and the N-terminal extension, in which conformational plasticity is predicted to occur to help the interaction with filamin, one of the CFTR cellular partners.

Holding the Nucleosome Together: A Quantitative Description of the DNA-Histone Interface in Solution.

The nucleosome is the fundamental unit of eukaryotic genome packaging in the chromatin. In this complex, the DNA wraps around eight histone proteins to form a superhelical double helix. The resulting bending, stronger than anything observed in free DNA, raises the question of how such a distortion is stabilized by the proteic and solvent environments. In this work, the DNA-histone interface in solution was exhaustively analyzed from nucleosome structures generated by molecular dynamics. An original Voronoi tessellation technique, measuring the topology of interacting elements without any empirical or subjective adjustment, was used to characterize the interface in terms of contact area and occurrence. Our results revealed an interface more robust than previously known, combining extensive, long-lived nonelectrostatic and electrostatic interactions between DNA and both structured and unstructured histone regions. Cation accumulation makes the proximity of juxtaposed DNA gyres in the superhelix possible by shielding the strong electrostatic repulsion of the charged phosphate groups. Overall, this study provides new insights on the nucleosome cohesion, explaining how DNA distortions can be maintained in a nucleoprotein complex.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics.

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.