Understanding the origins of bacterial resistance to aminoglycosides through molecular dynamics mutational study of the ribosomal A-site

A computational workflow to determine drug candidates alternative to aminoglycosides targeting the decoding center of E. coli ribosome

The global antibiotic resistance problem necessitates fast and effective approaches to finding novel inhibitors to treat bacterial infections. In this study, we propose a computational workflow to identify plausible high-affinity compounds from FDA-approved, investigational, and experimental libraries for the decoding center on the small subunit 30S of the E. coli ribosome. The workflow basically consists of two molecular docking calculations on the intact 30S, followed by molecular dynamics (MD) simulations coupled with MM-GBSA calculations on a truncated ribosome structure. The parameters used in the molecular docking suits, Glide and AutoDock Vina, as well as in the MD simulations with Desmond were carefully adjusted to obtain expected interactions for the ligand-rRNA complexes. A filtering procedure was followed, considering a fingerprint based on aminoglycoside's binding site on the 30S to obtain seven hit compounds either with different clinical usages or aminoglycoside derivatives under investigation, suggested for in vitro studies. The detailed workflow developed in this study promises an effective and fast approach for the estimation of binding free energies of large protein-RNA and ligand complexes.

Unveiling Synergistic Potency: Exploring Butyrolactone I to Enhance Gentamicin Efficacy against Methicillin-Resistant Staphylococcus aureus (MRSA) Strain USA300

Staphylococcus aureus, including MRSA strains, poses significant health risks, imposing a significant disease burden and mortality. We investigate butyrolactone I (BL-1), a marine-derived metabolite from Aspergillus terreus, enhancing aminoglycoside efficacy against MRSA. A promising synergy is observed with BL-1 and various aminoglycosides, marked by low fractional inhibitory concentration indexes (FICIs < 0.5). Comprehensive studies utilizing USA300 MRSA and gentamicin reveal a remarkable one-fourth reduction in minimum inhibitory concentration (MIC) with 20 μg/mL BL-1. A relative abundance assay indicates that BL-1 enhances gentamicin uptake while restraining extracellular presence, involving intricate transmembrane signaling and molecular interactions. RNA-Seq analysis yielded an unexpected revelation, unveiling a distinctive gene expression profile and distinguishing it from other treatment approaches. Furthermore, meticulous analyses validated the extensive perturbations induced by BL-1 exposure, affecting diverse biological functions, encompassing glycolysis, amino acid metabolisms, substance transmembrane transport, and virulence generation. These valuable insights inspired further confirmation of bacterial virulence and the modulation of membrane permeability resulting from BL-1 treatment. Phenotypic validations corroborated our observations, revealing reduced membrane permeability and hemolytic toxicity, albeit demanding a deeper comprehension of the intricate interplay underlying these actions. Our study contributes crucial mechanistic insights to the development of therapeutic strategies against this notorious pathogen and the judicious employment of aminoglycosides. Additionally, it elucidates marine-derived metabolites' ecological and functional roles, exemplified by fungal quorum sensing signals. These compounds could give producers a competitive edge, inhibiting microorganism proliferation and suggesting novel approaches for combating resistant pathogens.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BACKGROUND

Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.

METHODS

We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.

RESULTS

We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.

CONCLUSIONS

The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.

GENERAL SIGNIFICANCE

We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Interactions of aminoglycoside antibiotics with rRNA

Aminoglycoside antibiotics are protein synthesis inhibitors applied to treat infections caused mainly by aerobic Gram-negative bacteria. Due to their adverse side effects they are last resort antibiotics typically used to combat pathogens resistant to other drugs. Aminoglycosides target ribosomes. We describe the interactions of aminoglycoside antibiotics containing a 2-deoxystreptamine (2-DOS) ring with 16S rRNA. We review the computational studies, with a focus on molecular dynamics (MD) simulations performed on RNA models mimicking the 2-DOS aminoglycoside binding site in the small ribosomal subunit. We also briefly discuss thermodynamics of interactions of these aminoglycosides with their 16S RNA target.

