Modeling interactions between rubidium atom and magnetometer cell wall molecules.

Magnetometer cell wall coat molecules play an important role in preserving the lifetime of pumped alkali metal atoms for use in magnetometers that are capable of measuring very small magnetic fields. The goal of this study is to help rationalize the design of the cell coat molecules. Rubidium-87 is studied in terms of its interaction with three template cell coat molecules: ethane, ethene, and methyltrichlorosilane (MeTS). Ab initio electronic structure methods are applied to investigate the effect that the coat molecules have on the 2S ground state and 2P excited state of 87Rb. We find that, from the ab initio results, the three template molecules have differing effects, with MeTS having the largest effect on the ground state and ethane or ethene having the largest effect on the non-degenerate excited states.

Cysteine-Selective Modification of Peptides and Proteins via Desulfurative C-C Bond Formation.

The site-selective modification of peptides and proteins facilitates the preparation of targeted therapeutic agents and tools to interrogate biochemical pathways. Among the numerous bioconjugation techniques developed to install groups of interest, those that generate C(sp3 )-C(sp3 ) bonds are significantly underrepresented despite affording proteolytically stable, biogenic linkages. Herein, a visible-light-mediated reaction is described that enables the site-selective modification of peptides and proteins via desulfurative C(sp3 )-C(sp3 ) bond formation. The reaction is rapid and high yielding in peptide systems, with comparable translation to proteins. Using this chemistry, a range of moieties is installed into model systems and an effective PTM-mimic is successfully integrated into a recombinantly expressed histone.

QuantumScents: Quantum-Mechanical Properties for 3.5k Olfactory Molecules.

Quantitative structure-odor relationships are critically important for studies related to the function of olfaction. Current literature data sets contain expert-labeled molecules but lack feature data. This paper introduces QuantumScents, a quantum mechanics augmented derivative of the Leffingwell data set. QuantumScents contains 3.5k structurally and chemically diverse molecules ranging from 2 to 30 heavy atoms (CNOS) and their corresponding 3D coordinates, total PBE0 energy, molecular dipole moment, and per-atom Hirshfeld charges, dipoles, and ratios. The authors demonstrate that Hirshfeld charges and ratios contain sufficient information to perform molecular classification by training a Message Passing Neural Network with chemprop (Heid, E.; et al. ChemRxiv, 2023, DOI: 10.26434/chemrxiv-2023-3zcfl) to predict scent labels. The QuantumScents data set is freely available on Zenodo along with the authors' code, example models, and data set generation workflow (https://zenodo.org/doi/10.5281/zenodo.8239853).

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARS-CoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. These results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes.

With a longer-term goal of addressing the comparative behavior of the aqueous halides F-, Cl-, Br-, and I- on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl-, Br-, and I-, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

Near-Ultraviolet Circular Dichroism and Two-Dimensional Spectroscopy of Polypeptides.

A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper insights into biophysical and simulation studies of protein dynamics and folding. In contrast to intense bands in the far-ultraviolet, near-UV bands are much weaker and have been challenging to compute theoretically. We report some advances in the accuracy of calculations in the near-UV, which were realised through the consideration of the vibrational structure of the electronic transitions of aromatic side chains.

Supercomputing Pipelines Search for Therapeutics Against COVID-19.

The urgent search for drugs to combat SARS-CoV-2 has included the use of supercomputers. The use of general-purpose graphical processing units (GPUs), massive parallelism, and new software for high-performance computing (HPC) has allowed researchers to search the vast chemical space of potential drugs faster than ever before. We developed a new drug discovery pipeline using the Summit supercomputer at Oak Ridge National Laboratory to help pioneer this effort, with new platforms that incorporate GPU-accelerated simulation and allow for the virtual screening of billions of potential drug compounds in days compared to weeks or months for their ability to inhibit SARS-COV-2 proteins. This effort will accelerate the process of developing drugs to combat the current COVID-19 pandemic and other diseases.

Protein Conformational States-A First Principles Bayesian Method.

Automated identification of protein conformational states from simulation of an ensemble of structures is a hard problem because it requires teaching a computer to recognize shapes. We adapt the naïve Bayes classifier from the machine learning community for use on atom-to-atom pairwise contacts. The result is an unsupervised learning algorithm that samples a 'distribution' over potential classification schemes. We apply the classifier to a series of test structures and one real protein, showing that it identifies the conformational transition with >95% accuracy in most cases. A nontrivial feature of our adaptation is a new connection to information entropy that allows us to vary the level of structural detail without spoiling the categorization. This is confirmed by comparing results as the number of atoms and time-samples are varied over 1.5 orders of magnitude. Further, the method's derivation from Bayesian analysis on the set of inter-atomic contacts makes it easy to understand and extend to more complex cases.

