Electrostatics of Prokaryotic Ribosome and Its Biological Implication

Single-residue effects on the behavior of a nascent polypeptide chain inside the ribosome exit tunnel

Nascent polypeptide chains (NCs) are extruded from the ribosome through an exit tunnel (ET) traversing the large ribosomal subunit. The ET's irregular and chemically complex wall allows for various NC-ET interactions. Translational arrest peptides (APs) bind in the ET to induce translational arrest, a property that can be exploited to study NC-ET interactions by Force Profile Analysis (FPA). We employed FPA and molecular dynamics (MD) simulations to investigate how individual residues placed in a glycine-serine repeat segment within an AP-stalled NC interact with the ET to exert a pulling force on the AP and release stalling. Our results indicate that large and hydrophobic residues generate a pulling force on the NC when placed ≳10 residues away from the peptidyl transfer center (PTC). Moreover, an asparagine placed 12 residues from the PTC makes a specific stabilizing interaction with the tip of ribosomal protein uL22 that reduces the pulling force on the NC, while a lysine or leucine residue in the same position increases the pulling force. Finally, the MD simulations suggest how the Mannheimia succiniproducens SecM AP interacts with the ET to promote translational stalling.

Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations

Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.

A simple geometrical model of the electrostatic environment around the catalytic center of the ribosome and its significance for the elongation cycle kinetics

The central function of the large subunit of the ribosome is to catalyze peptide bond formation. This biochemical reaction is conducted at the peptidyl transferase center (PTC). Experimental evidence shows that the catalytic activity is affected by the electrostatic environment around the peptidyl transferase center. Here, we set up a minimal geometrical model fitting the available x-ray solved structures of the ribonucleic cavity around the catalytic center of the large subunit of the ribosome. The purpose of this phenomenological model is to estimate quantitatively the electrostatic potential and electric field that are experienced during the peptidyl transfer reaction. At least two reasons motivate the need for developing this quantification. First, we inquire whether the electric field in this particular catalytic environment, made only of nucleic acids, is of the same order of magnitude as the one prevailing in catalytic centers of the proteic enzymes counterparts. Second, the protein synthesis rate is dependent on the nature of the amino acid sequentially incorporated in the nascent chain. The activation energy of the catalytic reaction and its detailed kinetics are shown to be dependent on the mechanical work exerted on the amino acids by the electric field, especially when one of the four charged amino acid residues (R, K, E, D) has previously been incorporated at the carboxy-terminal end of the peptidyl-tRNA. Physical values of the electric field provide quantitative knowledge of mechanical work, activation energy and rate of the peptide bond formation catalyzed by the ribosome. We show that our theoretical calculations are consistent with two independent sets of previously published experimental results. Experimental results for E.coli in the minimal case of the dipeptide bond formation when puromycin is used as the final amino acid acceptor strongly support our theoretically derived reaction time courses. Experimental Ribo-Seq results on E. coli and S. cerevisiae comparing the residence time distribution of ribosomes upon specific codons are also well accounted for by our theoretical calculations. The statistical queueing time theory was used to model the ribosome residence time per codon during nascent protein elongation and applied for the interpretation of the Ribo-Seq data. The hypo-exponential distribution fits the residence time observed distribution of the ribosome on a codon. An educated deconvolution of this distribution is used to estimate the rates of each elongation step in a codon specific manner. Our interpretation of all these results sheds light on the functional role of the electrostatic profile around the PTC and its impact on the ribosome elongation cycle.

Biofilm characteristics and transcriptomic profiling of Acinetobacter johnsonii defines signatures for planktonic and biofilm cells

Most bacteria in the natural environment have a biofilm mode of life, which is intrinsically tolerant to antibiotics. While until now, the knowledge of biofilm formation by Acinetobacter johnsonii is not well understood. In this study, the characteristics and the effect of a sub-inhibitory concentration of antibiotic on A. johnsonii biofilm and planktonic cells were determined. We discovered a positive relationship between biofilm formation and tetracycline resistance, and biofilms rapidly evolve resistance to tetracycline they are treated with. Persister cells commonly exist in both planktonic and biofilm cells, with a higher frequency in the latter. Further transcriptomic analysis speculates that the overexpression of multidrug resistance genes and stress genes were mainly answered to sub lethal concentration of tetracycline in planktonic cells, and the lower metabolic levels after biofilm formation result in high resistance level of biofilm cells to tetracycline. Altogether, these data suggest that A. johnsonii can adjust its phenotype when grown as biofilm and change its metabolism under antibiotic stress, and provide implications for subsequent biofilm control.