Investigation of Ribosomes Using Molecular Dynamics Simulation Methods

The ribosome as a complex molecular machine undergoes significant conformational changes while synthesizing a protein molecule. Molecular dynamics simulations have been used as complementary approaches to X-ray crystallography and cryoelectron microscopy, as well as biochemical methods, to answer many questions that modern structural methods leave unsolved. In this review, we demonstrate that all-atom modeling of ribosome molecular dynamics is particularly useful in describing the process of tRNA translocation, atomic details of behavior of nascent peptides, antibiotics, and other small molecules in the ribosomal tunnel, and the putative mechanism of allosteric signal transmission to functional sites of the ribosome.

Antibiotics: where to throw the spanner in the ribosomal machinery?

The sheer molecular scale of the ribosome is intimidating to the traditional drug designer. By analyzing the ribosome as a series of 12 key target sites, this review seeks to make the ribosome ligand design process more manageable. Analysis of recently evaluated ribosomal structures, particularly those with bound antibiotics, indicates where the ligand target sites are located. This review employs current research data to map antibiotic binding across the ribosome. A number of neighboring ligand-binding sites are often contiguous and can be combined. Ligands that bind in close proximity can be combined into hybrid structures. The different ways antibiotics disrupt ribosomal function are also discussed. Antibiotics tend to inhibit conformational changes that are essential to the ribosomal mechanism.

Drug-target networks in aminoglycoside resistance: hierarchy of priority in structural drug design

Drug-target network analysis for advancing next-generation aminoglycoside therapies that combat antibiotic resistant infections.

Conformational dynamics of bacterial and human cytoplasmic models of the ribosomal A-site

The aminoacyl-tRNA binding site (A-site) is located in helix 44 of small ribosomal subunit. The mobile adenines 1492 and 1493 (Escherichia coli numbering), forming the A-site bulge, act as a functional switch that ensures mRNA decoding accuracy. Structural data on the oligonucleotide models mimicking the ribosomal A-site with sequences corresponding to bacterial and human cytoplasmic sites confirm that this RNA motif forms also without the ribosome context. We performed all-atom molecular dynamics simulations of these crystallographic A-site models to compare their conformational properties. We found that the human A-site bulge is more internally flexible than the bacterial one and has different base pairing preferences, which result in the overall different shapes of these bulges and cation density distributions. Also, in the human A-site model we observed repetitive destacking of A1492, while A1493 was more stably paired than in the bacterial variant. Based on the dynamics of the A-sites we suggest why aminoglycoside antibiotics, which target the bacterial A-site, have lower binding affinities and anti-translational activities toward the human variant.

Interplay of the bacterial ribosomal A-site, S12 protein mutations and paromomycin binding: a molecular dynamics study

The conformational properties of the aminoacyl-tRNA binding site (A-site), and its surroundings in the Escherichia coli 30S ribosomal subunit, are of great relevance in designing antibacterial agents. The 30S subunit A-site is near ribosomal protein S12, which neighbors helices h27 and H69; this latter helix, of the 50S subunit, is a functionally important component of an intersubunit bridge. Experimental work has shown that specific point mutations in S12 (K42A, R53A) yield hyper-accurate ribosomes, which in turn confers resistance to the antibiotic 'paromomycin' (even when this aminoglycoside is bound to the A-site). Suspecting that these effects can be elucidated in terms of the local atomic interactions and detailed dynamics of this region of the bacterial ribosome, we have used molecular dynamics simulations to explore the motion of a fragment of the E. coli ribosome, including the A-site. We found that the ribosomal regions surrounding the A-site modify the conformational space of the flexible A-site adenines 1492/93. Specifically, we found that A-site mobility is affected by stacking interactions between adenines A1493 and A1913, and by contacts between A1492 and a flexible side-chain (K43) from the S12 protein. In addition, our simulations reveal possible indirect pathways by which the R53A and K42A mutations in S12 are coupled to the dynamical properties of the A-site. Our work extends what is known about the atomistic dynamics of the A-site, and suggests possible links between the biological effects of hyper-accurate mutations in the S12 protein and conformational properties of the ribosome; the implications for S12 dynamics help elucidate how the miscoding effects of paromomycin may be evaded in antibiotic-resistant mutants of the bacterial ribosome.

