The power stroke driven by ATP binding in CFTR as studied by molecular dynamics simulations

Single cohesin molecules generate force by two distinct mechanisms

Spatial organization of DNA is facilitated by cohesin protein complexes that move on DNA and extrude DNA loops. How cohesin works mechanistically as a molecular machine is poorly understood. Here, we measure mechanical forces generated by conformational changes in single cohesin molecules. We show that bending of SMC coiled coils is driven by random thermal fluctuations leading to a ~32 nm head-hinge displacement that resists forces up to 1 pN; ATPase head engagement occurs in a single step of ~10 nm and is driven by an ATP dependent head-head movement, resisting forces up to 15 pN. Our molecular dynamic simulations show that the energy of head engagement can be stored in a mechanically strained conformation of NIPBL and released during disengagement. These findings reveal how single cohesin molecules generate force by two distinct mechanisms. We present a model, which proposes how this ability may power different aspects of cohesin-DNA interaction.

The molecular evolution of function in the CFTR chloride channel

The ATP-binding cassette (ABC) transporter superfamily includes many proteins of clinical relevance, with genes expressed in all domains of life. Although most members use the energy of ATP binding and hydrolysis to accomplish the active import or export of various substrates across membranes, the cystic fibrosis transmembrane conductance regulator (CFTR) is the only known animal ABC transporter that functions primarily as an ion channel. Defects in CFTR, which is closely related to ABCC subfamily members that bear function as bona fide transporters, underlie the lethal genetic disease cystic fibrosis. This article seeks to integrate structural, functional, and genomic data to begin to answer the critical question of how the function of CFTR evolved to exhibit regulated channel activity. We highlight several examples wherein preexisting features in ABCC transporters were functionally leveraged as is, or altered by molecular evolution, to ultimately support channel function. This includes features that may underlie (1) construction of an anionic channel pore from an anionic substrate transport pathway, (2) establishment and tuning of phosphoregulation, and (3) optimization of channel function by specialized ligand–channel interactions. We also discuss how divergence and conservation may help elucidate the pharmacology of important CFTR modulators.

Structural dynamics of ABC transporters: molecular simulation studies

The biological activities of living organisms involve various inputs and outputs. The ATP-driven substances (biomolecules) responsible for these kinds of activities through membrane (i.e. uptake and efflux of substrates) include ATP-binding cassette (ABC) transporters, some of which play important roles in multidrug resistance. The basic architecture of ABC transporters comprises transmembrane domains (TMDs) and nucleotide-binding domains (NBDs). The functional dynamics (substrate transport) of ABC transporters are realized by concerted motions, such as NBD dimerization, mechanical transmission via coupling helices (CHs), and the translocation of substrates through TMDs, which are induced by the binding and/or hydrolysis of ATP molecules and substrates. In this mini-review, we briefly discuss recent progresses in the structural dynamics as revealed by molecular simulation studies at all-atom (AA), coarse-grained (CG), and quantum mechanics/molecular mechanics (QM/MM) levels.

Recent Strategic Advances in CFTR Drug Discovery: An Overview

Cystic fibrosis transmembrane conductance regulator (CFTR)-rescuing drugs have already transformed cystic fibrosis (CF) from a fatal disease to a treatable chronic condition. However, new-generation drugs able to bind CFTR with higher specificity/affinity and to exert stronger therapeutic benefits and fewer side effects are still awaited. Computational methods and biosensors have become indispensable tools in the process of drug discovery for many important human pathologies. Instead, they have been used only piecemeal in CF so far, calling for their appropriate integration with well-tried CF biochemical and cell-based models to speed up the discovery of new CFTR-rescuing drugs. This review will give an overview of the available structures and computational models of CFTR and of the biosensors, biochemical and cell-based assays already used in CF-oriented studies. It will also give the reader some insights about how to integrate these tools as to improve the efficiency of the drug discovery process targeted to CFTR.

Development of an ATP force field for coarse-grained simulation of ATPases and its application to the maltose transporter

