Assembly and Analysis of Cell-Scale Membrane Envelopes.

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

A membrane-embedded pathway delivers general anesthetics to two interacting binding sites in the Gloeobacter violaceus ion channel.

General anesthetics exert their effects on the central nervous system by acting on ion channels, most notably pentameric ligand-gated ion channels. Although numerous studies have focused on pentameric ligand-gated ion channels, the details of anesthetic binding and channel modulation are still debated. A better understanding of the anesthetic mechanism of action is necessary for the development of safer and more efficacious drugs. Herein, we present a computational study identifying two anesthetic binding sites in the transmembrane domain of the Gloeobacter violaceus ligand-gated ion channel (GLIC) channel, characterize the putative binding pathway, and observe structural changes associated with channel function. Molecular simulations of desflurane reveal a binding pathway to GLIC via a membrane-embedded tunnel using an intrasubunit protein lumen as the conduit, an observation that explains the Meyer-Overton hypothesis, or why the lipophilicity of an anesthetic and its potency are generally proportional. Moreover, employing high concentrations of ligand led to the identification of a second transmembrane site (TM2) that inhibits dissociation of anesthetic from the TM1 site and is consistent with the high concentrations of anesthetics required to achieve clinical effects. Finally, asymmetric binding patterns of anesthetic to the channel were found to promote an iris-like conformational change that constricts and dehydrates the ion pore, creating a 13.5 kcal/mol barrier to ion translocation. Together with previous studies, the simulations presented herein demonstrate a novel anesthetic binding site in GLIC that is accessed through a membrane-embedded tunnel and interacts with a previously known site, resulting in conformational changes that produce a non-conductive state of the channel.

A "cross-stitched" peptide with improved helicity and proteolytic stability.

A new computational approach to obtain quantitative energy profiles for helix folding was used in the design of orthogonal hydrocarbon and lactam bicyclic peptides. The proteolytically stable, "cross-stitched" peptide SRC2-BCP1 shows nanomolar affinity for estrogen receptor α and X-ray crystallography confirms a helical binding pose.

The cellular membrane as a mediator for small molecule interaction with membrane proteins.

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Flexibility Coexists with Shape-Persistence in Cyanostar Macrocycles.

Shape-persistent macrocycles are attractive functional targets for synthesis, molecular recognition, and hierarchical self-assembly. Such macrocycles are noncollapsible and geometrically well-defined, and they are traditionally characterized by having repeat units and low conformational flexibility. Here, we find it necessary to refine these ideas in the face of highly flexible yet shape-persistent macrocycles. A molecule is shape-persistent if it has a small change in shape when perturbed by external stimuli (e.g., heat, light, and redox chemistry). In support of this idea, we provide the first examination of the relationships between a macrocycle's shape persistence, its conformational space, and the resulting functions. We do this with a star-shaped macrocycle called cyanostar that is flexible as well as being shape-persistent. We employed molecular dynamics (MD), density functional theory (DFT), and NMR experiments. Considering a thermal bath as a stimulus, we found a single macrocycle has 332 accessible conformers with olefins undergoing rapid interconversion by up-down and in-out motions on short time scales (0.2 ns). These many interconverting conformations classify single cyanostars as flexible. To determine and confirm that cyanostars are shape-persistent, we show that they have a high 87% shape similarity across these conformations. To further test the idea, we use the binding of diglyme to the single macrocycle as guest-induced stimulation. This guest has almost no effect on the conformational space. However, formation of a 2:1 sandwich complex involving two macrocycles enhances rigidity and dramatically shifts the conformer distribution toward perfect bowls. Overall, the present study expands the scope of shape-persistent macrocycles to include flexible macrocycles if, and only if, their conformers have similar shapes.

Stapled Peptides with γ-Methylated Hydrocarbon Chains for the Estrogen Receptor/Coactivator Interaction.

"Stapled" peptides are typically designed to replace two non-interacting residues with a constraining, olefinic staple. To mimic interacting leucine and isoleucine residues, we have created new amino acids that incorporate a methyl group in the γ-position of the stapling amino acid S5. We have incorporated them into a sequence derived from steroid receptor coactivator 2, which interacts with estrogen receptor α. The best peptide (IC50 =89 nm) replaces isoleucine 689 with an S-γ-methyl stapled amino acid, and has significantly higher affinity than unsubstituted peptides (390 and 760 nm). Through X-ray crystallography and molecular dynamics studies, we show that the conformation taken up by the S-γ-methyl peptide minimizes the syn-pentane interactions between the α- and γ-methyl groups.

