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.

A nearby transiting rocky exoplanet that is suitable for atmospheric investigation

Spectroscopy of transiting exoplanets can be used to investigate their atmospheric properties and habitability. Combining radial velocity (RV) and transit data provides additional information on exoplanet physical properties. We detect a transiting rocky planet with an orbital period of 1.467 days around the nearby red dwarf star Gliese 486. The planet Gliese 486 b is 2.81 Earth masses and 1.31 Earth radii, with uncertainties of 5%, as determined from RV data and photometric light curves. The host star is at a distance of ~8.1 parsecs, has a J-band magnitude of ~7.2, and is observable from both hemispheres of Earth. On the basis of these properties and the planet's short orbital period and high equilibrium temperature, we show that this terrestrial planet is suitable for emission and transit spectroscopy.

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.

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.

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.

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.