Conformational dynamics in a truncated epidermal growth factor receptor ectodomain

Structural studies have revealed two forms of the monomeric epidermal growth factor receptor (EGFR) ectodomain: a compact (tethered) form stabilized by interdomain interactions and an extended (untethered) form in the presence of ligand. An important question is whether the ligand induces a conformational transition from a tethered to untethered form or whether there is a preexisting conformational equilibrium between tethered and untethered states. To distinguish between these two possibilities, we investigated a truncated receptor, EGFR501 (spanning residues 1-501), that contains the minimal elements required for high-affinity ligand binding in solution. Conformational transitions and dynamics were inferred by means of fluorescence from five internal tryptophan residues that are located within or close to the ligand-binding domains of EGFR501. A preexisting conformational equilibrium between tethered and untethered states in EGFR501 was deduced from (1) the nonlinear Arrhenius temperature dependence of fluorescence and (2) fluorescence polarization showing independently mobile domains. In contrast, the ligand-EGFR501 complex revealed a linear Arrhenius temperature dependence of fluorescence and increased fluorescence polarization due to a lack of significant interdomain motions. The data suggest that the role of the ligand is to trap the EGFR501 in the untethered state that is transiently formed in solution through a preexisting conformational equilibrium.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

Structure-Based Simulations of the Translocation Mechanism of the Hepatitis C Virus NS3 Helicase along Single-Stranded Nucleic Acid

The NS3 helicase of Hepatitis C virus is an ATP-fueled molecular motor that can translocate along single-stranded (ss) nucleic acid, and unwind double-stranded nucleic acids. It makes a promising antiviral target and an important prototype system for helicase research. Despite recent progress, the detailed mechanism of NS3 helicase remains unknown. In this study, we have combined coarse-grained (CG) and atomistic simulations to probe the translocation mechanism of NS3 helicase along ssDNA. At the residue level of detail, our CG simulations have captured functionally important interdomain motions of NS3 helicase and reproduced single-base translocation of NS3 helicase along ssDNA in the 3′–5′ direction, which is in good agreement with experimental data and the inchworm model. By combining the CG simulations with residue-specific perturbations to protein-DNA interactions, we have identified a number of key residues important to the translocation machinery that agree with previous structural and mutational studies. Additionally, our atomistic simulations with targeted molecular dynamics have corroborated the findings of CG simulations and further revealed key protein-DNA hydrogen bonds that break/form during the transitions. This study offers, to our knowledge, the most detailed and realistic simulations of translocation mechanism of NS3 helicase. The simulation protocol established in this study will be useful for designing inhibitors that target the translocation machinery of NS3 helicase, and for simulations of a variety of nucleic-acid-based molecular motors.     

Hierarchical Modeling of Activation Mechanisms in the ABL and EGFR Kinase Domains: Thermodynamic and Mechanistic Catalysts of Kinase Activation by Cancer Mutations

Structural and functional studies of the ABL and EGFR kinase domains have recently suggested a common mechanism of activation by cancer-causing mutations. However, dynamics and mechanistic aspects of kinase activation by cancer mutations that stimulate conformational transitions and thermodynamic stabilization of the constitutively active kinase form remain elusive. We present a large-scale computational investigation of activation mechanisms in the ABL and EGFR kinase domains by a panel of clinically important cancer mutants ABL-T315I, ABL-L387M, EGFR-T790M, and EGFR-L858R. We have also simulated the activating effect of the gatekeeper mutation on conformational dynamics and allosteric interactions in functional states of the ABL-SH2-SH3 regulatory complexes. A comprehensive analysis was conducted using a hierarchy of computational approaches that included homology modeling, molecular dynamics simulations, protein stability analysis, targeted molecular dynamics, and molecular docking. Collectively, the results of this study have revealed thermodynamic and mechanistic catalysts of kinase activation by major cancer-causing mutations in the ABL and EGFR kinase domains. By using multiple crystallographic states of ABL and EGFR, computer simulations have allowed one to map dynamics of conformational fluctuations and transitions in the normal (wild-type) and oncogenic kinase forms. A proposed multi-stage mechanistic model of activation involves a series of cooperative transitions between different conformational states, including assembly of the hydrophobic spine, the formation of the Src-like intermediate structure, and a cooperative breakage and formation of characteristic salt bridges, which signify transition to the active kinase form. We suggest that molecular mechanisms of activation by cancer mutations could mimic the activation process of the normal kinase, yet exploiting conserved structural catalysts to accelerate a conformational transition and the enhanced stabilization of the active kinase form. The results of this study reconcile current experimental data with insights from theoretical approaches, pointing to general mechanistic aspects of activating transitions in protein kinases.     

Continuous Interdomain Orientation Distributions Reveal Components of Binding Thermodynamics