Macromolecular Crystallography for Synthetic Abiological Molecules: Combining xMDFF and PHENIX for Structure Determination of Cyanostar Macrocycles.

Crystal structure determination has long provided insight into structure and bonding of small molecules. When those same small molecules are designed to come together in multimolecular assemblies, such as in coordination cages, supramolecular architectures and organic-based frameworks, their crystallographic characteristics closely resemble biological macromolecules. This resemblance suggests that biomacromolecular refinement approaches be used for structure determination of abiological molecular complexes that arise in an aggregate state. Following this suggestion we investigated the crystal structure of a pentagonal macrocycle, cyanostar, by means of biological structure analysis methods and compared results to traditional small molecule methods. Cyanostar presents difficulties seen in supramolecular crystallography including whole molecule disorder and highly flexible solvent molecules sitting in macrocyclic and intermolecule void spaces. We used the force-field assisted refinement method, molecular dynamics flexible fitting algorithm for X-ray crystallography (xMDFF), along with tools from the macromolecular structure determination suite PHENIX. We found that a standard implementation of PHENIX, namely one without xMDFF, either fails to produce a solution by molecular replacement alone or produces an inaccurate structure when using generic geometry restraints, even at a very high diffraction data resolution of 0.84 Å. The problems disappear when taking advantage of xMDFF, which applies an optimized force field to realign molecular models during phasing by providing accurate restraints. The structure determination for this model system shows excellent agreement with the small-molecule methods. Therefore, the joint xMDFF-PHENIX refinement protocol provides a new strategy that uses macromolecule methods for structure determination of small molecules and their assemblies.

Ligand Accessibility and Bioactivity of a Hormone-Dendrimer Conjugate Depend on pH and pH History.

Estrogen conjugates with a polyamidoamine (PAMAM) dendrimer have shown remarkably selective regulation of the nongenomic actions of estrogens in target cells. In response to pH changes, however, these estrogen-dendrimer conjugates (EDCs) display a major morphological transition that alters the accessibility of the estrogen ligands that compromises the bioactivity of the EDC. A sharp break in dynamic behavior near pH 7 occurs for three different ligands on the surface of a PAMAM-G6 dendrimer: a fluorophore (tetramethylrhodamine [TMR]) and two estrogens (17α-ethynylestradiol and diphenolic acid). Collisional quenching and time-resolved fluorescence anisotropy experiments with TMR-PAMAM revealed high ligand shielding above pH 7 and low shielding below pH 7. Furthermore, when the pH was cycled from 8.5 (conditions of ligand-PAMAM conjugation) to 4.5 (e.g., endosome/lysosome) and through 6.5 (e.g., hypoxic environment) back to pH 8.5, the 17α-ethynylestradiol- and diphenolic acid-PAMAM conjugates experienced a dramatic, irreversible loss in cell stimulatory activity; dynamic NMR studies indicated that the hormonal ligands had become occluded within the more hydrophobic core of the PAMAM dendrimer. Thus, the active state of these estrogen-dendrimer conjugates appears to be metastable. This pH-dependent irreversible masking of activity is of considerable relevance to the design of drug conjugates with amine-bearing PAMAM dendrimers.

Neural-Network Scoring Functions Identify Structurally Novel Estrogen-Receptor Ligands.

The magnitude of the investment required to bring a drug to the market hinders medical progress, requiring hundreds of millions of dollars and years of research and development. Any innovation that improves the efficiency of the drug-discovery process has the potential to accelerate the delivery of new treatments to countless patients in need. "Virtual screening," wherein molecules are first tested in silico in order to prioritize compounds for subsequent experimental testing, is one such innovation. Although the traditional scoring functions used in virtual screens have proven useful, improved accuracy requires novel approaches. In the current work, we use the estrogen receptor to demonstrate that neural networks are adept at identifying structurally novel small molecules that bind to a selected drug target, ultimately allowing experimentalists to test fewer compounds in the earliest stages of lead identification while obtaining higher hit rates. We describe 39 novel estrogen-receptor ligands identified in silico with experimentally determined Ki values ranging from 460 nM to 20 µM, presented here for the first time.

