Genome-wide association study across European and African American ancestries identifies a SNP in DNMT3B contributing to nicotine dependence.

Cigarette smoking is a leading cause of preventable mortality worldwide. Nicotine dependence, which reduces the likelihood of quitting smoking, is a heritable trait with firmly established associations with sequence variants in nicotine acetylcholine receptor genes and at other loci. To search for additional loci, we conducted a genome-wide association study (GWAS) meta-analysis of nicotine dependence, totaling 38,602 smokers (28,677 Europeans/European Americans and 9925 African Americans) across 15 studies. In this largest-ever GWAS meta-analysis for nicotine dependence and the largest-ever cross-ancestry GWAS meta-analysis for any smoking phenotype, we reconfirmed the well-known CHRNA5-CHRNA3-CHRNB4 genes and further yielded a novel association in the DNA methyltransferase gene DNMT3B. The intronic DNMT3B rs910083-C allele (frequency=44-77%) was associated with increased risk of nicotine dependence at P=3.7 × 10-8 (odds ratio (OR)=1.06 and 95% confidence interval (CI)=1.04-1.07 for severe vs mild dependence). The association was independently confirmed in the UK Biobank (N=48,931) using heavy vs never smoking as a proxy phenotype (P=3.6 × 10-4, OR=1.05, and 95% CI=1.02-1.08). Rs910083-C is also associated with increased risk of squamous cell lung carcinoma in the International Lung Cancer Consortium (N=60,586, meta-analysis P=0.0095, OR=1.05, and 95% CI=1.01-1.09). Moreover, rs910083-C was implicated as a cis-methylation quantitative trait locus (QTL) variant associated with higher DNMT3B methylation in fetal brain (N=166, P=2.3 × 10-26) and a cis-expression QTL variant associated with higher DNMT3B expression in adult cerebellum from the Genotype-Tissue Expression project (N=103, P=3.0 × 10-6) and the independent Brain eQTL Almanac (N=134, P=0.028). This novel DNMT3B cis-acting QTL variant highlights the importance of genetically influenced regulation in brain on the risks of nicotine dependence, heavy smoking and consequent lung cancer.

A rare missense mutation in CHRNA4 associates with smoking behavior and its consequences.

Using Icelandic whole-genome sequence data and an imputation approach we searched for rare sequence variants in CHRNA4 and tested them for association with nicotine dependence. We show that carriers of a rare missense variant (allele frequency=0.24%) within CHRNA4, encoding an R336C substitution, have greater risk of nicotine addiction than non-carriers as assessed by the Fagerstrom Test for Nicotine Dependence (P=1.2 × 10(-4)). The variant also confers risk of several serious smoking-related diseases previously shown to be associated with the D398N substitution in CHRNA5. We observed odds ratios (ORs) of 1.7-2.3 for lung cancer (LC; P=4.0 × 10(-4)), chronic obstructive pulmonary disease (COPD; P=9.3 × 10(-4)), peripheral artery disease (PAD; P=0.090) and abdominal aortic aneurysms (AAAs; P=0.12), and the variant associates strongly with the early-onset forms of LC (OR=4.49, P=2.2 × 10(-4)), COPD (OR=3.22, P=2.9 × 10(-4)), PAD (OR=3.47, P=9.2 × 10(-3)) and AAA (OR=6.44, P=6.3 × 10(-3)). Joint analysis of the four smoking-related diseases reveals significant association (P=6.8 × 10(-5)), particularly for early-onset cases (P=2.1 × 10(-7)). Our results are in agreement with functional studies showing that the human α4ß2 isoform of the channel containing R336C has less sensitivity for its agonists than the wild-type form following nicotine incubation.

Expanding the range of ZNF804A variants conferring risk of psychosis.