Thermodynamics of co-translational folding and ribosome-nascent chain interactions

Proteins can begin the conformational search for their native structure in parallel with biosynthesis on the ribosome, in a process termed co-translational folding. In contrast to the reversible folding of isolated domains, as a nascent chain emerges from the ribosome exit tunnel during translation the free energy landscape it explores also evolves as a function of chain length. While this presents a substantially more complex measurement problem, this review will outline the progress that has been made recently in understanding, quantitatively, the process by which a nascent chain attains its full native stability, as well as the mechanisms through which interactions with the nearby ribosome surface can perturb or modulate this process.

Influence of Substitutions in the Binding Motif of Proline-Rich Antimicrobial Peptide ARV-1502 on 70S Ribosome Binding and Antimicrobial Activity

Proline-rich antimicrobial peptides (PrAMPs) are promising candidates to treat bacterial infections. The designer peptide ARV-1502 exhibits strong antimicrobial effects against Enterobacteriaceae both in vitro and in vivo. Since the inhibitory effects of ARV-1502 reported for the 70 kDa heat-shock protein DnaK do not fully explain the antimicrobial activity of its 176 substituted analogs, we further studied their effect on the bacterial 70S ribosome of Escherichia coli, a known target of PrAMPs. ARV-1502 analogues, substituted in positions 3, 4, and 8 to 12 (underlined) of the binding motif D³KPRPYLPRP¹² with aspartic acid, lysine, serine, phenylalanine or leucine, were tested in a competitive fluorescence polarization (FP) binding screening assay using 5(6)-carboxyfluorescein-labeled (Cf-) ARV-1502 and the 70S ribosome isolated from E. coli BW25113. While their effect on ribosomal protein expression was studied for green fluorescent protein (GFP) in a cell-free expression system (in vitro translation), the importance of known PrAMP transporters SbmA and MdtM was investigated using E. coli BW25113 and the corresponding knockout mutants. The dissociation constant (K_d) of 201 ± 16 nmol/L obtained for Cf-ARV-1502 suggests strong binding to the E. coli 70S ribosome. An inhibitory binding assay indicated that the binding site overlaps with those of other PrAMPs including Onc112 and pyrrhocoricin as well as the non-peptidic antibiotics erythromycin and chloramphenicol. All these drugs and drug candidates bind to the exit-tunnel of the 70S ribosome. Substitutions of the C-terminal fragment of the binding motif YLPRP reduced binding. At the same time, inhibition of GFP expression increased with net peptide charge. Interestingly, the MIC values of wild-type and ΔsbmA and ΔmdtM knockout mutants indicated that substitutions in the ribosomal binding motif altered also the bacterial uptake, which was generally improved by incorporation of hydrophobic residues. In conclusion, most substituted ARV-1502 analogs bound weaker to the 70S ribosome than ARV-1502 underlining the importance of the YLPRP binding motif. The weaker ribosomal binding correlated well with decreased antimicrobial activity in vitro. Substituted ARV-1502 analogs with a higher level of hydrophobicity or positive net charge improved the ribosome binding, inhibition of translation, and bacterial uptake.

Direct Measurement of Interhelical DNA Repulsion and Attraction by Quantitative Cross-Linking

To better understand the forces that mediate nucleic acid compaction in biology, we developed the disulfide cross-linking approach xHEED (X-linking of Helices to measure Electrostatic Effects at Distance) to measure the distance-dependent encounter frequency of two DNA helices in solution. Using xHEED, we determined the distance that the electrostatic potential extends from DNA helices, the dependence of this distance on ionic conditions, and the magnitude of repulsion when two helices approach one another. Across all conditions tested, the potential falls to that of the bulk solution within 15 Å of the major groove surface. For separations of ∼30 Å, we measured a repulsion of 1.8 kcal/mol in low monovalent ion concentration (30 mM Na+), with higher Na+ concentrations ameliorating this repulsion, and 2 M Na+ or 100 mM Mg2+ eliminating it. Strikingly, we found full screening at very low Co3+ concentrations and net attraction at higher concentrations, without the higher-order DNA condensation that typically complicates studies of helical attraction. Our measurements define the relevant distances for electrostatic interactions of nucleic-acid helices in biology and introduce a new method to propel further understanding of how these forces impact biological processes.

Nascent chain dynamics and ribosome interactions within folded ribosome-nascent chain complexes observed by NMR spectroscopy

