Inside the capsid: Revealing viral genome organization through multiscale simulations.

In a recent issue of Nature, Coshic et al. employ a computational multiscale approach to package the complete HK97 viral genome into its capsid. They find both good agreement with experimental observations and shed new light on the heterogeneity of genome structures and the mechanism by which they package.

Functional diversity among cardiolipin binding sites on the mitochondrial ADP/ATP carrier.

Lipid-protein interactions play a multitude of essential roles in membrane homeostasis. Mitochondrial membranes have a unique lipid-protein environment that ensures bioenergetic efficiency. Cardiolipin (CL), the signature mitochondrial lipid, plays multiple roles in promoting oxidative phosphorylation (OXPHOS). In the inner mitochondrial membrane, the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) exchanges ADP and ATP, enabling OXPHOS. AAC/ANT contains three tightly bound CLs, and these interactions are evolutionarily conserved. Here, we investigated the role of these buried CLs in AAC/ANT using a combination of biochemical approaches, native mass spectrometry, and molecular dynamics simulations. We introduced negatively charged mutations into each CL-binding site of yeast Aac2 and established experimentally that the mutations disrupted the CL interactions. While all mutations destabilized Aac2 tertiary structure, transport activity was impaired in a binding site-specific manner. Additionally, we determined that a disease-associated missense mutation in one CL-binding site in human ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

The α C helix is a central regulator of PKR activation.

Protein kinase R (PKR) functions in the eukaryotic innate immune system as a first-line defense against viral infections. PKR binds viral dsRNA, leading to autophosphorylation and activation. In its active state, PKR can phosphorylate its primary substrate, eIF2 α , which blocks initiation of translation in the infected cell. It has been established that PKR activation occurs when the kinase domain dimerizes in a back-to-back configuration. However, the mechanism by which dimerization leads to enzymatic activation is not fully understood. Here, we investigate the structural mechanistic basis and energy landscape for PKR activation, with a focus on the α C helix - a kinase activation and signal integration hub - using all-atom equilibrium and enhanced sampling molecular dynamics simulations. By employing window-exchange umbrella sampling, we compute free energy profiles of activation which show that back-to-back dimerization stabilizes a catalytically competent conformation of PKR. Key hydrophobic residues in the homodimer interface contribute to stabilization of the α C helix in an active conformation and the position of its glutamate residue. Using linear mutual information analysis, we analyze allosteric communication connecting the protomers' N-lobes and the α C helix dimer interface with the α C helix.

Bacterial diversity in agricultural drainage ditches shifts with increasing urea-N concentrations.

Urea-based fertilizers applied to crop fields can enter the surface waters of adjacent agricultural drainage ditches and contribute to the nitrogen (N) loading in nearby watersheds. Management practices applied in drainage ditches promote N removal by the bacterial communities, but little is known about the impacts of excess urea fertilizer from crop fields on the bacterial diversity in these ditches. In 2017, sediments from drainage ditches next to corn and soybean fields were sampled to determine if fertilizer application and high urea-N concentrations alters bacterial diversity and urease gene abundances. A mesocosm experiment was paired with a field study to determine which bacterial groups respond to high urea-N concentrations. The bacterial diversity in the ditch next to corn fields was significantly different from the other site. The bacterial orders of Rhizobiales, Bacteroidales, Acidobacteriales, Burkholderiales, and Anaerolineales were most abundant in the ditch next to corn and increased after the addition of urea-N (0.5 mg N L-1) during the mesocosm experiment. The results of our study suggests that urea-N concentrations >0.07 mg N L-1, which are higher than concentrations associated with downstream harmful algal blooms, can lead to shifts in the bacterial communities of agricultural drainage ditches.

Curvature sensing lipid dynamics in a mitochondrial inner membrane model.

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations.

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

High-Pressure Gravimetric Measurements for Binary Gas Adsorption Equilibria and Comparisons with Ideal Adsorbed Solution Theory (IAST).

