A combined molecular dynamics and rapid kinetics approach to identify conserved three-dimensional communication networks in elongation factor Tu

The C-terminal helix of ribosomal P stalk recognizes a hydrophobic groove of elongation factor 2 in a novel fashion

Archaea and eukaryotes have ribosomal P stalks composed of anchor protein P0 and aP1 homodimers (archaea) or P1•P2 heterodimers (eukaryotes). These P stalks recruit translational GTPases to the GTPase-associated center in ribosomes to provide energy during translation. The C-terminus of the P stalk is known to selectively recognize GTPases. Here we investigated the interaction between the P stalk and elongation factor 2 by determining the structures of Pyrococcus horikoshii EF-2 (PhoEF-2) in the Apo-form, GDP-form, GMPPCP-form (GTP-form), and GMPPCP-form bound with 11 C-terminal residues of P1 (P1C11). Helical structured P1C11 binds to a hydrophobic groove between domain G and subdomain G′ of PhoEF-2, where is completely different from that of aEF-1α in terms of both position and sequence, implying that such interaction characteristic may be requested by how GTPases perform their functions on the ribosome. Combining PhoEF-2 P1-binding assays with a structural comparison of current PhoEF-2 structures and molecular dynamics model of a P1C11-bound GDP form, the conformational changes of the P1C11-binding groove in each form suggest that in response to the translation process, the groove has three states: closed, open, and release for recruiting and releasing GTPases.

Elongation Factor Tu's Nucleotide Binding Is Governed by a Thermodynamic Landscape Unique among Bacterial Translation Factors

Molecular switches such as GTPases are powerful devices turning "on" or "off" biomolecular processes at the core of critical biological pathways. To develop molecular switches de novo, an intimate understanding of how they function is required. Here we investigate the thermodynamic parameters that define the nucleotide-dependent switch mechanism of elongation factor (EF) Tu as a prototypical molecular switch. EF-Tu alternates between GTP- and GDP-bound conformations during its functional cycle, representing the "on" and "off" states, respectively. We report for the first time that the activation barriers for nucleotide association are the same for both nucleotides, suggesting a guanosine nucleoside or ribose-first mechanism for nucleotide association. Additionally, molecular dynamics (MD) simulations indicate that enthalpic stabilization of GDP binding compared to GTP binding originates in the backbone hydrogen bonding network of EF-Tu. In contrast, binding of GTP to EF-Tu is entropically driven by the liberation of bound water during the GDP- to GTP-bound transition. GDP binding to the apo conformation of EF-Tu is both enthalpically and entropically favored, a feature unique among translational GTPases. This indicates that the apo conformation does not resemble the GDP-bound state. Finally, we show that antibiotics and single amino acid substitutions can be used to target specific structural elements in EF-Tu to redesign the thermodynamic landscape. These findings demonstrate how, through evolution, EF-Tu has fine-tuned the structural and dynamic features that define nucleotide binding, providing insight into how altering these properties could be exploited for protein engineering.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

Trade-off and flexibility in the dynamic regulation of the cullin-RING ubiquitin ligase repertoire

Cullin-RING ubiquitin ligases (CRLs) catalyze the ubiquitylation of substrates many of which are degraded by the 26S proteasome. Their modular architecture enables recognition of numerous substrates via exchangeable substrate receptors that competitively bind to a cullin scaffold with high affinity. Due to the plasticity of these interactions there is ongoing uncertainty how cells maintain a flexible CRL repertoire in view of changing substrate loads. Based on a series of in vivo and in vitro studies, different groups proposed that the exchange of substrate receptors is mediated by a protein exchange factor named Cand1. Here, we have performed mathematical modeling to provide a quantitative underpinning of this hypothesis. First we show that the exchange activity of Cand1 necessarily leads to a trade-off between high ligase activity and fast receptor exchange. Supported by measurements we argue that this trade-off yields an optimal Cand1 concentration in cells where the time scale for substrate degradation becomes minimal. In a second step we show through simulations that (i) substrates bias the CRL repertoire leading to preferential assembly of ligases for which substrates are available and (ii) differences in binding affinities or substrate receptor abundances create a temporal hierarchy for the degradation of substrates. Finally, we compare the Cand1-mediated exchange cycle with an alternative architecture lacking Cand1 which indicates superiority of a system with exchange factor if substrate receptors bind substrates and the cullin scaffold in a random order. Together, our results provide general constraints for the operating regimes of molecular exchange systems and suggest that Cand1 endows the CRL network with the properties of an “on demand” system allowing cells to dynamically adjust their CRL repertoire to fluctuating substrate abundances.

