Ion Transport, Selectivity, and Electronic Polarization in Fluoride Channels

Fluoride channels (Fluc) export toxic F- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the transport mechanism at the molecular level has remained elusive. Herein we have applied constant-pH molecular dynamics and free energy sampling methods to investigate fluoride transfer through a Fluc protein from Escherichia coli. We find that fluoride is facile to transfer in its charged form, i.e., F-, by traversing through a non-bonded network. The extraordinary F- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F- transfer rate calculated using an electronically polarizable force field is significantly more accurate compared to the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F- - Fluc interactions.

Key Factors Governing Initial Stages of Lipid Droplet Formation

Lipid droplets (LDs) are neutral lipid storage organelles surrounded by a phospholipid (PL)monolayer. LD biogenesis from the endoplasmic reticulum (ER) is driven by phase separation of neutral lipids, overcoming surface tension and membrane deformation. However, the core biophysics of the initial steps of LD formation remain relatively poorly understood. Here, we use a tunable, phenomenological coarse-grained (CG) model to study triacylglycerol (TG) nucleation in a bilayer membrane. We show that PL rigidity has a strong influence on TG lensing and membrane remodeling: When membrane rigidity increases, TG clusters remain more planar with high anisotropy but a minor degree of phase nucleation. This finding is confirmed by free energy sampling simulations that calculate the potential of mean force (PMF) as a function of the degree of nucleation and anisotropy. We also show that asymmetric tension, controlled by the number of PLs on each membrane leaflet, determines the budding direction. A TG lens buds in the direction of the monolayer containing excess PLs to allow for better PL coverage of TG, consistent with reported experiments. Finally, two governing mechanisms of the LD growth, Ostwald ripening and merging, are observed. Taken together, this study characterizes the interplay between two thermodynamic quantities during the initial LD phases, the TG bulk free energy and membrane remodeling energy.

Identifying the key steps determining the selectivity of toluene methylation with methanol over HZSM-5

Methylation of toluene with methanol to produce p-xylene has been investigated for decades, but the origin of selectivity is still under debate. Here we report computational studies based on ab initio molecular dynamics simulations and free energy sampling methods to identify the key steps determining the selectivity. The steps of toluene methylation to protonated-xylene, deprotonation of protonated-xylenes, and diffusion of xylene in HZSM-5 channels are compared. We find the pathways of formation for protonated p-/m-xylenes have similar free energy barriers. Meanwhile, the methylation is found rate-determining, thus the probability to generate p-/m-xylenes at the active site are similar. We then find that the diffusion for m-xylene along the zigzag channel is more difficult than its isomerization to p-xylene, which in turn further promotes the selectivity of p-xylene formation. These insights obtained at the molecular level are crucial for further development of high-performance zeolite catalysts for toluene methylation.

Key computational findings reveal proton transfer as driving the functional cycle in the phosphate transporter PiPT

Phosphate is an indispensable metabolite in a wide variety of cells and is involved in nucleotide and lipid synthesis, signaling, and chemical energy storage. Proton-coupled phosphate transporters within the major facilitator family are crucial for phosphate uptake in plants and fungi. Similar proton-coupled phosphate transporters have been found in different protozoan parasites that cause human diseases, in breast cancer cells with elevated phosphate demand, in osteoclast-like cells during bone reabsorption, and in human intestinal Caco2BBE cells for phosphate homeostasis. However, the mechanism of proton-driven phosphate transport remains unclear. Here, we demonstrate in a eukaryotic, high-affinity phosphate transporter from Piriformospora indica (PiPT) that deprotonation of aspartate 324 (D324) triggers phosphate release. Quantum mechanics/molecular mechanics molecular dynamics simulations combined with free energy sampling have been employed here to identify the proton transport pathways from D324 upon the transition from the occluded structure to the inward open structure and phosphate release. The computational insights so gained are then corroborated by studies of D45N and D45E amino acid substitutions via mutagenesis experiments. Our findings confirm the function of the structurally predicted cytosolic proton exit tunnel and suggest insights into the role of the titratable phosphate substrate.

Physical Characterization of Triolein and Implications for Its Role in Lipid Droplet Biogenesis

Lipid droplets (LDs) are neutral lipid storing organelles surrounded by a phospholipid (PL) monolayer. At present, how LDs are formed in the endoplasmic reticulum (ER) bilayer is poorly understood. In this study, we present a revised triolein (TG) model, the main constituent of the LD core, and characterize its properties in a bilayer membrane to demonstrate the implications of its behavior in LD biogenesis. In all-atom (AA) bilayer simulations, TG resides at the surface, adopting PL-like conformations (denoted in this work as SURF-TG). Free energy sampling simulation results estimate the barrier for TG relocating from the bilayer surface to the bilayer center to be ~2 kcal/mol in the absence of an oil lens. Conical SURF-TG is able to modulate membrane properties by increasing PL ordering, decreasing bending modulus, and creating local negative curvature. The other conical lipid, dioleoyl-glycerol (DAG), also reduces the membrane bending modulus and populates the negative curvature regions. A phenomenological coarse-grained (CG) model is also developed to observe larger scale SURF-TG-mediated membrane deformation. The CG simulations confirm that TG nucleates between the bilayer leaflets at a critical concentration when SURF-TG is evenly distributed. However, when one monolayer contains more SURF-TG, the membrane bends toward the other leaflet. The central conclusion of this study is that SURF-TG is a negative curvature inducer, as well as a membrane modulator. To this end, a model has proposed in which the accumulation of SURF-TG in the luminal leaflet bends the ER bilayer toward the cytosolic side, followed by TG nucleation.

Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale

ABSTRACTThe sarcoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca2+, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca2+. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca2+-binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics (MS-RMD) method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca2+ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.SIGNIFICANCEMultiscale reactive molecular dynamics combined with free energy sampling was applied to study proton transport through a transient water pore connecting the Ca2+-binding site to the lumen in SERCA. This is the first computational study of this large biomolecular system that treats the hydrated excess proton and its transport through water structures and amino acids explicitly. When also correctly accounting for the hydration fluctuations of the pore, it is found that a transiently hydrated channel can transport protons on a microsecond timescale. These results quantitatively support the hypothesis of the proton intake into the sarcoplasm via SERCA, in addition to the well-known proton pumping by SERCA to the cytoplasm along with Ca2+ transport.

Quantum Accurate Prediction of Plutonium–Plutonium Dihydride Phase Equilibrium Using a Spin-Lattice Model

<div> <div> <div> <p>Plutonium-based materials are vital for use as nuclear fuels and as portable power sources for space vehicles. However, elucidating their sensitivity to hydriding corrosion represents an extreme challenge due to the toxicity of Pu as well as its anomalous magnetic properties. In this work, we develop a spin-lattice model of plutonium–plutonium dihydride (Pu–PuH2) phase equilibrium that retains the accuracy of density functional theory (DFT) while yielding many orders of magnitude improvement in computational efficiency. Using Monte Carlo and free energy sampling algorithms, we compute a number of Pu–PuH2 equilibrium properties that are difficult to probe experimentally, including equilibrium pressures and phase compositions at room temperature and the PuH2 heat of formation. Our method will have particular impact on these types of materials studies, where there is a strong need for computationally efficient approaches to bridge time and length scale gaps between quantum calculations and experiments. </p> </div> </div> </div>