MDFF_NM: Improved Molecular Dynamics Flexible Fitting into Cryo-EM Density Maps with a Multireplica Normal Mode-Based Search.

Molecular Dynamics Flexible Fitting (MDFF) is a widely used tool to refine high-resolution structures into cryo-EM density maps. Despite many successful applications, MDFF is still limited by its high computational cost, overfitting, accuracy, and performance issues due to entrapment within wrong local minima. Modern ensemble-based MDFF tools have generated promising results in the past decade. In line with these studies, we present MDFF_NM, a stochastic hybrid flexible fitting algorithm combining Normal Mode Analysis (NMA) and simulation-based flexible fitting. Initial tests reveal that, besides accelerating the fitting process, MDFF_NM increases the diversity of fitting routes leading to the target, uncovering ensembles of conformations in closer agreement with experimental data. The potential integration of MDFF_NM with other existing methods and integrative modeling pipelines is also discussed.

Human Kir2.1 Potassium Channel: Protocols for Cryo-EM Data Processing and Molecular Dynamics Simulations.

Kir channels are potassium (K+) channels responsible for the mechanism of inward rectification, which plays a fundamental role in maintaining the resting membrane potential. There are seven Kir subfamilies, and their opening and closing mechanism is regulated by different regulatory factors. Genetically inherited defects in Kir channels are responsible for several rare human diseases, and for most of them, there are currently no effective therapeutic treatments. High-resolution structural information is not available for several members within the Kir subfamilies. Recently, our group achieved a significant breakthrough by utilizing cryo-EM single-particle analysis to elucidate the first structure of the human Kir2.1 channel. We present here the data processing protocol of the cryo-EM data of the human Kir2.1 channel, which is applicable to the structural determination of other ion channels by cryo-EM single-particle analysis. We also introduce a protocol designed to assess the structural heterogeneity within the cryo-EM data, allowing for the identification of other possible protein structure conformations present in the collected data. Moreover, we present a protocol for conducting all-atom molecular dynamics (MD) simulations for K+ channels, which can be incorporated into various membrane models to simulate different environments. We also propose some methods for analyzing the MD simulations, with a particular emphasis on assessing the local mobility of protein residues.

Overview of Membrane Protein Sample Preparation for Single-Particle Cryo-Electron Microscopy Analysis.

Single-particle cryo-electron microscopy (cryo-EM SPA) has recently emerged as an exceptionally well-suited technique for determining the structure of membrane proteins (MPs). Indeed, in recent years, huge increase in the number of MPs solved via cryo-EM SPA at a resolution better than 3.0 Å in the Protein Data Bank (PDB) has been observed. However, sample preparation remains a significant challenge in the field. Here, we evaluated the MPs solved using cryo-EM SPA deposited in the PDB in the last two years at a resolution below 3.0 Å. The most critical parameters for sample preparation are as follows: (i) the surfactant used for protein extraction from the membrane, (ii) the surfactant, amphiphiles, nanodiscs or other molecules present in the vitrification step, (iii) the vitrification method employed, and (iv) the type of grids used. The aim is not to provide a definitive answer on the optimal sample conditions for cryo-EM SPA of MPs but rather assess the current trends in the MP structural biology community towards obtaining high-resolution cryo-EM structures.

Cryo-electron microscopy unveils unique structural features of the human Kir2.1 channel.

We present the first structure of the human Kir2.1 channel containing both transmembrane domain (TMD) and cytoplasmic domain (CTD). Kir2.1 channels are strongly inward-rectifying potassium channels that play a key role in maintaining resting membrane potential. Their gating is modulated by phosphatidylinositol 4,5-bisphosphate (PIP2). Genetically inherited defects in Kir2.1 channels are responsible for several rare human diseases, including Andersen's syndrome. The structural analysis (cryo-electron microscopy), surface plasmon resonance, and electrophysiological experiments revealed a well-connected network of interactions between the PIP2-binding site and the G-loop through residues R312 and H221. In addition, molecular dynamics simulations and normal mode analysis showed the intrinsic tendency of the CTD to tether to the TMD and a movement of the secondary anionic binding site to the membrane even without PIP2. Our results revealed structural features unique to human Kir2.1 and provided insights into the connection between G-loop and gating and the pathological mechanisms associated with this channel.

Unexpected Gating Behaviour of an Engineered Potassium Channel Kir.