A trio of genome-wide association studies recently reported sequence variants at three loci to be significantly associated with schizophrenia. No sequence polymorphism had been unequivocally (P<5 × 10(-8)) associated with schizophrenia earlier. However, one variant, rs1344706[T], had come very close. This polymorphism, located in an intron of ZNF804A, was reported to associate with schizophrenia with a P-value of 1.6 × 10(-7), and with psychosis (schizophrenia plus bipolar disorder) with a P-value of 1.0 × 10(-8). In this study, using 5164 schizophrenia cases and 20,709 controls, we replicated the association with schizophrenia (odds ratio OR = 1.08, P = 0.0029) and, by adding bipolar disorder patients, we also confirmed the association with psychosis (added N = 609, OR = 1.09, P = 0.00065). Furthermore, as it has been proposed that variants such as rs1344706[T]-common and with low relative risk-may also serve to identify regions harboring less common, higher-risk susceptibility alleles, we searched ZNF804A for large copy number variants (CNVs) in 4235 psychosis patients, 1173 patients with other psychiatric disorders and 39,481 controls. We identified two CNVs including at least part of ZNF804A in psychosis patients and no ZNF804A CNVs in controls (P = 0.013 for association with psychosis). In addition, we found a ZNF804A CNV in an anxiety patient (P = 0.0016 for association with the larger set of psychiatric disorders).

Brain age prediction using deep learning uncovers associated sequence variants.

Machine learning algorithms can be trained to estimate age from brain structural MRI. The difference between an individual's predicted and chronological age, predicted age difference (PAD), is a phenotype of relevance to aging and brain disease. Here, we present a new deep learning approach to predict brain age from a T1-weighted MRI. The method was trained on a dataset of healthy Icelanders and tested on two datasets, IXI and UK Biobank, utilizing transfer learning to improve accuracy on new sites. A genome-wide association study (GWAS) of PAD in the UK Biobank data (discovery set: [Formula: see text], replication set: [Formula: see text]) yielded two sequence variants, rs1452628-T ([Formula: see text], [Formula: see text]) and rs2435204-G ([Formula: see text], [Formula: see text]). The former is near KCNK2 and correlates with reduced sulcal width, whereas the latter correlates with reduced white matter surface area and tags a well-known inversion at 17q21.31 (H2).

Transient channel-opening in bacteriorhodopsin: an EPR study.

Active translocation of ions across membranes requires alternating access of the ion binding site inside the pump to the two membrane surfaces. Proton translocation by bacteriorhodopsin (bR), the light-driven proton pump in Halobacterium salinarium, involves this kind of a change in the accessibility of the centrally located retinal Schiff base. This key event in bR's photocycle ensures that proton release occurs to the extracellular side and proton uptake from the cytoplasmic side. To study the role of protein conformational changes in this reprotonation switch, spin labels were attached to pairs of engineered cysteine residues in the cytoplasmic interhelical loops of bR. Light-induced changes in the distance between a spin label on the EF interhelical loop and a label on either the AB or the CD interhelical loop were observed, and the changes were monitored following photoactivation with time-resolved electron paramagnetic resonance (EPR) spectroscopy. Both distances increase transiently by about 5 A during the photocycle. This opening occurs between proton release and uptake, and may be the conformational switch that changes the accessibility of the retinal Schiff base to the cytoplasmic surface after proton release to the extracellular side.

Genome-wide meta-analysis reveals common splice site acceptor variant in CHRNA4 associated with nicotine dependence.

We conducted a 1000 Genomes-imputed genome-wide association study (GWAS) meta-analysis for nicotine dependence, defined by the Fagerström Test for Nicotine Dependence in 17 074 ever smokers from five European-ancestry samples. We followed up novel variants in 7469 ever smokers from five independent European-ancestry samples. We identified genome-wide significant association in the alpha-4 nicotinic receptor subunit (CHRNA4) gene on chromosome 20q13: lowest P=8.0 × 10(-9) across all the samples for rs2273500-C (frequency=0.15; odds ratio=1.12 and 95% confidence interval=1.08-1.17 for severe vs mild dependence). rs2273500-C, a splice site acceptor variant resulting in an alternate CHRNA4 transcript predicted to be targeted for nonsense-mediated decay, was associated with decreased CHRNA4 expression in physiologically normal human brains (lowest P=7.3 × 10(-4)). Importantly, rs2273500-C was associated with increased lung cancer risk (N=28 998, odds ratio=1.06 and 95% confidence interval=1.00-1.12), likely through its effect on smoking, as rs2273500-C was no longer associated with lung cancer after adjustment for smoking. Using criteria for smoking behavior that encompass more than the single 'cigarettes per day' item, we identified a common CHRNA4 variant with important regulatory properties that contributes to nicotine dependence and smoking-related consequences.

