Atomistic View of the Conformational Activation of Src Kinase Using the String Method with Swarms-of-Trajectories

The inactive-to-active conformational transition of the catalytic domain of human c-Src tyrosine kinase is characterized using the string method with swarms-of-trajectories with all-atom explicit solvent molecular dynamics simulations. The activation process occurs in two main steps in which the activation loop (A-loop) opens first, followed by the rotation of the αC helix. The computed potential of mean force energy along the activation pathway displays a local minimum, which allows the identification of an intermediate state. These results show that the string method with swarms-of-trajectories is an effective technique to characterize complex and slow conformational transitions in large biomolecular systems.     

A conformational transition in the myosin VI converter contributes to the variable step size

Myosin VI (MVI) is a dimeric molecular motor that translocates backwards on actin filaments with a surprisingly large and variable step size, given its short lever arm. A recent x-ray structure of MVI indicates that the large step size can be explained in part by a novel conformation of the converter subdomain in the prepowerstroke state, in which a 53-residue insert, unique to MVI, reorients the lever arm nearly parallel to the actin filament. To determine whether the existence of the novel converter conformation could contribute to the step-size variability, we used a path-based free-energy simulation tool, the string method, to show that there is a small free-energy difference between the novel converter conformation and the conventional conformation found in other myosins. This result suggests that MVI can bind to actin with the converter in either conformation. Models of MVI/MV chimeric dimers show that the variability in the tilting angle of the lever arm that results from the two converter conformations can lead to step-size variations of ∼12 nm. These variations, in combination with other proposed mechanisms, could explain the experimentally determined step-size variability of ∼25 nm for wild-type MVI. Mutations to test the findings by experiment are suggested.     

Locking the Active Conformation of c-Src Kinase through the Phosphorylation of the Activation Loop

Molecular dynamics umbrella sampling simulations are used to compare the relative stability of the active conformation of the catalytic domain of c-Src kinase while the tyrosine 416 in the activation loop (A-loop) is either unphosphorylated or phosphorylated. When the A-loop is unphosphorylated, there is considerable flexibility of the kinase. While the active conformation of the kinase is not forbidden and can be visited transiently, it is not the predominant state. This is consistent with the view that c-Src displays some catalytic activity even when the A-loop is unphosphorylated. In contrast, phosphorylation of the A-loop contributes to stabilize several structural features that are critical for catalysis, such as the hydrophobic regulatory spine, the HRD motif, and the electrostatic switch. In summary, the free-energy landscape calculations demonstrate that phosphorylation of tyrosine 416 in the A-loop essentially “locks” the kinase into its catalytically competent conformation.     

All-atom folding of the three-helix HIV accessory protein with an adaptive parallel tempering method

All-atom protein structure prediction from the amino acid sequence alone remains an important goal of biophysical chemistry. Recent progress in force field development and validation suggests that the PFF01 free-energy force field correctly predicts the native conformation of various helical proteins as the global optimum of its free-energy surface. Reproducible protein structure prediction requires the availability of efficient optimization methods to locate the global minima of such complex potentials. Here we investigate an adapted version of the parallel tempering method as an efficient parallel stochastic optimization method for protein structure prediction. Using this approach we report the reproducible all-atom folding of the three-helix 40 amino acid HIV accessory protein from random conformations to within 2.4 A backbone RMS deviation from the experimental structure with modest computational resources.     

Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

The βα-repeat class of proteins, represented by the (βα)₈ barrel and the α/β/α sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of βα-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Gō-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY.     

c-Src binds to the cancer drug imatinib with an inactive Abl/c-Kit conformation and a distributed thermodynamic penalty

The cancer drug imatinib inhibits the tyrosine kinases c-Abl, c-Kit, and the PDGF receptor. Imatinib is less effective against c-Src, which is difficult to understand because residues interacting with imatinib in crystal structures of Abl and c-Kit are conserved in c-Src. The crystal structure of the c-Src kinase domain in complex with imatinib closely resembles that of Abl*imatinib and c-Kit*imatinib, and differs significantly from the inactive "Src/CDK" conformation of the Src family kinases. Attempts to increase the affinity of c-Src for imatinib by swapping residues with the corresponding residues in Abl have not been successful, suggesting that the thermodynamic penalty for adoption of the imatinib-binding conformation by c-Src is distributed over a broad region of the structure. Two mutations that are expected to destabilize the inactive Src/CDK conformation increase drug sensitivity 15-fold, suggesting that the free-energy balance between different inactive states is a key to imatinib binding.     