The flexibility of biological macromolecules is an important structural determinant of function. Unfortunately, the correlations between different motional modes are poorly captured by discrete ensemble representations. Here, we present new ways to both represent and visualize correlated interdomain motions. Interdomain motions are determined directly from residual dipolar couplings, represented as a continuous conformational distribution, and visualized using the disk-on-sphere representation. Using the disk-on-sphere representation, features of interdomain motions, including correlations, are intuitively visualized. The representation works especially well for multidomain systems with broad conformational distributions.This analysis also can be extended to multiple probability density modes, using a Bingham mixture model. We use this new paradigm to study the interdomain motions of staphylococcal protein A, which is a key virulence factor contributing to the pathogenicity of Staphylococcus aureus. We capture the smooth transitions between important states and demonstrate the utility of continuous distribution functions for computing the reorientational components of binding thermodynamics. Such insights allow for the dissection of the dynamic structural components of functionally important intermolecular interactions.     

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.     

Differential Modulation of Functional Dynamics and Allosteric Interactions in the Hsp90-Cochaperone Complexes with p23 and Aha1: A Computational Study

Allosteric interactions of the molecular chaperone Hsp90 with a large cohort of cochaperones and client proteins allow for molecular communication and event coupling in signal transduction networks. The integration of cochaperones into the Hsp90 system is driven by the regulatory mechanisms that modulate the progression of the ATPase cycle and control the recruitment of the Hsp90 clientele. In this work, we report the results of computational modeling of allosteric regulation in the Hsp90 complexes with the cochaperones p23 and Aha1. By integrating protein docking, biophysical simulations, modeling of allosteric communications, protein structure network analysis and the energy landscape theory we have investigated dynamics and stability of the Hsp90-p23 and Hsp90-Aha1 interactions in direct comparison with the extensive body of structural and functional experiments. The results have revealed that functional dynamics and allosteric interactions of Hsp90 can be selectively modulated by these cochaperones via specific targeting of the regulatory hinge regions that could restrict collective motions and stabilize specific chaperone conformations. The protein structure network parameters have quantified the effects of cochaperones on conformational stability of the Hsp90 complexes and identified dynamically stable communities of residues that can contribute to the strengthening of allosteric interactions. According to our results, p23-mediated changes in the Hsp90 interactions may provide “molecular brakes” that could slow down an efficient transmission of the inter-domain allosteric signals, consistent with the functional role of p23 in partially inhibiting the ATPase cycle. Unlike p23, Aha1-mediated acceleration of the Hsp90-ATPase cycle may be achieved via modulation of the equilibrium motions that facilitate allosteric changes favoring a closed dimerized form of Hsp90. The results of our study have shown that Aha1 and p23 can modulate the Hsp90-ATPase activity and direct the chaperone cycle by exerting the precise control over structural stability, global movements and allosteric communications in Hsp90.     

KCTD5 is endowed with large,functionally relevant,interdomain motions

The KCTD family is an emerging class of proteins that are involved in important biological processes whose biochemical and structural properties are rather poorly characterized or even completely undefined. We here used KCTD5, the only member of the family with a known three-dimensional structure, to gain insights into the intrinsic structural stability of the C-terminal domain (CTD) and into the mutual dynamic interplay between the two domains of the protein. Molecular dynamics (MD) simulations indicate that in the simulation timescale (120 ns), the pentameric assembly of the CTD is endowed with a significant intrinsic stability. Moreover, MD analyses also led to the identification of exposed β-strand residues. Being these regions intrinsically sticky, they could be involved in the substrate recognition. More importantly, simulations conducted on the full-length protein provide interesting information of the relative motions between the BTB domain and the CTD of the protein. Indeed, the dissection of the overall motion of the protein is indicative of a large interdomain twisting associated with limited bending movements. Notably, MD data indicate that the entire interdomain motion is pivoted by a single residue (Ser150) of the hinge region that connects the domains. The functional relevance of these motions was evaluated in the context of the functional macromolecular machinery in which KCTD5 is involved. This analysis indicates that the interdomain twisting motion here characterized may be important for the correct positioning of the substrate to be ubiquitinated with respect to the other factors of the ubiquitination machinery.     

Glycosylation effects on the structure and dynamics of a full-length Cel7A cellulase

