Stability Prediction for Mutations in the Cytosolic Domains of Cystic Fibrosis Transmembrane Conductance Regulator.

Cystic Fibrosis (CF) is caused by mutations to the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) chloride channel. CFTR is composed of two membrane spanning domains, two cytosolic nucleotide-binding domains (NBD1 and NBD2) and a largely unstructured R-domain. Multiple CF-causing mutations reside in the NBDs and some are known to compromise the stability of these domains. The ability to predict the effect of mutations on the stability of the cytosolic domains of CFTR and to shed light on the mechanisms by which they exert their effect is therefore important in CF research. With this in mind, we have predicted the effect on domain stability of 59 mutations in NBD1 and NBD2 using 15 different algorithms and evaluated their performances via comparison to experimental data using several metrics including the correct classification rate (CCR), and the squared Pearson correlation (R2) and Spearman's correlation (ρ) calculated between the experimental ΔTm values and the computationally predicted ΔΔG values. Overall, the best results were obtained with FoldX and Rosetta. For NBD1 (35 mutations), FoldX provided R2 and ρ values of 0.64 and -0.71, respectively, with an 86% correct classification rate (CCR). For NBD2 (24 mutations), FoldX R2, ρ, and CCR were 0.51, -0.73, and 75%, respectively. Application of the Rosetta high-resolution protocol (Rosetta_hrp) to NBD1 yielded R2, ρ, and CCR of 0.64, -0.75, and 69%, respectively, and for NBD2 yielded R2, ρ, and CCR of 0.29, -0.27, and 50%, respectively. The corresponding numbers for the Rosetta's low-resolution protocol (Rosetta_lrp) were R2 = 0.47, ρ = -0.69, and CCR = 69% for NBD1 and R2 = 0.27, ρ = -0.24, and CCR = 63% for NBD2. For NBD1, both algorithms suggest that destabilizing mutations suffer from destabilizing vdW clashes, whereas stabilizing mutations benefit from favorable H-bond interactions. Two triple consensus approaches based on FoldX, Rosetta_lpr, and Rosetta_hpr were attempted using either "majority-voting" or "all-voting". The all-voting consensus outperformed the individual predictors, albeit on a smaller data set. In summary, our results suggest that the effect of mutations on the stability of CFTR's NBDs could be largely predicted. Since NBDs are common to all ABC transporters, these results may find use in predicting the effect and mechanism of the action of multiple disease-causing mutations in other proteins.

MeCP2 epigenetically regulates alpha-smooth muscle actin in human lung fibroblasts.

BACKGROUND: A critical feature for fibroblasts differentiation into myofibroblasts is the expression of alpha-smooth muscle actin (α-SMA) during the tissue injury and repair process. The epigenetic mechanism, DNA methylation, is involved in regulating α-SMA expression. It is not clear how methyl-CpG-binding protein 2 (MeCP2) interacts with CpG-rich region in α-SMA, and if the CpG methylation status would affect MeCP2 binding and regulation of α-SMA expression. METHODS: The association of MeCP2 with α-SMA CpG rich region were examined by chromatin immunoprecipitation (ChIP) assays in primary fibroblasts from idiopathic pulmonary fibrosis (IPF) and non-IPF control individuals, and in the lung fibroblasts treated with profibrotic cytokine transforming growth factor ß1 (TGF-ß1). The regulation of α-SMA by MeCP2 was examined by knocking down MeCP2 with small interfering RNA (siRNA). To explore the effects of the DNA methylation status of the CpG rich region on α-SMA expression, the cells were treated with DNA methyltransferase inhibitor, 5'-azacytidine (5'-aza). The expression of α-SMA was examined by Western blot and quantitative polymerase chain reaction, the association with MeCP2 was assessed by ChIP assays, and the methylation status was checked by bisulfate sequencing. RESULTS: The human lung fibroblasts with increased α-SMA showed an enriched association of MeCP2, while knockdown MeCP2 by siRNA reduced α-SMA upregulation by TGF-ß1. The 5'-Aza-treated cells have decreased α-SMA expression with reduced MeCP2 association. However, bisulfite sequencing revealed that most CpG sites are unmethylated despite the different expression levels of α-SMA after being treated by TGF-ß1 or 5'-aza. CONCLUSION: Our data indicate that the methyl-binding protein MeCP2 is critical for α-SMA expression in human lung myofibroblast, and the DNA methylation status at the CpG rich region of α-SMA is not a determinative factor for its inducible expression.