In this study, we investigated the dynamics and functional characteristics of the KirBac3.1 S129R, a mutated bacterial potassium channel for which the inner pore-lining helix (TM2) was engineered so that the bundle crossing is trapped in an open conformation. The structure of this channel has been previously determined at high atomic resolution. We explored the dynamical characteristics of this open state channel using an in silico method MDeNM that combines molecular dynamics simulations and normal modes. We captured the global and local motions at the mutation level and compared these data with HDX-MS experiments. MDeNM provided also an estimation of the probability of the different opening states that are in agreement with our electrophysiological experiments. In the S129R mutant, the Arg129 mutation releases the two constriction points in the channel that existed in the wild type but interestingly creates another restriction point.

Biochemistry, structure, and cellular internalization of a four nanobody-bearing Fc dimer.

VHH stands for the variable regions of heavy chain only of camelid IgGs. The VHH family forms a set of interesting proteins derived from antibodies that maintain their capacity to recognize the antigen, despite their relatively small molecular weight (in the 12,000 Da range). Continuing our exploration of the possibilities of those molecules, we chose to design alternative molecules with maintained antigen recognition, but enhanced capacity, by fusing four VHH with one Fc, the fragment crystallizable region of antibodies. In doing so, we aimed at having a molecule with superior quantitative antigen recognition (×4) while maintaining its size below the 110 kDa. In the present paper, we described the building of those molecules that we coined VHH2 -Fc-VHH2 . The structure of VHH2 -Fc-VHH2 in complex with HER2 antigen was determined using electronic microscopy and modeling. The molecule is shown to bind four HER2 proteins at the end of its flexible arms. VHH2 -Fc-VHH2 also shows an internalization capacity via HER2 receptor superior to the reference anti-HER2 monoclonal antibody, Herceptin®, and to a simple fusion of two VHH with one Fc (VHH2 -Fc). This new type of molecules, VHH2 -Fc-VHH2 , could be an interesting addition to the therapeutic arsenal with multiple applications, from diagnostic to therapy.

Structural modeling of a novel membrane-bound globin-coupled sensor in Geobacter sulfurreducens.

Globin-coupled sensors (GCS) usually consist of three domains: a sensor/globin, a linker, and a transmitter domain. The globin domain (GD), activated by ligand binding and/or redox change, induces an intramolecular signal transduction resulting in a response of the transmitter domain. Depending on the nature of the transmitter domain, GCSs can have different activities and functions, including adenylate and di-guanylate cyclase, histidine kinase activity, aerotaxis and/or oxygen sensing function. The gram-negative delta-proteobacterium Geobacter sulfurreducens expresses a protein with a GD covalently linked to a four transmembrane domain, classified, by sequence similarity, as GCS (GsGCS). While its GD is fully characterized, not so its transmembrane domain, which is rarely found in the globin superfamily. In the present work, GsGCS was characterized spectroscopically and by native ion mobility-mass spectrometry in combination with cryo-electron microscopy. Although lacking high resolution, the oligomeric state and the electron density map were valuable for further rational modeling of the full-length GsGCS structure. This model demonstrates that GsGCS forms a transmembrane domain-driven tetramer with minimal contact between the GDs and with the heme groups oriented outward. This organization makes an intramolecular signal transduction less likely. Our results, including the auto-oxidation rate and redox potential, suggest a potential role for GsGCS as redox sensor or in a membrane-bound e-/H+ transfer. As such, GsGCS might act as a player in connecting energy production to the oxidation of organic compounds and metal reduction. Database searches indicate that GDs linked to a four or seven helices transmembrane domain occur more frequently than expected.

MicroRNA Biophysically Modulates Cardiac Action Potential by Direct Binding to Ion Channel.