Rapid parameterization of small molecules using the Force Field Toolkit.

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics simulations to many biological systems of interest, especially in fields such as drug discovery. Although the release of generalized versions of common classical force fields, for example, General Amber Force Field and CHARMM General Force Field, have posited guidelines for parameterization of small molecules, many technical challenges remain that have hampered their wide-scale extension. The Force Field Toolkit (ffTK), described herein, minimizes common barriers to ligand parameterization through algorithm and method development, automation of tedious and error-prone tasks, and graphical user interface design. Distributed as a VMD plugin, ffTK facilitates the traversal of a clear and organized workflow resulting in a complete set of CHARMM-compatible parameters. A variety of tools are provided to generate quantum mechanical target data, setup multidimensional optimization routines, and analyze parameter performance. Parameters developed for a small test set of molecules using ffTK were comparable to existing CGenFF parameters in their ability to reproduce experimentally measured values for pure-solvent properties (<15% error from experiment) and free energy of solvation (±0.5 kcal/mol from experiment).

1,25-Dihydroxyvitamin D3 acts directly on the T lymphocyte vitamin D receptor to inhibit experimental autoimmune encephalomyelitis.

Multiple sclerosis (MS) is an incurable autoimmune neurodegenerative disease. Environmental factors may be key to MS prevention and treatment. MS prevalence and severity decrease with increasing sunlight exposure and vitamin D(3) supplies, supporting our hypothesis that the sunlight-dependent hormone, 1,25-dihydroxyvitamin D(3) (1,25-(OH)(2) D(3) ), inhibits autoimmune T-cell responses in MS. Moreover, 1,25-(OH)(2) D(3) inhibits and reverses experimental autoimmune encephalomyelitis (EAE), an MS model. Here, we investigated whether 1,25-(OH)(2) D(3) inhibits EAE via the vitamin D receptor (VDR) in T lymphocytes. Using bone marrow chimeric mice with a disrupted VDR only in radio-sensitive hematopoietic cells or radio-resistant non-hematopoietic cells, we found that hematopoietic cell VDR function was necessary for 1,25-(OH)(2) D(3) to inhibit EAE. Furthermore, conditional targeting experiments showed that VDR function in T cells was necessary. Neither 1,25-(OH)(2) D(3) nor T-cell-specific VDR targeting influenced CD4(+) Foxp3(+) T-cell proportions in the periphery or the CNS in these studies. These data support a model wherein 1,25-(OH)(2) D(3) acts directly on pathogenic CD4(+) T cells to inhibit EAE.

Defining the Energetic Basis for a Conformational Switch Mediating Ligand-Independent Activation of Mutant Estrogen Receptors in Breast Cancer.

Although most primary estrogen receptor (ER)-positive breast cancers respond well to endocrine therapies, many relapse later as metastatic disease due to endocrine therapy resistance. Over one third of these are associated with mutations in the ligand-binding domain (LBD) that activate the receptor independent of ligand. We have used an array of advanced computational techniques rooted in molecular dynamics simulations, in concert with and validated by experiments, to characterize the molecular mechanisms by which specific acquired somatic point mutations give rise to ER constitutive activation. By comparing structural and energetic features of constitutively active mutants and ligand-bound forms of ER-LBD with unliganded wild-type (WT) ER, we characterize a spring force originating from strain in the Helix 11-12 loop of WT-ER, opposing folding of Helix 12 into the active conformation and keeping WT-ER off and disordered, with the ligand-binding pocket open for rapid ligand binding. We quantify ways in which this spring force is abrogated by activating mutations that latch (Y537S) or relax (D538G) the folded form of the loop, enabling formation of the active conformation without ligand binding. We also identify a new ligand-mediated hydrogen-bonding network that stabilizes the active, ligand-bound conformation of WT-ER LBD, and similarly stabilizes the active conformation of the ER mutants in the hormone-free state. IMPLICATIONS: Our investigations provide deep insight into the energetic basis for the structural mechanisms of receptor activation through mutation, exemplified here with ER in endocrine-resistant metastatic breast cancers, with potential application to other dysregulated receptor signaling due to driver mutations.

Altered BAFF-receptor signaling and additional modifier loci contribute to systemic autoimmunity in A/WySnJ mice.