The folding of many proteins can begin during biosynthesis on the ribosome and can be modulated by the ribosome itself. Such perturbations are generally believed to be mediated through interactions between the nascent chain and the ribosome surface, but despite recent progress in characterising interactions of unfolded states with the ribosome, and their impact on the initiation of co-translational folding, a complete quantitative analysis of interactions across both folded and unfolded states of a nascent chain has yet to be realised. Here we apply solution-state NMR spectroscopy to measure transverse proton relaxation rates for methyl groups in folded ribosome-nascent chain complexes of the FLN5 filamin domain. We observe substantial increases in relaxation rates for the nascent chain relative to the isolated domain, which can be related to changes in effective rotational correlation times using measurements of relaxation and cross-correlated relaxation in the isolated domain. Using this approach, we can identify interactions between the nascent chain and the ribosome surface, driven predominantly by electrostatics, and by measuring the change in these interactions as the subsequent FLN6 domain emerges, we may deduce their impact on the free energy landscapes associated with the co-translational folding process.

Optimal design of adaptively sampled NMR experiments for measurement of methyl group dynamics with application to a ribosome-nascent chain complex

NMR measurements of cross-correlated nuclear spin relaxation provide powerful probes of polypeptide dynamics and rotational diffusion, free from contributions due to chemical exchange or interactions with external spins. Here, we report on the development of a sensitivity-optimized pulse sequence for the analysis of the differential relaxation of transitions within isolated ¹³CH₃ spin systems, in order to characterise rotational diffusion and side chain order through the product S²τ_c. We describe the application of optimal design theory to implement a real-time 'on-the-fly' adaptive sampling scheme that maximizes the accuracy of the measured parameters. The increase in sensitivity obtained using this approach enables quantitative measurements of rotational diffusion within folded states of translationally-arrested ribosome-nascent chain complexes of the FLN5 filamin domain, and can be used to place strong limits on interactions between the domain and the ribosome surface.

Electrostatic features for nucleocapsid proteins of SARS-CoV and SARS-CoV-2

COVID-19 is increasingly affecting human health and global economy. Understanding the fundamental mechanisms of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is highly demanded to develop treatments for COVID-19. SARS-CoV and SARS-CoV-2 share 92.06% identity in their N protein RBDs' sequences, which results in very similar structures. However, the SARS-CoV-2 is more easily to spread. Utilizing multi-scale computational approaches, this work studied the fundamental mechanisms of the nucleocapsid (N) proteins of SARS-CoV and SARS-CoV-2, including their stabilities and binding strengths with RNAs at different pH values. Electrostatic potential on the surfaces of N proteins show that both the N proteins of SARS-CoV and SARS-CoV-2 have dominantly positive potential to attract RNAs. The binding forces between SARS-CoV N protein and RNAs at different distances are similar to that of SARS-CoV-2, both in directions and magnitudes. The electric filed lines between N proteins and RNAs are also similar for both SARS-CoV and SARS-CoV-2. The folding energy and binding energy dependence on pH revealed that the best environment for N proteins to perform their functions with RNAs is the weak acidic environment.

Spike Proteins of SARS-CoV and SARS-CoV-2 Utilize Different Mechanisms to Bind With Human ACE2

The ongoing outbreak of COVID-19 has been a serious threat to human health worldwide. The virus SARS-CoV-2 initiates its infection to the human body via the interaction of its spike (S) protein with the human Angiotensin-Converting Enzyme 2 (ACE2) of the host cells. Therefore, understanding the fundamental mechanisms of how SARS-CoV-2 S protein receptor binding domain (RBD) binds to ACE2 is highly demanded for developing treatments for COVID-19. Here we implemented multi-scale computational approaches to study the binding mechanisms of human ACE2 and S proteins of both SARS-CoV and SARS-CoV-2. Electrostatic features, including electrostatic potential, electric field lines, and electrostatic forces of SARS-CoV and SARS-CoV-2 were calculated and compared in detail. The results demonstrate that SARS-CoV and SARS-CoV-2 S proteins are both attractive to ACE2 by electrostatic forces even at different distances. However, the residues contributing to the electrostatic features are quite different due to the mutations between SARS-CoV S protein and SARS-CoV-2 S protein. Such differences are analyzed comprehensively. Compared to SARS-CoV, the SARS-CoV-2 binds with ACE2 using a more robust strategy: The electric field line related residues are distributed quite differently, which results in a more robust binding strategy of SARS-CoV-2. Also, SARS-CoV-2 has a higher electric field line density than that of SARS-CoV, which indicates stronger interaction between SARS-CoV-2 and ACE2, compared to that of SARS-CoV. Key residues involved in salt bridges and hydrogen bonds are identified in this study, which may help the future drug design against COVID-19.