Measurements of gas mixture adsorption equilibria at high pressures are important for assessing actual adsorbent selectivities but are often out of reach, given the challenging nature of the required experiments. Here, we report a high-pressure gravimetric binary gas adsorption equilibrium measurement system based on simultaneous gas density and mixture adsorption measurements in a single gas cell coupled to a magnetic-suspension balance. Compared to traditional techniques which rely on analytical measurements of gas composition, this approach does not require any sampling. Adsorption measurements of two gas mixtures (0.500 N2 + 0.500 CH4 and 0.400 N2 + 0.600 CO2, mole fraction) on a commercially available molecular sieve (NaY, sodium molecular sieve type Y) were carried out in the temperature range 282 to 325 K with a pressure up to 10 MPa. A prediction method for the gas mixture adsorption equilibria in a closed system using the ideal adsorbed solution theory (IAST) model was used to compare the experimental results. For binary mixtures of components with similar adsorption capacities (here N2 and CH4), the system can measure the adsorption equilibria at pressures higher than 1.0 MPa and the result agrees well with the IAST model prediction. For two gases with very different adsorption capacities, the uncertainty in the adsorption equilibrium measurement is much larger. The dominant uncertainty source is the gas density measurement, whose uncertainty could potentially be cut to half if the current titanium sinker is replaced with a sinker made of single-crystal silicon and with a larger volume.

A cage-specific hydrate equilibrium model for robust predictions of industrially-relevant mixtures.

Understanding the thermophysical properties and phase behaviour of gas hydrates is essential for industrial applications ranging from energy transport and storage, CO2 capture and sequestration, to gas production from hydrates found on the seabed. Current tools for predicting hydrate equilibrium boundaries typically use van der Waals-Platteeuw-type models which are over-parameterised containing terms with limited physical basis. Here we present a new model for hydrate equilibrium calculations with 40% fewer parameters than existing tools but with equivalent accuracy, including for multicomponent gas mixtures and/or systems with thermodynamic inhibitors. By eliminating multi-layered shells from the model's conceptual basis and focusing on Kihara potential parameters for guest-water interactions specific to each hydrate cavity type, this new model provides insight into the physical chemistry governing hydrate thermodynamics. The model retains the improved description of the empty lattice developed recently by Hielscher et al. but couples the hydrate model with a Cubic-Plus-Association Equation of State (CPA-EOS) to describe fluid mixtures with many more components including inhibitors such as methanol and mono-ethylene glycol used by industry. An extensive database of over 4000 data points was used to train and evaluate the new model and compare its performance against existing tools. The absolute average deviation in temperature (AADT) achieved with the new model is 0.92 K for multicomponent gas mixtures, compared with 1.00 K for the widely-known model of Ballard and Sloan, and 0.86 K for the CPA-hydrates model implemented in the MultiFlash 7.0 software package. With fewer, more physically justified parameters, this new cage-specific model provides a robust basis for improved hydrate equilibrium predictions particularly for industrially-important, multi-component mixtures containing thermodynamic inhibitors.

Natural sources and controlling factors of urea-nitrogen concentrations in agricultural drainage ditches.

Agricultural drainage ditches accumulate high urea-nitrogen (N) concentrations even in the absence of urea fertilizer applications to adjacent crop fields. The accumulated urea, and other bioavailable forms of dissolved organic nitrogen (DON), can be flushed downstream during substantial rainfall events altering downstream water quality and phytoplankton communities. Sources of urea-N supporting its accumulation in agricultural drainage ditches are poorly understood. A ditch flooding event was simulated using mesocosms with N treatment solutions and monitored for changes in N concentrations, physicochemical properties, dissolved organic matter (DOM) composition, and N cycling enzymes. N concentrations were also monitored in field ditches after two rainfall events. Urea-N concentrations were higher with DON enrichment, but the treatment effects were temporary. The DOM released from the mesocosm sediments was dominated by terrestrial-derived, high molecular weight material. The lack of microbial-derived DOM and evidence from the bacterial gene abundances in the mesocosms suggests that urea-N accumulation after rainfall may not be associated with fresh biological inputs. The urea-N concentrations after spring rainfall and flooding with DON substrates indicated the urea from fertilizers may only temporarily affect urea-N concentrations in drainage ditches. Because urea-N concentrations increased with a high degree of DOM humification, sources of urea may derive from the slow decomposition of complex DOM structures. This study provides further insights of sources contributing to high urea-N concentrations and the types of DOM released from drainage ditches to nearby surface waters after hydrological events.