The biological functions of ATPases, such as myosin, kinesin, and ABC transporter, are due to large conformational motions driven by energy obtained from ATP. Elucidation of the mechanisms underlying these ATP-driven movements is one of the greatest challenges in computational chemistry. It has been shown that the MARTINI coarse-grained method is a promising tool for the investigation of large conformational motions in various proteins. However, this method has not yet been applied to ATPases because of the lack of a force field for the ATP molecule. Here, we developed force field parameters for the ATP molecule and conducted simulations using these parameters for the subunits (MalK₂ ) and the full-length structure (MalFGK₂ -E) of a maltose transporter. It was found for both targets that the dimerization of the nucleotide binding domains (NBDs) is induced upon ATP binding. Moreover, for the full-length transporter, the conformational transition from the pre-translocation state to the outward-facing state was observed and was accompanied by an initial transport motion of the substrate. It is expected that coarse-grained simulations utilizing the parameters for the ATP molecule developed here will serve as a powerful tool for investigating other ATPases as well. © 2019 Wiley Periodicals, Inc.

ATP-Dependent Signaling in Simulations of a Revised Model of Cystic Fibrosis Transmembrane Conductance Regulator (CFTR)

Cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ATP-binding cassette (ABC) transporter superfamily that has uniquely evolved to function as a chloride channel. It binds and hydrolyzes ATP at its nucleotide binding domains to form a pore providing a diffusive pathway within its transmembrane domains. CFTR is the only known protein from the ABC superfamily with channel activity, and its dysfunction causes the disease cystic fibrosis. While much is known about the functional aspects of CFTR, significant gaps remain, such as the structure-function relationship underlying signaling of ATP binding. In the present work, we refined an existing homology model using an intermediate-resolution (9 Å) published cryo-electron microscopy map. The newly derived models have been simulated in equilibrium molecular dynamics simulations for a total of 2.5 μs in multiple ATP-occupancy states. Putative conformational movements connecting ATP binding with pore formation are elucidated and quantified. Additionally, new interdomain interactions between E543, K968, and K1292 have been identified and confirmed experimentally; these interactions may be relevant for signaling ATP binding and hydrolysis to the transmembrane domains and induction of pore opening.

Molecular dynamics simulation study on the structural instability of the most common cystic fibrosis-associated mutant ΔF508-CFTR

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel that belongs to the ATP binding cassette protein superfamily. Deletion of phenylalanine at position 508 (ΔF508) is the most common CF-associated mutation and is present in nearly 90% of CF patients. Currently, atomistic level studies are insufficient for understanding the mechanism by which the deletion of a single amino acid causes greatly reduced folding as well as trafficking and gating defects. To clarify this mechanism, we first constructed an atomic model of the inward-facing ΔF508-CFTR and performed allatom molecular dynamics (MD) simulations of the protein in a membrane environment. All of the computational methodologies used are based on those developed in our previous study for wild-type CFTR. Two important findings were obtained. First, consistent with several previous computational results, the deletion of F508 causes a disruption of a hydrophobic cluster located at the interface between the nucleotide binding domain 1 (NBD1) and intracellular loop 4 (ICL4). This exerts unfavorable influences on the correlated motion between ICLs and transmembrane domains (TMDs), likely resulting in gating defects. Second, the F508 deletion affected the NBD1-NBD2 interface via allosteric communication originating from the correlated motion between NBDs and ICLs. As a result, several unusual inter-residue interactions are caused at the NBD1-NBD2 interface. In other words, their correct dimerization is impaired. This study provided insight into the atomic-level details of structural and dynamics changes caused by the ΔF508 mutation and thus provides good insight for drug design.

Transmembrane helical interactions in the CFTR channel pore