Photolysis of rhodopsin results in deprotonation of its retinal Schiff's base prior to formation of metarhodopsin II.

Absorption changes following photolysis of bovine rhodopsin in mildly sonicated membrane suspensions are monitored at 25 degrees C. Difference spectra collected at 17 times between 1 microsecond and 75 ms following excitation are analyzed globally using singular value decomposition and non-linear least-squares fitting techniques. The results are not consistent with the simple scheme: Lumirhodopsin-->Metarhodopsin I<-->Metarhodopsin II, but indicate that an intermediate with a deprotonated Schiff's base is formed nearly simultaneously with metarhodopsin I upon the decay of Lumirhodopsin.

A common biological basis of obesity and nicotine addiction.

Smoking influences body weight such that smokers weigh less than non-smokers and smoking cessation often leads to weight increase. The relationship between body weight and smoking is partly explained by the effect of nicotine on appetite and metabolism. However, the brain reward system is involved in the control of the intake of both food and tobacco. We evaluated the effect of single-nucleotide polymorphisms (SNPs) affecting body mass index (BMI) on smoking behavior, and tested the 32 SNPs identified in a meta-analysis for association with two smoking phenotypes, smoking initiation (SI) and the number of cigarettes smoked per day (CPD) in an Icelandic sample (N=34,216 smokers). Combined according to their effect on BMI, the SNPs correlate with both SI (r=0.019, P=0.00054) and CPD (r=0.032, P=8.0 × 10(-7)). These findings replicate in a second large data set (N=127,274, thereof 76,242 smokers) for both SI (P=1.2 × 10(-5)) and CPD (P=9.3 × 10(-5)). Notably, the variant most strongly associated with BMI (rs1558902-A in FTO) did not associate with smoking behavior. The association with smoking behavior is not due to the effect of the SNPs on BMI. Our results strongly point to a common biological basis of the regulation of our appetite for tobacco and food, and thus the vulnerability to nicotine addiction and obesity.

A limiting law for the electrostatics of the binding of polypeptides to phospholipid bilayers.

Two cysteine-substituted variants of a peptide derived from the first 25 residues of the presequence for subunit IV of cytochrome c oxidase were synthesized and modified with a nitroxide spin label. The equilibrium partitioning of these spin-labeled peptides into negatively charged phospholipid vesicles was studied with electron paramagnetic resonance (EPR) to investigate the binding energetics. It is found that the binding equilibrium constant is an explicit function of a unique variable, the membrane surface potential psi in the Gouy-Chapman-Stern theory. Moreover, at low psi (< 0.5RT/F) the binding equilibrium is described by the linear dependence of the transfer free energy delta G(el) on psi with a slope equal to the full formal charge of the peptides. However, the partition constant levels off at higher psi, suggesting departure from the ideal limiting behavior.

Topology of an amphiphilic mitochondrial signal sequence in the membrane-inserted state: a spin labeling study.

To investigate the interaction of the presequence of the precursor of yeast cytochrome C oxidase subunit IV (COX IV) with phospholipid membranes, a series of single- and double-cysteine-substituted peptide variants derived from the 25-residue NH2-terminal presequence has been synthesized and modified with nitroxide spin labels. The immersion depth, orientation, and secondary structure of the peptide in the POPC bilayer containing 10 mol % POPG were determined using electron paramagnetic resonance (EPR) spectroscopy. EPR saturation analysis of singly labeled variants reveals that the nitroxides attached to the NH2-terminal region of the peptide insert into the acyl chain region of the bilayer, approximately 13 A deep from the membrane surface. EPR line shape analysis of doubly labeled variants indicates that the peptide predominantly exists as an extended conformation, with little secondary structure. The experimental results, together with the energetic consideration of peptide-bilayer interactions, suggest that the presequence is located near the interface between the head group region and the acyl chain region, such that the hydrophobic side chains are solvated by the acyl chains and the charged side chains extended toward the polar environment at the bilayer surface.

The membrane affinities of the aliphatic amino acid side chains in an alpha-helical context are independent of membrane immersion depth.