Structural stability of purified human CFTR is systematically improved by mutations in nucleotide binding domain 1.

The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is an ABC transporter containing two transmembrane domains forming a chloride ion channel, and two nucleotide binding domains (NBD1 and NBD2). CFTR has presented a formidable challenge to obtain monodisperse, biophysically stable protein. Here we report a comprehensive study comparing effects of single and multiple NBD1 mutations on stability of both the NBD1 domain alone and on purified full length human CFTR. Single mutations S492P, A534P, I539T acted additively, and when combined with M470V, S495P, and R555K cumulatively yielded an NBD1 with highly improved structural stability. Strategic combinations of these mutations strongly stabilized the domain to attain a calorimetric Tm > 70 °C. Replica exchange molecular dynamics simulations on the most stable 6SS-NBD1 variant implicated fluctuations, electrostatic interactions and side chain packing as potential contributors to improved stability. Progressive stabilization of NBD1 directly correlated with enhanced structural stability of full-length CFTR protein. Thermal unfolding of the stabilized CFTR mutants, monitored by changes in intrinsic fluorescence, demonstrated that Tm could be shifted as high as 67.4 °C in 6SS-CFTR, more than 20 °C higher than wild-type. H1402S, an NBD2 mutation, conferred CFTR with additional thermal stability, possibly by stabilizing an NBD-dimerized conformation. CFTR variants with NBD1-stabilizing mutations were expressed at the cell surface in mammalian cells, exhibited ATPase and channel activity, and retained these functions to higher temperatures. The capability to produce enzymatically active CFTR with improved structural stability amenable to biophysical and structural studies will advance mechanistic investigations and future cystic fibrosis drug development.

Stabilization of a nucleotide-binding domain of the cystic fibrosis transmembrane conductance regulator yields insight into disease-causing mutations.

Characterization of the second nucleotide-binding domain (NBD2) of the cystic fibrosis transmembrane conductance regulator (CFTR) has lagged behind research into the NBD1 domain, in part because NBD1 contains the F508del mutation, which is the dominant cause of cystic fibrosis. Research on NBD2 has also been hampered by the overall instability of the domain and the difficulty of producing reagents. Nonetheless, multiple disease-causing mutations reside in NBD2, and the domain is critical for CFTR function, because channel gating involves NBD1/NBD2 dimerization, and NBD2 contains the catalytically active ATPase site in CFTR. Recognizing the paucity of structural and biophysical data on NBD2, here we have defined a bioinformatics-based method for manually identifying stabilizing substitutions in NBD2, and we used an iterative process of screening single substitutions against thermal melting points to both produce minimally mutated stable constructs and individually characterize mutations. We present a range of stable constructs with minimal mutations to help inform further research on NBD2. We have used this stabilized background to study the effects of NBD2 mutations identified in cystic fibrosis (CF) patients, demonstrating that mutants such as N1303K and G1349D are characterized by lower stability, as shown previously for some NBD1 mutations, suggesting a potential role for NBD2 instability in the pathology of CF.

Interactions and cooperativity between P-glycoprotein structural domains determined by thermal unfolding provides insights into its solution structure and function.

Structural changes in mouse P-glycoprotein (Pgp) induced by thermal unfolding were studied by differential scanning calorimetry (DSC), circular dichroism and fluorescence spectroscopy to gain insight into the solution conformation(s) of this ABC transporter that may not be apparent from current crystal structures. DSC of reconstituted Pgp showed two thermal unfolding transitions in the absence of MgATP, suggesting that each transition involved the cooperative unfolding of two or more interacting structural domains. A low calorimetric unfolding enthalpy and minimal structural changes were observed, which are hallmarks of the thermal unfolding of α-helical membrane proteins, because generally only the extramembranous regions undergo significant unfolding. Nucleotide binding increased the unfolding temperature of both transitions to the same extent, suggesting that one nucleotide binding domain (NBD) unfolds with each transition. Combined with the results from the two isolated NBDs, we propose that each DSC transition represents the cooperative unfolding of one NBD and the two contacting intracellular loops. Further, the presence of two transitions in both apo and MgATP bound wild-type Pgp suggests the NBD-dimeric conformation is transient, and that Pgp resides predominantly in the crystallographically observed inward-facing conformation with NBDs separated, even under conditions supporting continuous MgATP hydrolysis. In contrast, DSC of the vanadate-trapped MgADP·Pgp complex and the MgATP-bound catalytically inactive mutant, E552A/E1197A, show an additional transition at much higher temperature, corresponding to the unfolding of the nucleotide-trapped NBD-dimeric outward-facing conformation. The collective results indicate a strong preference for an NBD dissociated, inward-facing conformation of Pgp.

