Search for: All records

HD 106906 is a planetary system that hosts a wide-orbit companion, as well as an eccentric and flat debris disk, which hold important constraints on its formation and subsequent evolution. The recent observations of the companion constrain its orbit to be eccentric and inclined relative to the plane of the debris disk. Here, we show that, in the presence of the inclined companion, the debris disk quickly (≲5 Myr) becomes warped and puffy. This suggests that the current configuration of the system is relatively recent. We explore the possibility that a recent close encounter with a free-floating planet could produce a companion with orbital parameters that agree with observations of HD 106906 b. We find that this scenario is able to recreate the structure of the debris disk while producing a companion in agreement with observation.

 

Cyanuric acid (CA), a triazine heterocycle, is extensively utilized for noncovalent self‐assembly. The association between poly(adenine) and CA into micron‐length fibers was a remarkable observation made by Sleiman and co‐workers, who proposed that adenine and CA adopt a hexameric rosette configuration in analogy with previously reported structures for CA assemblies. However, recent experimental observations from the Krishnamurthy group led to a reevaluation of the hexameric rosette model, wherein they have proposed a hydrogen‐bonded helicene model as an alternative. Our molecular dynamics simulations show that the hexad model is indeed unlikely and that this novel noncovalent helicene geometry, where the adenine and CA bases form an extended helical hydrogen‐bond network across the system, is a more probable structural motif. The existence of noncovalent helicene compounds may have wide‐ranging applications in DNA nanotechnology and helicene chemistry.

 

The many-body expansion (MBE) is promising for the efficient, parallel computation of lattice energies in organic crystals. Very high accuracy should be achievable by employing coupled-cluster singles, doubles, and perturbative triples at the complete basis set limit [CCSD(T)/CBS] for the dimers, trimers, and potentially tetramers resulting from the MBE, but such a brute-force approach seems impractical for crystals of all but the smallest molecules. Here, we investigate hybrid or multi-level approaches that employ CCSD(T)/CBS only for the closest dimers and trimers and utilize much faster methods like Møller–Plesset perturbation theory (MP2) for more distant dimers and trimers. For trimers, MP2 is supplemented with the Axilrod–Teller–Muto (ATM) model of three-body dispersion. MP2(+ATM) is shown to be a very effective replacement for CCSD(T)/CBS for all but the closest dimers and trimers. A limited investigation of tetramers using CCSD(T)/CBS suggests that the four-body contribution is entirely negligible. The large set of CCSD(T)/CBS dimer and trimer data should be valuable in benchmarking approximate methods for molecular crystals and allows us to see that a literature estimate of the core-valence contribution of the closest dimers to the lattice energy using just MP2 was overbinding by 0.5 kJ mol−1, and an estimate of the three-body contribution from the closest trimers using the T0 approximation in local CCSD(T) was underbinding by 0.7 kJ mol−1. Our CCSD(T)/CBS best estimate of the 0 K lattice energy is −54.01 kJ mol−1, compared to an estimated experimental value of −55.3 ± 2.2 kJ mol−1. 

Abstract The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer’s receptor binding domain (RBD) switching from a "down” (closed) to an "up” (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μ s of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody’s epitope is continuously exposed, explaining its high efficacy. 

ABSTRACT TonB-dependent transporters (TDTs) are essential proteins for metal acquisition, an important step in the growth and pathogenesis of many pathogens, including Neisseria gonorrhoeae , the causative agent of gonorrhea. There is currently no available vaccine for gonorrhea; TDTs are being investigated as vaccine candidates because they are highly conserved and expressed in vivo . Transferrin binding protein A (TbpA) is an essential virulence factor in the initiation of experimental infection in human males and functions by acquiring iron upon binding to host transferrin (human transferrin [hTf]). The loop 3 helix (L3H) is a helix finger that inserts into the hTf C-lobe and is required for hTf binding and subsequent iron acquisition. This study identified and characterized the first TbpA single-point substitutions resulting in significantly decreased hTf binding and iron acquisition, suggesting that the helix structure is more important than charge for hTf binding and utilization. The tbpA D355P Δ tbpB and tbpA A356P Δ tbpB mutants demonstrated significantly reduced hTf binding and impaired iron uptake from Fe-loaded hTf; however, only the tbpA A356P Δ tbpB mutant was able to grow when hTf was the sole source of iron. The expression of tbpB was able to restore function in all tbpA mutants. These results implicate both D355 and A356 in the key binding, extraction, and uptake functions of gonococcal TbpA. 

Routinely assessing the stability of molecular crystals with high accuracy remains an open challenge in the computational sciences. The many-body expansion decomposes computation of the crystal lattice energy into an embarrassingly parallel collection of computations over molecular dimers, trimers, and so forth, making quantum chemistry techniques tractable for many crystals of small organic molecules. By examining the range-dependence of different types of energetic contributions to the crystal lattice energy, we can glean qualitative understanding of solid-state intermolecular interactions as well as practical, exploitable reductions in the number of computations required for accurate energies. Here, we assess the range-dependent character of two-body interactions of 24 small organic molecular crystals by using the physically interpretable components from symmetry-adapted perturbation theory (electrostatics, exchange-repulsion, induction/polarization, and London dispersion). We also examine correlations between the convergence rates of electrostatics and London dispersion terms with molecular dipole moments and polarizabilities, to provide guidance for estimating convergence rates in other molecular crystals.