Mutations in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene affect CFTR protein biogenesis or its function as a chloride channel, resulting in dysregulation of epithelial fluid transport in the lung, pancreas and other organs in cystic fibrosis (CF). Development of pharmaceutical strategies to treat CF requires understanding of the mechanisms underlying channel function. However, incomplete 3D structural information on the unique ABC ion channel, CFTR, hinders elucidation of its functional mechanism and correction of cystic fibrosis causing mutants. Several CFTR homology models have been developed using bacterial ABC transporters as templates but these have low sequence similarity to CFTR and are not ion channels. Here, we refine an earlier model in an outward (OWF) and develop an inward (IWF) facing model employing an integrated experimental-molecular dynamics simulation (200 ns) approach. Our IWF structure agrees well with a recently solved cryo-EM structure of a CFTR IWF state. We utilize cysteine cross-linking to verify positions and orientations of residues within trans-membrane helices (TMHs) of the OWF conformation and to reconstruct a physiologically relevant pore structure. Comparison of pore profiles of the two conformations reveal a radius sufficient to permit passage of hydrated Cl- ions in the OWF but not the IWF model. To identify structural determinants that distinguish the two conformations and possible rearrangements of TMHs within them responsible for channel gating, we perform cross-linking by bifunctional reagents of multiple predicted pairs of cysteines in TMH 6 and 12 and 6 and 9. To determine whether the effects of cross-linking on gating observed are the result of switching of the channel from open to close state, we also treat the same residue pairs with monofunctional reagents in separate experiments. Both types of reagents prevent ion currents indicating that pore blockage is primarily responsible.

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.

Maltose-binding protein effectively stabilizes the partially closed conformation of the ATP-binding cassette transporter MalFGK2

Maltose transporter MalFGK₂is a type-I importer in the ATP-binding cassette (ABC) transporter superfamily.

Conformational Dynamics and Protein-Substrate Interaction of ABC Transporter BtuCD at the Occluded State Revealed by Molecular Dynamics Simulations

ATP-binding cassette (ABC) transporters are ubiquitous in all three kingdoms of life and are implicated in many clinically relevant physiological processes. They couple the energy released by ATP hydrolysis to facilitate substrate translocation across cell membranes. The crystal structures of type II ABC importers have revealed their unique transmembrane domain architecture consisting of 10 transmembrane helices and their structurally conserved nucleotide-binding domains among all ABC transporters. However, molecular details of the interactions between the importers and their substrate remain largely elusive. Taking vitamin B₁₂ importer BtuCD as an exemplar of type II importers, we investigated the dynamics of its occluded state and the detailed protein-substrate interactions using molecular dynamics simulation. Our trajectories show that the importer accommodates the substrate through a nonspecific binding mode as the substrate undergoes evident vertical and tilt motions inside the translocation cavity. Extensive hydrogen bond and hydrophobic interactions were observed between the substrate and the importer; however, most of these interactions are weak, with <38% occurrence. The presence of substrate leads to enlargement of the translocation cavity, especially at its cytoplasmic end, which may activate cytoplasmic regions and probably facilitate the transportation. The perturbations caused by periplasmic binding protein and nucleotides were also investigated. The study provides deeper insight into the translocation mechanism of BtuCD.

Structural Dynamics of the Heterodimeric ABC Transporter TM287/288 Induced by ATP and Substrate Binding

TM287/288 is a heterodimeric ATP-binding cassette (ABC) transporter, which harnesses the energy of ATP binding and hydrolysis at the nucleotide-binding domains (NBDs) to transport a wide variety of molecules through the transmembrane domains (TMDs) by alternating inward- and outward-facing conformations. Here, we conducted multiple 100 ns molecular dynamics simulations of TM287/288 in different ATP- and substrate-bound states to elucidate the effects of ATP and substrate binding. As a result, the binding of two ATP molecules to the NBDs induced the formation of the consensus ATP-binding pocket (ABP2) or the NBD dimerization, whereas these processes did not occur in the presence of a single ATP molecule or when the protein was in its apo state. Moreover, binding of the substrate to the TMDs enhanced the formation of ABP2 through allosteric TMD-NBD communication. Furthermore, in the apo state, α-helical subdomains of the NBDs approached each other, acquiring a conformation with core half-pockets exposed to the solvent, appropriate for ATP binding. We propose a "core-exposed" model for this novel conformation found in the apo state of ABC transporters. These findings provide important insights into the structural dynamics of ABC transporters.

Functional characteristics of L1156F-CFTR associated with alcoholic chronic pancreatitis in Japanese

