The human RNA polymerase I structure reveals an HMG-like docking domain specific to metazoans

We characterize the human RNA polymerase I by evolutionary biochemistry and cryo-EM revealing a built-in structural domain that apparently serves as transcription factor–binding platform in metazoans.

Transcription of the ribosomal RNA precursor by RNA polymerase (Pol) I is a major determinant of cellular growth, and dysregulation is observed in many cancer types. Here, we present the purification of human Pol I from cells carrying a genomic GFP fusion on the largest subunit allowing the structural and functional analysis of the enzyme across species. In contrast to yeast, human Pol I carries a single-subunit stalk, and in vitro transcription indicates a reduced proofreading activity. Determination of the human Pol I cryo-EM reconstruction in a close-to-native state rationalizes the effects of disease-associated mutations and uncovers an additional domain that is built into the sequence of Pol I subunit RPA1. This “dock II” domain resembles a truncated HMG box incapable of DNA binding which may serve as a downstream transcription factor–binding platform in metazoans. Biochemical analysis, in situ modelling, and ChIP data indicate that Topoisomerase 2a can be recruited to Pol I via the domain and cooperates with the HMG box domain–containing factor UBF. These adaptations of the metazoan Pol I transcription system may allow efficient release of positive DNA supercoils accumulating downstream of the transcription bubble.

Allosteric control of Ubp6 and the proteasome via a bidirectional switch

The proteasome recognizes ubiquitinated proteins and can also edit ubiquitin marks, allowing substrates to be rejected based on ubiquitin chain topology. In yeast, editing is mediated by deubiquitinating enzyme Ubp6. The proteasome activates Ubp6, whereas Ubp6 inhibits the proteasome through deubiquitination and a noncatalytic effect. Here, we report cryo-EM structures of the proteasome bound to Ubp6, based on which we identify mutants in Ubp6 and proteasome subunit Rpt1 that abrogate Ubp6 activation. The Ubp6 mutations define a conserved region that we term the ILR element. The ILR is found within the BL1 loop, which obstructs the catalytic groove in free Ubp6. Rpt1-ILR interaction opens the groove by rearranging not only BL1 but also a previously undescribed network of three interconnected active-site-blocking loops. Ubp6 activation and noncatalytic proteasome inhibition are linked in that they are eliminated by the same mutations. Ubp6 and ubiquitin together drive proteasomes into a unique conformation associated with proteasome inhibition. Thus, a multicomponent allosteric switch exerts simultaneous control over both Ubp6 and the proteasome.

Structural basis for VIPP1 oligomerization and maintenance of thylakoid membrane integrity

Vesicle-inducing protein in plastids 1 (VIPP1) is essential for the biogenesis and maintenance of thylakoid membranes, which transform light into life. However, it is unknown how VIPP1 performs its vital membrane-remodeling functions. Here, we use cryo-electron microscopy to determine structures of cyanobacterial VIPP1 rings, revealing how VIPP1 monomers flex and interweave to form basket-like assemblies of different symmetries. Three VIPP1 monomers together coordinate a non-canonical nucleotide binding pocket on one end of the ring. Inside the ring's lumen, amphipathic helices from each monomer align to form large hydrophobic columns, enabling VIPP1 to bind and curve membranes. In vivo mutations in these hydrophobic surfaces cause extreme thylakoid swelling under high light, indicating an essential role of VIPP1 lipid binding in resisting stress-induced damage. Using cryo-correlative light and electron microscopy (cryo-CLEM), we observe oligomeric VIPP1 coats encapsulating membrane tubules within the Chlamydomonas chloroplast. Our work provides a structural foundation for understanding how VIPP1 directs thylakoid biogenesis and maintenance.

Time-resolved spectroscopic and electrophysiological data reveal insights in the gating mechanism of anion channelrhodopsin