Cullin-RING ubiquitin ligases (CRLs) are multisubunit protein complexes where exchangeable substrate receptors (SRs) assemble on a cullin scaffold to mediate ubiquitylation and subsequent degradation of a large variety of substrates. In humans there are hundreds of different CRLs having potentially thousands of substrates. Due to the high affinity of cullin-SR interactions, it has long been a mystery how cells would maintain flexibility to sample the entire SR repertoire in order to match fluctuating substrate loads. Recent experiments indicate that the exchange of different SRs is mediated by a novel protein exchange factor (Cand1). However, the proposed biochemical function of Cand1 as a promoter of CRL activity remained difficult to reconcile with previous reports of Cand1 acting as an inhibitor of CRL activity in vitro. Here we show that these two findings are not contradictory, but that the exchange activity of Cand1 necessarily leads to a trade-off between high ligase activity and fast receptor exchange which leads us to predict an optimal Cand1 concentration and a temporal hierarchy for substrate degradation. Our results support the view that Cand1 endows the CRL network with the flexibility of an “on demand” system where relative CRL abundances are dictated by substrate availability.

Tetracycline does not directly inhibit the function of bacterial elongation factor Tu

Understanding the molecular mechanism of antibiotics that are currently in use is important for the development of new antimicrobials. The tetracyclines, discovered in the 1940s, are a well-established class of antibiotics that still have a role in treating microbial infections in humans. It is generally accepted that the main target of their action is the ribosome. The estimated affinity for tetracycline binding to the ribosome is relatively low compared to the actual potency of the drug in vivo. Therefore, additional inhibitory effects of tetracycline on the translation machinery have been discussed. Structural evidence suggests that tetracycline inhibits the function of the essential bacterial GTPase Elongation Factor (EF)-Tu through interaction with the bound nucleotide. Based on this, tetracycline has been predicted to impede the nucleotide-binding properties of EF-Tu. However, detailed kinetic studies addressing the effect of tetracycline on nucleotide binding have been prevented by the fluorescence properties of the antibiotic. Here, we report a fluorescence-based kinetic assay that minimizes the effect of tetracycline autofluorescence, enabling the detailed kinetic analysis of the nucleotide-binding properties of Escherichia coli EF-Tu. Furthermore, using physiologically relevant conditions, we demonstrate that tetracycline does not affect EF-Tu’s intrinsic or ribosome-stimulated GTPase activity, nor the stability of the EF-Tu•GTP•Phe-tRNA^Phe complex. We therefore provide clear evidence that tetracycline does not directly impede the function of EF-Tu.

The C-terminal Helix of Pseudomonas aeruginosa Elongation Factor Ts Tunes EF-Tu Dynamics to Modulate Nucleotide Exchange

Little is known about the conservation of critical kinetic parameters and the mechanistic strategies of elongation factor (EF) Ts-catalyzed nucleotide exchange in EF-Tu in bacteria and particularly in clinically relevant pathogens. EF-Tu from the clinically relevant pathogen Pseudomonas aeruginosa shares over 84% sequence identity with the corresponding elongation factor from Escherichia coli Interestingly, the functionally closely linked EF-Ts only shares 55% sequence identity. To identify any differences in the nucleotide binding properties, as well as in the EF-Ts-mediated nucleotide exchange reaction, we performed a comparative rapid kinetics and mutagenesis analysis of the nucleotide exchange mechanism for both the E. coli and P. aeruginosa systems, identifying helix 13 of EF-Ts as a previously unnoticed regulatory element in the nucleotide exchange mechanism with species-specific elements. Our findings support the base side-first entry of the nucleotide into the binding pocket of the EF-Tu·EF-Ts binary complex, followed by displacement of helix 13 and rapid binding of the phosphate side of the nucleotide, ultimately leading to the release of EF-Ts.

