An Implicit Membrane Generalized Born Theory for the Study of Structure, Stability, and Interactions of Membrane Proteins

This is the published version. Copyright 2003 by Elsevier.Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available

Implicit solvation based on generalized Born theory in different dielectric environments

Copyright 2004 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1631258.In this paper we are investigating the effect of the dielectric environment on atomic Born radii used in generalized Born (GB) methods. Motivated by the Kirkwood expression for the reaction field of a single off-center charge in a spherical cavity, we are proposing extended formalisms for the calculation of Born radii as a function of external and internal dielectric constants. We demonstrate that reaction field energies calculated from environmentally dependent Born radii lead to much improved agreement with Poisson–Boltzmann solutions for low dielectric external environments, such as biological membranes or organic solvent, compared to previous methods where the calculation of Born radii does not depend on the environment. We also examine how this new approach can be applied for the calculation of transfer free energies from vacuum to a given external dielectric for a system with an internal dielectric larger than one. This has not been possible with standard GB theory but is relevant when scoring minimized or average structures with implicit solvent

Membrane Assembly of Simple Helix Homo-Oligomers Studied via Molecular Dynamics Simulations

This is the publisher's version. Copyright 2007 by Elsevier.The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins—glycophorin A, the M2 proton channel, and phospholamban—using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state

Deprotonation by Dehydration: The Origin of Ammonium Sensing in the AmtB Channel

The AmtB channel passively allows the transport of NH(4) (+) across the membranes of bacteria via a “gas” NH(3) intermediate and is related by homology (sequentially, structurally, and functionally) to many forms of Rh protein (both erythroid and nonerythroid) found in animals and humans. New structural information on this channel has inspired computational studies aimed at clarifying various aspects of NH(4) (+) recruitment and binding in the periplasm, as well as its deprotonation. However, precise mechanisms for these events are still unknown, and, so far, explanations for subsequent NH(3) translocation and reprotonation at the cytoplasmic end of the channel have not been rigorously addressed. We employ molecular dynamics simulations and free energy methods on a full AmtB trimer system in membrane and bathed in electrolyte. Combining the potential of mean force for NH(4) (+)/NH(3) translocation with data from thermodynamic integration calculations allows us to find the apparent pK(a) of NH(4) (+) as a function of the transport axis. Our calculations reveal the specific sites at which its deprotonation (at the periplasmic end) and reprotonation (at the cytoplasmic end) occurs. Contrary to most hypotheses, which ascribe a proton-accepting role to various periplasmic or luminal residues of the channel, our results suggest that the most plausible proton donor/acceptor at either of these sites is water. Free-energetic analysis not only verifies crystallographically determined binding sites for NH(4) (+) and NH(3) along the transport axis, but also reveals a previously undetermined binding site for NH(4) (+) at the cytoplasmic end of the channel. Analysis of dynamics and the free energies of all possible loading states for NH(3) inside the channel also reveal that hydrophobic pressure and the free-energetic profile provided by the pore lumen drives this species toward the cytoplasm for protonation just before reaching the newly discovered site

Predicting extreme p K a shifts in staphylococcal nuclease mutants with constant pH molecular dynamics

Accurate computational methods of determining protein and nucleic acid p K a values are vital to understanding pH‐dependent processes in biological systems. In this article, we use the recently developed method constant pH molecular dynamics (CPHMD) to explore the calculation of highly perturbed p K a values in variants of staphylococcal nuclease (SNase). Simulations were performed using the replica exchange (REX) protocol for improved conformational sampling with eight temperature windows, and yielded converged proton populations in a total sampling time of 4 ns. Our REX‐CPHMD simulations resulted in calculated p K a values with an average unsigned error (AUE) of 0.75 pK units for the acidic residues in Δ + PHS, a hyperstable variant of SNase. For highly p K a ‐perturbed SNase mutants with known crystal structures, our calculations yielded an AUE of 1.5 pK units and for those mutants based on modeled structures an AUE of 1.4 pK units was found. Although a systematic underestimate of pK shifts was observed in most of the cases for the highly perturbed pK mutants, correlations between conformational rearrangement and plasticity associated with the mutation and error in p K a prediction was not evident in the data. This study further extends the scope of electrostatic environments explored using the REX‐CPHMD methodology and suggests that it is a reliable tool for rapidly characterizing ionizable amino acids within proteins even when modeled structures are employed. Proteins 2011; © 2011 Wiley‐Liss, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/88038/1/23195_ftp.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/88038/2/PROT_23195_sm_SuppInfo.pd

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields

Single transcriptional and translational preQ1 riboswitches adopt similar pre-folded ensembles that follow distinct folding pathways into the same ligand-bound structure

