Transformations, Lineage Comparisons, and Analysis of Down-to-Up Protomer States of Variants of the SARS-CoV-2 Prefusion Spike Protein, Including the UK Variant B.1.1.7.

Monitoring and strategic response to variants in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represent a considerable challenge in the current pandemic and for future viral outbreaks. Mutations/deletions of the virion's prefusion Spike protein may have significant impact on vaccines and therapeutics that utilize this key structural protein in their mitigation strategies. In this study, we have demonstrated how dominant energetic landscape mappings ("glue points") based on ab inito all-atom force fields coupled with phylogenetic sequence alignment information can identify key residue mutations and deletions associated with variants. We also found several examples of excellent homology of stabilizing residue glue points across the lineages of betacoronavirus Spike proteins that we have called "sequence homologous glue points." SARS-CoV-2 demonstrates the least number of stabilizing glue points associated with interchain interactions among Down-state protomers across lineages. Additionally, we computationally studied variants among the trimeric Spike protein of SARS-CoV-2 using all-atom molecular dynamics to ascertain structural and energetic changes among variants. We examined both a theoretically based triple mutant and the UK or B.1.1.7 variant. For the theoretical triple mutant, we demonstrated through alanine substitutions that three key residues could cause the transition of Down-to-Up protomer states, where the transition is characterized by the "arm" length of the receptor-binding domain (RBD) rather than the hinge angle. For the B.1.1.7 variant, we demonstrated the critical importance of mutations D614G and N501Y on the structure and binding, respectively, of the Spike protein. We note that these same two key mutations are also found in the South African B.1.351 variant. IMPORTANCE Viral variants represent a major challenge to monitoring viral outbreaks and formulating strategic health care responses. Variants represent transmitting viruses that have specific mutations and deletions associated with their genome. In the case of SARS-CoV-2 and other related viruses (betacoronaviruses), many of these mutations and deletions are associated with the Spike protein that the virus uses to infect cells. Here, we have analyzed both SARS-CoV-2 variants and related viruses, such as Middle Eastern respiratory syndrome coronavirus (MERS-CoV), in order to understand not only differences, but also key similarities between them. Understanding similarities can be as important as differences in determining key functional features of a class of viruses, such as the betacoronaviruses. We have used both phylogenetic analysis, which traces genetic similarities and differences, along with independent biophysics analysis, which adds function or behavior, in order to determine possible functional differences and hence possible transmission and infection differences among variants and lineages.

Transformations, Lineage Comparisons, and Analysis of Down to Up Protomer States of Variants of the SARS-CoV-2 Prefusion Spike Protein Including the UK Variant B.1.1.7.

Monitoring and strategic response to variants in SARS-CoV-2 represents a considerable challenge in the current pandemic, as well as potentially future viral outbreaks of similar magnitude. In particular mutations and deletions involving the virion's prefusion Spike protein have significant potential impact on vaccines and therapeutics that utilize this key structural viral protein in their mitigation strategies. In this study, we have demonstrated how dominant energetic landscape mappings ("glue points") coupled with sequence alignment information can potentially identify or flag key residue mutations and deletions associated with variants. Surprisingly, we also found excellent homology of stabilizing residue glue points across the lineage of ß coronavirus Spike proteins, and we have termed this as "sequence homologous glue points". In general, these flagged residue mutations and/or deletions are then computationally studied in detail using all-atom biocomputational molecular dynamics over approximately one microsecond in order to ascertain structural and energetic changes in the Spike protein associated variants. Specifically, we examined both a theoretically-based triple mutant and the so-called UK or B.1.1.7 variant. For the theoretical triple mutant, we demonstrated through Alanine mutations, which help "unglue" key residue-residue interactions, that these three key stabilizing residues could cause the transition of Down to Up protomer states, where the Up protomer state allows binding of the prefusion Spike protein to hACE2 host cell receptors, whereas the Down state is believed inaccessible. Thus, we are able to demonstrate the importance of glue point residue identification in the overall stability of the prefusion Spike protein. For the B.1.1.7 variant, we demonstrated the critical importance of D614G and N5017 on the structure and binding, respectively, of the Spike protein. Notably, we had previously identified D614 as a key glue point in the inter-protomer stabilization of the Spike protein prior to the emergence of its mutation. The mutant D614G is a structure breaking Glycine mutation demonstrating a relatively more distal Down state RBD and a more stable conformation in general. In addition, we demonstrate that the mutation N501Y may significantly increase the Spike protein binding to hACE2 cell receptors through its interaction with Y41 of hACE2 forming a potentially strong hydrophobic residue binding pair. We note that these two key mutations, D614G and N501Y, are also found in the so-called South African (SA; B.1.351) variant of SARS-CoV-2. Future studies along these lines are, therefore, aimed at mapping glue points to residue mutations and deletions of associated prefusion Spike protein variants in order to help identify and analyze possible "variants of interest" and optimize efforts aimed at the mitigation of this current and future virions.