Understanding, predicting, and designing the binding of peptides and proteins to bilayers require quantifying the intrinsic propensities of individual amino acid residues to bind membranes as a function of structural context and bilayer depth. A host-guest study was performed using the peptide host named helix5 in order to determine the membrane affinities of the aliphatic side chains both in an alpha-helical context and as a function of bilayer depth. Use of the alpha-helical host with a constrained geometry allowed the placement of guest sites at three different depths in bilayers and minimized secondary structural changes due to guest substitutions. Circular dichroism and electron paramagnetic resonance (EPR) were used to characterize the aqueous and bilayer-bound structures of the peptide variants. EPR was also used to measure the bilayer-water partition constants of the peptide variants, and the Delta DeltaGtr values (relative to Gly) of the aliphatic amino acid side chains were subsequently calculated. Surprisingly, the DeltaDeltaGtr values did not significantly vary as a function of the guest site depth in bilayers. In addition, the Delta DeltaGtr values determined in an alpha-helical context are reduced to approximately two-thirds of Delta DeltaGtr values determined in other studies for the bilayer-water and octanol-water partitioning of amino acid side chains in extended and unstructured hosts. Both the relative reduction in Delta DeltaGtr values in the context of an alpha-helical host and the invariance of Delta DeltaGtr values with respect to bilayer depth are consistent with the membrane affinities of the aliphatic residues being largely determined by the classical hydrophobic effect.

Temperature dependence of polypeptide partitioning between water and phospholipid bilayers.

Various thermodynamic forces (e.g., the hydrophobic effect, electrostatic interactions, peptide immobilization, peptide conformational changes, "bilayer effects," and van der Waals dispersion forces) can participate in the transfer of polypeptides from aqueous solution into lipid bilayers. To investigate the contributions of these forces to peptide-membrane thermodynamics, we have studied the temperature dependence of the water-bilayer partitioning of 4 polypeptides derived from the first 25 amino acid residues in the presequence of subunit IV of yeast cytochrome c oxidase (Cox IVp) using electron paramagnetic resonance spectroscopy. The partitioning of the Cox IVp peptides into phospholipid bilayers increases as the temperature is increased from 3 to 40 degrees C. The contribution of bilayer surface expansion to the temperature-dependent partitioning is estimated to be relatively small and to contribute minimally to the increased bilayer binding of the peptides with increasing temperature. Thermodynamic analysis of the data shows that the transfer of the peptides from water into bilayers at 298 K is driven by the entropic term (-T delta Str) with values ranging from -6.7 to -10 kcal mol-1, opposed by the enthalpic term (delta Htr) by approximately 4 kcal mol-1, and accompanied by a change in heat capacity (delta Cp) ranging from -117 to -208 cal K-1 mol-1. Our results indicate that while a variety of forces do, in fact, contribute to the transfer free energies (delta Gtr), the major driving force for the water-to-bilayer transfer is the hydrophobic effect.

Effects of temperature on rhodopsin photointermediates from lumirhodopsin to metarhodopsin II.

Absorbance changes following the photolysis of mildly sonicated membrane suspensions of bovine rhodopsin are monitored using multichannel detection at 15, 20, 25, 30, and 35 degrees C. Difference spectra collected with microsecond time resolution are analyzed by singular value decomposition and multiexponential fitting. Several kinetic schemes are tested using methods that compare the observed rates and associated spectral amplitudes to the eigenvalues and eigenvectors of kinetic matrices. The time evolution of the spectra is more complex than can be accounted for by the traditional lumi-->metarhodopsin I<-->metarhodopsin II scheme. Above 25 degrees C, the formation of metarhodopsin II is achieved without a large transient accumulation of metarhodopsin I. Within the framework of first-order kinetics, the observations are explained by simple kinetic schemes that lead to the formation of a deprotonated Schiff's base species temporally distinct from metarhodopsin II directly upon the decay of lumirhodopsin.

Two modes of ligand binding in maltose-binding protein of Escherichia coli. Electron paramagnetic resonance study of ligand-induced global conformational changes by site-directed spin labeling.