Helical order in tarantula thick filaments requires the "closed" conformation of the myosin head

Myosin heads are helically ordered on the thick filament surface in relaxed muscle. In mammalian and avian filaments this helical arrangement is dependent on temperature and it has been suggested that helical order is related to ATP hydrolysis by the heads. To test this hypothesis, we have used electron microscopy and image analysis to study the ability and temperature dependence of analogs of ATP and ADP.Pi to induce helical order in tarantula thick filaments. ATP or analogs were added to rigor myofibrils or purified thick filaments at 22 degrees C and 4 degrees C and the samples negatively stained. The ADP.Pi analogs ADP.AlF4 and ADP.Vi, and the ATP analogs ADP.BeFx, AMPPNP and ATPgammaNH2, all induced helical order in tarantula thick filaments, independent of temperature. In the absence of nucleotide, or in the presence of ADP or the ATP analog, ATPgammaS, there was no helical ordering. According to crystallographic and tryptophan fluorescence studies, all of these analogs, except ATPgammaS and ADP, induce the "closed" conformation of the myosin head (in which the gamma phosphate pocket is closed). We suggest that helical order requires the closed conformation of the myosin head but is not dependent on the hydrolysis of ATP.     

Acetonitrile-induced conformational transitions in poly-l-lysine

The effect of acetonitrile on the random coil, α-helix and β-sheet conformations induced in poly-

Full-size image

-lysine is studied. It is found that acetonitrile at higher concentrations transforms the backbone of polylysine from a random coil to a helical conformation. Addition of acetonitrile to polylysine (pH 11.5) in the α-helix conformation, induces conformational changes in two stages. At concentrations below 60% v/v, acetonitrile stabilizes the helical conformation and at higher concentrations (>70% v/v), it destabilizes the helix. β-sheet→α-helix→random coil conformational transitions are found to occur when polylysine in the heat-induced conformation is titrated with acetonitrile. The possible mechanism(s) of action of acetonitrile in inducing these structural transitions is discussed.     

Thermodynamic measurements of bilayer insertion of a single transmembrane helix chaperoned by fluorinated surfactants

Accurate determination of the free energy of transfer of a helical segment from an aqueous into a transmembrane (TM) conformation is essential for understanding and predicting the folding and stability of membrane proteins. Until recently, direct thermodynamically sound measurements of free energy of insertion of hydrophobic TM peptides were impossible due to peptide aggregation outside the lipid bilayer. Here, we overcome this problem by using fluorinated surfactants that are capable of preventing aggregation but, unlike detergents, do not themselves interact with the bilayer. We have applied the fluorescence correlation spectroscopy methodology to study surfactant-chaperoned insertion into preformed POPC (palmitoyloleoylphosphatidylcholine) vesicles of the two well-studied dye-labeled TM peptides of different lengths: WALP23 and WALP27. Extrapolation of the apparent free-energy values measured in the presence of surfactants to a zero surfactant concentration yielded free-energy values of -9.0±0.1 and -10.0±0.1 kcal/mol for insertion of WALP23 and WALP27, respectively. Circular dichroism measurements confirmed helical structure of peptides in lipid bilayer, in the presence of surfactants, and in aqueous mixtures of organic solvents. From a combination of thermodynamic and conformational measurements, we conclude that the partitioning of a four-residue L-A-L-A segment in the context of a continuous helical conformation from an aqueous environment into the hydrocarbon core of the membrane has a favorable free energy of 1 kcal/mol. Our measurements, combined with the predictions of two independent experimental hydrophobicity scales, indicate that the per-residue cost of transfer of the helical backbone from water to the hydrocarbon core of the lipid bilayer is unfavorable and is equal to +2.13±0.17 kcal/mol.     

Active and inactive orientations of the transmembrane and cytosolic domains of the erythropoietin receptor dimer