Static all-atom energetic mappings of the SARS-Cov-2 spike protein and dynamic stability analysis of "Up" versus "Down" protomer states.

The SARS-CoV-2 virion responsible for the current world-wide pandemic COVID-19 has a characteristic Spike protein (S) on its surface that embellishes both a prefusion state and fusion state. The prefusion Spike protein (S) is a large trimeric protein where each protomer may be in a so-called Up state or Down state, depending on the configuration of its receptor binding domain (RBD) within its distal, prefusion S1 domain. The Up state is believed to allow binding of the virion to ACE-2 receptors on human epithelial cells, whereas the Down state is believed to be relatively inactive or reduced in its binding behavior. We have performed detailed all-atom, dominant energy landscape mappings for noncovalent interactions (charge, partial charge, and van der Waals) of the SARS-CoV-2 Spike protein in its static prefusion state based on two recent and independent experimental structure publications. We included both interchain interactions and intrachain (domain) interactions in our mappings in order to determine any telling differences (different so-called "glue" points) between residues in the Up and Down state protomers. The S2 proximal, fusion domain demonstrated no appreciable energetic differences between Up and Down protomers, including interchain as well as each protomer's intrachain, S1-S2 interactions. However, the S1 domain interactions across neighboring protomers, which include the RBD-NTD cross chain interactions, showed significant energetic differences between Up-Down and Down-Down neighboring protomers. This included, for example, a key RBD residue ARG357 in the Up-Down interaction and a three residue sequence ALA520-PRO521-ALA522, associated with a turn structure in the RBD of the Up state protomer, acting as a stabilizing interaction with the NTD of its neighbor protomer. Additionally, our intra chain dominant energy mappings within each protomer, identified a significant "glue" point or possible "latch" for the Down state protomer between the S1 subdomain, SD1, and the RBD domain of the same protomer that was completely missing in the Up state protomer analysis. Ironically, this dominant energetic interaction in the Down state protomer involved the backbone atoms of the same three residue sequence ALA520-PRO521-ALA522 of the RBD with the amino acid R-group of GLN564 in the SD1 domain. Thus, this same three residue sequence acts as a stabilizer of the RBD in the Up conformation through its interactions with its neighboring NTD chain and a kind of latch in the Down state conformation through its interactions with its own SD1 domain. The dominant interaction energy residues identified here are also conserved across reported variations of SARS-CoV-2, as well as the closely related virions SARS-Cov and the bat corona virus RatG13. We conducted preliminary molecular dynamics simulations across 0.1 µ seconds to see if this latch provided structural stability and indeed found that a single point mutation (Q564G) resulted in the latch releasing transforming the protomer from the Down to the Up state conformation. Full trimeric Spike protein studies of the same mutation across all protomers, however, did not exhibit latch release demonstrating the critical importance of interchain interactions across the S1 domain, including RBD-NTD neighboring chain interactions. Therapies aimed at disrupting these noncovalent interactions could be a viable route for the physico-chemical mitigation of this deadly virion.

Static All-Atom Energetic Mappings of the SARS-Cov-2 Spike Protein with Potential Latch Identification of the Down State Protomer.

