Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

EmrE is a small, homodimeric membrane transporter that exploits the established electrochemical proton gradient across the Escherichia coli inner membrane to export toxic polyaromatic cations, prototypical of the wider small-multidrug resistance transporter family. While prior studies have established many fundamental aspects of the specificity and rate of substrate transport in EmrE, low resolution of available structures has hampered identification of the transport coupling mechanism. Here we present a complete, refined atomic structure of EmrE optimized against available cryo-electron microscopy (cryo-EM) data to delineate the critical interactions by which EmrE regulates its conformation during the transport process. With the model, we conduct molecular dynamics simulations of the transporter in explicit membranes to probe EmrE dynamics under different substrate loading and conformational states, representing different intermediates in the transport cycle. The refined model is stable under extended simulation. The water dynamics in simulation indicate that the hydrogen-bonding networks around a pair of solvent-exposed glutamate residues (E14) depend on the loading state of EmrE. One specific hydrogen bond from a tyrosine (Y60) on one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. The hydrogen bond provided by Y60 lowers the pKa of one glutamate relative to the other, suggesting both glutamates should be protonated for the hydrogen bond to break and a substrate-free transition to take place. These findings establish the molecular mechanism for the coupling between proton transfer reactions and protein conformation in this proton-coupled secondary transporter.

Abstract not provided.

Abstract not provided.

Electrochemical double-layer capacitances of charged carbon nanotube (CNT) forests with tetraethyl ammonium tetrafluoro borate electrolyte in propylene carbonate are studied on the basis of molecular dynamics simulation. Direct molecular simulation of the filling of pore spaces of the forest is feasible even with realistic, small CNT spacings. The numerical solution of the Poisson equation based on the extracted average charge densities then yields a regular experimental dependence on the width of the pore spaces, in contrast to the anomalous pattern observed in experiments on other carbon materials and also in simulations on planar slot-like pores. The capacitances obtained have realistic magnitudes but are insensitive to electric potential differences between the electrodes in this model. This agrees with previous calculations on CNT forest supercapacitors, but not with experiments which have suggested electrochemical doping for these systems. Those phenomena remain for further theory/modeling work.

The role that van der Waals (vdW) attractive forces play in the hydration and association of atomic hydrophobic solutes such as argon (Ar) in water is reanalyzed using the local molecular field (LMF) theory of those interactions. In this problem, solute vdW attractive forces can reduce or mask hydrophobic interactions as measured by contact peak heights of the ArAr correlation function compared to reference results for purely repulsive core solutes. Nevertheless, both systems exhibit a characteristic hydrophobic inverse temperature behavior in which hydrophobic association becomes stronger with increasing temperature through a moderate temperature range. The new theoretical approximation obtained here is remarkably simple and faithful to the statistical mechanical LMF assessment of the necessary force balance. Our results extend and significantly revise approximations made in a recent application of the LMF approach to this problem and, unexpectedly, support a theory of nearly 40 years ago.

Laying a basis for molecularly specific theory for the mobilities of ions in solutions of practical interest, we report a broad survey of velocity autocorrelation functions (VACFs) of Li+ and PF6- ions in water, ethylene carbonate, propylene carbonate, and acetonitrile solutions. We extract the memory function, γ(t), which characterizes the random forces governing the mobilities of ions. We provide comparisons controlling for the effects of electrolyte concentration and ion-pairing, van der Waals attractive interactions, and solvent molecular characteristics. For the heavier ion (PF6-), velocity relaxations are all similar: negative tail relaxations for the VACF and a clear second relaxation for γt, observed previously also for other molecular ions and with n-pentanol as the solvent. For the light Li+ ion, short time-scale oscillatory behavior masks simple, longer time-scale relaxation of γt. But the corresponding analysis of the solventberg Li+H2O4 does conform to the standard picture set by all the PF6- results.

Abstract not provided.

Progress in understanding liquid ethylene carbonate (EC) and propylene carbonate (PC) on the basis of molecular simulation, emphasizing simple models of interatomic forces, is reviewed. Results on the bulk liquids are examined from the perspective of anticipated applications to materials for electrical energy storage devices. Preliminary results on electrochemical double-layer capacitors based on carbon nanotube forests and on model solid-electrolyte interphase (SEI) layers of lithium ion batteries are considered as examples. The basic results discussed suggest that an empirically parameterized, non-polarizable force field can reproduce experimental structural, thermodynamic, and dielectric properties of EC and PC liquids with acceptable accuracy. More sophisticated force fields might include molecular polarizability and Buckingham-model description of inter-atomic overlap repulsions as extensions to Lennard-Jones models of van der Waals interactions. Simple approaches should be similarly successful also for applications to organic molecular ions in EC/PC solutions, but the important case of Li+ deserves special attention because of the particularly strong interactions of that small ion with neighboring solvent molecules. To treat the Li+ ions in liquid EC/PC solutions, we identify interaction models defined by empirically scaled partial charges for ion-solvent interactions. The empirical adjustments use more basic inputs, electronic structure calculations and ab initio molecular dynamics simulations, and also experimental results on Li+ thermodynamics and transport in EC/PC solutions. Application of such models to the mechanism of Li+ transport in glassy SEI models emphasizes the advantage of long time-scale molecular dynamics studies of these non-equilibrium materials.

Abstract not provided.

