Theoretical study of the interaction between the antibiotic linezolid and the active site of the 50S ribosomal subunit of the bacterium Haloarcula marismortui.

First antibiotic in the oxazolidinone class, linezolid fights gram-positive multiresistant bacteria by inhibiting protein synthesis through its interaction with the 50S subunit of the functional bacterial ribosome. For its antimicrobial action, it is necessary that its chiral carbon located in the oxazolidinone ring is in the S-conformation. Computational calculation at time-dependent density functional theory methodology, ultraviolet-visible (UV-Vis), and electronic circular dichroism spectra was obtained for noncomplexed and complexed forms of linezolid to verify the possible chirality of nitrogen atom in the acetamide group of the molecule. The molecular system has two chiral centers. So, there are now four possible configurations: RR, RS, SR, and SS. For a better understanding of the system, the electronic spectra at the PBE0/6-311++G(3df,2p) level of theory were obtained. The complexed form was obtained from the crystallographic data of the ribosome, containing the S-linezolid molecular system. The computational results obtained for the electronic properties are in good agreement with the experimental crystallographic data and available theoretical results.

Structure and folding of four putative kink turns identified in structured RNA species in a test of structural prediction rules.

k-Turns are widespread key architectural elements that occur in many classes of RNA molecules. We have shown previously that their folding properties (whether or not they fold into their tightly kinked structure on addition of metal ions) and conformation depend on their local sequence, and we have elucidated a series of rules for prediction of these properties from sequence. In this work, we have expanded the rules for prediction of folding properties, and then applied the full set to predict the folding and conformation of four probable k-turns we have identified amongst 224 structured RNA species found in bacterial intergenenic regions by the Breaker lab (1). We have analyzed the ion-dependence of folding of the four k-turns using fluorescence resonance energy transfer, and determined the conformation of two of them using X-ray crystallography. We find that the experimental data fully conform to both the predicted folding and conformational properties. We conclude that our folding rules are robust, and can be applied to new k-turns of unknown characteristics with confidence.

Experimental study of proteome halophilicity using nanoDSF: a proof of concept.

Adaption to environmental conditions is reflected by protein adaptation. In particular, proteins of extremophiles display distinctive traits ensuring functional, structural and dynamical properties under permanently extreme physical and chemical conditions. While it has mostly been studied with approaches focusing on specific proteins, biophysical approaches have also confirmed this link between environmental and protein adaptation at the more complex and diverse scale of the proteome. However, studies of this type remain challenging and often require large amounts of biological material. We report here the use of nanoDSF as a tool to study proteome stability and solubility in cell lysates of the model halophilic archaeon Haloarcula marismortui. Notably, our results show that, as with single halophilic protein studies, proteome stability was correlated to the concentration of NaCl or KCl under which the cells were lysed and hence the proteome exposed. This work highlights that adaptation to environmental conditions can be experimentally observed at the scale of the proteome. Still, we show that the biochemical properties of single halophilic proteins can only be partially extrapolated to the whole proteome.

Characterization of Regulatory Elements of L11 and L1 Operons in Thermophilic Bacteria and Archaea.

Ribosomal protein L1 is a conserved two-domain protein that is involved in formation of the L1 stalk of the large ribosomal subunit. When there are no free binding sites available on the ribosomal 23S RNA, the protein binds to the specific site on the mRNA of its own operon (L11 operon in bacteria and L1 operon in archaea) preventing translation. Here we show that the regulatory properties of the r-protein L1 and its domain I are conserved in the thermophilic bacteria Thermus and Thermotoga and in the halophilic archaeon Haloarcula marismortui. At the same time the revealed features of the operon regulation in thermophilic bacteria suggest presence of two regulatory regions.

Prebiotic Iron Originates the Peptidyl Transfer Origin.

The ribosome is responsible for protein synthesis in all living organisms. It is best known to exist around 3.5-3.7 Ga whereat life on Earth inhabited anoxic environment with abundant soluble irons. The RNAs and proteins are the two biopolymers that constitute the ribosome. However, both proteins and RNAs require metal cations to fold and to function. There are four Mg-microcluster (Mg2+-µc) structures conserved in core of large subunit, and the 23S ribosomal RNA (rRNA) was shown to catalyze electron transfer in an anoxic environment in the presence of Fe2+. The Mg2+-µc features two idiosyncratic Mg2+ ions that are chelated and bridged by a common phosphate group and along with that, the adjacent residues of RNA backbone together forming ten-membered chelation ring(s). Here, we utilized four rRNA fragments of the large subunit 23S rRNA of Haloarcula marismortui, that includes the residues that form the four Mg2+-µc's. These four rRNA fragments are shown competent to assemble with Mg2+. Our results show that when these rRNA fragments fold or assembly in the presence of Fe2+ under anoxic conditions, each Fe2+-microcluster can catalyze electron transfer. We propose that Fe2+-microclusters of the ribosome, which use Fe2+ as a cofactor to regulate electron transfer, are pivotal and primordial and may be an origin in evolution of the ribosome.