The SARS-Cov-2 virion responsible for the current world-wide pandemic Covid-19 has a characteristic Spike protein (S) on its surface that embellishes both a prefusion state and fusion state. The prefusion Spike protein (S) is a large trimeric protein where each protomer may be in a so-called Up state or Down state, depending on the configuration of its receptor binding domain (RBD). The Up state is believed to allow binding of the virion to ACE-2 receptors on human epithelial cells, whereas the Down state is believed to be relatively inactive or reduced in its binding behavior. We have performed detailed all-atom, dominant energy landscape mappings for noncovalent interactions (charge, partial charge, and van der Waals) of the SARS-Cov-2 Spike protein in its static prefusion state based on recent structural information. We included both interchain interactions and intrachain (domain) interactions in our mappings in order to determine any telling differences (different so-called "glue" points) between residues in the Up and Down state protomers. In general, the S2 or fusion machinery domain of S is relatively rigid with strong noncovalent interactions facilitated by helical secondary structures, whereas the S1 domain, which contains the RBD and N-terminal domain (NTD), is relatively more flexible and characterized by beta strand structural motifs. The S2 domain demonstrated no appreciable energetic differences between Up and Down protomers, including interchain as well as each protomer's intrachain, S1-S2 interactions. However, the S1 domain interactions across neighboring protomers, which include the RBD-NTD cross chain interactions, showed significant energetic differences between Up-Down and Down-Down neighboring protomers. Surprisingly, the Up-Down, RBD-NTD interactions were overall stronger and more numerous than the Down-Down cross chain interactions, including the appearance of the three residue sequence ALA520-PRO521-ALA522 associated with a turn structure in the RBD of the Up state protomer. Additionally, our intrachain dominant energy mappings within each protomer, identified a significant "glue" point or possible "latch" for the Down state protomer between the S1 subdomain, SD1, and the RBD domain of the same protomer that was completely missing in the Up state protomer analysis. Ironically, this dominant energetic interaction in the Down state protomer involved the backbone atoms of the same three residue sequence ALA520-PRO521-ALA522 of the RBD with the R-group of GLN564 in the SD1 domain. Thus, this same three residue sequence acts as a stabilizer of the RBD in the Up conformation through its interactions with its neighboring NTD chain and a kind of latch in the Down state conformation through its interactions with its own SD1 domain. The dominant interaction energy residues identified here are also conserved across reported variations of SARS-Cov-2, as well as the closely related virions SARS-Cov and the bat corona virus RatG13. To help verify the potential latch for the Down state protomer, we conducted some preliminary molecular dynamic simulations that effectively turn off this specific latch glue point via a single point mutation of GLN564. Interestingly, the single point mutation lead to the latch releasing in less than a few nanoseconds, but the latch remained fixed in the wild state protomer for up to 0.1 microseconds that were simulated. Many more detailed studies are needed to understand the dynamics of the Up and Down states of the Spike protein, including the stabilizing chain-chain interactions and the mechanisms of transition from Down to Up state protomers. Nonetheless, static dominant energy landscape mappings and preliminary molecular dynamic studies given here may represent a useful starting point for more detailed dynamic analyses and hopefully an improved understanding of the structure-function relationship of this highly complex protein associated with COVID-19.

An Algorithmic Approach for Detecting Bolides with the Geostationary Lightning Mapper.

The Geostationary Lightning Mapper (GLM) instrument onboard the GOES 16 and 17 satellites can be used to detect bolides in the atmosphere. This capacity is unique because GLM provides semi-global, continuous coverage and releases its measurements publicly. Here, six filters are developed that are aggregated into an automatic algorithm to extract bolide signatures from the GLM level 2 data product. The filters exploit unique bolide characteristics to distinguish bolide signatures from lightning and other noise. Typical lightning and bolide signatures are introduced and the filter functions are presented. The filter performance is assessed on 144845 GLM L2 files (equivalent to 34 days-worth of data) and the algorithm selected 2252 filtered files (corresponding to a pass rate of 1.44%) with bolide-similar signatures. The challenge of identifying frequent but small, decimeter-sized bolide signatures is discussed as GLM reaches its resolution limit for these meteors. The effectiveness of the algorithm is demonstrated by its ability to extract confirmed and new bolide discoveries. We provide discovery numbers for November 2018 when seven likely bolides were discovered of which four are confirmed by secondary observations. The Cuban meteor on Feb 1st 2019 serves as an additional example to demonstrate the algorithms capability and the first light curve as well as correct ground track was available within 8.5 hours based on GLM data for this event. The combination of the automatic bolide extraction algorithm with GLM can provide a wealth of new measurements of bolides in Earth's atmosphere to enhance the study of asteroids and meteors.