Here, we study quasi-chemical theory (QCT) for the free energies of divalent alkaline earth ions (Ba 2+, Sr 2+, Ca 2+, Mg 2+) in water, emphasizing that: (a) interactions between metal ions and proximal water molecules are as strong as traditional chemical effects; (b) QCT builds directly from accessible electronic structure calculations but rests on fully elaborated molecular statistical thermodynamics; (c) QCT offers choices of convenience in identifying coordination numbers for analysis. We investigate utilisation of direct QCT with inner-shell conditioning (Formula presented.), alternative to the traditional nλ=0 conditioning motivated by a generalised van der Waals view. The alternative (Formula presented.) works well: deleterious non-Gaussian effects of van der Waals repulsive interactions are not serious, and the alternative conditioning improves the convenience of QCT calculations. Comparison between ab initio and force field molecular dynamics (AIMD and FFMD) with standard models suggests that FFMD likely exaggerates the anharmonicity in the thermal motion of inner-shell ion-water clusters. Together with the general encouraging support for the harmonic approximations implied by the (Formula presented.) conditioning, that observation helps explain the remarkable success of the cluster-based QCT solution free energies, which do not require assessment of all inner-shell occupancies by simulation.

Teixobactin (Txb) is a recently discovered antibiotic against Gram-positive bacteria that induces no detectable resistance. The bactericidal mechanism is believed to be the inhibition of cell wall biosynthesis by Txb binding to lipid II and lipid III. Txb binding specificity likely arises from targeting of the shared lipid component, the pyrophosphate moiety. Despite synthesis and functional assessment of numerous chemical analogs of Txb, and consequent identification of the Txb pharmacophore, the detailed structural information of Txb-substrate binding is still lacking. Here, we use molecular modeling and microsecond-scale molecular dynamics simulations to capture the formation of Txb-lipid II complexes at a membrane surface. Two dominant binding conformations were observed, both showing characteristic lipid II phosphate binding by the Txb backbone amides near the C-terminal cyclodepsipeptide (d-Thr8-Ile11) ring. Additionally, binding by Txb also involved the side chain hydroxyl group of Ser7, as well as a secondary phosphate binding provided by the side chain of l-allo-enduracididine. Interestingly, those conformations differ by swapping two groups of hydrogen bond donors that coordinate the two phosphate moieties of lipid II, resulting in opposite orientations of lipid II binding. In addition, residues d-allo-Ile5 and Ile6 serve as the membrane anchors in both Txb conformations, regardless of the detailed phosphate binding interactions near the cyclodepsipeptide ring. The role of hydrophobic residues in Txb activity is primarily for its membrane insertion, and subsidiarily to provide non-polar interactions with the lipid II tail. Based on the Txb-lipid II interactions captured in their complexes, as well as their partitioning depths into the membrane, we propose that the bactericidal mechanism of Txb is to arrest cell wall synthesis by selectively inhibiting the transglycosylation of peptidoglycan, while possibly leaving the transpeptidation step unaffected. The observed "pyrophosphate caging" mechanism of lipid II inhibition appears to be similar to some lantibiotics, but different from that of vancomycin or bacitracin.

This work demonstrates that the ionic selectivity and ionic conductivity of nanoporous membranes can be controlled independently via layer-by-layer (LbL) deposition of polyelectrolytes and subsequent selective cross-linking of these polymer layers. LbL deposition offers a scalable, inexpensive method to tune the ion transport properties of nanoporous membranes by sequentially dip coating layers of cationic polyethyleneimine and anionic poly(acrylic acid) onto polycarbonate membranes. The cationic and anionic polymers are self-assembled through electrostatic and hydrogen bonding interactions and are chemically crosslinked to both change the charge distribution and improve the intermolecular integrity of the deposited films. Both the thickness of the deposited coating and the use of chemical cross-linking agents influence charge transport properties significantly. Increased polyelectrolyte thickness increases the selectivity for cationic transport through the membranes while adding polyelectrolyte films decreases the ionic conductivity compared to an uncoated membrane. Once the nanopores are filled, no additional decrease in conductivity is observed with increasing film thickness and, upon cross-linking, a portion of the lost conductivity is recovered. The cross-linking agent also influences the ionic selectivity of the resulting polyelectrolyte membranes. Increased selectivity for cationic transport occurs when using glutaraldehyde as the cross-linking agent, as expected due to the selective cross-linking of primary amines that decreases the net positive charge. Together, these results inform deposition of chemically robust, highly conductive, ion-selective membranes onto inexpensive porous supports for applications ranging from energy storage to water purification.

Abstract not provided.

We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Enzymes that degrade specific small molecules could save lives by neutralizing threats from chemical agents in the blood or environment, or by starving pathogenic cells, but promiscuous interactions with other molecules typically limit their effectiveness by blocking the enzyme active site. An obvious solution would be to re-engineer the enzyme to enhance catalytic fidelity, but lack of understanding about how enzymes discriminate between molecules remains a formidable challenge to this approach. Our recent work in collaboration with the University of Texas (UT) suggested a new approach and a model system for understanding enzyme specificity. Asparaginase enzymes catalyze degradation of asparagine, which forms the basis of a medical treatment. Com- petition by the abundant and chemically similar molecule, glutamine, interferes with asparagine decomposition, thus hindering enzyme efficacy. Asparaginase is advantageous as a model degra- dation enzyme because variants that demonstrate different binding affinities and catalytic rates can be compared. Here, we leveraged Sandia and the University of Maryland's strengths in molecu- lar simulation, and UT experimental expertise in asparaginase modification and functional assays, to understand asparaginase specificity. Our results advanced a new hypothesis about asparagi- nase catalytic mechanism that explains for the first time why proximity between the substrate's alpha-carboxyl and carboxamide is absolutely required for catalysis. Based on those insights, we developed the first mutant (Q59L) asparaginase from E. coli that lacks activity toward glutamine. We used that mutant to show that glutaminase activity is required to kill cancer cells that have asparagine synthetase enzymes (ASNS), but not ASNS-negative cancer cells.