Structure-based engineering of Tor complexes reveals that two types of yeast TORC1 produce distinct phenotypes.

Certain proteins assemble into diverse complex states, each having a distinct and unique function in the cell. Target of rapamycin (Tor) complex 1 (TORC1) plays a central role in signalling pathways that allow cells to respond to the environment, including nutritional status signalling. TORC1 is widely recognised for its association with various diseases. The budding yeast Saccharomyces cerevisiae has two types of TORC1, Tor1-containing TORC1 and Tor2-containing TORC1, which comprise different constituent proteins but are considered to have the same function. Here, we computationally modelled the relevant complex structures and then, based on the structures, rationally engineered a Tor2 mutant that could form Tor complex 2 (TORC2) but not TORC1, resulting in a redesign of the complex states. Functional analysis of the Tor2 mutant revealed that the two types of TORC1 induce different phenotypes, with changes observed in rapamycin, caffeine and pH dependencies of cell growth, as well as in replicative and chronological lifespan. These findings uncovered by a general approach with huge potential - model structure-based engineering - are expected to provide further insights into various fields such as molecular evolution and lifespan.

Design of complicated all-α protein structures.

A wide range of de novo protein structure designs have been achieved, but the complexity of naturally occurring protein structures is still far beyond these designs. Here, to expand the diversity and complexity of de novo designed protein structures, we sought to develop a method for designing 'difficult-to-describe' α-helical protein structures composed of irregularly aligned α-helices like globins. Backbone structure libraries consisting of a myriad of α-helical structures with five or six helices were generated by combining 18 helix-loop-helix motifs and canonical α-helices, and five distinct topologies were selected for de novo design. The designs were found to be monomeric with high thermal stability in solution and fold into the target topologies with atomic accuracy. This study demonstrated that complicated α-helical proteins are created using typical building blocks. The method we developed will enable us to explore the universe of protein structures for designing novel functional proteins.

Design of allosteric sites into rotary motor V1-ATPase by restoring lost function of pseudo-active sites.

Allostery produces concerted functions of protein complexes by orchestrating the cooperative work between the constituent subunits. Here we describe an approach to create artificial allosteric sites in protein complexes. Certain protein complexes contain subunits with pseudo-active sites, which are believed to have lost functions during evolution. Our hypothesis is that allosteric sites in such protein complexes can be created by restoring the lost functions of pseudo-active sites. We used computational design to restore the lost ATP-binding ability of the pseudo-active site in the B subunit of a rotary molecular motor, V1-ATPase. Single-molecule experiments with X-ray crystallography analyses revealed that binding of ATP to the designed allosteric site boosts this V1's activity compared with the wild-type, and the rotation rate can be tuned by modulating ATP's binding affinity. Pseudo-active sites are widespread in nature, and our approach shows promise as a means of programming allosteric control over concerted functions of protein complexes.

State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein.

G-protein coupled receptors (GPCRs) are known for their low stability and large conformational changes upon transitions between multiple states. A widely used method for stabilizing these receptors is to make chimeric receptors by fusing soluble proteins (i.e., fusion partner proteins) into the intracellular loop 3 (ICL3) connecting the transmembrane helices 5 and 6 (TM5 and TM6). However, this fusion approach requires experimental trial and error to identify appropriate soluble proteins, residue positions, and linker lengths for making the fusion. Moreover, this approach has not provided state-targeting stabilization of GPCRs. Here, to rationally stabilize a class A GPCR, adenosine A2A receptor (A2AR) in a target state, we carried out the custom-made de novo design of α-helical fusion partner proteins, which can fix the conformation of TM5 and TM6 to that in an inactive state of A2AR through straight helical connections without any kinks or intervening loops. The chimeric A2AR fused with one of the designs (FiX1) exhibited increased thermal stability. Moreover, compared with the wild type, the binding affinity of the chimera against the agonist NECA was significantly decreased, whereas that against the inverse agonist ZM241385 was similar, indicating that the inactive state was selectively stabilized. Our strategy contributes to the rational state-targeting stabilization of GPCRs.

Robust folding of a de novo designed ideal protein even with most of the core mutated to valine.

Protein design provides a stringent test for our understanding of protein folding. We previously described principles for designing ideal protein structures stabilized by consistent local and nonlocal interactions, based on a set of rules relating local backbone structures to tertiary packing motifs. The principles have made possible the design of protein structures having various topologies with high thermal stability. Whereas nonlocal interactions such as tight hydrophobic core packing have traditionally been considered to be crucial for protein folding and stability, the rules proposed by our previous studies suggest the importance of local backbone structures to protein folding. In this study, we investigated the robustness of folding of de novo designed proteins to the reduction of the hydrophobic core, by extensive mutation of large hydrophobic residues (Leu, Ile) to smaller ones (Val) for one of the designs. Surprisingly, even after 10 Leu and Ile residues were mutated to Val, this mutant with the core mostly filled with Val was found to not be in a molten globule state and fold into the same backbone structure as the original design, with high stability. These results indicate the importance of local backbone structures to the folding ability and high thermal stability of designed proteins and suggest a method for engineering thermally stabilized natural proteins.

ShadowR: a novel chromoprotein with reduced non-specific binding and improved expression in living cells.