Local and Global Motions Underlying Antibiotic Binding in Bacterial Ribosome.

The bacterial ribosome is one of the most important targets in the treatment of infectious diseases. As antibiotic resistance in bacteria poses a growing threat, a significant amount of effort is concentrated on exploring new drug-binding sites where testable predictions are of significance. Here, we study the dynamics of a ribosomal complex and 67 small and large subunits of the ribosomal crystal structures (64 antibiotic-bound, 3 antibiotic-free) from Deinococcus radiodurans, Escherichia coli, Haloarcula marismortui, and Thermus thermophilus by the Gaussian network model. Interestingly, a network of nucleotides coupled in high-frequency fluctuations reveals known antibiotic-binding sites. These sites are seen to locate at the interface of dynamic domains that have an intrinsic dynamic capacity to interfere with functional globular motions. The nucleotides and the residues fluctuating in the fast and slow modes of motion thus have promise for plausible antibiotic-binding and allosteric sites that can alter antibiotic binding and resistance. Overall, the present analysis brings a new dynamic perspective to the long-discussed link between small-molecule binding and large conformational changes of the supramolecule.

Verification of the Stabilized Protein Design Based on the Prediction of Intrinsically Disordered Regions: Ribosomal Proteins L1.

In our previous papers, we proposed the idea that programs predicting intrinsically disordered regions in amino acid sequences can be used for finding weakened sites in proteins. The regions predicted by such programs are suitable targets for the introduction of protein-stabilizing mutations. However, for each specific protein, it remains unclear what determines protein stabilization - the amino acid sequence (and accordingly, prediction of weakened sites) or the 3D structure. To answer this question, it is necessary to study two proteins with similar structures but different amino acid sequences and, consequently, different predictions of weakened regions. By introducing identical mutations into identical elements of the two proteins, we will be able to reveal whether predictions of the weakened sites or the 3D protein structure are the key factors in the protein stability increase. Here, we have chosen ribosomal proteins L1 from the halophilic archaeon Haloarcula marismortui (HmaL1) and extremophilic bacterium Aquifex aeolicus (AaeL1). These proteins are identical in their structure but different in amino acid sequences. A disulfide bond introduced into the region predicted as the structured one in AaeL1 did not lead to the increase in the protein melting temperature. At the same time, a disulfide bond introduced into the same region in HmaL1 that was predicted as a weakened one, resulted in the increase in the protein melting temperature by approximately 10°C.

Stationary Phase EPR Spectroscopy for Monitoring Membrane Protein Refolding by Conformational Response.

Protein production remains a major bottleneck in membrane protein structural biology. In many cases, large-scale recombinant protein expression is either unfeasible or impossible, driving structural biologists to explore new production avenues. Several membrane proteins have been successfully refolded from solubilized E. coli inclusion bodies. In recent years, a structure of the G-protein-coupled receptor CXCR1 was obtained using refolded material from E. coli inclusion bodies. However, aggregation during the refolding process is a common difficulty, which is often addressed by immobilization of the protein onto a solid support. Most spectroscopic methods are incompatible with these light-scattering matrices, which renders automated buffer exchange to screen refolding conditions impossible. This work explores a potential approach to overcome this problem by utilizing site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy of protein bound to standard, commercially available Ni-NTA agarose resin. With this approach, the correct protein fold is determined by activity, which is inferred from a protein conformational response to a known stimulant. EPR spectra at each state of the refolding workflow of spin-labeled Haloarcula marismortui bacteriorhodopsin-I (HmbRI) are obtained, and refolded fractions of HmbRI with this platform are quantitated using both protein from inclusion bodies and denatured recombinant protein from E. coli membranes. The stimulant used for HmbRI is visible light. The solid support allows for multiple refolding trials through buffer exchanges, and the EPR spectra are collected on the order of seconds under ambient conditions.

Topology independent comparison of RNA 3D structures using the CLICK algorithm.

