Converting structural information into an allosteric-energy-based picture for elongation factor Tu activation by the ribosome

RAC1P29S is a spontaneously activating cancer-associated GTPase

RAC1 is a small, Ras-related GTPase that was recently reported to harbor a recurrent UV-induced signature mutation in melanoma, resulting in substitution of P29 to serine (RAC1(P29S)), ranking this the third most frequently occurring gain-of-function mutation in melanoma. Although the Ras family GTPases are mutated in about 30% of all cancers, mutations in the Rho family GTPases have rarely been observed. In this study, we demonstrate that unlike oncogenic Ras proteins, which are primarily activated by mutations that eliminate GTPase activity, the activated melanoma RAC1(P29S) protein maintains intrinsic GTP hydrolysis and is spontaneously activated by substantially increased inherent GDP/GTP nucleotide exchange. Determination and comparison of crystal structures for activated RAC1 GTPases suggest that RAC1(F28L)--a known spontaneously activated RAC1 mutant--and RAC1(P29S) are self-activated in distinct fashions. Moreover, the mechanism of RAC1(P29S) and RAC1(F28L) activation differs from the common oncogenic mutations found in Ras-like GTPases that abrogate GTP hydrolysis. The melanoma RAC1(P29S) gain-of-function point mutation therefore represents a previously undescribed class of cancer-related GTPase activity.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Hydrogen tunneling links protein dynamics to enzyme catalysis

The relationship between protein dynamics and function is a subject of considerable contemporary interest. Although protein motions are frequently observed during ligand binding and release steps, the contribution of protein motions to the catalysis of bond making/breaking processes is more difficult to probe and verify. Here, we show how the quantum mechanical hydrogen tunneling associated with enzymatic C-H bond cleavage provides a unique window into the necessity of protein dynamics for achieving optimal catalysis. Experimental findings support a hierarchy of thermodynamically equilibrated motions that control the H-donor and -acceptor distance and active-site electrostatics, creating an ensemble of conformations suitable for H-tunneling. A possible extension of this view to methyl transfer and other catalyzed reactions is also presented. The impact of understanding these dynamics on the conceptual framework for enzyme activity, inhibitor/drug design, and biomimetic catalyst design is likely to be substantial.

Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping

The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

Bridging the gap between ribosome structure and biochemistry by mechanistic computations

The wealth of structural and biochemical data now available for protein synthesis on the ribosome presents major new challenges for computational biochemistry. Apart from technical difficulties in modeling ribosome systems, the complexity of the overall translation cycle with a multitude of different kinetic steps presents a formidable problem for computational efforts where we have only seen the beginning. However, a range of methodologies including molecular dynamics simulations, free energy calculations, molecular docking and quantum chemical approaches have already been put to work with promising results. In particular, the combined efforts of structural biology, biochemistry, kinetics and computational modeling can lead towards a quantitative structure-based description of translation.

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.

A master switch couples Mg²⁺-assisted catalysis to domain motion in B. stearothermophilus tryptophanyl-tRNA Synthetase

We demonstrate how tryptophanyl-tRNA synthetase uses conformation-dependent Mg(2+) activation to couple catalysis of tryptophan activation to specific, functional domain movements. Rate acceleration by Mg(2+) requires ∼-6.0 kcal/mol in protein⋅Mg(2+) interaction energy, none of which arises from the active site. A highly cooperative interaction between Mg(2+) and four residues from a remote, conserved motif that mediates the shear of domain movement (1) destabilizes the pretransition state conformation, thereby (2) inducing the Mg(2+) to stabilize the transition state for k(cat) by ∼-5.0 kcal/mol. Cooperative, long-range conformational effects on the metal therefore convert an inactive Mg(2+) coordination into one that can stabilize the transition state if, and only if, domain motion occurs. Transient, conformation-dependent Mg(2+) activation, analogous to the escapement in mechanical clocks, explains vectorial coupling.

The busiest of all ribosomal assistants: elongation factor Tu

