J. P. Mornon

Full-open and closed CFTR channels, with lateral tunnels from the cytoplasm and an alternative position of the F508 region, as revealed by molecular dynamics

In absence of experimental 3D structures, several homology models, based on ABC exporter 3D structures, have provided significant insights into the molecular mechanisms underlying the function of the cystic fibrosis transmembrane conductance regulator (CFTR) protein, a chloride channel whose defects are associated with cystic fibrosis (CF). Until now, these models, however, did not furnished much insights into the continuous way that ions could follow from the cytosol to the extracellular milieu in the open form of the channel. Here, we have built a refined model of CFTR, based on the outward-facing Sav1866 experimental 3D structure and integrating the evolutionary and structural information available today. Molecular dynamics simulations revealed significant conformational changes, resulting in a full-open channel, accessible from the cytosol through lateral tunnels displayed in the long intracellular loops (ICLs). At the same time, the region of nucleotide-binding domain 1 in contact with one of the ICLs and carrying amino acid F508, the deletion of which is the most common CF-causing mutation, was found to adopt an alternative but stable position. Then, in a second step, this first stable full-open conformation evolved toward another stable state, in which only a limited displacement of the upper part of the transmembrane helices leads to a closure of the channel, in a conformation very close to that adopted by the Atm1 ABC exporter, in an inward-facing conformation. These models, supported by experimental data, provide significant new insights into the CFTR structure–function relationships and into the possible impact of CF-causing mutations.

MANOSK: a graphics program for analyzing and modeling molecular structure and functions

Abstract A program written for the Evans and Sutherland PS300 graphic displays is described in this paper. Called MANOSK, it provides a flexible and powerful environment for displaying, manipulating and analyzing small molecules and macromolecules from atomic coordinates. Translations and rotations of up to four independent molecules are available from the dials, and screen-relative orientations of the molecules are used in all geometrical and energetical calculations. A gradual progression of functions complexity makes the program easy to use for occasional works and efficient for more sophisticated studies.

Searching for geometric molecular shape complementarity using bidimensional surface profiles

The study presented herein is a bidimensional approach to the complementarity of two molecular surfaces. From two chosen sections we have established a methodology of generating the optimal matching of two shapes. Our approach consists in describing two molecular surface sections by a shape vector (the angular profile), in finding their matching patterns by comparison of the two profiles, and in optimizing the relative locations of the two sections in two-dimensional space, using rotations and translations defined by geometric characteristics. The set of optimal configurations are successively displayed on a screen. Satisfying results have been obtained for the matching of the complex kallikreine A-trypsin pancreatic bovin 2. This efficient method could be used as a preprocessing for a tridimensional shape complementarity approach between two molecular surfaces.

Molecular modelling and molecular dynamics of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) protein is a member of the ATP-binding cassette (ABC) transporter superfamily that functions as an ATP-gated channel. Considerable progress has been made over the last years in the understanding of the molecular basis of the CFTR functions, as well as dysfunctions causing the common genetic disease cystic fibrosis (CF). This review provides a global overview of the theoretical studies that have been performed so far, especially molecular modelling and molecular dynamics (MD) simulations. A special emphasis is placed on the CFTR-specific evolution of an ABC transporter framework towards a channel function, as well as on the understanding of the effects of disease-causing mutations and their specific modulation. This in silico work should help structure-based drug discovery and design, with a view to develop CFTR-specific pharmacotherapeutic approaches for the treatment of CF in the context of precision medicine.

Order in Disorder as Observed by the “Hydrophobic Cluster Analysis” of Protein Sequences

Hydrophobic cluster analysis (HCA) is an original approach for protein sequence analysis, which provides access to the foldable repertoire of the protein universe, including yet unannotated protein segments (“dark proteome”). Foldable segments correspond to ordered regions, as well as to intrinsically disordered regions (IDRs) undergoing disorder to order transitions. In this review, how HCA can be used to give insight into this last category of foldable segments is illustrated, with examples matching known 3D structures. After reviewing the HCA principles, examples of short foldable segments are given, which often contain short linear motifs, typically matching hydrophobic clusters. These segments become ordered upon contact with partners, with secondary structure preferences generally corresponding to those observed in the 3D structures within the complexes. Such small foldable segments are sometimes larger than the segments of known 3D structures, including flanking hydrophobic clusters that may be critical for interaction specificity or regulation, as well as intervening sequences allowing fuzziness. Cases of larger conditionally disordered domains are also presented, with lower density in hydrophobic clusters than well‐folded globular domains or with exposed hydrophobic patches, which are stabilized by interaction with partners.

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.