RNA molecules are attractive therapeutic targets because non-coding RNA molecules have increasingly been found to play key regulatory roles in the cell. Comparing and classifying RNA 3D structures yields unique insights into RNA evolution and function. With the rapid increase in the number of atomic-resolution RNA structures, it is crucial to have effective tools to classify RNA structures and to investigate them for structural similarities at different resolutions. We previously developed the algorithm CLICK to superimpose a pair of protein 3D structures by clique matching and 3D least squares fitting. In this study, we extend and optimize the CLICK algorithm to superimpose pairs of RNA 3D structures and RNA-protein complexes, independent of the associated topologies. Benchmarking Rclick on four different datasets showed that it is either comparable to or better than other structural alignment methods in terms of the extent of structural overlaps. Rclick also recognizes conformational changes between RNA structures and produces complementary alignments to maximize the extent of detectable similarity. Applying Rclick to study Ribonuclease III protein correctly aligned the RNA binding sites of RNAse III with its substrate. Rclick can be further extended to identify ligand-binding pockets in RNA. A web server is developed at http://mspc.bii.a-star.edu.sg/minhn/rclick.html.

Synthesis and Biological Evaluation of Pyrimidine-oxazolidin-2-arylimino Hybrid Molecules as Antibacterial Agents.

Pyrimidine-1,3-oxazolidin-2-arylimino hybrids have been synthesized as a new class of antibacterial agents. The synthetic approach exploits a Cu(II)-catalyzed intramolecular halkoxyhalogenation of alkynyl ureas, followed by a Suzuki coupling reaction with 2,4-dimethoxypyrimidin-5-boronic acid. Biological screenings revealed that most of the compounds showed moderate to good activity against two Gram-positive (B. subtilis, S. aureus) and three Gram-negative (P. aeruginosa, S. typhi, K. pneumonia) pathogenic strains. A molecular docking study, performed in the crystal structure of 50S ribosomal unit of Haloarcula marismortui, indicated that pyrimidine-oxazolidin-2-arylimino hybrids 8c and 8h exhibited a high binding affinity (-9.65 and -10.74 kcal/mol), which was in agreement with their good antibacterial activity. The obtained results suggest that the combination of pyrimidine and oxazolidone moieties can be considered as a valid basis to develop new further modifications towards more efficacious antibacterial compounds.

Identification and codon reading properties of 5-cyanomethyl uridine, a new modified nucleoside found in the anticodon wobble position of mutant haloarchaeal isoleucine tRNAs.

Most archaea and bacteria use a modified C in the anticodon wobble position of isoleucine tRNA to base pair with A but not with G of the mRNA. This allows the tRNA to read the isoleucine codon AUA without also reading the methionine codon AUG. To understand why a modified C, and not U or modified U, is used to base pair with A, we mutated the C34 in the anticodon of Haloarcula marismortui isoleucine tRNA (tRNA2(Ile)) to U, expressed the mutant tRNA in Haloferax volcanii, and purified and analyzed the tRNA. Ribosome binding experiments show that although the wild-type tRNA2(Ile) binds exclusively to the isoleucine codon AUA, the mutant tRNA binds not only to AUA but also to AUU, another isoleucine codon, and to AUG, a methionine codon. The G34 to U mutant in the anticodon of another H. marismortui isoleucine tRNA species showed similar codon binding properties. Binding of the mutant tRNA to AUG could lead to misreading of the AUG codon and insertion of isoleucine in place of methionine. This result would explain why most archaea and bacteria do not normally use U or a modified U in the anticodon wobble position of isoleucine tRNA for reading the codon AUA. Biochemical and mass spectrometric analyses of the mutant tRNAs have led to the discovery of a new modified nucleoside, 5-cyanomethyl U in the anticodon wobble position of the mutant tRNAs. 5-Cyanomethyl U is present in total tRNAs from euryarchaea but not in crenarchaea, eubacteria, or eukaryotes.

The plasticity of a structural motif in RNA: structural polymorphism of a kink turn as a function of its environment.

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2' at the -1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.

The molecular recognition of kink-turn structure by the L7Ae class of proteins.

L7Ae is a member of a protein family that binds kink-turns (k-turns) in many functional RNA species. We have solved the X-ray crystal structure of the near-consensus sequence Kt-7 of Haloarcula marismortui bound by Archaeoglobus fulgidus L7Ae at 2.3-Å resolution. We also present a structure of Kt-7 in the absence of bound protein at 2.2-Å resolution. As a result, we can describe a general mode of recognition of k-turn structure by the L7Ae family proteins. The protein makes interactions in the widened major groove on the outer face of the k-turn. Two regions of the protein are involved. One is an α-helix that enters the major groove of the NC helix, making both nonspecific backbone interactions and specific interactions with the guanine nucleobases of the conserved G • A pairs. A hydrophobic loop makes close contact with the L1 and L2 bases, and a glutamate side chain hydrogen bonds with L1. Taken together, these interactions are highly selective for the structure of the k-turn and suggest how conformational selection of the folded k-turn occurs.