Membrane protein stability can be compromised by detergent interactions with the extramembranous soluble domains.

Detergent interaction with extramembranous soluble domains (ESDs) is not commonly considered an important determinant of integral membrane protein (IMP) behavior during purification and crystallization, even though ESDs contribute to the stability of many IMPs. Here we demonstrate that some generally nondenaturing detergents critically destabilize a model ESD, the first nucleotide-binding domain (NBD1) from the human cystic fibrosis transmembrane conductance regulator (CFTR), a model IMP. Notably, the detergents show equivalent trends in their influence on the stability of isolated NBD1 and full-length CFTR. We used differential scanning calorimetry (DSC) and circular dichroism (CD) spectroscopy to monitor changes in NBD1 stability and secondary structure, respectively, during titration with a series of detergents. Their effective harshness in these assays mirrors that widely accepted for their interaction with IMPs, i.e., anionic > zwitterionic > nonionic. It is noteworthy that including lipids or nonionic detergents is shown to mitigate detergent harshness, as will limiting contact time. We infer three thermodynamic mechanisms from the observed thermal destabilization by monomer or micelle: (i) binding to the unfolded state with no change in the native structure (all detergent classes); (ii) native state binding that alters thermodynamic properties and perhaps conformation (nonionic detergents); and (iii) detergent binding that directly leads to denaturation of the native state (anionic and zwitterionic). These results demonstrate that the accepted model for the harshness of detergents applies to their interaction with an ESD. It is concluded that destabilization of extramembranous soluble domains by specific detergents will influence the stability of some IMPs during purification.

SAR studies for a new class of antibacterial NAD biosynthesis inhibitors.

A new lead class of antibacterial drug-like NAD synthetase (NADs) inhibitors was previously identified from a virtual screening study. Here a solution-phase synthetic library of 76 compounds, analogs of the urea-sulfonamide 5838, was synthesized in parallel to explore SAR on the sulfonamide aryl group. All library members were tested for enzyme inhibition against NADs and nicotinic acid mononucleotide adenylyltransferase (NaMNAT), the last two enzymes in the biosynthesis of NAD, and for growth inhibition in a Bacillus anthracis antibacterial assay. Most compounds that inhibited bacterial growth also showed inhibition against one of the enzymes tested. While only modest enhancements in the enzyme inhibition potency against NADs were observed, of significance was the observation that the antibacterial urea-sulfonamides more consistently inhibited NaMNAT.

Purification, crystallization and preliminary X-ray crystallographic analysis of the nucleocapsid protein of Bunyamwera virus.

Bunyamwera virus (BUNV) is the prototypic member of the Bunyaviridae family of segmented negative-sense RNA viruses. The BUNV nucleocapsid protein has been cloned and expressed in Escherichia coli. The purified protein has been crystallized and a complete data set has been collected to 3.3 angstroms resolution at a synchrotron source. Crystals of the nucleocapsid protein belong to space group C2, with unit-cell parameters a = 384.7, b = 89.8, c = 89.2 angstroms, beta = 94.4 degrees . Self-rotation function analysis of the X-ray diffraction data has provided insight into the oligomeric state of the protein as well as the orientation of the oligomers in the asymmetric unit. The structure determination of the protein is ongoing.

Structural analysis of Pseudomonas 1-aminocyclopropane-1-carboxylate deaminase complexes: insight into the mechanism of a unique pyridoxal-5'-phosphate dependent cyclopropane ring-opening reaction.

