Protocols for Molecular Modeling with Rosetta3 and RosettaScripts.

Previously, we published an article providing an overview of the Rosetta suite of biomacromolecular modeling software and a series of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987-2998]. The overwhelming positive response to this publication we received motivates us to here share the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking, small molecule docking, and protein design. This updated and expanded set of tutorials is needed, as since 2010 Rosetta has been fully redesigned into an object-oriented protein modeling program Rosetta3. Notable improvements include a substantially improved energy function, an XML-like language termed "RosettaScripts" for flexibly specifying modeling task, new analysis tools, the addition of the TopologyBroker to control conformational sampling, and support for multiple templates in comparative modeling. Rosetta's ability to model systems with symmetric proteins, membrane proteins, noncanonical amino acids, and RNA has also been greatly expanded and improved.

A Conserved Hydrophobic Core in Gαi1 Regulates G Protein Activation and Release from Activated Receptor.

G protein-coupled receptor-mediated heterotrimeric G protein activation is a major mode of signal transduction in the cell. Previously, we and other groups reported that the α5 helix of Gαi1, especially the hydrophobic interactions in this region, plays a key role during nucleotide release and G protein activation. To further investigate the effect of this hydrophobic core, we disrupted it in Gαi1 by inserting 4 alanine amino acids into the α5 helix between residues Gln(333) and Phe(334) (Ins4A). This extends the length of the α5 helix without disturbing the ß6-α5 loop interactions. This mutant has high basal nucleotide exchange activity yet no receptor-mediated activation of nucleotide exchange. By using structural approaches, we show that this mutant loses critical hydrophobic interactions, leading to significant rearrangements of side chain residues His(57), Phe(189), Phe(191), and Phe(336); it also disturbs the rotation of the α5 helix and the π-π interaction between His(57) and Phe(189) In addition, the insertion mutant abolishes G protein release from the activated receptor after nucleotide binding. Our biochemical and computational data indicate that the interactions between α5, α1, and ß2-ß3 are not only vital for GDP release during G protein activation, but they are also necessary for proper GTP binding (or GDP rebinding). Thus, our studies suggest that this hydrophobic interface is critical for accurate rearrangement of the α5 helix for G protein release from the receptor after GTP binding.

A conserved phenylalanine as a relay between the α5 helix and the GDP binding region of heterotrimeric Gi protein α subunit.

G protein activation by G protein-coupled receptors is one of the critical steps for many cellular signal transduction pathways. Previously, we and other groups reported that the α5 helix in the G protein α subunit plays a major role during this activation process. However, the precise signaling pathway between the α5 helix and the guanosine diphosphate (GDP) binding pocket remains elusive. Here, using structural, biochemical, and computational techniques, we probed different residues around the α5 helix for their role in signaling. Our data showed that perturbing the Phe-336 residue disturbs hydrophobic interactions with the ß2-ß3 strands and α1 helix, leading to high basal nucleotide exchange. However, mutations in ß strands ß5 and ß6 do not perturb G protein activation. We have highlighted critical residues that leverage Phe-336 as a relay. Conformational changes are transmitted starting from Phe-336 via ß2-ß3/α1 to Switch I and the phosphate binding loop, decreasing the stability of the GDP binding pocket and triggering nucleotide release. When the α1 and α5 helices were cross-linked, inhibiting the receptor-mediated displacement of the C-terminal α5 helix, mutation of Phe-336 still leads to high basal exchange rates. This suggests that unlike receptor-mediated activation, helix 5 rotation and translocation are not necessary for GDP release from the α subunit. Rather, destabilization of the backdoor region of the Gα subunit is sufficient for triggering the activation process.