Role of S-turn2 in the structure, dynamics, and function of mitochondrial ribosomal A-site. A bioinformatics and molecular dynamics simulation study

The mRNA decoding site (A-site) in the small ribosomal subunit controls fidelity of the translation process. Here, using molecular dynamics simulations and bioinformatic analyses, we investigated the structural dynamics of the human mitochondrial A-site (native and A1490G mutant) and compared it with the dynamics of the bacterial A-site. We detected and characterized a specific RNA backbone configuration, S-turn2, which occurs in the human mitochondrial but not in the bacterial A-site. Mitochondrial and bacterial A-sites show different propensities to form S-turn2 that may be caused by different base-pairing patterns of the flanking nucleotides. Also, the S-turn2 structural stability observed in the simulations supports higher accuracy and lower speed of mRNA decoding in mitochondria in comparison with bacteria. In the mitochondrial A-site, we observed collective movement of stacked nucleotides A1408·C1409·C1410, which may explain the known differences in aminoglycoside antibiotic binding affinities toward the studied A-site variants.

Flipping of the ribosomal A-site adenines provides a basis for tRNA selection

Ribosomes control the missense error rate of ~10(-4) during translation though quantitative contributions of individual mechanistic steps of the conformational changes yet to be fully determined. Biochemical and biophysical studies led to a qualitative tRNA selection model in which ribosomal A-site residues A1492 and A1493 (A1492/3) flip out in response to cognate tRNA binding, promoting the subsequent reactions, but not in the case of near-cognate or non-cognate tRNA. However, this model was recently questioned by X-ray structures revealing conformations of extrahelical A1492/3 and domain closure of the decoding center in both cognate and near-cognate tRNA bound ribosome complexes, suggesting that the non-specific flipping of A1492/3 has no active role in tRNA selection. We explore this question by carrying out molecular dynamics simulations, aided with fluorescence and NMR experiments, to probe the free energy cost of extrahelical flipping of 1492/3 and the strain energy associated with domain conformational change. Our rigorous calculations demonstrate that the A1492/3 flipping is indeed a specific response to the binding of cognate tRNA, contributing 3kcal/mol to the specificity of tRNA selection. Furthermore, the different A-minor interactions in cognate and near-cognate complexes propagate into the conformational strain and contribute another 4kcal/mol in domain closure. The recent structure of ribosome with features of extrahelical A1492/3 and closed domain in near-cognate complex is reconciled by possible tautomerization of the wobble base pair in mRNA-tRNA. These results quantitatively rationalize other independent experimental observations and explain the ribosomal discrimination mechanism of selecting cognate versus near-cognate tRNA.

Interactions of amikacin with the RNA model of the ribosomal A-site: computational, spectroscopic and calorimetric studies