Insertion of Dengue E into lipid bilayers studied by neutron reflectivity and molecular dynamics simulations.

The envelope (E) protein of Dengue virus rearranges to a trimeric hairpin to mediate fusion of the viral and target membranes, which is essential for infectivity. Insertion of E into the target membrane serves to anchor E and possibly also to disrupt local order within the membrane. Both aspects are likely to be affected by the depth of insertion, orientation of the trimer with respect to the membrane normal, and the interactions that form between trimer and membrane. In the present work, we resolved the depth of insertion, the tilt angle, and the fundamental interactions for the soluble portion of Dengue E trimers (sE) associated with planar lipid bilayer membranes of various combinations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-rac-glycerol (POPG), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and cholesterol (CHOL) by neutron reflectivity (NR) and by molecular dynamics (MD) simulations. The results show that the tip of E containing the fusion loop (FL) is located at the interface of the headgroups and acyl chains of the outer leaflet of the lipid bilayers, in good agreement with prior predictions. The results also indicate that E tilts with respect to the membrane normal upon insertion, promoted by either the anionic lipid POPG or CHOL. The simulations show that tilting of the protein correlates with hydrogen bond formation between lysines and arginines located on the sides of the trimer close to the tip (K246, K247, and R73) and nearby lipid headgroups. These hydrogen bonds provide a major contribution to the membrane anchoring and may help to destabilize the target membrane.

Extension of Kirkwood-Buff theory to the canonical ensemble.

Kirkwood-Buff (KB) integrals are notoriously difficult to converge from a canonical simulation because they require estimating the grand-canonical radial distribution. The same essential difficulty is encountered when attempting to estimate the direct correlation function of Ornstein-Zernike theory by inverting the pair correlation functions. We present a new theory that applies to the entire, finite, simulation volume, so that no cutoff issues arise at all. The theory gives the direct correlation function for closed systems, while smoothness of the direct correlation function in reciprocal space allows calculating canonical KB integrals via a well-posed extrapolation to the origin. The present analysis method represents an improvement over previous work because it makes use of the entire simulation volume and its convergence can be accelerated using known properties of the direct correlation function. Using known interaction energy functions can make this extrapolation near perfect accuracy in the low-density case. Because finite size effects are stronger in the canonical than in the grand-canonical ensemble, we state ensemble correction formulas for the chemical potential and the KB coefficients. The new theory is illustrated with both analytical and simulation results on the 1D Ising model and a supercritical Lennard-Jones fluid. For the latter, the finite-size corrections are shown to be small.

Ultrafast photodissociation dynamics of 1,4-diiodobenzene.

The photodissociation dynamics of 1,4-diiodobenzene is investigated using ultrafast time-resolved photoelectron spectroscopy. Following excitation by laser pulses at 271 nm, the excited-state dynamics is probed by resonance-enhanced multiphoton ionization with 405 nm probe pulses. A progression of Rydberg states, which come into resonance sequentially, provide a fingerprint of the dissociation dynamics of the molecule. The initial excitation decays with a lifetime of 33 ± 4 fs, in good agreement with a previous study. The spectrum is interpreted by reference to ab initio calculations at the CASPT2(18,14) level, including spin-orbit coupling. We propose that both the 5B1 and 6B1 states are excited initially, and based on the calculations, we identify diabatic spin-orbit coupled states corresponding to the main dissociation pathways.

Method for measuring the unbinding energy of strongly-bound membrane-associated proteins.

We describe a new method to measure the activation energy for unbinding (enthalpy ΔH*u and free energy ΔG*u) of a strongly-bound membrane-associated protein from a lipid membrane. It is based on measuring the rate of release of a liposome-bound protein during centrifugation on a sucrose gradient as a function of time and temperature. The method is used to determine ΔH*u and ΔG*u for the soluble dengue virus envelope protein (sE) strongly bound to 80:20 POPC:POPG liposomes at pH5.5. ΔH*u is determined from the Arrhenius equation whereas ΔG*u is determined by fitting the data to a model based on mean first passage time for escape from a potential well. The binding free energy ΔGb of sE was also measured at the same pH for the initial, predominantly reversible, phase of binding to a 70:30 PC:PG lipid bilayer. The unbinding free energy (20±3kcal/mol, 20% PG) was found to be roughly three times the binding energy per monomer, (7.8±0.3kcal/mol for 30% PG, or est. 7.0kcal/mol for 20% PG). This is consistent with data showing that free sE is a monomer at pH5.5, but assembles into trimers after associating with membranes. This new method to determine unbinding energies should be useful to understand better the complex interactions of integral monotopic proteins and strongly-bound peripheral membrane proteins with lipid membranes.