Conserved cardiolipin-mitochondrial ADP/ATP carrier interactions assume distinct structural and functional roles that are clinically relevant.

The mitochondrial phospholipid cardiolipin (CL) promotes bioenergetics via oxidative phosphorylation (OXPHOS). Three tightly bound CLs are evolutionarily conserved in the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) which resides in the inner mitochondrial membrane and exchanges ADP and ATP to enable OXPHOS. Here, we investigated the role of these buried CLs in the carrier using yeast Aac2 as a model. We introduced negatively charged mutations into each CL-binding site of Aac2 to disrupt the CL interactions via electrostatic repulsion. While all mutations disturbing the CL-protein interaction destabilized Aac2 monomeric structure, transport activity was impaired in a pocket-specific manner. Finally, we determined that a disease-associated missense mutation in one CL-binding site in ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

Non-equilibrium virus particle dynamics: Microsecond MD simulations of the complete Flock House virus capsid under different conditions.

Flock House virus (FHV) is an animal virus and considered a model system for non-enveloped viruses. It has a small, icosahedral capsid (T=3) and a bipartite positive-sense RNA genome. We present an extensive study of the FHV capsid dynamics from all-atom molecular dynamics simulations of the complete capsid. The simulations explore different biologically relevant conditions (neutral/low pH, with/without RNA in the capsid) using the CHARMM force field. The results show that low pH destabilizes the capsid, causing radial expansion, and RNA stabilizes the capsid. The finding of low pH destabilization is biologically relevant because the capsid is exposed to low pH in the endosome, where conformational changes occur leading to genome release. We also observe structural changes at the fivefold and twofold symmetry axes that likely relate to the externalization of membrane active γ peptides through the fivefold vertex and extrusion of RNA at the twofold axis. Simulations using the Amber force field at neutral pH are also performed and display similar characteristics to the CHARMM simulations.

Aqueous Solid Formation Kinetics in High-Pressure Methane at Trace Water Concentrations.

Natural gas containing trace amounts of water is frequently liquefied at conditions where aqueous solids are thermodynamically stable. However, no data are available to describe the kinetics of aqueous solid formation at these conditions. Here, we present experimental measurements of both solid formation kinetics and solid-fluid equilibrium for trace concentrations of (12 ± 0.7) ppm water in methane using a stirred, high-pressure apparatus and visual microscopy. Along isochoric pathways with cooling rates around 1 K·min-1, micron-scale aqueous solids were observed to form at subcoolings of (0.3-8.6) K, relative to an average equilibrium melting temperature of (253 ± 1.9) K at (8.9 ± 0.08) MPa; these data are consistent with predicted methane hydrate dissociation conditions within the uncertainty of both the experiment and model. The 36 measured formation events were used to construct a cumulative formation probability distribution, which was then fitted with a model from Classical Nucleation Theory, enabling the extraction of kinetic and thermodynamic nucleation parameters. While the resulting nucleation parameter values were comparable to those published for methane hydrate formation in bulk-water systems, the observed growth kinetics were distinctly different with only a small percentage of the water in the system converting into micron-scale solids over the experimental time scale. These results may help explain how cryogenic heat exchangers in liquefied natural gas facilities can operate for long periods without blockages forming despite being at very high subcoolings for aqueous solids.

Variations in Bacterial Communities and Antibiotic Resistance Genes Across Diverse Recycled and Surface Water Irrigation Sources in the Mid-Atlantic and Southwest United States: A CONSERVE Two-Year Field Study.