Binding of ligands to the maltose-binding protein (MBP) of Escherichia coli often causes a global conformational change involving the closure of its two lobes. We have introduced a cysteine residue onto each of these lobes by site-directed mutagenesis and modified these residues with spin labels. Using EPR spectroscopy, we examined the changes, caused by the ligand binding, in distance between the two spin labels, hence between the two lobes. The binding of both maltose and maltotetraose induced a considerable closure of the N- and C-terminal lobes of MBP. Little closure occurred upon the binding of maltotetraitol or beta-cyclodextrin. Previous study by fluorescence and UV differential absorbance spectroscopy (Hall, J. A., Gehring, K., and Nikaido, H. (1997) J. Biol. Chem. 272, 17605-17609) showed that maltose and a large portion of maltotetraose bound to MBP via one mode (R mode or "end-on" mode), which is physiologically active and leads to the subsequent transport of the ligands across the cytoplasmic membrane. In contrast, maltotetraitol and beta-cyclodextrin bound to MBP via a different mode (B mode or "middle" mode), which is physiologically inactive. The present work suggests that the B mode is nonproductive because ligands binding in this manner prevent the closure of the two domains of MBP, and, as a result, the resulting ligand-MBP complex is incapable of interacting properly with the inner membrane-associated transporter complex.

Direct determination of the membrane affinities of individual amino acids.

Amino acids have distinct lipid bilayer affinities which influence the insertion and topology of membrane-bound polypeptides and proteins. To measure membrane affinities, 14 uncharged amino acids were introduced individually at a guest site in a 25-residue peptide derived from the membrane-binding presequence of yeast cytochrome c oxidase, and the peptides were labeled with a nitroxide spin-label. The free energies of transfer from phospholipid bilayers to water (delta delta Gbilayer) were determined directly by examination of partitioning into phospholipid bilayers using electron paramagnetic resonance. The delta delta Gbilayer values are in agreement with hydrophobicities assessed from 1-octanol-water partitioning of N-acetyl amino acid amides [Fauchere, J.-L., & Pliska, V. (1983) Eur. J. Med. Chem. 18, 369-375; Eisenberg, D., & McLachlan, A. (1986) Nature 319, 199-203] and quantitatively demonstrate the role of the hydrophobic effect in membrane-protein interactions.

De novo design of a peptide which partitions between water and phospholipid bilayers as a monomeric alpha-helix.

To dissect the determinants of protein insertion into membranes, we designed a model peptide which partitions between water and phospholipid bilayers as an alpha-helical monomer. We used a simplex method to optimize the 'a, d hydrophobicity' and 'e, g charge' of a series of five peptides, where 'abcdefg' correspond to the positions in two turns of an alpha-helix. Circular dichroism and analytical ultra-centrifugation experiments showed that the final peptide (helix5) is monomeric and has an alpha-helix content of approximately 89% at 0 degrees C in aqueous solution. In the presence of large unilamellar vesicles (LUVs), helix5 partitions between the aqueous and membranous phases with a partition constant well suited for measurements by electron paramagnetic resonance (EPR) spectroscopy. EPR power saturation experiments with a cysteine-scanning strategy showed that the alpha-helicity of helix5 is conserved upon binding to LUVs and that the alpha-helix binds parallel to the membrane surface with the central axis approximately 5 A below the lipid phosphate groups. Helix5 should be a useful model peptide for studies aimed at dissecting the determinants of the membrane binding of alpha-helices. The simplex-based strategy may be useful in the rational design of proteins when desired structural or partitioning properties cannot be selected or screened from libraries.

Photolysis intermediates of human rhodopsin.

Photochemical studies were conducted on human rhodopsin at 20 degrees C to characterize the intermediates which precede the formation of metarhodopsin II, the trigger for the enzyme cascade mechanism of visual transduction. Human rhodopsin was prepared from eyes which had previously been used for corneal donations. Time resolved absorption spectra collected from 10(-8) to 10(-6) s after photolysis of human rhodopsin in detergent suspensions displayed biexponential decay kinetics. The apparent lifetimes obtained from the data are 65 +/- 20 and 292 +/- 25 ns, almost a factor of 2 slower than the corresponding rates in bovine rhodopsin. The spectra can be fit well using a model in which human bathorhodopsin decays toward equilibrium with a blue-shifted intermediate (BSI) which then decays to lumirhodopsin. Spectra and kinetic rate constants were determined for all these intermediates using a global analysis which showed that the spectra of the human intermediates are remarkably similar to bovine intermediates. Microscopic rate constants derived from this model are 7.4 x 10(6) s-1 for bathorhodopsin decay and 7.5 x 10(6) s-1 and 4.6 x 10(6) s-1 for the forward and reverse reactions of BSI, respectively. Decay of lumirhodopsin to later intermediates was studied from 10(-6) to 10(-1) s after photolysis of rhodopsin in human disk membrane suspensions. The human metarhodopsin I in equilibrium metarhodopsin II equilibrium appears to be more forward shifted than in comparable bovine studies.