Although cystic fibrosis is rare in Japanese, measurement of sweat Cl(-) has suggested mild dysfunction of cystic fibrosis transmembrane conductance regulator (CFTR) in some patients with chronic pancreatitis. In the present study, we have investigated the association of CFTR variants and chronic pancreatitis in Japanese and the functional characteristics of a Japanese- and pancreatitis-specific CFTR variant, L1156F. Seventy patients with alcoholic chronic pancreatitis, 18 patients with idiopathic chronic pancreatitis, and 180 normal subjects participated. All exons and their boundaries and promoter region of the CFTR gene were sequenced. Human embryonic kidney-293 cells were transfected with three CFTR variants (M470V, L1156F, and M470V+L1156F), and the protein expression was examined. Xenopus laevis oocytes were injected with the CFTR variants, and bicarbonate (HCO3 (-)) transport activity was examined. CFPAC-1 cells were transfected with the CFTR variants and Cl(-)/HCO3 (-) exchange activity was examined. Six variants (E217G, I556V, M470V, L1156F, Q1352H, and R1453W) were identified in the coding region of the CFTR gene. Cystic fibrosis-causing mutations were not found. The allele frequencies of L1156F and Q1352H in alcoholic chronic pancreatitis (5.0 and 7.9%) were significantly (P < 0.01) higher than those in normal subjects (0.6 and 1.9%). L1156F was linked with a worldwide CFTR variant, M470V. Combination of M470V and L1156F significantly reduced CFTR expression to ∼60%, impaired CFTR-mediated HCO3 (-)/Cl(-) transport activity to 50-60%, and impaired CFTR-coupled Cl(-)/HCO3 (-) exchange activity to 20-30%. The data suggest that the Japanese-specific CFTR variant L1156F causes mild dysfunction of CFTR and increases the risk of alcoholic chronic pancreatitis in Japanese.

Analysis of the Free Energy Landscapes for the Opening-Closing Dynamics of the Maltose Transporter ATPase MalK2 Using Enhanced-Sampling Molecular Dynamics Simulation

Protein dynamics are considered significant for many physiological processes, such as metabolism, biomolecular recognition, and the regulation of several vital cellular processes. Due to their flexibility, proteins may stay in different substates with or without the existence of the cognate substrates. To describe these phenomena, two models have been proposed: the "induced fit" and the "conformational selection" mechanisms. In this study, we used MalK2, the subunits that mainly include the nucleotide-binding domains (NBDs) of the maltose transporter from Escherichia coli, as a target to understand the NBD dimerization mechanism. Accelerated and conventional molecular dynamics have been performed. The results revealed that Mg-ATP binding to MalK2 led to a significant change in the free energy profile and thus stabilized the closed conformation. On the contrary, when Mg-ATP was removed, the open conformation would be favored. The fact that ligand binding induces a drastic free energy change leads to a significant inference: MalK2 dimerization would occur through the induced-fit mechanism rather than the conformational selection mechanism. This study sheds new light on the NBD dimerization mechanism and would be of wide applicability to other ABC transporters.

The Nucleotide-Free State of the Multidrug Resistance ABC Transporter LmrA: Sulfhydryl Cross-Linking Supports a Constant Contact, Head-to-Tail Configuration of the Nucleotide-Binding Domains

ABC transporters are integral membrane pumps that are responsible for the import or export of a diverse range of molecules across cell membranes. ABC transporters have been implicated in many phenomena of medical importance, including cystic fibrosis and multidrug resistance in humans. The molecular architecture of ABC transporters comprises two transmembrane domains and two ATP-binding cassettes, or nucleotide-binding domains (NBDs), which are highly conserved and contain motifs that are crucial to ATP binding and hydrolysis. Despite the improved clarity of recent structural, biophysical, and biochemical data, the seemingly simple process of ATP binding and hydrolysis remains controversial, with a major unresolved issue being whether the NBD protomers separate during the catalytic cycle. Here chemical cross-linking data is presented for the bacterial ABC multidrug resistance (MDR) transporter LmrA. These indicate that in the absence of nucleotide or substrate, the NBDs come into contact to a significant extent, even at 4°C, where ATPase activity is abrogated. The data are clearly not in accord with an inward-closed conformation akin to that observed in a crystal structure of V. cholerae MsbA. Rather, they suggest a head-to-tail configuration ‘sandwich’ dimer similar to that observed in crystal structures of nucleotide-bound ABC NBDs. We argue the data are more readily reconciled with the notion that the NBDs are in proximity while undergoing intra-domain motions, than with an NBD ‘Switch’ mechanism in which the NBD monomers separate in between ATP hydrolysis cycles.

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.

ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force

ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.