Channelrhodopsins are widely used in optogenetic applications. High photocurrents and low current inactivation levels are desirable. Two parallel photocycles evoked by different retinal conformations cause cation-conducting channelrhodopsin-2 (CrChR2) inactivation: one with efficient conductivity; one with low conductivity. Given the longer half-life of the low conducting photocycle intermediates, which accumulate under continuous illumination, resulting in a largely reduced photocurrent. Here, we demonstrate that for channelrhodopsin-1 of the cryptophyte Guillardia theta (GtACR1), the highly conducting C = N-anti-photocycle was the sole operating cycle using time-resolved step-scan FTIR spectroscopy. The correlation between our spectroscopic measurements and previously reported electrophysiological data provides insights into molecular gating mechanisms and their role in the characteristic high photocurrents. The mechanistic importance of the central constriction site amino acid Glu-68 is also shown. We propose that canceling out the poorly conducting photocycle avoids the inactivation observed in CrChR2, and anticipate that this discovery will advance the development of optimized optogenetic tools.

The Ras dimer structure

Oncogenic mutated Ras is a key player in cancer, but despite intense and expensive approaches its catalytic center seems undruggable. The Ras dimer interface is a possible alternative drug target. Dimerization at the membrane affects cell growth signal transduction. In vivo studies indicate that preventing dimerization of oncogenic mutated Ras inhibits uncontrolled cell growth. Conventional computational drug-screening approaches require a precise atomic dimer model as input to successfully access drug candidates. However, the proposed dimer structural models are controversial. Here, we provide a clear-cut experimentally validated N-Ras dimer structural model. We incorporated unnatural amino acids into Ras to enable the binding of labels at multiple positions via click chemistry. This labeling allowed the determination of multiple distances of the membrane-bound Ras-dimer measured by fluorescence and electron paramagnetic resonance spectroscopy. In combination with protein-protein docking and biomolecular simulations, we identified key residues for dimerization. Site-directed mutations of these residues prevent dimer formation in our experiments, proving our dimer model to be correct. The presented dimer structure enables computational drug-screening studies exploiting the Ras dimer interface as an alternative drug target.

Structural insights into photosystem II assembly

Biogenesis of photosystem II (PSII), nature's water-splitting catalyst, is assisted by auxiliary proteins that form transient complexes with PSII components to facilitate stepwise assembly events. Using cryo-electron microscopy, we solved the structure of such a PSII assembly intermediate from Thermosynechococcus elongatus at 2.94 Å resolution. It contains three assembly factors (Psb27, Psb28 and Psb34) and provides detailed insights into their molecular function. Binding of Psb28 induces large conformational changes at the PSII acceptor side, which distort the binding pocket of the mobile quinone (Q_B) and replace the bicarbonate ligand of non-haem iron with glutamate, a structural motif found in reaction centres of non-oxygenic photosynthetic bacteria. These results reveal mechanisms that protect PSII from damage during biogenesis until water splitting is activated. Our structure further demonstrates how the PSII active site is prepared for the incorporation of the Mn₄CaO₅ cluster, which performs the unique water-splitting reaction.

Expanded Coverage of the 26S Proteasome Conformational Landscape Reveals Mechanisms of Peptidase Gating

The proteasome is the central protease for intracellular protein breakdown. Coordinated binding and hydrolysis of ATP by the six proteasomal ATPase subunits induces conformational changes that drive the unfolding and translocation of substrates into the proteolytic 20S core particle for degradation. Here, we combine genetic and biochemical approaches with cryo-electron microscopy and integrative modeling to dissect the relationship between individual nucleotide binding events and proteasome conformational dynamics. We demonstrate unique impacts of ATP binding by individual ATPases on the proteasome conformational distribution and report two conformational states of the proteasome suggestive of a rotary ATP hydrolysis mechanism. These structures, coupled with functional analyses, reveal key roles for the ATPases Rpt1 and Rpt6 in gating substrate entry into the core particle. This deepened knowledge of proteasome conformational dynamics reveals key elements of intersubunit communication within the proteasome and clarifies the regulation of substrate entry into the proteolytic chamber.

Lamprey Parapinopsin ("UVLamP"): a Bistable UV-Sensitive Optogenetic Switch for Ultrafast Control of GPCR Pathways

