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.

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.

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.

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.

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)."

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.

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.

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.

Adamantyl Antiestrogens with Novel Side Chains Reveal a Spectrum of Activities in Suppressing Estrogen Receptor Mediated Activities in Breast Cancer Cells.

To search for new antiestrogens more effective in treating breast cancers, we explored alternatives to the acrylic acid side chain used in many antiestrogens. To facilitate our search, we used a simple adamantyl ligand core that by avoiding stereochemical issues enabled rapid synthesis of acrylate ketone, ester, and amide analogs. All compounds were high affinity estrogen receptor α (ERα) ligands but displayed a range of efficacies and potencies as antiproliferative and ERα-downregulating agents. There were large differences in activity between compounds having minor structural changes, but antiproliferative and ERα-downregulating efficacies generally paralleled one another. Some compounds with side chain polar groups had particularly high affinities. The secondary carboxamides had the best cellular activities, and the 3-hydroxypropylamide was as efficacious as fulvestrant in suppressing cell proliferation and gene expression. This study has produced structurally novel antiestrogens based on a simple adamantyl core structure with acrylate side chains optimized for cellular antagonist activity.

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.

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.

Adaptive Multilevel Splitting Method for Molecular Dynamics Calculation of Benzamidine-Trypsin Dissociation Time.

Adaptive multilevel splitting (AMS) is a rare event sampling method that requires minimal parameter tuning and allows unbiased sampling of transition pathways of a given rare event. Previous simulation studies have verified the efficiency and accuracy of AMS in the calculation of transition times for simple systems in both Monte Carlo and molecular dynamics (MD) simulations. Now, AMS is applied for the first time to an MD simulation of protein-ligand dissociation, representing a leap in complexity from the previous test cases. Of interest is the dissociation rate, which is typically too low to be accessible to conventional MD. The present study joins other recent efforts to develop advanced sampling techniques in MD to calculate dissociation rates, which are gaining importance in the pharmaceutical field as indicators of drug efficacy. The system investigated here, benzamidine bound to trypsin, is an example common to many of these efforts. The AMS estimate of the dissociation rate was found to be (2.6 ± 2.4) × 10(2) s(-1), which compares well with the experimental value.

Design of pathway preferential estrogens that provide beneficial metabolic and vascular effects without stimulating reproductive tissues.

There is great medical need for estrogens with favorable pharmacological profiles that support desirable activities for menopausal women, such as metabolic and vascular protection, but that lack stimulatory activities on the breast and uterus. We report the development of structurally novel estrogens that preferentially activate a subset of estrogen receptor (ER) signaling pathways and result in favorable target tissue-selective activity. Through a process of structural alteration of estrogenic ligands that was designed to preserve their essential chemical and physical features but greatly reduced their binding affinity for ERs, we obtained "pathway preferential estrogens" (PaPEs), which interacted with ERs to activate the extranuclear-initiated signaling pathway preferentially over the nuclear-initiated pathway. PaPEs elicited a pattern of gene regulation and cellular and biological processes that did not stimulate reproductive and mammary tissues or breast cancer cells. However, in ovariectomized mice, PaPEs triggered beneficial responses both in metabolic tissues (adipose tissue and liver) that reduced body weight gain and fat accumulation and in the vasculature that accelerated repair of endothelial damage. This process of designed ligand structure alteration represents a novel approach to develop ligands that shift the balance in ER-mediated extranuclear and nuclear pathways to obtain tissue-selective, non-nuclear PaPEs, which may be beneficial for postmenopausal hormone replacement. The approach may also have broad applicability for other members of the nuclear hormone receptor superfamily.

All-Atom Molecular Dynamics of Virus Capsids as Drug Targets.

Virus capsids are protein shells that package the viral genome. Although their morphology and biological functions can vary markedly, capsids often play critical roles in regulating viral infection pathways. A detailed knowledge of virus capsids, including their dynamic structure, interactions with cellular factors, and the specific roles that they play in the replication cycle, is imperative for the development of antiviral therapeutics. The following Perspective introduces an emerging area of computational biology that focuses on the dynamics of virus capsids and capsid-protein assemblies, with particular emphasis on the effects of small-molecule drug binding on capsid structure, stability, and allosteric pathways. When performed at chemical detail, molecular dynamics simulations can reveal subtle changes in virus capsids induced by drug molecules a fraction of their size. Here, the current challenges of performing all-atom capsid-drug simulations are discussed, along with an outlook on the applicability of virus capsid simulations to reveal novel drug targets.

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.

Estrogen receptor alpha somatic mutations Y537S and D538G confer breast cancer endocrine resistance by stabilizing the activating function-2 binding conformation.

Somatic mutations in the estrogen receptor alpha (ERα) gene (ESR1), especially Y537S and D538G, have been linked to acquired resistance to endocrine therapies. Cell-based studies demonstrated that these mutants confer ERα constitutive activity and antiestrogen resistance and suggest that ligand-binding domain dysfunction leads to endocrine therapy resistance. Here, we integrate biophysical and structural biology data to reveal how these mutations lead to a constitutively active and antiestrogen-resistant ERα. We show that these mutant ERs recruit coactivator in the absence of hormone while their affinities for estrogen agonist (estradiol) and antagonist (4-hydroxytamoxifen) are reduced. Further, they confer antiestrogen resistance by altering the conformational dynamics of the loop connecting Helix 11 and Helix 12 in the ligand-binding domain of ERα, which leads to a stabilized agonist state and an altered antagonist state that resists inhibition.

Preparation of o-Fluorophenols from Nonaromatic Precursors: Mechanistic Considerations for Adaptation to Fluorine-18 Radiolabeling.

The preparation of fluorine-18 labeled o-fluorophenols at high specific activity is challenging and requires use of [(18)F]fluoride ion as the radioisotope source. As a novel, alternative approach, we found that treatment of α-diazocyclohexenones with Selectfluor and Et3N·3HF followed by HF elimination and tautomerization afforded o-fluorophenols regioselectively and rapidly. To adapt this chemistry to (18)F radiolabeling, using bromine electrophiles in place of Selectfluor gave the o-fluorophenol via an α-bromo-α-fluoroketone intermediate in lower but still reasonable yields.