Fungi cellulases are used to degrade cellulose-containing biomass for bioethanol production. Industrial cellulases such as Cel7A from Trichoderma reesei (TrCel7A) are critical in this process. Thus, the understanding of structure and dynamics is crucial for engineering variants with improved cellulolytic activity. This cellulase consists of two domains connected by a flexible and highly glycosylated linker. However, the linker flexibility has hindered the determination of Cel7A complete structure. Herein, based on atomic and sparse data, we applied integrative modelling to build a model of the complete enzyme structure. Next, through simulations, we studied the glycosylation effects on the structure and dynamics of a solubilized TrCel7A. Essential dynamics analysis showed that O-glycosylation in the linker led to the stabilization of protein overall dynamics. O-linked glycans seem to restrict protein dihedral angles distribution in this region, selecting more elongated conformations. Besides the reduced flexibility, functional interdomain motions occurred in a more concerted way in the glycosylated system. In contrast, in the absence of glycosylation, we observed vast conformational plasticity with the functional domains frequently collapsing. We report here evidence that targeting Cel7A linker flexibility by point mutations including modification of glycosylation sites could be a promising design strategy to improve cellulase activity.     

Mechanisms of intramolecular communication in a hyperthermophilic acylaminoacyl peptidase: a molecular dynamics investigation

Protein dynamics and the underlying networks of intramolecular interactions and communicating residues within the three-dimensional (3D) structure are known to influence protein function and stability, as well as to modulate conformational changes and allostery. Acylaminoacyl peptidase (AAP) subfamily of enzymes belongs to a unique class of serine proteases, the prolyl oligopeptidase (POP) family, which has not been thoroughly investigated yet. POPs have a characteristic multidomain three-dimensional architecture with the active site at the interface of the C-terminal catalytic domain and a β-propeller domain, whose N-terminal region acts as a bridge to the hydrolase domain. In the present contribution, protein dynamics signatures of a hyperthermophilic acylaminoacyl peptidase (AAP) of the prolyl oligopeptidase (POP) family, as well as of a deletion variant and alanine mutants (I12A, V13A, V16A, L19A, I20A) are reported. In particular, we aimed at identifying crucial residues for long range communications to the catalytic site or promoting the conformational changes to switch from closed to open ApAAP conformations. Our investigation shows that the N-terminal α1-helix mediates structural intramolecular communication to the catalytic site, concurring to the maintenance of a proper functional architecture of the catalytic triad. Main determinants of the effects induced by α1-helix are a subset of hydrophobic residues (V16, L19 and I20). Moreover, a subset of residues characterized by relevant interaction networks or coupled motions have been identified, which are likely to modulate the conformational properties at the interdomain interface.     

Solution NMR and computer simulation studies of active site loop motion in triosephosphate isomerase

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.     

Allosteric Regulation of the Hsp90 Dynamics and Stability by Client Recruiter Cochaperones: Protein Structure Network Modeling

The fundamental role of the Hsp90 chaperone in supporting functional activity of diverse protein clients is anchored by specific cochaperones. A family of immune sensing client proteins is delivered to the Hsp90 system with the aid of cochaperones Sgt1 and Rar1 that act cooperatively with Hsp90 to form allosterically regulated dynamic complexes. In this work, functional dynamics and protein structure network modeling are combined to dissect molecular mechanisms of Hsp90 regulation by the client recruiter cochaperones. Dynamic signatures of the Hsp90-cochaperone complexes are manifested in differential modulation of the conformational mobility in the Hsp90 lid motif. Consistent with the experiments, we have determined that targeted reorganization of the lid dynamics is a unifying characteristic of the client recruiter cochaperones. Protein network analysis of the essential conformational space of the Hsp90-cochaperone motions has identified structurally stable interaction communities, interfacial hubs and key mediating residues of allosteric communication pathways that act concertedly with the shifts in conformational equilibrium. The results have shown that client recruiter cochaperones can orchestrate global changes in the dynamics and stability of the interaction networks that could enhance the ATPase activity and assist in the client recruitment. The network analysis has recapitulated a broad range of structural and mutagenesis experiments, particularly clarifying the elusive role of Rar1 as a regulator of the Hsp90 interactions and a stability enhancer of the Hsp90-cochaperone complexes. Small-world organization of the interaction networks in the Hsp90 regulatory complexes gives rise to a strong correspondence between highly connected local interfacial hubs, global mediator residues of allosteric interactions and key functional hot spots of the Hsp90 activity. We have found that cochaperone-induced conformational changes in Hsp90 may be determined by specific interaction networks that can inhibit or promote progression of the ATPase cycle and thus control the recruitment of client proteins.     

The influence of interdomain interactions on the intradomain motions in yeast phosphoglycerate kinase: a molecular dynamics study

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations—in which interdomain interactions cannot occur—the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Flexibility in a drug transport accessory protein: molecular dynamics simulations of MexA

Drug resistance in gram-negative bacteria may be conferred via efflux through a tripartite complex of an inner membrane pump, an outer membrane pore, and a periplasmic adaptor protein. These are AcrB, TolC, and AcrA, respectively, in Escherichia coli. In Pseudomonas aerugonisa, their homologs are MexB, OprM, and MexA. Defining the interdomain dynamics of the adaptor protein is essential to understanding the mechanism of complex formation. Extended (25 ns) molecular dynamics simulations of MexA have been performed to determine such interdomain dynamics. Analysis of conformational drift demonstrates substantial motions of the three domains of MexA relative to one another. Principal components analysis reveals a hinge-bending motion and rotation of the alpha-helical hairpin relative to the other domains to be the two dominant motions. These two motions provide an element of considerable flexibility which is likely to be exploited in the adaptor function of MexA.     

Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