Binding of erythropoietin to the erythropoietin receptor (EpoR) extracellular domain orients the transmembrane (TM) and cytosolic regions of the receptor dimer into an unknown activated conformation. By replacing the EpoR extracellular domain with a dimeric coiled coil, we engineered TM EpoR fusion proteins where the helical TM domains were constrained into seven possible relative orientations. We identify one dimeric TM conformation that imparts full activity to the cytosolic domain of the receptor and signals via JAK2, STAT proteins, and MAP kinase, one partially active orientation that preferentially activates MAP kinase, and one conformation corresponding to the inactive receptor. The active and inactive conformations were independently identified by computational searches for low-energy TM dimeric structures. We propose a specific EpoR-activated interface and suggest its use for structural and signaling studies.     

A Catalytically Disabled Double Mutant of Src Tyrosine Kinase Can Be Stabilized into an Active-Like Conformation

Tyrosine kinases are enzymes playing a critical role in cellular signaling. Molecular dynamics umbrella sampling potential of mean force computations are used to quantify the impact of activating and inactivating mutations of c-Src kinase. The potential of mean force computations predict that a specific double mutant can stabilize c-Src kinase into an active-like conformation while disabling the binding of ATP in the catalytic active site. The active-like conformational equilibrium of this catalytically dead kinase is affected by a hydrophobic unit that connects to the hydrophobic spine network via the C-helix. The αC-helix plays a crucial role in integrating the hydrophobic residues, making it a hub for allosteric regulation of kinase activity and the active conformation. The computational free-energy landscapes reported here illustrate novel design principles focusing on the important role of the hydrophobic spines. The relative stability of the spines could be exploited in future efforts to artificially engineer active-like but catalytically dead forms of protein kinases.     

Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Conformation and assembly of polypeptide scaffolds in templating the synthesis of silica: an example of a polylysine macromolecular "switch"

Although the role of polycationic macromolecules in catalyzing the synthesis of silica structures is well established, detailed understanding of the mechanisms behind the production of silica structures of controlled morphologies remains unclear. In this study, we have used both poly-L-lysine (PLL) and/or poly-D-lysine (PDL) for silica synthesis to investigate mechanisms controlling inorganic morphologies. The formation of both spherical silica particles and hexagonal plates was observed. The formation of hexagonal plates was suggested, via circular dichroic spectroscopy (CD), to result from the assembly of helical polylysine molecules. We confirm that the formation of PLL helices is a prerequisite to the hexagonal silica synthesis. In addition, we present for the first time that the handedness of the helicity of the macromolecule does not affect the formation of hexagonal silica. We also show, by using two different silica precursors, that the precursor does not have a direct effect on the formation of hexagonal silica plates. Furthermore, when polylysine helices were converted to beta-sheet structure, only silica particles were obtained, thus suggesting that the adoption of a helical conformation by PLL is required for the formation of hexagonally organized silica. These results demonstrate that the change in polylysine conformation can act as a "switch" in silica structure formation and suggest the potential for controlling morphologies and structures of inorganic materials via control of the conformation of soft macromolecular templates.     

Free-energy landscape of the villin headpiece in an all-atom force field

We investigate the landscape of the internal free-energy of the 36 amino acid villin headpiece with a modified basin hopping method in the all-atom force field PFF01, which was previously used to predictively fold several helical proteins with atomic resolution. We identify near native conformations of the protein as the global optimum of the force field. More than half of the twenty best simulations started from random initial conditions converge to the folding funnel of the native conformation, but several competing low-energy metastable conformations were observed. From 76,000 independently generated conformations we derived a decoy tree which illustrates the topological structure of the entire low-energy part of the free-energy landscape and characterizes the ensemble of metastable conformations. These emerge as similar in secondary content, but differ in tertiary arrangement.     

Free-energy profiles of membrane insertion of the M2 transmembrane peptide from influenza A virus

The insertion of the M2 transmembrane peptide from influenza A virus into a membrane has been studied with molecular-dynamics simulations. This system is modeled by an atomically detailed peptide interacting with a continuum representation of a membrane bilayer in aqueous solution. We performed replica-exchange molecular-dynamics simulations with umbrella-sampling techniques to characterize the probability distribution and conformation preference of the peptide in the solution, at the membrane interface, and in the membrane. The minimum in the calculated free-energy surface of peptide insertion corresponds to a fully inserted, helical peptide spanning the membrane. The free-energy profile also shows that there is a significant barrier for the peptide to enter into this minimum in a nonhelical conformation. The sequence of the peptide is such that hydrophilic amino acid residues at the ends of the otherwise primarily hydrophobic peptide create a trapped, U-shaped conformation with the hydrophilic residues associated with the aqueous phase and the hydrophobic residues embedded in the membrane. Analysis of the free energy shows that the barrier to insertion is largely enthalpic in nature, whereas the membrane-spanning global minimum is favored by entropy.     