During translation, the nucleic acid language employed by genes is translated into the amino acid language used by proteins. The translator is the ribosome, while the dictionary employed is known as the genetic code. The genetic information is presented to the ribosome in the form of a mRNA, and tRNAs connect the two languages. Translation takes place in three steps: initiation, elongation, and termination. After a protein has been synthesized, the components of the translation apparatus are recycled. During each phase of translation, the ribosome collaborates with specific translation factors, which secure a proper balance between speed and fidelity. Notably, initiation, termination, and ribosomal recycling occur only once per protein produced during normal translation, while the elongation step is repeated a large number of times, corresponding to the number of amino acids constituting the protein of interest. In bacteria, elongation factor Tu plays a central role during the selection of the correct amino acids throughout the elongation phase of translation. Elongation factor Tu is the main subject of this review.

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.

Examining the promiscuous phosphatase activity of Pseudomonas aeruginosa arylsulfatase: a comparison to analogous phosphatases

Pseudomonas aeruginosa arylsulfatase (PAS) is a bacterial sulfatase capable of hydrolyzing a range of sulfate esters. Recently, it has been demonstrated to also show very high proficiency for phosphate ester hydrolysis. Such proficient catalytic promiscuity is significant, as promiscuity has been suggested to play an important role in enzyme evolution. Additionally, a comparative study of the hydrolyses of the p-nitrophenyl phosphate and sulfate monoesters in aqueous solution has demonstrated that despite superficial similarities, the two reactions proceed through markedly different transition states with very different solvation effects, indicating that the requirements for the efficient catalysis of the two reactions by an enzyme will also be very different (and yet they are both catalyzed by the same active site). This work explores the promiscuous phosphomonoesterase activity of PAS. Specifically, we have investigated the identity of the most likely base for the initial activation of the unusual formylglycine hydrate nucleophile (which is common to many sulfatases), and demonstrate that a concerted substrate-as-base mechanism is fully consistent with the experimentally observed data. This is very similar to other related systems, and suggests that, as far as the phosphomonoesterase activity of PAS is concerned, the sulfatase behaves like a "classical" phosphatase, despite the fact that such a mechanism is unlikely to be available to the native substrate (based on pK(a) considerations and studies of model systems). Understanding such catalytic versatility can be used to design novel artificial enzymes that are far more proficient than the current generation of designer enzymes.

Quality control of mRNA decoding on the bacterial ribosome

The ribosome is a major player in providing accurate gene expression in the cell. The fidelity of substrate selection is tightly controlled throughout the translation process, including the initiation, elongation, and termination phases. Although each phase of translation involves different players, that is, translation factors and tRNAs, the general principles of selection appear surprisingly similar for very different substrates. At essentially every step of translation, differences in complex stabilities as well as induced fit are sources of selectivity. A view starts to emerge of how the ribosome uses local and global conformational switches to govern induced-fit mechanisms that ensure fidelity. This review describes the mechanisms of tRNA and mRNA selection at all phases of protein synthesis in bacteria.

Electrostatic origin of the mechanochemical rotary mechanism and the catalytic dwell of F1-ATPase

Understanding the nature of energy transduction in life processes requires a quantitative description of the energetics of the conversion of ATP to ADP by ATPases. Previous attempts to do so have provided an interesting insight but could not account for the rotary mechanism by a nonphenomenological structure/energy description. In particular it has been very challenging to account for the observations of the 80° and 40° rotational substates, without any prior information about such states in the simulation procedure. Here we use a coarse-grained model of F1-ATPase and generate, without the adjustment of phenomenological parameters, a structure-based free energy landscape that reproduces the energetics of the mechanochemical process. It is found that the landscape along the relevant rotary path is determined by the electrostatic free energy and not by steric effects. Furthermore, the generated surface and the corresponding Langevin dynamics simulations identify a hidden conformational barrier that provides a new fundamental interpretation of the catalytic dwell and illuminate the nature of the energy conversion process.