Systematic analysis of domain motions in proteins from conformational change: New results on citrate synthase and T4 lysozyme

Methods developed originally to analyze domain motions from simulation [Proteins 27:425–437, 1997] are adapted and extended for the analysis of X-ray conformers and for proteins with more than two domains. The method can be applied as an automatic procedure to any case where more than one conformation is available. The basis of the methodology is that domains can be recognized from the difference in the parameters governing their quasi-rigid body motion, and in particular their rotation vectors. A clustering algorithm is used to determine clusters of rotation vectors corresponding to main-chain segments that form possible dynamic domains. Domains are accepted for further analysis on the basis of a ratio of interdomain to intradomain fluctuation, and Chasles' theorem is used to determine interdomain screw axes. Finally residues involved in the interdomain motion are identified. The methodology is tested on citrate synthase and the M6I mutant of T4 lysozyme. In both cases new aspects to their conformational change are revealed, as are individual residues intimately involved in their dynamics. For citrate synthase the beta sheet is identified to be part of the hinging mechanism. In the case of T4 lysozyme, one of the four transitions in the pathway from the closed to the open conformation, furnished four dynamic domains rather than the expected two. This result indicates that the number of dynamic domains a protein possesses may not be a constant of the motion. Proteins 30:144–154, 1998. © 1998 Wiley-Liss, Inc.     

Effects of cyclization on conformational dynamics and binding properties of Lys48-linked di-ubiquitin

In solution, Lys48-linked di-ubiquitin exists in dynamic equilibrium between closed and open conformations. To understand the effect of interdomain motion in polyubiquitin chains on their ability to bind ligands, we cyclized di-ubiquitin by cross-linking the free C terminus of the proximal ubiquitin with the side chain of residue 48 in the distal ubiquitin, using a chemical cross-linker, 1,6-Hexane-bis-vinylsulfone. Our NMR studies confirm that the cyclization affects conformational dynamics in di-ubiquitin by restricting opening of the interface and shifting the conformational equilibrium toward closed conformations. The cyclization, however, did not rigidly lock di-ubiquitin in a single closed conformation: The chain undergoes slow exchange between at least two closed conformations, characterized by interdomain contacts involving the same hydrophobic patch residues (Leu8-Ile44-Val70) as in the uncyclized di-ubiquitin. Lowering the pH changes the relative populations of these conformations, but in contrast with the uncyclized di-ubiquitin, does not lead to opening of the interface. This restriction of domain motions inhibits direct access of protein molecules to the hydrophobic patch residues located at the very center of the interdomain interface in di-ubiquitin, although the residual motions are sufficient to allow access of small molecules to the interface. This renders di-ubiquitin unable to bind protein molecules (e.g., UBA2 domain) in the normal manner, and thus could interfere with Ub(2) recognition by various downstream effectors. These results emphasize the importance of the opening/closing domain motions for the recognition and function of di-ubiquitin and possibly longer polyubiquitin chains.     

The influence of amino acid side-chains on α-helix stability: S-peptide analogues and related ribonucleases S′

In order to determine the influence of amino acid side-chains on α-helix stability, in relation to the protein folding process, the coil-helix transitions of some synthetic [Orn 10]-S-peptide analogues, containing, in position 8, Phe, Tyr, Ile, Ala, cpGly² and Gly, were investigated by the technique of circular dichroism under two different sets of conditions. First, the transitions of the Speptide analogues in water/trifluoroethanol mixtures were recorded. From the pattern of the transitions and from the ellipticity values in 97% trifluoroethanol, the following increasing order of amino acids as α-helix formers was found: Gly < Tyr ≤ Phe < cpGly < Ala < Ile. This finding indicates that the conformational parameters (Chou & Fasman, 1974) of the residues in position 8 play an important but not exclusive role in α-helix stability, since the hydrophobicity of the side-chain (Nozaki & Tanford, 1971) of residue 8 exerts a strong influence. From the second approach, studying the capability of the S-peptide analogues to bind to S-protein, the following increasing order was found: (Gly, Ala) < Ile < cpGly < Tyr < Phe. This result reveals that the conformational parameters of the residues in position 8 play no role, whereas their hydrophobic character and side-chain interactions with surrounding residues in the S-protein portion are the determining binding factors. This finding explains the reason for the Phe8 invariance in RNAase A during evolution, and furnishes evidence for the relevant role of long-range interactions in the protein folding process.