The molecular dynamics of possible inhibitors for SARS-CoV-2.

The novel coronavirus SARS-CoV-2, responsible for the present COVID-19 global pandemic, is known to bind to the angiotensin converting enzyme-2 (ACE2) receptor in human cells. A possible treatment of COVID-19 could involve blocking ACE2 and/or disabling the spike protein on the virus. Here, molecular dynamics simulations were performed to test the binding affinities of nine candidate compounds. Of these, three drugs showed significant therapeutic potential that warrant further investigation: SN35563, a ketamine ester analogue, was found to bind strongly to the ACE2 receptor but weakly within the spike receptor-binding domain (RBD); in contrast, arbidol and hydroxychloroquine bound preferentially with the spike RBD rather than ACE2. A fourth drug, remdesivir, bound approximately equally to both the ACE2 and viral spike RBD, thus potentially increasing risk of viral infection by bringing the spike protein into closer proximity to the ACE2 receptor. We suggest more experimental investigations to test that SN35563-in combination with arbidol or hydroxychloroquine-might act synergistically to block viral cell entry by providing therapeutic blockade of the host ACE2 simultaneous with reduction of viral spike receptor-binding; and that this combination therapy would allow the use of smaller doses of each drug.Communicated by Ramaswamy H. Sarma.

Polysaccharide Structures in the Outer Mucilage of Arabidopsis Seeds Visualized by AFM.

Evidence is presented that the polysaccharide rhamnogalacturonan I (RGI) can be biosynthesized in remarkably organized branched configurations and surprisingly long versions and can self-assemble into a plethora of structures. AFM imaging has been applied to study the outer mucilage obtained from wild-type (WT) and mutant (bxl1-3 and cesa5-1) Arabidopsis thaliana seeds. For WT mucilage, ordered, multichain structures of the polysaccharide RGI were observed, with a helical twist visible in favorable circumstances. Molecular dynamics (MD) simulations demonstrated the stability of several possible multichain complexes and the possibility of twisted fibril formation. For bxl1-3 seeds, the imaged polymers clearly showed the presence of side chains. These were surprisingly regular and well organized with an average length of ∼100 nm and a spacing of ∼50 nm. The heights of the side chains imaged were suggestive of single polysaccharide chains, while the backbone was on average 4 times this height and showed regular height variations along its length consistent with models of multichain fibrils examined in MD. Finally, in mucilage extracts from cesa5-1 seeds, a minor population of chains in excess of 30 µm long was observed.

On the electrophoretic mobilities of partially charged oligosaccharides as a function of charge patterning and degree of polymerization.

Fully or partially charged oligosaccharide molecules play a key role in many areas of biology, where their fine structures are crucial in determining their functionality. However, the separation of specific charged oligosaccharides from similar moieties that typically coexist in extracted samples, even for those that are unbranched, and in cases where each saccharide moiety can only carry a single charge or not, is far from trivial. Typically such molecules are characterized by a degree of polymerization n and a number m (and distribution) of charged residues, and must be separated from a plethora of similar species possessing different combinations of n and m. Furthermore, the separation of the possible n!/m!(n-m)! isomers of each species of fixed n and m is a formidable challenge to analytical chemists. Herein, we report the results of molecular dynamics simulations that have been performed in order to calculate the free solution electrophoretic mobilities of galacturonides and charged oligosaccharides derived from digests of the important plant cell-wall polysaccharide pectin. The simulations are compared with an experiment and are found to correctly predict the loss of resolution of fully charged species above a critical degree of polymerization n and the ionic strength dependence of the electrophoretic mobilities of different partially charged oligosaccharides. It is expected that having a predictive tool for the calculation of the electrophoretic mobilities of differently charged oligosaccharide species in hand will allow experimental conditions that optimize the resolution of particular species to be ascertained and understood.

Molecular Dynamics Simulations Illuminate the Role of Counterion Condensation in the Electrophoretic Transport of Homogalacturonans.

Homogalacturonans (HGs) are polysaccharide copolymers of galacturonic acid and its methylesterified counterpart. The inter- and intramolecular distributions of the methylesterifed residues are vital behavior-determining characteristics of a sample's structure, and much experimental effort has been directed to their measurement. While many techniques are able to measure the sample-averaged degree of methylesterification (DM), the measurement of inter- and intramolecular charge distributions are challenging. Here, molecular dynamics (MD) simulations are used to calculate the electrophoretic mobilities of HGs that have different amounts and distributions of charges placed along the backbone. The simulations are shown to capture experimental results well, even for low-DM samples that possess high charge densities. In addition, they illuminate the role that local counterion condensation can play in the determination of the electrophoretic mobility of heterogeneous blocky polyelectrolytes that cannot be adequately described by a single chain-averaged charge spacing.

Ab initio parameterization of an all-atom polarizable and dissociable force field for water.

A novel all-atom, dissociative, and polarizable force field for water is presented. The force field is parameterized based on forces, stresses, and energies obtained form ab initio calculations of liquid water at ambient conditions. The accuracy of the force field is tested by calculating structural and dynamical properties of liquid water and the energetics of small water clusters. The transferability of the force field to dissociated states is studied by considering the solvation of a proton and the ionization of water at extreme conditions of pressure and temperature. In the case of the solvated proton, the force field properly describes the presence of both Eigen and Zundel configurations. In the case of the pressure-induced ice VIII/ice X transition and the temperature-induced transition to a superionic phase, the force field is found to describe accurately the proton symmetrization and the melting of the proton sublattice, respectively.