Structures, Properties, and Dynamics of Intermediates in eEF2-Diphthamide Biosynthesis

The eukaryotic translation Elongation Factor 2 (eEF2) is an essential enzyme in protein synthesis. eEF2 contains a unique modification of a histidine (His699 in yeast; HIS) into diphthamide (DTA), obtained via 3-amino-3-carboxypropyl (ACP) and diphthine (DTI) intermediates in the biosynthetic pathway. This essential and unique modification is also vulnerable, in that it can be efficiently targeted by NAD(+)-dependent ADP-ribosylase toxins, such as diphtheria toxin (DT). However, none of the intermediates in the biosynthesis path is equally vulnerable against the toxins. This study aims to address the different susceptibility of DTA and its precursors against bacterial toxins. We have herein undertaken a detailed in silico study of the structural features and dynamic motion of different His699 intermediates along the diphthamide synthesis pathway (HIS, ACP, DTI, DTA). The study points out that DTA forms a strong hydrogen bond with an asparagine which might explain the ADP-ribosylation mechanism caused by the diphtheria toxin (DT). Finally, in silico mutagenesis studies were performed on the DTA modified protein, in order to hamper the formation of such a hydrogen bond. The results indicate that the mutant structure might in fact be less susceptible to attack by DT and thereby behave similarly to DTI in this respect.

A conserved P-loop anchor limits the structural dynamics that mediate nucleotide dissociation in EF-Tu

The phosphate-binding loop (P-loop) is a conserved sequence motif found in mononucleotide-binding proteins. Little is known about the structural dynamics of this region and its contribution to the observed nucleotide binding properties. Understanding the underlying design principles is of great interest for biomolecular engineering applications. We have used rapid-kinetics measurements in vitro and molecular dynamics (MD) simulations in silico to investigate the relationship between GTP-binding properties and P-loop structural dynamics in the universally conserved Elongation Factor (EF) Tu. Analysis of wild type EF-Tu and variants with substitutions at positions in or adjacent to the P-loop revealed a correlation between P-loop flexibility and the entropy of activation for GTP dissociation. The same variants demonstrate more backbone flexibility in two N-terminal amino acids of the P-loop during force-induced EF-Tu · GTP dissociation in Steered Molecular Dynamics simulations. Amino acids Gly18 and His19 are involved in stabilizing the P-loop backbone via interactions with the adjacent helix C. We propose that these P-loop/helix C interactions function as a conserved P-loop anchoring module and identify the presence of P-loop anchors within several GTPases and ATPases suggesting their evolutionary conservation.

An arginine-aspartate network in the active site of bacterial TruB is critical for catalyzing pseudouridine formation

Pseudouridine synthases introduce the most common RNA modification and likely use the same catalytic mechanism. Besides a catalytic aspartate residue, the contributions of other residues for catalysis of pseudouridine formation are poorly understood. Here, we have tested the role of a conserved basic residue in the active site for catalysis using the bacterial pseudouridine synthase TruB targeting U55 in tRNAs. Substitution of arginine 181 with lysine results in a 2500-fold reduction of TruB's catalytic rate without affecting tRNA binding. Furthermore, we analyzed the function of a second-shell aspartate residue (D90) that is conserved in all TruB enzymes and interacts with C56 of tRNA. Site-directed mutagenesis, biochemical and kinetic studies reveal that this residue is not critical for substrate binding but influences catalysis significantly as replacement of D90 with glutamate or asparagine reduces the catalytic rate 30- and 50-fold, respectively. In agreement with molecular dynamics simulations of TruB wild type and TruB D90N, we propose an electrostatic network composed of the catalytic aspartate (D48), R181 and D90 that is important for catalysis by fine-tuning the D48-R181 interaction. Conserved, negatively charged residues similar to D90 are found in a number of pseudouridine synthases, suggesting that this might be a general mechanism.

Elongation factor Ts directly facilitates the formation and disassembly of the Escherichia coli elongation factor Tu·GTP·aminoacyl-tRNA ternary complex

BACKGROUND

Aminoacyl-tRNA (aa-tRNA) enters the ribosome in a ternary complex with the G-protein elongation factor Tu (EF-Tu) and GTP.

RESULTS

