Assessing the Performance of the Nonbonded Mg2+ Models in a Two-Metal-Dependent Ribonuclease

SLC24A-mediated calcium exchange as an indispensable component of the diatom cell density-driven signaling pathway

Diatom bloom is characterized by a rapid increase of population density. Perception of population density and physiological responses can significantly influence their survival strategies, subsequently impacting bloom fate. The population density itself can serve as a signal, which is perceived through chemical signals or chlorophyll fluorescence signals triggered by high cell density, and their intracellular signaling mechanisms remain to be elucidated. In this study, we focused on the model diatom, Phaeodactylum tricornutum, and designed an orthogonal experiment involving varying cell densities and light conditions, to stimulate the release of chemical signals and light-induced chlorophyll fluorescence signals. Utilizing RNA-Seq and Weighted Gene Co-expression Network Analysis, we identified four gene clusters displaying density-dependent expression patterns. Within these, a potential hub gene, PtSLC24A, encoding a Na+/Ca2+ exchanger, was identified. Based on molecular genetics, cellular physiology, computational structural biology, and in situ oceanic data, we propose a potential intracellular signaling mechanism related to cell density in marine diatoms using Ca2+: upon sensing population density signals mediated by chemical cues, the membrane-bound PtSLC24A facilitates the efflux of Ca2+ to maintain specific intracellular calcium levels, allowing the transduction of intracellular density signals, subsequently regulating physiological responses, including cell apoptosis, ultimately affecting algal blooms fate. These findings shed light on the calcium-mediated intracellular signaling mechanism of marine diatoms to changing population densities, and enhances our understanding of diatom bloom dynamics and their ecological implications.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Molecular Mechanism of D1135E-Induced Discriminated CRISPR-Cas9 PAM Recognition

The off-target effects of Streptococcus pyogenes Cas9 (SpCas9) pose a significant challenge to harness it as a therapeutical approach. Two major factors can result in SpCas9 off-targeting: tolerance to target DNA-guide RNA (gRNA) mismatch and less stringent recognition of protospacer adjacent motif (PAM) flanking the target DNA. Despite the abundance of engineered SpCas9-gRNA variants with improved sensitivity to target DNA-gRNA mismatch, studies focusing on enhancing SpCas9 PAM recognition stringency are quite few. A recent pioneering study identified a D1135E variant of SpCas9 that exhibits much-reduced editing activity at the noncanonical NAG/NGA PAM sites while preserving robust on-target activity at the canonical NGG-flanking sites (N is any nucleobase). Herein, we aim to clarify the molecular mechanism by which this single D1135E mutation confers on SpCas9 enhanced specificity for PAM recognition by molecular dynamics simulations. The results suggest that the variant maintains the base-specific recognition for the canonical NGG PAM via four hydrogen bonds, akin to that in the wild type (WT) SpCas9. While the noncanonical NAG PAM is engaged to the two PAM-interacting arginine residues (i.e., R1333 and R1335) in WT SpCas9 via two to three hydrogen bonds, the D1135E variant prefers to establish two hydrogen bonds with the PAM bases, accounting for its minimal editing activity on the off-target sites with an NAG PAM. The impaired NAG recognition by D1135E SpCas9 results from the PAM duplex displacement such that the hydrogen bond of R1333 to the second PAM base is disfavored. We further propose a mechanistic model to delineate how the mutation perturbs the noncanonical PAM recognition. We anticipate that the mechanistic knowledge could be leveraged for continuous optimization of SpCas9 PAM recognition specificity toward high-precision demanding applications.

Rational Engineering of CRISPR-Cas9 Nuclease to Attenuate Position-Dependent Off-Target Effects