Systemic lupus erythematosus pathology reflects autoantibody-mediated damage due to a failure of B-lymphocyte tolerance. We previously reported that B-lymphopenic A/WySnJ mice develop a lupus-like syndrome and linked this syndrome to the B-cell maturation defect-1 (Bcmd-1) mutant allele of the B-cell-activating factor belonging to the TNF family-receptor (Baffr) gene. Here, we further evaluate the genetic basis for autoimmunity in A/WySnJ mice. We produced B6.Bcmd-1 and AW.Baffr(-/-) congenic mice (N5), and compared them with B6.Baffr(-/-) and A/WySnJ mice with respect to B-lymphocyte development. Bcmd-1-expressing mice had more B cells with greater maturity than Baffr(-/-) mice regardless of genetic background, indicating that Bcmd-1 encodes a partially functional BAFF-R. We also compared these mice for lupus phenotypes to determine whether Bcmd-1 is necessary and sufficient for disease, or whether the Baffr(-/-) (-) allele can also cause autoimmunity. The Baffr(-/-) allele did not lead to autoimmunity on either genetic background. In contrast, the Bcmd-1 allele was necessary and sufficient for development of low levels of IgM autoantibodies in B6.Bcmd-1 mice. However, Bcmd-1 plus unidentified A/WySnJ modifier genes were necessary for development of IgG autoantibodies and renal pathology. We propose that in A/WySnJ mice an excess of BAFF per B cell rescues self-reactive B cells through a partially functional BAFF-R in a B-lymphopenic environment.

A mutant form of ERα associated with estrogen insensitivity affects the coupling between ligand binding and coactivator recruitment.

A homozygous missense mutation in the gene encoding the estrogen receptor α (ERα) was previously identified in a female patient with estrogen insensitivity syndrome. We investigated the molecular features underlying the impaired transcriptional response of this mutant (ERα-Q375H) and four other missense mutations at this position designed to query alternative mechanisms. The identity of residue 375 greatly affected the sensitivity of the receptor to agonists without changing the ligand binding affinity. Instead, the mutations caused changes in the affinity of coactivator binding and alterations in the balance of coactivator and corepressor recruitment. Comparisons among the transcriptional regulatory responses of these six ERα genotypes to a set of ER agonists showed that both steric and electrostatic factors contributed to the functional deficits in gene regulatory activity of the mutant ERα proteins. ERα-coregulator peptide binding in vitro and RIME (rapid immunoprecipitation mass spectrometry of endogenous) analysis in cells showed that the degree of functional impairment paralleled changes in receptor-coregulator binding interactions. These findings uncover coupling between ligand binding and coregulator recruitment that affects the potency rather than the efficacy of the receptor response without substantially altering ligand binding affinity. This highlights a molecular mechanism for estrogen insensitivity syndrome involving mutations that perturb a bidirectional allosteric coupling between ligand binding and coregulator binding that determines receptor transcriptional output.

The Role of Micronutrients in Graft-VS.-Host Disease: Immunomodulatory Effects of Vitamins A and D.

Graft-vs.-host disease (GVHD) remains a major obstacle to the success of allogeneic hematopoietic stem cell transplantation (HSCT). GVHD occurs because donor T cells in the allograft recognize the genetically disparate host as foreign and attack the transplant recipient's tissues. While genetic incompatibility between donor and recipient is the primary determinant for the extent of alloimmune response, GVHD incidence and severity are also influenced by non-genetic factors. Recent advances in immunology establish that environmental factors, including dietary micronutrients, contribute significantly to modulating various immune responses and may influence the susceptibility to autoimmune and inflammatory diseases of experimental animals and humans. Emerging clinical and preclinical evidence indicates that certain micronutrients may participate in regulating GVHD risk after allogeneic HSCT. In this review, we summarize recent advances in our understanding with respect to the potential role of micronutrients in the pathogenesis of acute and chronic GVHD, focusing on vitamins A and D.

Publisher Correction: Structural underpinnings of oestrogen receptor mutations in endocrine therapy resistance.

The originally published article contained an error in the legend of supplementary figure 1. A figure permission line was left off. The correct figure permission line has now been added to the HTML and PDF versions of the article, stating that "Data shown in (B) and (C) of this figure were originally published in Jeyakumar, M., Carlson, K. E., Gunther, J. R. & Katzenellenbogen, J. A. Exploration of dimensions of estrogen potency: parsing ligand binding and coactivator binding affinities. J. Biol. Chem. 286, 12971-12982, (2011) (c) the American Society for Biochemistry and Molecular Biology (Ref. 53)."

Structural underpinnings of oestrogen receptor mutations in endocrine therapy resistance.