PLANETARY CANDIDATES OBSERVED BY Kepler. VIII. A FULLY AUTOMATED CATALOG WITH MEASURED COMPLETENESS AND RELIABILITY BASED ON DATA RELEASE 25.

We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new in this catalog and include two new candidates in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten new high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs) found by the Kepler Pipeline (Twicken et al. 2016). Because of this automation, we were also able to vet simulated data sets and therefore measure how well the Robovetter separates those TCEs caused by noise from those caused by low signal-to-noise transits. Because of these measurements we fully expect that this catalog can be used to accurately calculate the frequency of planets out to Kepler's detection limit, which includes temperate, super-Earth size planets around GK dwarf stars in our Galaxy. This paper discusses the Robovetter and the metrics it uses to decide which TCEs are called planet candidates in the DR25 KOI catalog. We also discuss the simulated transits, simulated systematic noise, and simulated astrophysical false positives created in order to characterize the properties of the final catalog. For orbital periods less than 100 d the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates found between 200 and 500 days, our measurements indicate that the Robovetter is 73.5% complete and 37.2% reliable across all searched stars (or 76.7% complete and 50.5% reliable when considering just the FGK dwarf stars). We describe how the measured completeness and reliability varies with period, signal-to-noise, number of transits, and stellar type. Also, we discuss a value called the disposition score which provides an easy way to select a more reliable, albeit less complete, sample of candidates. The entire KOI catalog, the transit fits using Markov chain Monte Carlo methods, and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive.

Exploring the SDSS Dataset with Linked Scatter Plots: I. EMP, CEMP, and CV Stars.

We present the results of a search for EMP, CEMP, and cataclysmic variable stars using a new exploration tool based on linked scatter plots (LSPs). Our approach is especially designed to work with very large spectrum data sets such as the SDSS, LAMOST, RAVE, and Gaia data sets and can be applied to stellar, galaxy, and quasar spectra. As a demonstration, we conduct a search for EMP, CEMP, and cataclysmic variable stars in the SDSS DR10 data set. We first created a 3326-dimensional phase space containing nearly 2 billion measures of the strengths of over 1600 spectral features in 569,738 SDSS stars. These measures capture essentially all the stellar atomic and molecular species visible at the resolution of SDSS spectra. We show how LSPs can be used to quickly isolate and examine interesting portions of this phase space. To illustrate, we use LSPs coupled with cuts in selected portions of phase space to extract EMP stars, C-rich EMP stars, and CV stars. We present identifications for 59 previously unrecognized candidate EMP stars and 11 previously unrecognized candidate CEMP stars. We also call attention to 2 candidate He II emission CV stars found by the LSP approach that have not yet been discussed in the literature.

Volume-rendering on a 3D hyperwall: A molecular visualization platform for research, education and outreach.

We present a unique platform for molecular visualization and design that uses novel subatomic feature detection software in tandem with 3D hyperwall visualization technology. We demonstrate the fleshing-out of pharmacophores in drug molecules, as well as reactive sites in catalysts, focusing on subatomic features. Topological analysis with picometer resolution, in conjunction with interactive volume-rendering of the Laplacian of the electronic charge density, leads to new insight into docking and catalysis. Visual data-mining is done efficiently and in parallel using a 4×4 3D hyperwall (a tiled array of 3D monitors driven independently by slave GPUs but displaying high-resolution, synchronized and functionally-related images). The visual texture of images for a wide variety of molecular systems are intuitive to experienced chemists but also appealing to neophytes, making the platform simultaneously useful as a tool for advanced research as well as for pedagogical and STEM education outreach purposes.

An Earth-sized planet in the habitable zone of a cool star.

The quest for Earth-like planets is a major focus of current exoplanet research. Although planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surfaces. We present the detection of Kepler-186f, a 1.11 ± 0.14 Earth-radius planet that is the outermost of five planets, all roughly Earth-sized, that transit a 0.47 ± 0.05 solar-radius star. The intensity and spectrum of the star's radiation place Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and water at its surface, then some of this water is likely to be in liquid form.

The same frequency of planets inside and outside open clusters of stars.