Amikacin is a 2-deoxystreptamine aminoglycoside antibiotic possessing a unique l-HABA (l-(-)-γ-amino-α-hydroxybutyric acid) group and applied in the treatment of hospital-acquired infections. Amikacin influences bacterial translation by binding to the decoding region of the small ribosomal subunit that overlaps with the binding site of aminoacylated-tRNA (A-site). Here, we have characterized thermodynamics of interactions of amikacin with a 27-mer RNA oligonucleotide mimicking the aminoglycoside binding site in the bacterial ribosome. We applied isothermal titration and differential scanning calorimetries, circular dichroism and thermal denaturation experiments, as well as computer simulations. Thermal denaturation studies have shown that amikacin affects only slightly the melting temperatures of the A-site mimicking RNA model suggesting a moderate stabilization of RNA by amikacin. Isothermal titration calorimetry gives the equilibrium dissociation constants for the binding reaction between amikacin and the A-site oligonucleotide in the micromolar range with a favorable enthalpic contribution. However, for amikacin we observe a positive entropic contribution to binding, contrary to other aminoglycosides, paromomycin and ribostamycin. Circular dichroism spectra suggest that the observed increase in entropy is not caused by structural changes of RNA because amikacin binding does not destabilize the helicity of the RNA model. To investigate the origins of this positive entropy change we performed all-atom molecular dynamics simulations in explicit solvent for the 27-mer RNA oligonucleotide mimicking one A-site and the crystal structure of an RNA duplex containing two A-sites. We observed that the diversity of the conformational states of the l-HABA group sampled in the simulations of the complex was larger than for the free amikacin in explicit water. Therefore, the larger flexibility of the l-HABA group in the bound form may contribute to an increase of entropy upon binding.

Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop

The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.

Comparing aminoglycoside binding sites in bacterial ribosomal RNA and aminoglycoside modifying enzymes

Aminoglycoside antibiotics are used against severe bacterial infections. They bind to the bacterial ribosomal RNA and interfere with the translation process. However, bacteria produce aminoglycoside modifying enzymes (AME) to resist aminoglycoside actions. AMEs form a variable group and yet they specifically recognize and efficiently bind aminoglycosides, which are also diverse in terms of total net charge and the number of pseudo-sugar rings. Here, we present the results of 25 molecular dynamics simulations of three AME representatives and aminoglycoside ribosomal RNA binding site, unliganded and complexed with an aminoglycoside, kanamycin A. A comparison of the aminoglycoside binding sites in these different receptors revealed that the enzymes efficiently mimic the nucleic acid environment of the ribosomal RNA binding cleft. Although internal dynamics of AMEs and their interaction patterns with aminoglycosides differ, the energetical analysis showed that the most favorable sites are virtually the same in the enzymes and RNA. The most copied interactions were of electrostatic nature, but stacking was also replicated in one AME:kanamycin complex. In addition, we found that some water-mediated interactions were very stable in the simulations of the complexes. We show that our simulations reproduce well findings from NMR or X-ray structural studies, as well as results from directed mutagenesis. The outcomes of our analyses provide new insight into aminoglycoside resistance mechanism that is related to the enzymatic modification of these drugs.

Novel insights of structure-based modeling for RNA-targeted drug discovery

Substantial progress in RNA biology highlights the importance of RNAs (e.g., microRNAs) in diseases and the potential of targeting RNAs for drug discovery. However, the lack of RNA-specific modeling techniques demands the development of new tools for RNA-targeted rational drug design. Herein, we implemented integrated approaches of accurate RNA modeling and virtual screening for RNA inhibitor discovery with the most comprehensive evaluation to date of five docking and 11 scoring methods. For the first time, statistical analysis was heavily employed to assess the significance of our predictions. We found that GOLD:GOLD Fitness and rDock:rDock_solv could accurately predict the RNA ligand poses, and ASP rescoring further improved the ranking of ligand binding poses. Due to the weak correlations (R(2) < 0.3) of existing scoring with experimental binding affinities, we implemented two new RNA-specific scoring functions, iMDLScore1 and iMDLScore2, and obtained better correlations with R(2) = 0.70 and 0.79, respectively. We also proposed a multistep virtual screening approach and demonstrated that rDock:rDock_solv together with iMDLScore2 rescoring obtained the best enrichment on the flexible RNA targets, whereas GOLD:GOLD Fitness combined with rDock_solv rescoring outperformed other methods for rigid RNAs. This study provided practical strategies for RNA modeling and offered new insights into RNA-small molecule interactions for drug discovery.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.