BACKGROUND: MicroRNAs (miRs) play critical roles in regulation of numerous biological events, including cardiac electrophysiology and arrhythmia, through a canonical RNA interference mechanism. It remains unknown whether endogenous miRs modulate physiologic homeostasis of the heart through noncanonical mechanisms. METHODS: We focused on the predominant miR of the heart (miR1) and investigated whether miR1 could physically bind with ion channels in cardiomyocytes by electrophoretic mobility shift assay, in situ proximity ligation assay, RNA pull down, and RNA immunoprecipitation assays. The functional modulations of cellular electrophysiology were evaluated by inside-out and whole-cell patch clamp. Mutagenesis of miR1 and the ion channel was used to understand the underlying mechanism. The effect on the heart ex vivo was demonstrated through investigating arrhythmia-associated human single nucleotide polymorphisms with miR1-deficient mice. RESULTS: We found that endogenous miR1 could physically bind with cardiac membrane proteins, including an inward-rectifier potassium channel Kir2.1. The miR1-Kir2.1 physical interaction was observed in mouse, guinea pig, canine, and human cardiomyocytes. miR1 quickly and significantly suppressed IK1 at sub-pmol/L concentration, which is close to endogenous miR expression level. Acute presence of miR1 depolarized resting membrane potential and prolonged final repolarization of the action potential in cardiomyocytes. We identified 3 miR1-binding residues on the C-terminus of Kir2.1. Mechanistically, miR1 binds to the pore-facing G-loop of Kir2.1 through the core sequence AAGAAG, which is outside its RNA interference seed region. This biophysical modulation is involved in the dysregulation of gain-of-function Kir2.1-M301K mutation in short QT or atrial fibrillation. We found that an arrhythmia-associated human single nucleotide polymorphism of miR1 (hSNP14A/G) specifically disrupts the biophysical modulation while retaining the RNA interference function. It is remarkable that miR1 but not hSNP14A/G relieved the hyperpolarized resting membrane potential in miR1-deficient cardiomyocytes, improved the conduction velocity, and eliminated the high inducibility of arrhythmia in miR1-deficient hearts ex vivo. CONCLUSIONS: Our study reveals a novel evolutionarily conserved biophysical action of endogenous miRs in modulating cardiac electrophysiology. Our discovery of miRs' biophysical modulation provides a more comprehensive understanding of ion channel dysregulation and may provide new insights into the pathogenesis of cardiac arrhythmias.

Integrative Study of the Structural and Dynamical Properties of a KirBac3.1 Mutant: Functional Implication of a Highly Conserved Tryptophan in the Transmembrane Domain.

ATP-sensitive potassium (K-ATP) channels are ubiquitously expressed on the plasma membrane of cells in several organs, including the heart, pancreas, and brain, and they govern a wide range of physiological processes. In pancreatic ß-cells, K-ATP channels composed of Kir6.2 and SUR1 play a key role in coupling blood glucose and insulin secretion. A tryptophan residue located at the cytosolic end of the transmembrane helix is highly conserved in eukaryote and prokaryote Kir channels. Any mutation on this amino acid causes a gain of function and neonatal diabetes mellitus. In this study, we have investigated the effect of mutation on this highly conserved residue on a KirBac channel (prokaryotic homolog of mammalian Kir6.2). We provide the crystal structure of the mutant KirBac3.1 W46R (equivalent to W68R in Kir6.2) and its conformational flexibility properties using HDX-MS. In addition, the detailed dynamical view of the mutant during the gating was investigated using the in silico method. Finally, functional assays have been performed. A comparison of important structural determinants for the gating mechanism between the wild type KirBac and the mutant W46R suggests interesting structural and dynamical clues and a mechanism of action of the mutation that leads to the gain of function.

New Structural insights into Kir channel gating from molecular simulations, HDX-MS and functional studies.

Inward rectifier potassium (Kir) channels play diverse and important roles in shaping action potentials in biological membranes. An increasing number of diseases are now known to be directly associated with abnormal Kir function. However, the gating of Kir still remains unknown. To increase our understanding of its gating mechanism, a dynamical view of the entire channel is essential. Here the gating activation was studied using a recent developped in silico method, MDeNM, which combines normal mode analysis and molecular dynamics simulations that showed for the very first time the importance of interrelated collective and localized conformational movements. In particular, we highlighted the role played by concerted movements of the different regions throughout the entire protein, such as the cytoplasmic and transmembrane domains and the slide helices. In addition, the HDX-MS analysis achieved in these studies provided a comprehensive and detailed view of the dynamics associated with open/closed transition of the Kir channel in coherence with the theoretical results. MDeNM gives access to the probability of the different opening states that are in agreement with our electrophysiological experiments. The investigations presented in this article are important to remedy dysfunctional channels and are of interest for designing new pharmacological compounds.

A New Strategy for Atomic Flexible Fitting in Cryo-EM Maps by Molecular Dynamics with Excited Normal Modes (MDeNM-EMfit).

Previous studies demonstrated the efficiency of the Molecular Dynamics with excited Normal Modes (MDeNM) method on the characterization of large structural changes at a low computational cost. We present here MDeNM-EMfit, an extension of the original method designed to the flexible fit of structures into cryo-EM maps. Here, instead of a uniform exploration of the collective motions described by normal modes, sampling is directed toward conformations with increased correlations with the experimental map. Future perspectives to improve the accuracy of fitting and speed of calculations are discussed in light of the results.

The CRACAM Robot: Two-Dimensional Crystallization of Membrane Protein.

