Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Modelling ligand exchange in metal complexes with machine learning potentials

Metal ions are irreplaceable in many areas of chemistry, including (bio)catalysis, self-assembly and charge transfer processes. Yet, modelling their structural and dynamic properties in diverse chemical environments remains challenging for both force fields and ab initio methods. Here, we introduce a strategy to train machine learning potentials (MLPs) using MACE, an equivariant message-passing neural network, for metal-ligand complexes in explicit solvents. We explore the structure and ligand exchange dynamics of Mg2+ in water and Pd2+ in acetonitrile as two illustrative model systems. The trained potentials accurately reproduce equilibrium structures of the complexes in solution, including different coordination numbers and geometries. Furthermore, the MLPs can model structural changes between metal ions and ligands in the first coordination shell, and reproduce the free energy barriers for the corresponding ligand exchange. The strategy presented here provides a computationally efficient approach to model metal ions in solution, paving the way for modelling larger and more diverse metal complexes relevant to biomolecules and supramolecular assemblies.

Multimodal binding and inhibition of bacterial ribosomes by the antimicrobial peptides Api137 and Api88

Proline-rich antimicrobial peptides (PrAMPs) inhibit bacterial protein biosynthesis by binding to the polypeptide exit tunnel (PET) near the peptidyl transferase center. Api137, an optimized derivative of honeybee PrAMP apidaecin, inhibits protein expression by trapping release factors (RFs), which interact with stop codons on ribosomes to terminate translation. This study uses cryo-EM, functional assays and molecular dynamic (MD) simulations to show that Api137 additionally occupies a second binding site near the exit of the PET and can repress translation independently of RF-trapping. Api88, a C-terminally amidated (-CONH2) analog of Api137 (-COOH), binds to the same sites, occupies a third binding pocket and interferes with the translation process presumably without RF-trapping. In conclusion, apidaecin-derived PrAMPs inhibit bacterial ribosomes by multimodal mechanisms caused by minor structural changes and thus represent a promising pool for drug development efforts.

RAPP-containing arrest peptides induce translational stalling by short circuiting the ribosomal peptidyltransferase activity

Arrest peptides containing RAPP (ArgAlaProPro) motifs have been discovered in both Gram-positive and Gram-negative bacteria, where they are thought to regulate expression of important protein localization machinery components. Here we determine cryo-EM structures of ribosomes stalled on RAPP arrest motifs in both Bacillus subtilis and Escherichia coli. Together with molecular dynamics simulations, our structures reveal that the RAPP motifs allow full accommodation of the A-site tRNA, but prevent the subsequent peptide bond from forming. Our data support a model where the RAP in the P-site interacts and stabilizes a single hydrogen atom on the Pro-tRNA in the A-site, thereby preventing an optimal geometry for the nucleophilic attack required for peptide bond formation to occur. This mechanism to short circuit the ribosomal peptidyltransferase activity is likely to operate for the majority of other RAPP-like arrest peptides found across diverse bacterial phylogenies.

The SecM arrest peptide traps a pre-peptide bond formation state of the ribosome

Nascent polypeptide chains can induce translational stalling to regulate gene expression. This is exemplified by the E. coli secretion monitor (SecM) arrest peptide that induces translational stalling to regulate expression of the downstream encoded SecA, an ATPase that co-operates with the SecYEG translocon to facilitate insertion of proteins into or through the cytoplasmic membrane. Here we present the structure of a ribosome stalled during translation of the full-length E. coli SecM arrest peptide at 2.0 Å resolution. The structure reveals that SecM arrests translation by stabilizing the Pro-tRNA in the A-site, but in a manner that prevents peptide bond formation with the SecM-peptidyl-tRNA in the P-site. By employing molecular dynamic simulations, we also provide insight into how a pulling force on the SecM nascent chain can relieve the SecM-mediated translation arrest. Collectively, the mechanisms determined here for SecM arrest and relief are also likely to be applicable for a variety of other arrest peptides that regulate components of the protein localization machinery identified across a wide range of bacteria lineages.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Differential Divalent Metal Binding by SpyCas9's RuvC Active Site Contributes to Nonspecific DNA Cleavage

To protect against mobile genetic elements (MGEs), some bacteria and archaea have clustered regularly interspaced short palindromic repeats-CRISPR associated (CRISPR-Cas) adaptive immune systems. CRISPR RNAs (crRNAs) bound to Cas nucleases hybridize to MGEs based on sequence complementarity to guide the nucleases to cleave the MGEs. This programmable DNA cleavage has been harnessed for gene editing. Safety concerns include off-target and guide RNA (gRNA)-free DNA cleavages, both of which are observed in the Cas nuclease commonly used for gene editing, Streptococcus pyogenes Cas9 (SpyCas9). We developed a SpyCas9 variant (SpyCas9H982A) devoid of gRNA-free DNA cleavage activity that is more selective for on-target cleavage. The H982A substitution in the metal-dependent RuvC active site reduces Mn2+-dependent gRNA-free DNA cleavage by ∼167-fold. Mechanistic molecular dynamics analysis shows that Mn2+, but not Mg2+, produces a gRNA-free DNA cleavage competent state that is disrupted by the H982A substitution. Our study demonstrates the feasibility of modulating cation:protein interactions to engineer safer gene editing tools.

Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch

Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.