EF-Tu·GTP·aa-tRNA ternary complex formation and decay rates are accelerated in the presence of the nucleotide exchange factor elongation factor Ts (EF-Ts).

CONCLUSION

EF-Ts directly facilitates the formation and disassociation of ternary complex.

SIGNIFICANCE

This system demonstrates a novel function of EF-Ts. Aminoacyl-tRNA enters the translating ribosome in a ternary complex with elongation factor Tu (EF-Tu) and GTP. Here, we describe bulk steady state and pre-steady state fluorescence methods that enabled us to quantitatively explore the kinetic features of Escherichia coli ternary complex formation and decay. The data obtained suggest that both processes are controlled by a nucleotide-dependent, rate-determining conformational change in EF-Tu. Unexpectedly, we found that this conformational change is accelerated by elongation factor Ts (EF-Ts), the guanosine nucleotide exchange factor for EF-Tu. Notably, EF-Ts attenuates the affinity of EF-Tu for GTP and destabilizes ternary complex in the presence of non-hydrolyzable GTP analogs. These results suggest that EF-Ts serves an unanticipated role in the cell of actively regulating the abundance and stability of ternary complex in a manner that contributes to rapid and faithful protein synthesis.

A computational investigation on the connection between dynamics properties of ribosomal proteins and ribosome assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins.

Ribosomes are complex cellular machines that synthesize new proteins in the cell. The accurate and efficient assembly of ribosomal proteins (r-proteins) and ribosomal RNA (rRNA) to form a functional ribosome is important for cell growth, metabolic reactions, and other cellular processes. Additionally, some antibacterial drugs are believed to target the bacterial ribosome during its construction. Hence, ribosomal assembly has been an active research topic for many years because understanding the assembly mechanisms can provide insight into protein/RNA recognitions important in many other cellular processes, as well as optimize the development of antibacterial therapeutics. Experimental studies thus far have provided still limited understanding about the assembly process. To further understand the assembly process, we have computationally studied the dynamic properties that r-proteins exhibit during assembly and the relationship between dynamics, physical properties, and binding propensity. We observe significant charged interactions between r-proteins and rRNA. We also detect a strong correlation between contact residues and their dynamic mobilities. Protein residues contacting with rRNA are observed to be more mobile in comparison with other residues. We also relate the location of the r-protein in the fully assembled ribosome to its susceptibility for large conformational changes prior to binding.

Characterization of the emergent properties of a synthetic quasi-cellular system

BACKGROUND

The process of solutes entrapment during liposomes formation is interesting for the investigation of the relationship between the formation of compartments and the distribution of molecules inside them; a relevant issue in the studies of the origin of life. Theoretically, when no interactions are supposed among the chemical species to be entrapped, the entrapment is described by a standard Poisson process. But very recent experimental findings show that, for small liposomes (100 nm diameter), the distribution of entrapped molecules is best described by a power-law function. This is of a great importance, as the two random processes give rise to two completely different scenarios. Here we present an in silico stochastic simulation of the encapsulation of a cell-free molecular translation system (the PURE system), obtained following two different entrapment models: a pure Poisson process, and a power-law. The protein synthesis inside the liposomes has been studied in both cases, with the aim to highlight experimental observables that could be measured to assess which model gives a better representation of the real process.

RESULTS

Firstly, a minimal model for in vitro protein synthesis, based on the PURE system molecular composition, has been formalized. Then, we have designed a reliable experimental simulation where stochastic factors affect the reaction course inside the compartment. To this end, 24 solutes, which represent the PURE system components, have been stochastically distributed among vesicles by following either a Poisson or a power-law distribution. The course of the protein synthesis within each vesicle has been consequently calculated, as a function of vesicle size. Our study can predict translation yield in a population of small liposomes down to the attoliter (10(-18) L) range. Our results show that the efficiency of protein synthesis peaks at approximately 3 · 10(-16) L (840 nm diam.) with a Poisson distribution of solutes, while a relative optimum is found at around 10(-17) L (275 nm diam.) for the power-law statistics.

CONCLUSIONS

Our simulation clearly shows that the wet-lab measurement of an effective protein synthesis at smaller volumes than 10(-17) L would rule out, according to our models, a Poisson distribution of solutes.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.