Mechanisms of CFTR functional variants that impair regulated bicarbonate permeation and increase risk for pancreatitis but not for cystic fibrosis

CFTR is a dynamically regulated anion channel. Intracellular WNK1-SPAK activation causes CFTR to change permeability and conductance characteristics from a chloride-preferring to bicarbonate-preferring channel through unknown mechanisms. Two severe CFTR mutations (CFTRsev) cause complete loss of CFTR function and result in cystic fibrosis (CF), a severe genetic disorder affecting sweat glands, nasal sinuses, lungs, pancreas, liver, intestines, and male reproductive system. We hypothesize that those CFTR mutations that disrupt the WNK1-SPAK activation mechanisms cause a selective, bicarbonate defect in channel function (CFTRBD) affecting organs that utilize CFTR for bicarbonate secretion (e.g. the pancreas, nasal sinus, vas deferens) but do not cause typical CF. To understand the structural and functional requirements of the CFTR bicarbonate-preferring channel, we (a) screened 984 well-phenotyped pancreatitis cases for candidate CFTRBD mutations from among 81 previously described CFTR variants; (b) conducted electrophysiology studies on clones of variants found in pancreatitis but not CF; (c) computationally constructed a new, complete structural model of CFTR for molecular dynamics simulation of wild-type and mutant variants; and (d) tested the newly defined CFTRBD variants for disease in non-pancreas organs utilizing CFTR for bicarbonate secretion. Nine variants (CFTR R74Q, R75Q, R117H, R170H, L967S, L997F, D1152H, S1235R, and D1270N) not associated with typical CF were associated with pancreatitis (OR 1.5, p = 0.002). Clones expressed in HEK 293T cells had normal chloride but not bicarbonate permeability and conductance with WNK1-SPAK activation. Molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation. Comparing pancreatitis patients and controls, CFTRBD increased risk for rhinosinusitis (OR 2.3, p<0.005) and male infertility (OR 395, p<<0.0001). WNK1-SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTRBD variants through multiple mechanisms. CFTRBD variants are associated with clinically significant disorders of the pancreas, sinuses, and male reproductive system.

Modeling the conformational changes underlying channel opening in CFTR

Mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator protein (CFTR) cause cystic fibrosis (CF), the most common life-shortening genetic disease among Caucasians. Although general features of the structure of CFTR have been predicted from homology models, the conformational changes that result in channel opening and closing have yet to be resolved. We created new closed- and open-state homology models of CFTR, and performed targeted molecular dynamics simulations of the conformational transitions in a channel opening event. The simulations predict a conformational wave that starts at the nucleotide binding domains and ends with the formation of an open conduction pathway. Changes in side-chain interactions are observed in all major domains of the protein, and experimental confirmation was obtained for a novel intra-protein salt bridge that breaks near the end of the transition. The models and simulation add to our understanding of the mechanism of ATP-dependent gating in this disease-relevant ion channel.

CFTR: effect of ICL2 and ICL4 amino acids in close spatial proximity on the current properties of the channel

BACKGROUND

CFTR is the only ABC transporter functioning as a chloride (Cl(-)) channel. We studied molecular determinants, which might distinguish CFTR from standard ABC transporters, and focused on the interface formed by the intracellular loops from the membrane spanning domains.

METHODS

Residues from ICL2 and ICL4 in close proximity were targeted, and their involvement in the functioning of CFTR was studied by whole cell patch clamp recording.

RESULTS

We identified 2 pairs of amino acids, at the extremity of the bundle formed by the four intracellular loops, whose mutation i) decreases the Cl(-) current of CFTR (couple E267-K1060) or ii) increases it with a change of the electrophysiological signature (couple S263-V1056).

CONCLUSIONS

These results highlight the critical role of these ICL residues in the assembly of the different domains and/or in the Cl(-) permeation pathway of CFTR.