Bacteriorhodopsin D85N: three spectroscopic species in equilibrium.

Ground-state absorbance measurements show that BR from Halobacterium halobium containing asparagine at residue 85 (D85N) exists as three distinct chromophoric states in equilibrium. In the pH range 6-12 the absorbance spectra of the three states are demonstrated to be similar to flash-induced spectral intermediates which comprise the latter portion of the wild-type BR photocycle. One of the states absorbs maximally at 405 nm, has a deprotonated Schiff base, and contains predominantly the 13-cis form of retinal, identifying it as a close homologue of the M intermediate in the BR photocycle. The other species possess absorbance maxima with correspondence to those of the wild-type N (570 nm) and O (615 nm) photointermediates. The retinal composition of the O-like form was found to be dominated by all-trans isomer. The pH dependence of the concentrations of the equilibrium species corresponds closely with the pH dependence of the M, N, and O photointermediates. These data support kinetic models which emphasize the role of back-reactions during the photocycle of bacteriorhodopsin. Energetic and spectral characterization of the D85N ground-state equilibrium supports its use as a model for elucidating molecular transitions comprising the latter portion of the BR photocycle.

Effects of detergent environments on the photocycle of purified monomeric bacteriorhodopsin.

Time-resolved difference spectra have been obtained for the photocycle of delipidated bacteriorhodopsin monomers (d-BR) in six different detergent micelle environments that were prepared by two new detergent-exchange techniques. A global kinetic analysis of the photocycle spectra for d-BR in each detergent environment was performed. Comparison of these results with those obtained for the photocycle of bacteriorhodopsin in purple membrane (PM) shows that there is one fewer kinetically distinguishable process for monomeric BR between the decay of the K intermediate and the rise of the M intermediate. Assuming a sequential pathway occurs in the photocycle, it appears that the equilibrium between the L and M intermediates is reached much more rapidly in the detergent micelles. This is attributed to a more direct interaction between Asp-85 and the proton on the nitrogen of the Schiff base of retinal for BR in the detergents. Equilibrium concentrations of late photocycle intermediates are also altered in detergents. The later steps of the photocycle, including the decay of the M intermediate, are slowed in detergents with rings in their hydrocarbon region. This is attributed to effects on conformational changes occurring during the decay of M and/or other later photocycle intermediates. The lifetime of dark adaptation of light-adapted d-BR in different detergent environments increases in environments where the lifetime of the M intermediate increases. These results suggest that the high percentage of either unsaturated or methyl-branched lipids in PM and the membranes of other retinal proteins may be important for their effective functioning.

Direct measurement of small ligand-induced conformational changes in the aspartate chemoreceptor using EPR.

Ligand-binding-induced conformational changes in the Salmonella typhimurium aspartate receptor were studied using spin-labeling electron paramagnetic resonance. Cysteine residues, introduced by site-directed mutagenesis at several positions in the aspartate receptor periplasmic domain, were used to attach covalently a thiol-specific spin label. The electron paramagnetic resonance spectra of these labeled proteins were obtained in the presence and absence of the ligand aspartate, and used to calculate the distance change between spin labels. The results support a model in which transmembrane signaling is executed by a combined movement of alpha helix 4 (which leads into transmembrane domain 2) relative to alpha helix 1 (connected to transmembrane domain 1), as well as a coming together of the two subunits. Ligand binding causes spin labels at position 39 and 179 (within one subunit) to move further from each other and spin labels at position 39 and 39' (between two subunits) to move closer to each other. Both of these changes are very small-less than 2.5 A. No similar changes were detected in any aspartate receptor samples solubilized in detergent, suggesting that the membrane is required for these conformational changes. This is the first case of physically measured ligand-induced changes in a full-length 1-2 transmembrane domain receptor, and the results suggest that very small ligand-induced movements can result in large effects on the activity of downstream proteins.