Optogenetics uses light‐sensitive proteins, so‐called optogenetic tools, for highly precise spatiotemporal control of cellular states and signals. The major limitations of such tools include the overlap of excitation spectra, phototoxicity, and lack of sensitivity. The protein characterized in this study, the Japanese lamprey parapinopsin, which we named UVLamP, is a promising optogenetic tool to overcome these limitations. Using a hybrid strategy combining molecular, cellular, electrophysiological, and computational methods we elucidated a structural model of the dark state and probed the optogenetic potential of UVLamP. Interestingly, it is the first described bistable vertebrate opsin that has a charged amino acid interacting with the Schiff base in the dark state, that has no relevance for its photoreaction. UVLamP is a bistable UV‐sensitive opsin that allows for precise and sustained optogenetic control of G protein‐coupled receptor (GPCR) pathways and can be switched on, but more importantly also off within milliseconds via lowintensity short light pulses. UVLamP exhibits an extremely narrow excitation spectrum in the UV range allowing for sustained activation of the G_i/o pathway with a millisecond UV light pulse. Its sustained pathway activation can be switched off, surprisingly also with a millisecond blue light pulse, minimizing phototoxicity. Thus, UVLamP serves as a minimally invasive, narrow‐bandwidth probe for controlling the G_i/o pathway, allowing for combinatorial use with multiple optogenetic tools or sensors. Because UVLamP activated G_i/o signals are generally inhibitory and decrease cellular activity, it has tremendous potential for health‐related applications such as relieving pain, blocking seizures, and delaying neurodegeneration.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Design of an Ultrafast G Protein Switch Based on a Mouse Melanopsin Variant

The primary goal of optogenetics is the light-controlled noninvasive and specific manipulation of various cellular processes. Herein, we present a hybrid strategy for targeted protein engineering combining computational techniques with electrophysiological and UV/visible spectroscopic experiments. We validated our concept for channelrhodopsin-2 and applied it to modify the less-well-studied vertebrate opsin melanopsin. Melanopsin is a promising optogenetic tool that functions as a selective molecular light switch for G protein-coupled receptor pathways. Thus, we constructed a model of the melanopsin G_q protein complex and predicted an absorption maximum shift of the Y211F variant. This variant displays a narrow blue-shifted action spectrum and twofold faster deactivation kinetics compared to wild-type melanopsin on G protein-coupled inward rectifying K⁺ (GIRK) channels in HEK293 cells. Furthermore, we verified the in vivo activity and optogenetic potential for the variant in mice. Thus, we propose that our developed concept will be generally applicable to designing optogenetic tools.

Monitoring transient events in infrared spectra using local mode analysis

Time-resolved Fourier transformed infrared (FTIR) spectroscopy of chemical reactions is highly sensitive to minimal spatiotemporal changes. Structural features are decoded and represented in a comprehensible manner by combining FTIR spectroscopy with biomolecular simulations. Local mode analysis (LMA) is a tool to connect molecular motion based on a quantum mechanics simulation with infrared (IR) spectral features and vice versa. Here, we present the python-based software tool of LMA and demonstrate the novel feature of LMA to extract transient structural details and identify the related IR spectra at the case example of malonaldehyde (MA). Deuterated MA exists in two almost equally populated tautomeric states separated by a low barrier for proton transfer so IR spectra represent a mixture of both states. By state-dependent LMA, we obtain pure spectra for each tautomeric state occurring within the quantum mechanics trajectory. By time-resolved LMA, we obtain a clear view of the transition between states in the spectrum. Through local mode decomposition and the band-pass filter, marker bands for each state are identified. Thus, LMA is beneficial to analyze the experimental spectra based on a mixture of states by determining the individual contributions to the spectrum and motion of each state.

In Situ Structure of Neuronal C9orf72 Poly-GA Aggregates Reveals Proteasome Recruitment