Structural conservation of antibiotic interaction with ribosomes

The ribosome is a major target for clinically used antibiotics, but multidrug resistant pathogenic bacteria are making our current arsenal of antimicrobials obsolete. Here we present cryo-electron-microscopy structures of 17 distinct compounds from six different antibiotic classes bound to the bacterial ribosome at resolutions ranging from 1.6 to 2.2 Å. The improved resolution enables a precise description of antibiotic-ribosome interactions, encompassing solvent networks that mediate multiple additional interactions between the drugs and their target. Our results reveal a high structural conservation in the binding mode between antibiotics with the same scaffold, including ordered water molecules. Water molecules are visualized within the antibiotic binding sites that are preordered, become ordered in the presence of the drug and that are physically displaced on drug binding. Insight into RNA-ligand interactions will facilitate development of new antimicrobial agents, as well as other RNA-targeting therapies.

Structural insights on ionizable Dlin-MC3-DMA lipids in DOPC layers by combining accurate atomistic force fields, molecular dynamics simulations and neutron reflectivity

Ionizable lipids such as the promising Dlin-MC3-DMA (MC3) are essential for the successful design of lipid nanoparticles (LNPs) as drug delivery agents. Combining molecular dynamics simulations with experimental data, such as neutron reflectivity experiments and other scattering techniques, is essential to provide insights into the internal structure of LNPs, which is not fully understood to date. However, the accuracy of the simulations relies on the choice of force field parameters and high-quality experimental data is indispensable to verify the parametrization. For MC3, different parameterizations in combination with the CHARMM and the Slipids force fields have recently emerged. Here, we complement the existing efforts by providing parameters for cationic and neutral MC3 compatible with the AMBER Lipid17 force field. Subsequently, we carefully assess the accuracy of the different force fields by providing a direct comparison to neutron reflectivity experiments of mixed lipid bilayers consisting of MC3 and DOPC at different pHs. At low pH (cationic MC3) and at high pH (neutral MC3) the newly developed MC3 parameters in combination with AMBER Lipid17 for DOPC give good agreement with the experiments. Overall, the agreement is similar compared to the Park-Im parameters for MC3 in combination with the CHARMM36 force field for DOPC. The Ermilova-Swenson MC3 parameters in combination with the Slipids force field underestimate the bilayer thickness. While the distribution of cationic MC3 is very similar, the different force fields for neutral MC3 reveal distinct differences ranging from strong accumulation in the membrane center (current MC3/AMBER Lipid17 DOPC), over mild accumulation (Park-Im MC3/CHARMM36 DOPC) to surface accumulation (Ermilova-Swenson MC3/Slipids DOPC). These pronounced differences highlight the importance of accurate force field parameters and their experimental validation.

Trans-palmitoleic acid, a dairy fat biomarker, stimulates insulin secretion and activates G protein-coupled receptors with a different mechanism from the cis isomer

Dietary trans-palmitoleic acid (trans 16:1n-7, tPOA), a biomarker for high-fat dairy product intake, has been associated with a lower risk of type 2 diabetes mellitus (T2DM) in some cross-sectional and prospective epidemiological studies. Here, we investigated the insulin secretion-promoting activity of tPOA and compared them with the effects evoked by the cis-POA isomer (cPOA), an endogenous lipokine biosynthesized in the liver and adipose tissue, and found in some natural food sources. The debate about the positive and negative relationships of those two POA isomers with metabolic risk factors and the underlying mechanisms is still going on. Therefore, we examined the potency of both POA isomers to potentiate insulin secretion in murine and human pancreatic β cell lines. We also investigated whether POA isomers activate G protein-coupled receptors proposed as potential targets for T2DM treatment. We show that tPOA and cPOA augment glucose-stimulated insulin secretion (GSIS) to a similar extent; however, their insulin secretagogue activity is associated with different signaling pathways. We also performed ligand docking and molecular dynamics simulations to predict the preferred orientation of POA isomers and the strength of association between those two fatty acids and GPR40, GPR55, GPR119, and GPR120 receptors. Overall, this study provides insight into the bioactivity of tPOA and cPOA toward selected GPCR functions, indicating them as targets responsible for the insulin secretagogue action of POA isomers. It reveals that both tPOA and cPOA may promote insulin secretion and subsequently regulate glucose homeostasis.

RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations

The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li⁺, Na⁺, K⁺, Cs⁺, Ca²⁺, Sr²⁺, and Ba²⁺, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.