Membrane proteins are key cellular components that perform essential functions. They are major therapeutic targets. Electron crystallography can provide structural experimental information at atomic scale for membrane proteins forming two-dimensional (2D) crystals. There are two different methods to produce 2D crystals of membrane proteins. (1) either directly in the bulk of the solution (2) or under a lipid monolayer at the air-water interface. This extra lipid monolayer helps to pre-orient the proteins in order to facilitate the growth of 2D crystals. We present here these two methods for 2D crystallization of membrane proteins implemented in a fully automated robot called CRACAM. These methods require small volume of low concentration of proteins, making it possible to explore more conditions with the same amount of protein. These automated methods outperform traditional 2D crystallization approaches in terms of accuracy, flexibility, and throughput.

Cryo-electron microscopy and X-ray crystallography: complementary approaches to structural biology and drug discovery.

The invention of the electron microscope has greatly enhanced the view scientists have of small structural details. Since its implementation, this technology has undergone considerable evolution and the resolution that can be obtained for biological objects has been extended. In addition, the latest generation of cryo-electron microscopes equipped with direct electron detectors and software for the automated collection of images, in combination with the use of advanced image-analysis methods, has dramatically improved the performance of this technique in terms of resolution. While calculating a sub-10 Šresolution structure was an accomplishment less than a decade ago, it is now common to generate structures at sub-5 Šresolution and even better. It is becoming possible to relatively quickly obtain high-resolution structures of biological molecules, in particular large ones (>500 kDa) which, in some cases, have resisted more conventional methods such as X-ray crystallography or nuclear magnetic resonance (NMR). Such newly resolved structures may, for the first time, shed light on the precise mechanisms that are essential for cellular physiological processes. The ability to attain atomic resolution may support the development of new drugs that target these proteins, allowing medicinal chemists to understand the intimacy of the relationship between their molecules and targets. In addition, recent developments in cryo-electron microscopy combined with image analysis can provide unique information on the conformational variability of macromolecular complexes. Conformational flexibility of macromolecular complexes can be investigated using cryo-electron microscopy and multiconformation reconstruction methods. However, the biochemical quality of the sample remains the major bottleneck to routine cryo-electron microscopy-based determination of structures at very high resolution.

Utilization of interferometric light microscopy for the rapid analysis of virus abundance in a river.

There is a constant need for direct counting of biotic nanoparticles such as viruses to unravel river functioning. We used, for the first time in freshwater, a new method based on interferometry differentiating viruses from other particles such as membrane vesicles. In the French Marne River, viruses represented between 42 and 72% of the particles. A spring monitoring in 2014 revealed their increase (2.1 × 107 to 2.1 × 108 mL-1) linked to an increase in algal biomass and diversity of bacterial plankton. Predicted virus size distributions were in agreement with transmission electron microscopy analysis suggesting a dominance of large viruses (≥60 nm).

[Cryo-microscopy, an alternative to the X-ray crystallography?]. / La cryo-microscopie, une alternative à la cristallographie aux rayons X ?

Recent technological advances have revolutionized the field of structural biologists. Specifically, dramatic progress related to the development of new electron microscopes and image capture (direct electron detection camera) and the provision of new image analysis software has led to a breakthrough in terms of resolution attained using cryo-electron transmission microscopy. It is thus possible to calculate relatively quickly high-resolution structures of biological molecules whom structural study still resists to more conventional methods such as X-ray diffraction or nuclear magnetic resonance (NMR). These structures thus obtained may also bring complementary structural information to those already described by other methods. Some of these new structures resolved through cryo-electron microscopy revealed for the first time the precise operation of essential mechanisms necessary for the good physiological process of a cell. The ability to solve these structures at atomic resolution detail is essential for the development of new drugs that target these proteins of therapeutic interest. Thanks to these advanced techniques that we summarize in this revew, biological and medical issues have now become accessible, whereas this approach was inconceivable only five yeras ago. ‡.

The FANCD2-FANCI complex is recruited to DNA interstrand crosslinks before monoubiquitination of FANCD2.

The Fanconi anaemia (FA) pathway is important for the repair of DNA interstrand crosslinks (ICL). The FANCD2-FANCI complex is central to the pathway, and localizes to ICLs dependent on its monoubiquitination. It has remained elusive whether the complex is recruited before or after the critical monoubiquitination. Here, we report the first structural insight into the human FANCD2-FANCI complex by obtaining the cryo-EM structure. The complex contains an inner cavity, large enough to accommodate a double-stranded DNA helix, as well as a protruding Tower domain. Disease-causing mutations in the Tower domain are observed in several FA patients. Our work reveals that recruitment of the complex to a stalled replication fork serves as the trigger for the activating monoubiquitination event. Taken together, our results uncover the mechanism of how the FANCD2-FANCI complex activates the FA pathway, and explains the underlying molecular defect in FA patients with mutations in the Tower domain.