Reduced availability of agricultural water has spurred increased interest in using recycled irrigation water for U.S. food crop production. However, there are significant knowledge gaps concerning the microbiological quality of these water sources. To address these gaps, we used 16S rRNA gene and metagenomic sequencing to characterize taxonomic and functional variations (e.g., antimicrobial resistance) in bacterial communities across diverse recycled and surface water irrigation sources. We collected 1 L water samples (n = 410) between 2016 and 2018 from the Mid-Atlantic (12 sites) and Southwest (10 sites) U.S. Samples were filtered, and DNA was extracted. The V3-V4 regions of the 16S rRNA gene were then PCR amplified and sequenced. Metagenomic sequencing was also performed to characterize antibiotic, metal, and biocide resistance genes. Bacterial alpha and beta diversities were significantly different (p < 0.001) across water types and seasons. Pathogenic bacteria, such as Salmonella enterica, Staphylococcus aureus, and Aeromonas hydrophilia were observed across sample types. The most common antibiotic resistance genes identified coded against macrolides/lincosamides/streptogramins, aminoglycosides, rifampin and elfamycins, and their read counts fluctuated across seasons. We also observed multi-metal and multi-biocide resistance across all water types. To our knowledge, this is the most comprehensive longitudinal study to date of U.S. recycled water and surface water used for irrigation. Our findings improve understanding of the potential differences in the risk of exposure to bacterial pathogens and antibiotic resistance genes originating from diverse irrigation water sources across seasons and U.S. regions.

Surface-Enhanced Raman Scattering Imaging of Cetylpyridinium Chloride Adsorption to a Solid Surface.

Surface active agents (surfactants) have found a variety of critical technological applications, from helping infant lungs breathe to fugitive dust control at industrial sites. Surfactant molecules adsorb to an interface and facilitate a decrease in the surface free energy (interfacial tension) between two immiscible phases. However, a limited number of methods (e.g., holography and fluorescence microscopy) achieved visualization of surfactant molecule distribution in multiphase systems qualitatively. To probe the efficacy and/or adsorption density of surfactants at such interfaces quantitatively, we demonstrate here a direct observation of surfactant adsorption by surface-enhanced Raman scattering (SERS). This work details the development of a research platform to study surfactant adsorption using Raman imaging. The imaging and analysis were successfully benchmarked against conventional interfacial tension measurements and thermodynamic theory employed to estimate surfactant adsorption at equilibrium. This in situ Raman-based experimental method provides a platform to interrogate structure-function relationships that inform the design process for new surfactant species.

Structure-activity relationships of mitochondria-targeted tetrapeptide pharmacological compounds.

Mitochondria play a central role in metabolic homeostasis, and dysfunction of this organelle underpins the etiology of many heritable and aging-related diseases. Tetrapeptides with alternating cationic and aromatic residues such as SS-31 (elamipretide) show promise as therapeutic compounds for mitochondrial disorders. In this study, we conducted a quantitative structure-activity analysis of three alternative tetrapeptide analogs, benchmarked against SS-31, that differ with respect to aromatic side chain composition and sequence register. We present the first structural models for this class of compounds, obtained with Nuclear Magnetic Resonance (NMR) and molecular dynamics approaches, showing that all analogs except for SS-31 form compact reverse turn conformations in the membrane-bound state. All peptide analogs bound cardiolipin-containing membranes, yet they had significant differences in equilibrium binding behavior and membrane interactions. Notably, analogs had markedly different effects on membrane surface charge, supporting a mechanism in which modulation of membrane electrostatics is a key feature of their mechanism of action. The peptides had no strict requirement for side chain composition or sequence register to permeate cells and target mitochondria in mammalian cell culture assays. All four peptides were pharmacologically active in serum withdrawal cell stress models yet showed significant differences in their abilities to restore mitochondrial membrane potential, preserve ATP content, and promote cell survival. Within our peptide set, the analog containing tryptophan side chains, SPN10, had the strongest impact on most membrane properties and showed greatest efficacy in cell culture studies. Taken together, these results show that side chain composition and register influence the activity of these mitochondria-targeted peptides, helping provide a framework for the rational design of next-generation therapeutics with enhanced potency.

Low-Field NMR Relaxation Analysis of High-Pressure Ethane Adsorption in Mesoporous Silicas.