Riboswitches are structural elements in the 50 untranslated regions of many bacterial messenger RNAs that regulate gene expression in response to changing metabolite concentrations by inhibition of either transcription or translation initiation. The preQ1 (7-aminomethyl-7-deazaguanine) riboswitch family comprises some of the smallest metabolite sensing RNAs found in nature. Once ligand-bound, the transcriptional Bacillus subtilis and translational Thermoanaerobacter tengcongensis preQ1 riboswitch aptamers are structurally similar RNA pseudoknots; yet, prior structural studies have characterized their ligand-free conformations as largely unfolded and folded, respectively. In contrast, through single molecule observation, we now show that, at nearphysiological Mg2+ concentration and pH, both ligand-free aptamers adopt similar pre-folded state ensembles that differ in their ligand-mediated folding. Structure-based Go¯ -model simulations of the two aptamers suggest that the ligand binds late (Bacillus subtilis) and early (Thermoanaerobacter tengcongensis) relative to pseudoknot folding, leading to the proposal that the principal distinction between the two riboswitches lies in their relative tendencies to fold via mechanisms of conformational selection and induced fit, respectively. These mechanistic insights are put to the test by rationally designing a single nucleotide swap distal from the ligand binding pocket that we find to predictably control the aptamers0 pre-folded states and their ligand binding affinities

Assessing the carcinogenic potential of low-dose exposures to chemical mixtures in the environment: the challenge ahead.

Lifestyle factors are responsible for a considerable portion of cancer incidence worldwide, but credible estimates from the World Health Organization and the International Agency for Research on Cancer (IARC) suggest that the fraction of cancers attributable to toxic environmental exposures is between 7% and 19%. To explore the hypothesis that low-dose exposures to mixtures of chemicals in the environment may be combining to contribute to environmental carcinogenesis, we reviewed 11 hallmark phenotypes of cancer, multiple priority target sites for disruption in each area and prototypical chemical disruptors for all targets, this included dose-response characterizations, evidence of low-dose effects and cross-hallmark effects for all targets and chemicals. In total, 85 examples of chemicals were reviewed for actions on key pathways/mechanisms related to carcinogenesis. Only 15% (13/85) were found to have evidence of a dose-response threshold, whereas 59% (50/85) exerted low-dose effects. No dose-response information was found for the remaining 26% (22/85). Our analysis suggests that the cumulative effects of individual (non-carcinogenic) chemicals acting on different pathways, and a variety of related systems, organs, tissues and cells could plausibly conspire to produce carcinogenic synergies. Additional basic research on carcinogenesis and research focused on low-dose effects of chemical mixtures needs to be rigorously pursued before the merits of this hypothesis can be further advanced. However, the structure of the World Health Organization International Programme on Chemical Safety 'Mode of Action' framework should be revisited as it has inherent weaknesses that are not fully aligned with our current understanding of cancer biology

TMEM106B is a genetic modifier of frontotemporal lobar degeneration with C9orf72 hexanucleotide repeat expansions

Hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72) have recently been linked to frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis, and may be the most common genetic cause of both neurodegenerative diseases. Genetic variants at TMEM106B influence risk for the most common neuropathological subtype of FTLD, characterized by inclusions of TAR DNA-binding protein of 43 kDa (FTLD-TDP). Previous reports have shown that TMEM106B is a genetic modifier of FTLD-TDP caused by progranulin (GRN) mutations, with the major (risk) allele of rs1990622 associating with earlier age at onset of disease. Here, we report that rs1990622 genotype affects age at death in a single-site discovery cohort of FTLD patients with C9orf72 expansions (n = 14), with the major allele correlated with later age at death (p = 0.024). We replicate this modifier effect in a 30-site international neuropathological cohort of FTLD-TDP patients with C9orf72 expansions (n = 75), again finding that the major allele associates with later age at death (p = 0.016), as well as later age at onset (p = 0.019). In contrast, TMEM106B genotype does not affect age at onset or death in 241 FTLD-TDP cases negative for GRN mutations or C9orf72 expansions. Thus, TMEM106B is a genetic modifier of FTLD with C9orf72 expansions. Intriguingly, the genotype that confers increased risk for developing FTLD-TDP (major, or T, allele of rs1990622) is associated with later age at onset and death in C9orf72 expansion carriers, providing an example of sign epistasis in human neurodegenerative disease

Physiological Correlates of Volunteering

We review research on physiological correlates of volunteering, a neglected but promising research field. Some of these correlates seem to be causal factors influencing volunteering. Volunteers tend to have better physical health, both self-reported and expert-assessed, better mental health, and perform better on cognitive tasks. Research thus far has rarely examined neurological, neurochemical, hormonal, and genetic correlates of volunteering to any significant extent, especially controlling for other factors as potential confounds. Evolutionary theory and behavioral genetic research suggest the importance of such physiological factors in humans. Basically, many aspects of social relationships and social activities have effects on health (e.g., Newman and Roberts 2013; Uchino 2004), as the widely used biopsychosocial (BPS) model suggests (Institute of Medicine 2001). Studies of formal volunteering (FV), charitable giving, and altruistic behavior suggest that physiological characteristics are related to volunteering, including specific genes (such as oxytocin receptor [OXTR] genes, Arginine vasopressin receptor [AVPR] genes, dopamine D4 receptor [DRD4] genes, and 5-HTTLPR). We recommend that future research on physiological factors be extended to non-Western populations, focusing specifically on volunteering, and differentiating between different forms and types of volunteering and civic participation