Molecular basis of endosomal-membrane association for the dengue virus envelope protein.

Dengue virus is coated by an icosahedral shell of 90 envelope protein dimers that convert to trimers at low pH and promote fusion of its membrane with the membrane of the host endosome. We provide the first estimates for the free energy barrier and minimum for two key steps in this process: host membrane bending and protein-membrane binding. Both are studied using complementary membrane elastic, continuum electrostatics and all-atom molecular dynamics simulations. The predicted host membrane bending required to form an initial fusion stalk presents a 22-30 kcal/mol free energy barrier according to a constrained membrane elastic model. Combined continuum and molecular dynamics results predict a 15 kcal/mol free energy decrease on binding of each trimer of dengue envelope protein to a membrane with 30% anionic phosphatidylglycerol lipid. The bending cost depends on the preferred curvature of the lipids composing the host membrane leaflets, while the free energy gained for protein binding depends on the surface charge density of the host membrane. The fusion loop of the envelope protein inserts exactly at the level of the interface between the membrane's hydrophobic and head-group regions. The methods used in this work provide a means for further characterization of the structures and free energies of protein-assisted membrane fusion.

The glutaminase activity of L-asparaginase is not required for anticancer activity against ASNS-negative cells.

L-Asparaginase (L-ASP) is a key component of therapy for acute lymphoblastic leukemia. Its mechanism of action, however, is still poorly understood, in part because of its dual asparaginase and glutaminase activities. Here, we show that L-ASP's glutaminase activity is not always required for the enzyme's anticancer effect. We first used molecular dynamics simulations of the clinically standard Escherichia coli L-ASP to predict what mutated forms could be engineered to retain activity against asparagine but not glutamine. Dynamic mapping of enzyme substrate contacts identified Q59 as a promising mutagenesis target for that purpose. Saturation mutagenesis followed by enzymatic screening identified Q59L as a variant that retains asparaginase activity but shows undetectable glutaminase activity. Unlike wild-type L-ASP, Q59L is inactive against cancer cells that express measurable asparagine synthetase (ASNS). Q59L is potently active, however, against ASNS-negative cells. Those observations indicate that the glutaminase activity of L-ASP is necessary for anticancer activity against ASNS-positive cell types but not ASNS-negative cell types. Because the clinical toxicity of L-ASP is thought to stem from its glutaminase activity, these findings suggest the hypothesis that glutaminase-negative variants of L-ASP would provide larger therapeutic indices than wild-type L-ASP for ASNS-negative cancers.

Can Dispersion Forces Govern Aromatic Stacking in an Organic Solvent?

Experimental support for the dominance of van der Waals dispersion forces in aromatic stacking interactions occurring in organic solution is surprisingly limited. The size-dependence of aromatic stacking in an organic solvent was examined. The interaction energy was found to vary by about 7.5 kJ mol(-1) on going from a phenyl-phenyl to an anthracene-pyrene stack. Strikingly, the experimental data were highly correlated with dispersion energies determined using symmetry-adapted perturbation theory (SAPT), while the induction, exchange, electrostatic, and solvation energy components correlated poorly. Both the experimental data and the SAPT-dispersion energies gave high-quality correlations with the change in solvent accessible area upon complexation. Thus, the size-dependence of aromatic stacking interactions is consistent with the dominance of van der Waals dispersion forces even in the presence of a competing polarizable solvent.

Real-space quadrature: a convenient, efficient representation for multipole expansions.

Multipoles are central to the theory and modeling of polarizable and nonpolarizable molecular electrostatics. This has made a representation in terms of point charges a highly sought after goal, since rotation of multipoles is a bottleneck in molecular dynamics implementations. All known point charge representations are orders of magnitude less efficient than spherical harmonics due to either using too many fixed charge locations or due to nonlinear fitting of fewer charge locations. We present the first complete solution to this problem-completely replacing spherical harmonic basis functions by a dramatically simpler set of weights associated to fixed, discrete points on a sphere. This representation is shown to be space optimal. It reduces the spherical harmonic decomposition of Poisson's operator to pairwise summations over the point set. As a corollary, we also shows exact quadrature-based formulas for contraction over trace-free supersymmetric 3D tensors. Moreover, multiplication of spherical harmonic basis functions translates to a direct product in this representation.

Spatiotemporal pH dynamics in concentration polarization near ion-selective membranes.