Most stars and their planets form in open clusters. Over 95 per cent of such clusters have stellar densities too low (less than a hundred stars per cubic parsec) to withstand internal and external dynamical stresses and fall apart within a few hundred million years. Older open clusters have survived by virtue of being richer and denser in stars (1,000 to 10,000 per cubic parsec) when they formed. Such clusters represent a stellar environment very different from the birthplace of the Sun and other planet-hosting field stars. So far more than 800 planets have been found around Sun-like stars in the field. The field planets are usually the size of Neptune or smaller. In contrast, only four planets have been found orbiting stars in open clusters, all with masses similar to or greater than that of Jupiter. Here we report observations of the transits of two Sun-like stars by planets smaller than Neptune in the billion-year-old open cluster NGC6811. This demonstrates that small planets can form and survive in a dense cluster environment, and implies that the frequency and properties of planets in open clusters are consistent with those of planets around field stars in the Galaxy.

Kepler-62: a five-planet system with planets of 1.4 and 1.6 Earth radii in the habitable zone.

We present the detection of five planets--Kepler-62b, c, d, e, and f--of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii (R⊕), orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4, and 267.3 days, respectively. The outermost planets, Kepler-62e and -62f, are super-Earth-size (1.25 R⊕ < planet radius ≤ 2.0 R⊕) planets in the habitable zone of their host star, respectively receiving 1.2 ± 0.2 times and 0.41 ± 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 billion years suggest that both planets could be solid, either with a rocky composition or composed of mostly solid water in their bulk.

A sub-Mercury-sized exoplanet.

Since the discovery of the first exoplanets, it has been known that other planetary systems can look quite unlike our own. Until fairly recently, we have been able to probe only the upper range of the planet size distribution, and, since last year, to detect planets that are the size of Earth or somewhat smaller. Hitherto, no planets have been found that are smaller than those we see in the Solar System. Here we report a planet significantly smaller than Mercury. This tiny planet is the innermost of three that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of the Moon, and highly irradiated surface, the planet, Kepler-37b, is probably rocky with no atmosphere or water, similar to Mercury.

Two Earth-sized planets orbiting Kepler-20.

Since the discovery of the first extrasolar giant planets around Sun-like stars, evolving observational capabilities have brought us closer to the detection of true Earth analogues. The size of an exoplanet can be determined when it periodically passes in front of (transits) its parent star, causing a decrease in starlight proportional to its radius. The smallest exoplanet hitherto discovered has a radius 1.42 times that of the Earth's radius (R(⊕)), and hence has 2.9 times its volume. Here we report the discovery of two planets, one Earth-sized (1.03R(⊕)) and the other smaller than the Earth (0.87R(⊕)), orbiting the star Kepler-20, which is already known to host three other, larger, transiting planets. The gravitational pull of the new planets on the parent star is too small to measure with current instrumentation. We apply a statistical method to show that the likelihood of the planetary interpretation of the transit signals is more than three orders of magnitude larger than that of the alternative hypothesis that the signals result from an eclipsing binary star. Theoretical considerations imply that these planets are rocky, with a composition of iron and silicate. The outer planet could have developed a thick water vapour atmosphere.

Intracellular localization of a group II chaperonin indicates a membrane-related function.

Chaperonins are protein complexes that are believed to function as part of a protein folding system in the cytoplasm of the cell. We observed, however, that the group II chaperonins known as rosettasomes in the hyperthermophilic archaeon Sulfolobus shibatae, are not cytoplasmic but membrane associated. This association was observed in cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C by using immunofluorescence microscopy and in thick sections of rapidly frozen cells grown at 76 degrees C by using immunogold electron microscopy. We observed that increased abundance of rosettasomes after heat shock correlated with decreased membrane permeability at lethal temperature (92 degrees C). This change in permeability was not seen in cells heat-shocked in the presence of the amino acid analogue azetidine 2-carboxylic acid, indicating functional protein synthesis influences permeability. Azetidine experiments also indicated that observed heat-induced changes in lipid composition in S. shibatae could not account for changes in membrane permeability. Rosettasomes purified from cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C bind to liposomes made from either the bipolar tetraether lipids of Sulfolobus or a variety of artificial lipid mixtures. The presence of rosettasomes did not significantly change the transition temperature of liposomes, as indicated by differential scanning calorimetry, or the proton permeability of liposomes, as indicated by pyranine fluorescence. We propose that these group II chaperonins function as a structural element in the natural membrane based on their intracellular location, the correlation between their functional abundance and membrane permeability, and their potential distribution on the membrane surface.