Oestrogen receptor-α (ERα), a key driver of breast cancer, normally requires oestrogen for activation. Mutations that constitutively activate ERα without the need for hormone binding are frequently found in endocrine-therapy-resistant breast cancer metastases and are associated with poor patient outcomes. The location of these mutations in the ER ligand-binding domain and their impact on receptor conformation suggest that they subvert distinct mechanisms that normally maintain the low basal state of wild-type ERα in the absence of hormone. Such mutations provide opportunities to probe fundamental issues underlying ligand-mediated control of ERα activity. Instructive contrasts between these ERα mutations and those that arise in the androgen receptor (AR) during anti-androgen treatment of prostate cancer highlight differences in how activation functions in ERs and AR control receptor activity, how hormonal pressures (deprivation versus antagonism) drive the selection of phenotypically different mutants, how altered protein conformations can reduce antagonist potency and how altered ligand-receptor contacts can invert the response that a receptor has to an agonist ligand versus an antagonist ligand. A deeper understanding of how ligand regulation of receptor conformation is linked to receptor function offers a conceptual framework for developing new anti-oestrogens that might be more effective in preventing and treating breast cancer.

A model of TH17-associated ileal hyperplasia that requires both IL-17A and IFNγ to generate self-tolerance and prevent colitis.

Homeostasis in the ileum, which is commonly disrupted in patients with Crohn's disease, involves ongoing immune responses. To study how homeostatic processes of the ileum impact CD4+T cell responses, we used TCR transgenic tools to breed mice that spontaneously produced CD4+T cells reactive to an antigen expressed in the ileum. At an early age, the ilea of these mice exhibit crypt hyperplasia and accumulate increased numbers of TH17 cells bearing non-transgenic clonotypes. Half of these mice subsequently developed colitis linked to broad mucosal infiltration by TH17 and TH1 cells expressing non-transgenic clonotypes, chronic wasting disease and loss of ileal crypt hyperplasia. By contrast, adult mice with normal growth continued to exhibit TH17-associated ileal crypt hyperplasia and additionally accumulated ileal-reactive Treg cells. Both IL-17A and IFNγ were protective, as their deficiency precluded ileal-reactive Treg accumulation and exacerbated colitic disease. IL-23R blockade prevented progression to colitis, whereas nTreg cell transfers prevented colitic disease, ileal crypt hyperplasia and ileal-reactive Treg accumulation. Thus, our studies identify an IL-17A and IFNγ-dependent homeostatic process that mobilizes ileal-reactive Treg cells and is disrupted by IL-23.

Antagonists for Constitutively Active Mutant Estrogen Receptors: Insights into the Roles of Antiestrogen-Core and Side-Chain.

A major risk for patients having estrogen receptor α (ERα)-positive breast cancer is the recurrence of drug-resistant metastases after initial successful treatment with endocrine therapies. Recent studies have implicated a number of activating mutations in the ligand-binding domain of ERα that stabilize the agonist conformation as a prominent mechanism for this acquired resistance. There are several critical gaps in our knowledge regarding the specific pharmacophore requirements of an antagonist that could effectively inhibit all or most of the different mutant ERs. To address this, we screened various chemotypes for blocking mutant ER-mediated transcriptional signaling and identified RU58668 as a model compound that contains structural elements that support potent ligand-induced inhibition of mutant ERs. We designed and synthesized a focused library of novel antagonists and probed how small and large perturbations in different ligand structural regions influenced inhibitory activity on individual mutant ERs in breast cancer cells. Effective inhibition derives from both nonpolar and moderately polar motifs in a multifunctional side chain of the antagonists, with the nature of the ligand core making important contributions by increasing the potency of ligands possessing similar types of side chains. Some of our new antagonists potently blocked the transcriptional activity of the three most common mutant ERs (L536R, Y537S, D538G) and inhibited mutant ER-mediated cell proliferation. Supported by our molecular modeling, these studies provide new insights into the role of specific components, involving both the ligand core and multifunctional side chain, in suppressing wild-type and mutant ER-mediated transcription and breast cancer cell proliferation.

The SERM/SERD bazedoxifene disrupts ESR1 helix 12 to overcome acquired hormone resistance in breast cancer cells.

Acquired resistance to endocrine therapy remains a significant clinical burden for breast cancer patients. Somatic mutations in the ESR1 (estrogen receptor alpha (ERα)) gene ligand-binding domain (LBD) represent a recognized mechanism of acquired resistance. Antiestrogens with improved efficacy versus tamoxifen might overcome the resistant phenotype in ER +breast cancers. Bazedoxifene (BZA) is a potent antiestrogen that is clinically approved for use in hormone replacement therapies. We found that BZA possesses improved inhibitory potency against the Y537S and D538G ERα mutants compared to tamoxifen and has additional inhibitory activity in combination with the CDK4/6 inhibitor palbociclib. In addition, comprehensive biophysical and structural biology studies show BZA's selective estrogen receptor degrading (SERD) properties that override the stabilizing effects of the Y537S and D538G ERα mutations.