Nanomechanics combined with HDX reveals allosteric drug binding sites of CFTR NBD1.

Cystic fibrosis (CF) is a frequent genetic disease in Caucasians that is caused by the deletion of F508 (ΔF508) in the nucleotide binding domain 1 (NBD1) of the CF transmembrane conductance regulator (CFTR). The ΔF508 compromises the folding energetics of the NBD1, as well as the folding of three other CFTR domains. Combination of FDA approved corrector molecules can efficiently but incompletely rescue the ΔF508-CFTR folding and stability defect. Thus, new pharmacophores that would reinstate the wild-type-like conformational stability of the ΔF508-NBD1 would be highly beneficial. The most prominent molecule, 5-bromoindole-3-acetic acid (BIA) that can thermally stabilize the NBD1 has low potency and efficacy. To gain insights into the NBD1 (un)folding dynamics and BIA binding site localization, we combined molecular dynamics (MD) simulations, atomic force spectroscopy (AFM) and hydrogen-deuterium exchange (HDX) experiments. We found that the NBD1 α-subdomain with three adjacent strands from the ß-subdomain plays an important role in early folding steps, when crucial non-native interactions are formed via residue F508. Our AFM and HDX experiments showed that BIA associates with this α-core region and increases the resistance of the ΔF508-NBD1 against mechanical unfolding, a phenomenon that could be exploited in future developments of folding correctors.

Ins and outs of AlphaFold2 transmembrane protein structure predictions.

Transmembrane (TM) proteins are major drug targets, but their structure determination, a prerequisite for rational drug design, remains challenging. Recently, the DeepMind's AlphaFold2 machine learning method greatly expanded the structural coverage of sequences with high accuracy. Since the employed algorithm did not take specific properties of TM proteins into account, the reliability of the generated TM structures should be assessed. Therefore, we quantitatively investigated the quality of structures at genome scales, at the level of ABC protein superfamily folds and for specific membrane proteins (e.g. dimer modeling and stability in molecular dynamics simulations). We tested template-free structure prediction with a challenging TM CASP14 target and several TM protein structures published after AlphaFold2 training. Our results suggest that AlphaFold2 performs well in the case of TM proteins and its neural network is not overfitted. We conclude that cautious applications of AlphaFold2 structural models will advance TM protein-associated studies at an unexpected level.

Using MemBlob to Analyze Transmembrane Regions Based on Cryo-EM Maps.

Transmembrane proteins include membrane channels, pores, and receptors and, as such, comprise an important part of the proteome, yet our knowledge about them is much less complete than about soluble, globular proteins. An important aspect of transmembrane protein structure is their exact position within the lipid bilayer, a feature hard to investigate experimentally at the atomic level. Here we describe MemBlob, a novel approach utilizing difference electron density maps obtained by cryo-EM studies of transmembrane proteins. The idea behind is that the nonprotein part of such maps carries information on the exact localization of the membrane mimetics used in the experiment and can be used to extract the positional information of the protein within the membrane. MemBlob uses a structural model of the protein and an experimental electron density map to provide an estimation of the surface residues interacting with the membrane.

MemBlob database and server for identifying transmembrane regions using cryo-EM maps.

SUMMARY: The identification of transmembrane helices in transmembrane proteins is crucial, not only to understand their mechanism of action but also to develop new therapies. While experimental data on the boundaries of membrane-embedded regions are sparse, this information is present in cryo-electron microscopy (cryo-EM) density maps and it has not been utilized yet for determining membrane regions. We developed a computational pipeline, where the inputs of a cryo-EM map, the corresponding atomistic structure, and the potential bilayer orientation determined by TMDET algorithm of a given protein result in an output defining the residues assigned to the bulk water phase, lipid interface and the lipid hydrophobic core. Based on this method, we built a database involving published cryo-EM protein structures and a server to be able to compute this data for newly obtained structures. AVAILABILITY AND IMPLEMENTATION: http://memblob.hegelab.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discovering the chloride pathway in the CFTR channel.

Cystic fibrosis (CF), a lethal monogenic disease, is caused by pathogenic variants of the CFTR chloride channel. The majority of CF mutations affect protein folding and stability leading overall to diminished apical anion conductance of epithelial cells. The recently published cryo-EM structures of full-length human and zebrafish CFTR provide a good model to gain insight into structure-function relationships of CFTR variants. Although, some of the structures were determined in the phosphorylated and ATP-bound active state, none of the static structures showed an open pathway for chloride permeation. Therefore, we performed molecular dynamics simulations to generate a conformational ensemble of the protein and used channel detecting algorithms to identify conformations with an opened channel. Our simulations indicate a main intracellular entry at TM4/6, a secondary pore at TM10/12, and a bottleneck region involving numerous amino acids from TM1, TM6, and TM12 in accordance with experiments. Since chloride ions entered the pathway in our equilibrium simulations, but did not traverse the bottleneck region, we performed metadynamics simulations, which revealed two possible exits. One of the chloride ions exits includes hydrophobic lipid tails that may explain the lipid-dependency of CFTR function. In summary, our in silico study provides a detailed description of a potential chloride channel pathway based on a recent cryo-EM structure and may help to understand the gating of the CFTR chloride channel, thus contributing to novel strategies to rescue dysfunctional mutants.

Quantitative comparison of ABC membrane protein type I exporter structures in a standardized way.

An increasing number of ABC membrane protein structures are determined by cryo-electron microscopy and X-ray crystallography, consequently identifying differences between their conformations has become an arising issue. Therefore, we propose to define standardized measures for ABC Type I exporter structure characterization. We set conformational vectors, conftors, which describe the relative orientation of domains and can highlight structural differences. In addition, continuum electrostatics calculations were performed to characterize the energetics of membrane insertion illuminating functionally crucial regions. In summary, the proposed metrics contribute to deeper understanding of ABC membrane proteins' structural features, structure validation, and analysis of movements observed in a molecular dynamics trajectory. Moreover, the concept of standardized metrics can be applied not only to ABC membrane protein structures (http://conftors.hegelab.org).