Structure and regulation of the human Nek2 centrosomal kinase

The dimeric Ser/Thr kinase Nek2 regulates centrosome cohesion and separation through phosphorylation of structural components of the centrosome, and aberrant regulation of Nek2 activity can lead to aneuploid defects characteristic of cancer cells. Mutational analysis of autophosphorylation sites within the kinase domain identified by mass spectrometry shows a complex pattern of positive and negative regulatory effects on kinase activity that are correlated with effects on centrosomal splitting efficiency in vivo. The 2.2-A resolution x-ray structure of the Nek2 kinase domain in complex with a pyrrole-indolinone inhibitor reveals an inhibitory helical motif within the activation loop. This helix presents a steric barrier to formation of the active enzyme and generates a surface that may be exploitable in the design of specific inhibitors that selectively target the inactive state. Comparison of this "auto-inhibitory" conformation with similar arrangements in cyclin-dependent kinase 2 and epidermal growth factor receptor kinase suggests a role for dimerization-dependent allosteric regulation that combines with autophosphorylation and protein phosphatase 1c phosphatase activity to generate the precise spatial and temporal control required for Nek2 function in centrosomal maturation.     

Innovative, multifunctional sequential oligopeptide carriers Soc(n)-I and SOC(n)-II: functions-technology-perspectives

An innovative type of multifunctional helicoid artificial carriers, formed by the repetitive Lys-Aib-Gly (SOC(n)-I) or by the Aib-Lys-Aib-Gly (SOC(n)-II), with structural rigidity and regularity were successfully applied in our laboratory for anchoring antigenic/immunogenic peptides. The carriers, designed to display a predetermined 3D structure, adopt the 3(10) helical conformation, while the attached peptides preserve their original "active" conformation. The constructed conjugates were used as substrates in solid phase immunoassays, as well as for generating potent and specific immune responses.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

Exploring the protein G helix free-energy surface by solute tempering metadynamics

The free-energy landscape of the alpha-helix of protein G is studied by means of metadynamics coupled with a solute tempering algorithm. Metadynamics allows to overcome large energy barriers, whereas solute tempering improves the sampling with an affordable computational effort. From the sampled free-energy surface we are able to reproduce a number of experimental observations, such as the fact that the lowest minimum corresponds to a globular conformation displaying some degree of beta-structure, that the helical state is metastable and involves only 65% of the chain. The calculations also show that the system populates consistently a pi-helix state and that the hydrophobic staple motif is present only in the free-energy minimum associated with the helices, and contributes to their stabilization. The use of metadynamics coupled with solute tempering results then particularly suitable to provide the thermodynamics of a short peptide, and its computational efficiency is promising to deal with larger proteins.     

Design, synthesis, and biophysical properties of a helical Abeta1-42 analog: Inhibition of fibrillogenesis and cytotoxicity

The aggregation of amyloid β-peptide (Aβ) into β-sheet-rich aggregates is a crucial step in the etiology of Alzheimer’s disease. Helical forms of Aβ have been suggested to be intermediates in the aggregation process of the peptide in aqueous phase, micelles and membranes. A stable helical Aβ analog would be useful to investigate the role of helical intermediates in fibrillization by Aβ. Here we designed a helical analog by simply cross-linking the Cys residues of A30C, G37C-Aβ1-42 with 1,6-bismaleimidohexane. The analog assumed a weak α-helical conformation in model membranes mimicking lipid raft microdomains of neuronal membranes under conditions in which the wild-type Aβ1-42 formed a β-sheet, indicating the cross-linking locally induced a helical conformation. Furthermore, addition of equimolar helical Aβ analog significantly reduced the amyloid formation and cytotoxicity by Aβ1-42. Thus, our helical Aβ1-42 is not only a model peptide to investigate the role of helical intermediates in fibrillization by Aβ, but also an inhibitor of Aβ-induced cytotoxicity.