Potassium stress growth characteristics and energetics in the haloarchaeon Haloarcula marismortui.

Growth characteristics surrounding halophilic archaeal organisms are extremely limited in the scientific literature, with studies tending toward observing changes in cellular generation times under growth conditions limited to changes in temperature and sodium chloride concentrations. Currently, knowledge of the ionic stress experienced by haloarchaeal species through an excess or depletion of other required ions is lacking at best. The halophilic archaeon, Haloarcula marismortui, was analyzed under extreme ionic stress conditions with a specific focus on induced potassium ion stress using growth curves and analysis of the intracellular ion concentrations. Generation times were determined under potassium chloride concentrations ranging from 8 to 720 mM, and also in the presence of the alternative monovalent cations of lithium, rubidium, and cesium under limiting potassium conditions. Intracellular ion concentrations, as determined by inductively coupled mass spectrometry (ICP-MS), indicate a minimum intracellular total ion requirement of 1.13 M while tolerating up to 2.43 M intracellular concentrations. The presence of intracellular rubidium and cesium indicates that monovalent ion transport is important for energy production. Comparison of eight archaeal genomes indicates an increased diversity of potassium transport complex subunits in the halophilic organisms. Analysis of the generation times, intracellular concentrations and genome survey shows Har. marismortui exhibits an ability to cope with monovalent cation concentration changes in its native environment and provides insight into the organisms ion transport capability and specificity.

Substrate promiscuity: AglB, the archaeal oligosaccharyltransferase, can process a variety of lipid-linked glycans.

Across evolution, N-glycosylation involves oligosaccharyltransferases that transfer lipid-linked glycans to selected Asn residues of target proteins. While these enzymes catalyze similar reactions in each domain, differences exist in terms of the chemical composition, length and degree of phosphorylation of the lipid glycan carrier, the sugar linking the glycan to the lipid carrier, and the composition and structure of the transferred glycan. To gain insight into how oligosaccharyltransferases cope with such substrate diversity, the present study analyzed the archaeal oligosaccharyltransferase AglB from four haloarchaeal species. Accordingly, it was shown that despite processing distinct lipid-linked glycans in their native hosts, AglB from Haloarcula marismortui, Halobacterium salinarum, and Haloferax mediterranei could readily replace their counterpart from Haloferax volcanii when introduced into Hfx. volcanii cells deleted of aglB. As the four enzymes show significant sequence and apparently structural homology, it appears that the functional similarity of the four AglB proteins reflects the relaxed substrate specificity of these enzymes. Such demonstration of AglB substrate promiscuity is important not only for better understanding of N-glycosylation in Archaea and elsewhere but also for efforts aimed at transforming Hfx. volcanii into a glycoengineering platform.

Haloarcula marismortui archaellin genes as ecoparalogs.

The genome of haloarchaeon Haloarcula marismortui contains two archaellin genes-flaA2 and flaB. Earlier we isolated and characterized two H. marismortui strains in that archaella consisting of FlaA2 archaellin (with a minor FlaB fraction) or of FlaB only. Both the FlaA2 and FlaB strains were motile and produced functional helical archaella. Thus, it may seem that the FlaA2 archaellin is redundant. In this study we investigated the biological roles of archaellin redundancy and demonstrated that FlaA2 archaellin is better adapted to more severe conditions of high temperature/low salinity, while FlaB has an advantage with increasing salinity. We used the thermodynamic data and bioinformatics sequence analysis to demonstrate that archaella formed by FlaA2 are more stable than those formed by FlaB. Our combined data indicate that the monomer FlaA2 archaellin is more flexible and leads to more compact and stable formation of filamentous structures. The difference in response to environmental stress indicates that FlaA2 and FlaB replace each other under different environmental conditions and can be considered as ecoparalogs.

Estimation of bacterioruberin by Raman spectroscopy during the growth of halophilic archaeon Haloarcula marismortui.