Structural and Functional Highlights of Vacuolar Soluble Protein 1 from Pathogen Trypanosoma brucei brucei.

Trypanosoma brucei (T. brucei) is responsible for the fatal human disease called African trypanosomiasis, or sleeping sickness. The causative parasite, Trypanosoma, encodes soluble versions of inorganic pyrophosphatases (PPase), also called vacuolar soluble proteins (VSPs), which are localized to its acidocalcisomes. The latter are acidic membrane-enclosed organelles rich in polyphosphate chains and divalent cations whose significance in these parasites remains unclear. We here report the crystal structure of T. brucei brucei acidocalcisomal PPases in a ternary complex with Mg(2+) and imidodiphosphate. The crystal structure reveals a novel structural architecture distinct from known class I PPases in its tetrameric oligomeric state in which a fused EF hand domain arranges around the catalytic PPase domain. This unprecedented assembly evident from TbbVSP1 crystal structure is further confirmed by SAXS and TEM data. SAXS data suggest structural flexibility in EF hand domains indicative of conformational plasticity within TbbVSP1.

Mapping of the epitopes of poliovirus type 2 in complex with antibodies.

The inactivated polio vaccine (IPV) contains poliovirus (PV) samples that belong to serotypes 1, 2 and 3. All three serotypes contain the D-antigen, which induces protective antibodies. The antigenic structure of PVs consists of at least four different antigenic sites and the D-antigen content represents the combined activity of multiple epitopes (Ferguson et al., 1993; Minor, 1990; Minor et al., 1986). The potency of IPV vaccines is determined by measuring the D-antigen content. Several ELISA methods have been developed using polyclonal or monoclonal antibodies (Mabs) in order to quantify the D-antigen content. Characterization of the epitopes recognized by the different Mabs is crucial to map the entire virus surface and ensure the presence of epitopes able to induce neutralizing antibodies. Using a new approach that we developed to study the interaction between monoclonal antibodies and poliovirus type 2, which combines cryo-electron microscopy, image analysis and X-ray crystallography along with identification of exposed amino acids, we have mapped in 3D the epitope sites recognized by three specific Fabs at the surface of poliovirus type 2 (PV2) and characterized precisely the antigenic sites for these Fabs.

Top-down mass spectrometry of intact membrane protein complexes reveals oligomeric state and sequence information in a single experiment.

Here we study the intact stoichiometry and top-down fragmentation behavior of three integral membrane proteins which were natively reconstituted into detergent micelles: the mechano-sensitive ion channel of large conductance (MscL), the Kirbac potassium channel and the p7 viroporin from the hepatitis C virus. By releasing the proteins under nondenaturing conditions inside the mass spectrometer, we obtained their oligomeric sizes. Increasing the ion activation (collision energy) causes unfolding and subsequent ejection of a highly charged monomer from the membrane protein complexes. Further increase of the ion activation then causes collision-induced dissociation (CID) of the ejected monomers, with fragments observed which were predominantly found to stem from membrane-embedded regions. These experiments show how in a single experiment, we can probe the relation between higher-order structure and protein sequence, by combining the native MS data with fragmentation obtained from top-down MS.

Structural studies of virus-antibody immune complexes (poliovirus type I): Characterization of the epitopes in 3D.

The inactivated polio vaccine (IPV) contains poliovirus (PVs) samples that belong to serotypes 1, 2 and 3. All three serotypes contain the D-antigen, which induces protective antibodies. The antigenic structure of PVs consists of at least four different antigenic sites and the D-antigen content represents the combined activity of multiple epitopes (Ferguson et al., 1993; Minor, 1990; Minor et al., 1986). The potency of IPV vaccines is determined by measuring the D-antigen content. Several ELISA methods have been developed using polyclonal or monoclonal antibodies (Mabs) in order to quantify the D-antigen content. Characterization of the epitopes recognized by the different Mabs is crucial to map the entire virus surface and ensure the presence of epitopes able to induce neutralizing antibodies. In a new approach, combining cryo-electron microscopy and image analysis with X-ray crystallography data available along with identification of exposed amino acids we have mapped in 3D the epitope sites recognized by five specific Fabs and one Mab and characterized precisely the antigenic sites for these Mabs. We propose this method to be used to map the entire "epitopic" surface of virus.