The RNA-guided Cas9 nuclease from Streptococcus pyogenes has become an important gene-editing tool. However, its intrinsic off-target activity is a major challenge for biomedical applications. Distinct from some reported engineering strategies that specifically target a single domain, we rationally introduced multiple amino acid substitutions across multiple domains in the enzyme to create potential high-fidelity variants, considering the Cas9 specificity is synergistically determined by various domains. We also exploited our previously derived atomic model of activated Cas9 complex structure for guiding new modifications. This approach has led to the identification of the HSC1.2 Cas9 variant with enhanced specificity for DNA cleavage. While the enhanced specificity associated with the HSC1.2 variant appeared to be position-dependent in the in vitro cleavage assays, the frequency of off-target DNA editing with this Cas9 variant is much less than that of the wild-type Cas9 in human cells. The potential mechanisms causing the observed position-dependent effect were investigated through molecular dynamics simulation. Our discoveries establish a solid foundation for leveraging structural and dynamic information to develop Cas9-like enzymes with high specificity in gene editing.

A multiscale approach to predict the binding mode of metallo beta-lactamase inhibitors

Antibiotic resistance is a major threat to global public health. β-lactamases, which catalyze breakdown of β-lactam antibiotics, are a principal cause. Metallo β-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all β-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems.

Open simultaneous saccharification and fermentation of l-lactic acid by complete utilization of sweet sorghum stalk: a water-saving process

A complete and efficient utilization of sweet sorghum stalk including sweet sorghum juice (SSJ) and sweet sorghum bagasse (SSB) was achieved via the open simultaneous saccharification and fermentation (SSF) of l-lactic acid. To simplify the pretreatment process and reduce water consumption, a combined hydrolysis approach was applied and the NaOH-pretreated liquor (SL) was utilized as a partial neutralizing agent. In order to further enhance the product titer, the acid hydrolysate of SSJ (SSJAH) was fed, and MgO was used as a neutralizing agent. A product titer of 94 g L-1 was obtained with a productivity of 1.55 g L-1 h-1, and the yield reached 98.31%. Totally, 274.79 g l-lactic acid was produced from 1 kg sweet sorghum stalk, and 83.22% water was saved compared with the previous study based on alkali pretreatment of SSB. This study provides an effective process for l-lactic acid biosynthesis from lignocellulosic substrates.

Comparing Nonbonded Metal Ion Models in the Divalent Cation Binding Protein PsaA

Divalent metal cations are essential for many biological processes; however, accurately modeling divalent metal ions has proved a significant challenge for molecular dynamics force fields. Here we show that the choice of ion model influences the observed dynamics in PsaA, a metal binding protein from Streptococcus pneumoniae. We conduct extensive unbiased simulations and free energy calculations of PsaA bound to its cognate ligand Mn²⁺ and inhibitory ligand Zn²⁺ using three nonbonded ion models: a 12-6 model, a 12-6-4 model, and a multisite model. The observed coordination geometries and metal binding dynamics are sensitive to the choice of ion model, with the most dramatic differences observed in free energy calculations of ion release. We show that the conformational ensemble of Mn-bound PsaA is more similar to the crystallographic metal bound open state. This work extends the current model of PsaA metal binding and provides a framework for the rationalization of experimentally determined metal binding behavior. Our findings support the use of the 12-6-4 ion model for further simulations of divalent cation binding proteins.

Structural and functional insights into the bona fide catalytic state of Streptococcus pyogenes Cas9 HNH nuclease domain

The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (SpyCas9), along with a programmable single-guide RNA (sgRNA), has been exploited as a significant genome-editing tool. Despite the recent advances in determining the SpyCas9 structures and DNA cleavage mechanism, the cleavage-competent conformation of the catalytic HNH nuclease domain of SpyCas9 remains largely elusive and debatable. By integrating computational and experimental approaches, we unveiled and validated the activated Cas9-sgRNA-DNA ternary complex in which the HNH domain is neatly poised for cleaving the target DNA strand. In this catalysis model, the HNH employs the catalytic triad of D839-H840-N863 for cleavage catalysis, rather than previously implicated D839-H840-D861, D837-D839-H840, or D839-H840-D861-N863. Our study contributes critical information to defining the catalytic conformation of the HNH domain and advances the knowledge about the conformational activation underlying Cas9-mediated DNA cleavage.