1-Aminocyclopropane-1-carboxylate (ACC) deaminase is a pyridoxal 5'-phosphate (PLP) dependent enzyme catalyzing the opening of the cyclopropane ring of ACC to give alpha-ketobutyric acid and ammonia as the products. This ring cleavage reaction is unusual because the substrate, ACC, contains no abstractable alpha-proton and the carboxyl group is retained in the product. How the reaction is initiated to generate an alpha-carbanionic intermediate, which is the common entry for most PLP-dependent reactions, is not obvious. To gain insight into this unusual ring-opening reaction, we have solved the crystal structures of ACC deaminase from Pseudomonas sp. ACP in complex with substrate ACC, an inhibitor, 1-aminocyclopropane-1-phosphonate (ACP), the product alpha-ketobutyrate, and two d-amino acids. Several notable observations of these structural studies include the following: (1) a typically elusive gem-diamine intermediate is trapped in the enzyme complex with ACC or ACP; (2) Tyr294 is in close proximity (3.0 A) to the pro-S methylene carbon of ACC in the gem-diamine complexes, implicating a direct role of this residue in the ring-opening reaction; (3) Tyr294 may also be responsible for the abstraction of the alpha-proton from d-amino acids, a prelude to the subsequent deamination reaction; (4) the steric hindrance precludes accessibility of active site functional groups to the l-amino acid substrates and may account for the stereospecificity of this enzyme toward d-amino acids. These structural data provide evidence favoring a mechanism in which the ring cleavage is induced by a nucleophilic attack at the pro-S beta-methylene carbon of ACC, with Tyr294 as the nucleophile. However, these observations are also consistent with an alternative mechanistic possibility in which the ring opening is acid-catalyzed and may be facilitated by charge relay through PLP, where Tyr294 functions as a general acid. The results of mutagenesis studies corroborated the assigned critical role for Tyr294 in the catalysis.

Ligand binding-induced conformational changes in riboflavin kinase: structural basis for the ordered mechanism.

Riboflavin kinase (RFK) is an essential enzyme catalyzing the phosphorylation of riboflavin (vitamin B(2)) in the presence of ATP and Mg(2+) to form the active cofactor FMN, which can be further converted to FAD. Previously, the crystal structures of RFKs from human and Schizosaccharomyces pombe have been determined in the apo form and in complex with MgADP. These structures revealed that RFK adopts a novel kinase fold and utilizes a unique nucleotide binding site. The structures of the flavin-bound RFK obtained by soaking pre-existing crystals were also reported. Because of crystal packing restraints, these flavin-bound RFK complexes adopt conformations nearly identical with that of corresponding flavin-free structures. Here we report the structure of human RFK cocrystallized with both MgADP and FMN. Drastic conformational changes associated with flavin binding are observed primarily at the so-called Flap I and Flap II loop regions. As a result, the bound FMN molecule now interacts with the enzyme extensively and is well-ordered. Residues from Flap II interact with Flap I and shield the FMN molecule from the solvent. The conformational changes in Flap I resulted in a new Mg(2+) coordination pattern in which a FMN phosphate oxygen and Asn36 side chain carbonyl are directly coordinating to the Mg(2+) ion. The proposed catalytic base Glu86 is well-positioned for activation of the O5' hydroxyl group of riboflavin for the phosphoryl transfer reaction. The structural data obtained so far on human and yeast RFK complexes provide a rationale for the ordered kinetic mechanism of RFK.

Crystal structure of human riboflavin kinase reveals a beta barrel fold and a novel active site arch.

Riboflavin kinase (RFK) is an essential enzyme catalyzing the phosphorylation of riboflavin (vitamin B(2)) to form FMN, an obligatory step in vitamin B(2) utilization and flavin cofactor synthesis. The structure of human RFK revealed a six-stranded antiparallel beta barrel core structurally similar to the riboflavin synthase/ferredoxin reductase FAD binding domain fold. The binding site of an intrinsically bound MgADP defines a novel nucleotide binding motif that encompasses a loop, a 3(10) helix, and a reverse turn followed by a short beta strand. This active site loop forms an arch with ATP and riboflavin binding at the opposite side and the phosphoryl transfer appears to occur through the hole underneath the arch. The invariant residues Asn36 and Glu86 are implicated in the catalysis.