Understanding the behaviour of short-chain hydrocarbons confined to porous solids informs the targeted extraction of natural resources from geological features, and underpins rational developments in separation, storage and catalytic conversion processes. Herein, we report the application of low-field (12.7 MHz) 1 H nuclear magnetic resonance (NMR) relaxation measurements to characterise ethane dynamics within mesoporous silica materials exhibiting mean pore diameters between 6 and 50 nm. Our measurements provide NMR-based adsorption isotherms within the range 25-50 bar and at ambient temperature, incorporating the ethane condensation point (40.7 bar at our experimental temperature of 23.6 °C). The quantitative nature of the acquired data is validated via a direct comparison of NMR-derived excess adsorption capacities with ex situ gravimetric ethane adsorption measurements, which are demonstrated to agree to within 0.2 mmol g-1 of the observed ethane capacity. NMR T2 relaxation time distributions are further demonstrated as a means to decouple interparticle and mesopore dominated adsorption phenomena, with unexpectedly rapid relaxation rates associated with interparticle ethane gas confirmed via a direct comparison with NMR self-diffusion analysis.

A favorable path to domain separation in the orange carotenoid protein.

The orange carotenoid protein (OCP) is responsible for nonphotochemical quenching (NPQ) in cyanobacteria, a defense mechanism against potentially damaging effects of excess light conditions. This soluble two-domain protein undergoes profound conformational changes upon photoactivation, involving translocation of the ketocarotenoid inside the cavity followed by domain separation. Domain separation is a critical step in the photocycle of OCP because it exposes the N-terminal domain (NTD) to perform quenching of the phycobilisomes. Many details regarding the mechanism and energetics of OCP domain separation remain unknown. In this work, we apply metadynamics to elucidate the protein rearrangements that lead to the active, domain-separated, form of OCP. We find that translocation of the ketocarotenoid canthaxanthin has a profound effect on the energetic landscape and that domain separation only becomes favorable following translocation. We further explore, characterize, and validate the free energy surface (FES) using equilibrium simulations initiated from different states on the FES. Through pathway optimization methods, we characterize the most probable path to domain separation and reveal the barriers along that pathway. We find that the free energy barriers are relatively small (<5 kcal/mol), but the overall estimated kinetic rate is consistent with experimental measurements (>1 ms). Overall, our results provide detailed information on the requirement for canthaxanthin translocation to precede domain separation and an energetically feasible pathway to dissociation.

Application of Raman Spectroscopy for Sorption Analysis of Functionalized Porous Materials.

Functionalized porous materials could play a key role in improving the efficiency of gas separation processes as required by applications such as carbon capture and storage (CCS) and across the hydrogen value chain. Due to the large number of different functionalizations, new experimental approaches are needed to determine if an adsorbent is suitable for a specific separation task. Here, it is shown for the first time that Raman spectroscopy is an efficient tool to characterize the adsorption capacity and selectivity of translucent functionalized porous materials at high pressures, whereby translucence is the precondition to study mass transport inside of a material. As a proof of function, the performance of three silica ionogels to separate an equimolar (hydrogen + carbon dioxide) gas mixture is determined by both accurate gravimetric sorption measurements and Raman spectroscopy, with the observed consistency establishing the latter as a novel measurement technique for the determination of adsorption capacity. These results encourage the use of the spectroscopic approach as a rapid screening method for translucent porous materials, particularly since only very small amounts of sample are required.

Low-Field Functional Group Resolved Nuclear Spin Relaxation in Mesoporous Silica.

Solid-fluid interactions underpin the efficacy of functional porous materials across a diverse array of chemical reaction and separation processes. However, detailed characterization of interfacial phenomena within such systems is hampered by their optically opaque nature. Motivated by the need to bridge this capability gap, we report low-magnetic-field two-dimensional (2D) 1H nuclear spin relaxation measurements as a noninvasive probe of adsorbate identity and interfacial dynamics, exploring the relaxation characteristics exhibited by liquid hydrocarbon adsorbates confined to a model mesoporous silica. For the first time, we demonstrate the capacity of this approach in distinguishing functional group-specific relaxation phenomena across a diverse range of alcohols and carboxylic acids employed as solvents, reagents, and liquid hydrogen carriers, with distinct relaxation responses assigned to the alkyl and hydroxyl moieties of each confined liquid. Uniquely, this relaxation behavior is shown to correlate with adsorbate acidity, with the observed relationship rationalized on the basis of surface-adsorbate proton-exchange dynamics. Our results demonstrate that nuclear spin relaxation provides a molecular-level perspective on sorbent/sorbate interactions, motivating the exploration of such measurements as a unique probe of adsorbate identity within optically opaque porous media.