Here we developed an orange light-absorbing chromoprotein named ShadowR as a novel acceptor for performing fluorescence lifetime imaging microscopy-based Förster resonance energy transfer (FLIM-FRET) measurement in living cells. ShadowR was generated by replacing hydrophobic amino acids located at the surface of the chromoprotein Ultramarine with hydrophilic amino acids in order to reduce non-specific interactions with cytosolic proteins. Similar to Ultramarine, ShadowR shows high absorption capacity and no fluorescence. However, it exhibits reduced non-specific binding to cytosolic proteins and is highly expressed in HeLa cells. Using tandem constructs and a LOVTRAP system, we showed that ShadowR can be used as a FRET acceptor in combination with donor mRuby2 or mScarlet in HeLa cells. Thus, ShadowR is a useful, novel FLIM-FRET acceptor.

QM/MM Geometry Optimization on Extensive Free-Energy Surfaces for Examination of Enzymatic Reactions and Design of Novel Functional Properties of Proteins.

Many remarkable molecular functions of proteins use their characteristic global and slow conformational dynamics through coupling of local chemical states in reaction centers with global conformational changes of proteins. To theoretically examine the functional processes of proteins in atomic detail, a methodology of quantum mechanical/molecular mechanical (QM/MM) free-energy geometry optimization is introduced. In the methodology, a geometry optimization of a local reaction center is performed with a quantum mechanical calculation on a free-energy surface constructed with conformational samples of the surrounding protein environment obtained by a molecular dynamics simulation with a molecular mechanics force field. Geometry optimizations on extensive free-energy surfaces by a QM/MM reweighting free-energy self-consistent field method designed to be variationally consistent and computationally efficient have enabled examinations of the multiscale molecular coupling of local chemical states with global protein conformational changes in functional processes and analysis and design of protein mutants with novel functional properties.

Theoretical study on excited states of bacteriochlorophyll a in solutions with density functional assessment.

The excited-state properties of bacteriochlorophyll (BChl) a in triethylamine, 1-propanol, and methanol are investigated with the time-dependent density functional theory by using the quantum mechanical and molecular mechanical reweighting free energy self-consistant field method. It is found that no prevalent density functionals can reproduce the experimental excited-state properties, i.e., the absorption and reorganization energies, of BChl a in the solutions. The parameter µ in the range-separated hybrid functional is therefore optimized to reproduce the differences of the absorption energies in the solutions. We examine the origin of the differences of the absorption energies in the solutions and find that sensitive balance between contributions of structural changes and solute-solvent interactions determines the differences. The accurate description of the excitation with the density functional with the adjusted parameter is therefore essential to the understanding of the excited-state properties of BChl a in proteins and also the mechanism of the photosynthetic systems.

Crucial role of protein flexibility in formation of a stable reaction transition state in an α-amylase catalysis.

Conformational flexibility of proteins provides enzymes with high catalytic activity. Although the conformational flexibility is known to be pivotal for the ligand binding and release, its role in the chemical reaction process of the reactive substrate remains unclear. We determined a transition state of an enzymatic reaction in a psychrophilic α-amylase by a hybrid molecular simulation that allows one to identify the optimal chemical state in an extensive conformational ensemble of protein. The molecular simulation uncovered that formation of the reaction transition state accompanies a large and slow movement of a loop adjacent to the catalytic site. Free energy calculations revealed that, although catalytic electrostatic potentials on the reactive moiety are formed by local and fast reorganization around the catalytic site, reorganization of the large and slow movement of the loop significantly contributes to reduction of the free energy barrier by stabilizing the local reorganization.

QM/MM Reweighting Free Energy SCF for Geometry Optimization on Extensive Free Energy Surface of Enzymatic Reaction.

We developed a quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimization method by which the geometry of a quantum chemically treated (QM) molecule is optimized on a free energy surface defined with thermal distribution of the surrounding molecular environment obtained by molecular dynamics simulation with a molecular mechanics (MM) force field. The method called QM/MM reweighting free energy self-consistent field combines a mean field theory of QM/MM free energy geometry optimization developed by Yamamoto (Yamamoto, T. J. Chem. Phys.2008, 129, 244104) with a reweighting scheme for updating the MM distribution introduced by Hu et al. (Hu, H., et al. J. Chem. Phys.2008, 128, 034105) and features high computational efficiency suitable for exploring the reaction free energy surface of extensive protein conformational space. The computational efficiency with improved treatment of a long-range electrostatic (ES) interaction using the Ewald summation technique permits one to take into account global conformational relaxation of an entire protein of an enzyme in the free energy geometry optimization of its reaction center. We applied the method to an enzymatic reaction of a substrate complex of psychrophilic α-amylase from Antarctic bacterium Pseudoalteromonas haloplanktis and succeeded in geometry optimizations of the reactant and the product of the catalytic reaction that involve large conformational changes of protein loops adjacent to the reaction center on time scales reaching sub-microseconds. We found that the adjacent loops in the reactant and the product form in different conformations and produce catalytic ES potentials on the reaction center.

Binding free energy calculations of adenosine deaminase inhibitor and the effect of methyl substitution in inhibitors.

The binding affinity of an inhibitor is often improved ten times or more by introducing a simple substituent, such as a methyl group or a chlorine atom. We have investigated this phenomenon in the case of adenosine deaminase (ADA) inhibitors using molecular dynamics (MD) simulations and binding free energy calculations, by the linear interaction energy (LIE) method. For MD simulations, the coordination bond parameters and partial charges of atoms around the zinc ion in ADA have been determined by referring to ab initio MO calculations. The calculated binding free energies for seven inhibitors agreed well with the experimental ones, with a maximum error of 1.2 kcal/mol. The effect of methyl substitution in inhibitor molecules was examined on the basis of MD trajectories. It is suggested that the increase in binding affinity is caused by both van der Waals stabilizations by amino acid residues in contact with the introduced methyl group and through favored overall interactions with surrounding residues in the binding pocket.