We present a detailed analysis of the transient pH dynamics for a weak, buffered electrolyte subject to voltage-driven transport through an ion-selective membrane. We show that pH fronts emanate from the concentration polarization zone next to the membrane and that these propagating fronts change the pH in the system several units from its equilibrium value. The analysis is based on a 1D model using the unsteady Poisson-Nernst-Planck equations with nonequilibrium chemistry and without assumptions of electroneutrality or asymptotically thin electric double layers. Nonequilibrium chemical effects, especially for water splitting, are shown to be important for the dynamical and spatiotemporal evolution of the pH fronts. Nonetheless, the model also shows that at steady state the assumption of chemical equilibrium can still lead to good approximations of the global pH distribution. Moreover, our model shows that the transport of the hydronium ion in the extended space charge region is governed by a balance between electromigration and water self-ionization. On the basis of this observation, we present a simple model showing that the net flux of the hydronium ion is proportional to the length of the extended space charge region and the water self-ionization rate. To demonstrate these effects in practice, we have adopted the experiment of Mai et al. (Mai, J.; Miller, H.; Hatch, A. V. Spatiotemporal Mapping of Concentration Polarization Induced pH Changes at Nanoconstrictions. ACS Nano 2012, 6, 10206) as a model problem, and by including the full chemistry and transport, we show that the present model can capture the experimentally observed pH fronts. Our model can, among other things, be used to predict and engineer pH dynamics, which can be essential to the performance of membrane-based systems for biochemical separation and analysis.

A combined resonance enhanced multiphoton ionization and ab initio study of the first absorption band of 1,2,4,5-tetrafluorobenzene, pentafluorobenzene, and hexafluorobenzene.

The resonance enhanced multiphoton ionization (REMPI) spectra of jet-cooled penta- and hexafluorobenzene when excited in the region λ(ex) = 265-253 nm of the first absorption band and observed only in the CF(+) mass channel is dominated by rotational structure in the A←X transition of CF. However, structure in the CF(+) channel for λ(ex) > 265 nm is not a continuation of this CF spectrum and is assigned to vibrational activity in two low-frequency modes of a distorted excited state of the parent molecule. The vibrational structure is assigned to the lowest ππ* state from a comparison with the equivalent spectrum of 1,2,4,5-tetrafluorobenzene. Ab initio calculations at the CIS level of theory of the ππ* state of the penta- and hexafluorobenzene reveal a much more distorted equilibrium geometry (C1 symmetry) than that of 1,2,4,5-tetrafluorobenzene. Long progressions observed in the λ(ex) > 265 nm REMPI spectra of C6HF5 and C6F6 are assigned to two very low frequency (~30 and 80 cm(-1)) modes. The role of the close-lying ππ* and πσ* states in determining the energy redistribution of the initially excited state by internal conversion is discussed. Both the fluorescent yield and the direct production of CF(X) are associated with transfer to the lower-lying πσ* state.

A facile one step route that introduces functionality to polymer powders for laser sintering.

Laser Sintering (LS) is a type of Additive Manufacturing (AM) exploiting laser processing of polymeric particles to produce 3D objects. Because of its ease of processability and thermo-physical properties, polyamide-12 (PA-12) represents ~95% of the polymeric materials used in LS. This constrains the functionality of the items produced, including limited available colours. Moreover, PA-12 objects tend to biofoul in wet environments. Therefore, a key challenge is to develop an inexpensive route to introduce desirable functionality to PA-12. We report a facile, clean, and scalable approach to modification of PA-12, exploiting supercritical carbon dioxide (scCO2) and free radical polymerizations to yield functionalised PA-12 materials. These can be easily printed using commercial apparatus. We demonstrate the potential by creating coloured PA-12 materials and show that the same approach can be utilized to create anti-biofouling objects. Our approach to functionalise materials could open significant new applications for AM.

Theoretical analysis of divalent cation effects on aptamer recognition of neurotransmitter targets.

Aptamer-based sensing of small molecules such as dopamine and serotonin in the brain, requires characterization of the specific aptamer sequences in solutions mimicking the in vivo environment with physiological ionic concentrations. In particular, divalent cations (Mg2+ and Ca2+) present in brain fluid, have been shown to affect the conformational dynamics of aptamers upon target recognition. Thus, for biosensors that transduce aptamer structure switching as the signal response, it is critical to interrogate the influence of divalent cations on each unique aptamer sequence. Herein, we demonstrate the potential of molecular dynamics (MD) simulations to predict the behaviour of dopamine and serotonin aptamers on sensor surfaces. The simulations enable molecular-level visualization of aptamer conformational changes that, in some cases, are significantly influenced by divalent cations. The correlations of theoretical simulations with experimental findings validate the potential for MD simulations to predict aptamer-specific behaviors on biosensors.