Halophilic archaea are interesting microorganisms that produce low biomass and metabolites, complicating their quantification. Raman spectroscopy (RS) is a powerful technique, which requires small samples, attractive for using in archaeal research. The objective of this work was the estimation of bacterioruberin content along with Haloarcula marismortui growth and their correlation with biomass concentration. RS was used to detect characteristic bands of bacterioruberin (vibrational modes C═CH, C─C, and C═C) in H. marismortui culture samples. The intensity of Raman spectra in bacterioruberin and the biomass concentration were adequately correlated. The highest production of bacterioruberin occurred at 60 h. RS is revealed as a reliable technique for the estimation of bacterioruberin in the biomass of H. marismortui, which could be considered as a promising qualitative and quantitative technique to assay metabolites in cell cultures.

Revisiting the Haloarcula marismortui 50S ribosomal subunit model.

The structure of the large ribosomal subunit from the halophilic archaeon Haloarcula marismortui (Hma) is the only crystal structure of an archaeal ribosomal particle that has been determined to date. However, the first model of the Hma 50S ribosomal subunit contained some gaps: the structures of functionally important mobile lateral protuberances were not visualized. Subsequently, some parts of the P (L12) stalk base were visualized at 3.0 Å resolution [Kavran & Steitz (2007), J. Mol. Biol. 371, 1047-1059]: the RNA-binding domain of r-protein P0 (L10), the C-terminal domain of L11 and helices 43 and 44 of the 23 S rRNA. Here, the 2.4 Å resolution electron-density map of the Hma 50S ribosomal subunit was revisited and approximately two-thirds of the P0 protein, residues 1-58 of the N-terminal domains of two P1 protein molecules, residues 130-156 of L11, the full-length r-protein LX, nucleotides 2137-2149 and 2226-2237 of the 23S rRNA helix H76 forming the L1 stalk, nucleotides 2339-2343 of the 23S rRNA (contacting L5 protein) and loops 29-34 and 108-128 of protein L5 could be visualized. Thus, this paper provides a supplemented version of the Hma 50S ribosomal subunit model.

Determinants of the species selectivity of oxazolidinone antibiotics targeting the large ribosomal subunit.

Oxazolidinone antibiotics bind to the highly conserved peptidyl transferase center in the ribosome. For developing selective antibiotics, a profound understanding of the selectivity determinants is required. We have performed for the first time technically challenging molecular dynamics simulations in combination with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculations of the oxazolidinones linezolid and radezolid bound to the large ribosomal subunits of the eubacterium Deinococcus radiodurans and the archaeon Haloarcula marismortui. A remarkably good agreement of the computed relative binding free energy with selectivity data available from experiment for linezolid is found. On an atomic level, the analyses reveal an intricate interplay of structural, energetic, and dynamic determinants of the species selectivity of oxazolidinone antibiotics: A structural decomposition of free energy components identifies influences that originate from first and second shell nucleotides of the binding sites and lead to (opposing) contributions from interaction energies, solvation, and entropic factors. These findings add another layer of complexity to the current knowledge on structure-activity relationships of oxazolidinones binding to the ribosome and suggest that selectivity analyses solely based on structural information and qualitative arguments on interactions may not reach far enough. The computational analyses presented here should be of sufficient accuracy to fill this gap.

Revisiting the structures of several antibiotics bound to the bacterial ribosome.

The increasing prevalence of antibiotic-resistant pathogens reinforces the need for structures of antibiotic-ribosome complexes that are accurate enough to enable the rational design of novel ribosome-targeting therapeutics. Structures of many antibiotics in complex with both archaeal and eubacterial ribosomes have been determined, yet discrepancies between several of these models have raised the question of whether these differences arise from species-specific variations or from experimental problems. Our structure of chloramphenicol in complex with the 70S ribosome from Thermus thermophilus suggests a model for chloramphenicol bound to the large subunit of the bacterial ribosome that is radically different from the prevailing model. Further, our structures of the macrolide antibiotics erythromycin and azithromycin in complex with a bacterial ribosome are indistinguishable from those determined of complexes with the 50S subunit of Haloarcula marismortui, but differ significantly from the models that have been published for 50S subunit complexes of the eubacterium Deinococcus radiodurans. Our structure of the antibiotic telithromycin bound to the T. thermophilus ribosome reveals a lactone ring with a conformation similar to that observed in the H. marismortui and D. radiodurans complexes. However, the alkyl-aryl moiety is oriented differently in all three organisms, and the contacts observed with the T. thermophilus ribosome are consistent with biochemical studies performed on the Escherichia coli ribosome. Thus, our results support a mode of macrolide binding that is largely conserved across species, suggesting that the quality and interpretation of electron density, rather than species specificity, may be responsible for many of the discrepancies between the models.