Protein aggregation and dysfunction of the ubiquitin-proteasome system are hallmarks of many neurodegenerative diseases. Here, we address the elusive link between these phenomena by employing cryo-electron tomography to dissect the molecular architecture of protein aggregates within intact neurons at high resolution. We focus on the poly-Gly-Ala (poly-GA) aggregates resulting from aberrant translation of an expanded GGGGCC repeat in C9orf72, the most common genetic cause of amyotrophic lateral sclerosis and frontotemporal dementia. We find that poly-GA aggregates consist of densely packed twisted ribbons that recruit numerous 26S proteasome complexes, while other macromolecules are largely excluded. Proximity to poly-GA ribbons stabilizes a transient substrate-processing conformation of the 26S proteasome, suggesting stalled degradation. Thus, poly-GA aggregates may compromise neuronal proteostasis by driving the accumulation and functional impairment of a large fraction of cellular proteasomes.

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Recognition of Poly-Ubiquitins by the Proteasome through Protein Refolding Guided by Electrostatic and Hydrophobic Interactions

Specificity of protein degradation by cellular proteasomes comes from tetra-ubiquitin recognition. We carry out molecular dynamics simulations to characterize how the ubiquitin receptor Rpn10 recognizes in the 26S proteasome K48-linked tetra-ubiquitin. In the binding pose, ubiquitin and Rpn10 interact primarily through hydrophobic patches. However, K48-linked tetra-ubiquitin mostly assumes a closed form in solution prior to binding, in which its hydrophobic patches are not exposed to solvent. Likewise, the hydrophobic ubiquitin interacting motifs (UIMs) of Rpn10 are mostly protected prior to binding. As a result, ubiquitin recognition in the proteasome requires refolding of both K48-linked tetra-ubiquitin and Rpn10. Simulations suggest that conserved complementary electrostatic patterns of Rpn10 and ubiquitins guide protein association (stage 1 in the recognition process), which induces refolding (stage 2), and then facilitates formation of hydrophobic contacts (stage 3). The simulations also explain why Rpn10 has a higher affinity for K48-linked tetra-ubiquitin than for mono-ubiquitin and K48-linked di- and tri-ubiquitins. Simulation results expand on the current view that the flexible arm of Rpn10 acts as an extended fragment of α-helices and flexible coils in the recognition process.

Catalysis of GTP hydrolysis by small GTPases at atomic detail by integration of X-ray crystallography, experimental, and theoretical IR spectroscopy

Small GTPases regulate key processes in cells. Malfunction of their GTPase reaction by mutations is involved in severe diseases. Here, we compare the GTPase reaction of the slower hydrolyzing GTPase Ran with Ras. By combination of time-resolved FTIR difference spectroscopy and QM/MM simulations we elucidate that the Mg(2+) coordination by the phosphate groups, which varies largely among the x-ray structures, is the same for Ran and Ras. A new x-ray structure of a Ran·RanBD1 complex with improved resolution confirmed this finding and revealed a general problem with the refinement of Mg(2+) in GTPases. The Mg(2+) coordination is not responsible for the much slower GTPase reaction of Ran. Instead, the location of the Tyr-39 side chain of Ran between the γ-phosphate and Gln-69 prevents the optimal positioning of the attacking water molecule by the Gln-69 relative to the γ-phosphate. This is confirmed in the RanY39A·RanBD1 crystal structure. The QM/MM simulations provide IR spectra of the catalytic center, which agree very nicely with the experimental ones. The combination of both methods can correlate spectra with structure at atomic detail. For example the FTIR difference spectra of RasA18T and RanT25A mutants show that spectral differences are mainly due to the hydrogen bond of Thr-25 to the α-phosphate in Ran. By integration of x-ray structure analysis, experimental, and theoretical IR spectroscopy the catalytic center of the x-ray structural models are further refined to sub-Å resolution, allowing an improved understanding of catalysis.

Detailed structure of the H2PO4(-)-guanosine diphosphate intermediate in Ras-GAP decoded from FTIR experiments by biomolecular simulations

Essential biochemical processes such as signal transduction, energy conversion, or substrate conversion depend on transient ligand binding. Thus, identifying the detailed structure and transient positioning of small ligands, and their stabilization by the surrounding protein, is of great importance. In this study, by decoding information from Fourier transform infrared (FTIR) spectra with biomolecular simulation methods, we identify the precise position and hydrogen network of a small compound, the guanosine diphosphate (GDP)-H(2)PO(4)(-) intermediate, in the surrounding protein-protein complex of Ras and its GTPase-activating protein, a central molecular switch in cellular signal transduction. We validate the simulated structure by comparing the calculated fingerprint vibrational modes of H(2)PO(4)(-) with those obtained from FTIR experiments. The new structural information, below the resolution of X-ray structural analysis, gives detailed insight into the catalytic mechanism.

N-Ras forms dimers at POPC membranes

Ras is a central regulator of cellular signaling pathways. It is mutated in 20-30% of human tumors. To perform its function, Ras has to be bound to a membrane by a posttranslationally attached lipid anchor. Surprisingly, we identified here dimerization of membrane anchored Ras by combining attenuated total reflectance Fourier transform infrared spectroscopy, biomolecular simulations, and Förster resonance energy transfer experiments. By analyzing x-ray structural models and molecular-dynamics simulations, we propose a dimerization interface between α-helices 4 and 5 and the loop between β2 and β3. This seems to explain why the residues D47, E49, R135, R161, and R164 of this interface are influencing Ras signaling in cellular physiological experiments, although they are not positioned in the catalytic site. Dimerization could catalyze nanoclustering, which is well accepted for membrane-bound Ras. The interface could provide a new target for a seemingly novel type of small molecule interfering with signal transduction in oncogenic Ras mutants.

The role of magnesium for geometry and charge in GTP hydrolysis, revealed by quantum mechanics/molecular mechanics simulations

The coordination of the magnesium ion in proteins by triphosphates plays an important role in catalytic hydrolysis of GTP or ATP, either in signal transduction or energy conversion. For example, in Ras the magnesium ion contributes to the catalysis of GTP hydrolysis. The cleavage of GTP to GDP and P(i) in Ras switches off cellular signaling. We analyzed GTP hydrolysis in water, Ras, and Ras·Ras-GTPase-activating protein using quantum mechanics/molecular mechanics simulations. By comparison of the theoretical IR-difference spectra for magnesium ion coordinated triphosphate to experimental ones, the simulations are validated. We elucidated thereby how the magnesium ion contributes to catalysis. It provides a temporary storage for the electrons taken from the triphosphate and it returns them after bond cleavage and P(i) release back to the diphosphate. Furthermore, the Ras·Mg(2+) complex forces the triphosphate into a stretched conformation in which the β- and γ-phosphates are coordinated in a bidentate manner. In this conformation, the triphosphate elongates the bond, which has to be cleaved during hydrolysis. Furthermore, the γ-phosphate adopts a more planar structure, driving the conformation of the molecule closer to the hydrolysis transition state. GTPase-activating protein enhances these changes in GTP conformation and charge distribution via the intruding arginine finger.

The specific vibrational modes of GTP in solution and bound to Ras: a detailed theoretical analysis by QM/MM simulations

The hydrolysis of guanosine triphosphate (GTP) in general, and especially by GTPases like the Ras protein, is in the focus of biological investigations. A huge amount of experimental data from Fourier-transformed infrared studies is currently available, and many vibrational bands of free GTP, GTP·Mg(2+), and Ras·GTP·Mg(2+) in solution have been assigned by isotopic labeling. In the Ras environment, bands between 800 cm(-1) and 1300 cm(-1) have already been assigned, but not those below 800 cm(-1). The combination of quantum and molecular mechanics (QM/MM) methods takes the quantum effects for selected relevant atoms into account. This provides structural details, vibrational frequencies and electron distributions of the region of interest. We therefore used MM and QM/MM simulations to investigate the normal vibrational modes of GTP, GTP·Mg(2+), and Ras·GTP·Mg(2+) in solution, and assigned the vibrational frequencies for each normal vibration mode. In this study, the quantum box contains the nucleoside and the Mg(2+). The comparison of calculated and experimental vibrational spectra provides a very good control for the quality of the calculations. Structurally, MM and QM/MM simulations reveal a stable tridentate coordination of the Mg(2+) by GTP in water, and a stable bidentate coordination by GTP in complex with Ras. For validation, we compare the calculated frequencies and isotopic shifts with the experimental results available in the range of 800 cm(-1) to 1300 cm(-1). For the first time we suggest band assignments of the vibrational modes below 800 cm(-1) by comparison of calculated and experimental spectra.