Protein folding and function: the N-terminal fragment in adenylate kinase

Structural determinants of the dominant conformational epitopes of phospholipase A2 receptor in primary membranous nephropathy

Anti-phospholipase A2 receptor autoantibody (PLA2R-Ab) plays a critical role in the pathogenesis of primary membranous nephropathy (PMN), an autoimmune kidney disease characterized by immune deposits in the glomerular subepithelial spaces and proteinuria. However, the mechanism of how PLA2R-Abs interact with the conformational epitope(s) of PLA2R has remained elusive. PLA2R is a single transmembrane helix receptor containing ten extracellular domains that begin with a CysR domain followed by a FnII and eight CTLD domains. Here, we examined the interactions of PLA2R-Ab with the full PLA2R protein, N-terminal domain truncations, and C-terminal domain deletions under either denaturing or physiological conditions. Our data demonstrate that the PLA2R-Abs against the dominant epitope (the N-terminal CysR-CTLD1 triple domain) possess weak cross-reactivities to the C-terminal domains beyond CTLD1. Moreover, both the CysR and CTLD1 domains are required to form a conformational epitope for PLA2R-Ab interaction, with FnII serving as a linker domain. Upon close examination, we also observed that patients with newly diagnosed PMN carry two populations of PLA2R-Abs in sera that react to the denatured CysR-CTLD3 (the PLA2R-Ab₁) and denatured CysR-CTLD1 (the PLA2R-Ab₂) domain complexes on Western blots, respectively. Furthermore, the PLA2R-Ab₁ appeared at an earlier time point than PLA2R-Ab₂ in patients, whereas the increased levels of PLA2R-Ab₂ coincided with the worsening of proteinuria. In summary, our data support that an integrated folding of the three PLA2R N-terminal domains, CysR, FnII, and CTLD1, is a prerequisite to forming the PLA2R conformational epitope and that the dominant epitope-reactive PLA2R-Ab₂ plays a critical role in PMN clinical progression.

Nature's Shortcut to Protein Folding

This Feature Article presents a view of the protein folding transition based on the hypothesis that Nature has built features within the sequences that enable a Shortcut to efficient folding. Nature's Shortcut is proposed to be the early establishment of a set of nonlocal weak contacts, constituting protein loops that significantly constrain regions of the collapsed disordered protein into a native-like low-resolution fluctuating topology of major sections of the backbone. Nature's establishment of this scaffold of nonlocal contacts is claimed to bypass what would otherwise be a nearly hopeless unaided search for the final three-dimensional structure in proteins longer than ∼100 amino acids. To support this main contention of the Feature Article, the loop hypothesis (LH) description of early folding events is experimentally tested with time-resolved Förster resonance energy transfer techniques for adenylate kinase, and the data are shown to be consistent with theoretical predictions from the sequential collapse model (SCM). The experimentally based LH and the theoretically founded SCM are argued to provide a unified picture of the role of nonlocal contacts as constituting Nature's Shortcut to protein folding. Importantly, the SCM is shown to reliably predict key nonlocal contacts utilizing only primary sequence information. This view on Nature's Shortcut is open to the protein community for further detailed assessment, including its practical consequences, by suitable application of advanced experimental and computational techniques.

Rational design of viscosity reducing mutants of a monoclonal antibody: hydrophobic versus electrostatic inter-molecular interactions

High viscosity of monoclonal antibody formulations at concentrations ≥100 mg/mL can impede their development as products suitable for subcutaneous delivery. The effects of hydrophobic and electrostatic intermolecular interactions on the solution behavior of MAB 1, which becomes unacceptably viscous at high concentrations, was studied by testing 5 single point mutants. The mutations were designed to reduce viscosity by disrupting either an aggregation prone region (APR), which also participates in 2 hydrophobic surface patches, or a negatively charged surface patch in the variable region. The disruption of an APR that lies at the interface of light and heavy chain variable domains, VH and VL, via L45K mutation destabilized MAB 1 and abolished antigen binding. However, mutation at the preceding residue (V44K), which also lies in the same APR, increased apparent solubility and reduced viscosity of MAB 1 without sacrificing antigen binding or thermal stability. Neutralizing the negatively charged surface patch (E59Y) also increased apparent solubility and reduced viscosity of MAB 1, but charge reversal at the same position (E59K/R) caused destabilization, decreased solubility and led to difficulties in sample manipulation that precluded their viscosity measurements at high concentrations. Both V44K and E59Y mutations showed similar increase in apparent solubility. However, the viscosity profile of E59Y was considerably better than that of the V44K, providing evidence that inter-molecular interactions in MAB 1 are electrostatically driven. In conclusion, neutralizing negatively charged surface patches may be more beneficial toward reducing viscosity of highly concentrated antibody solutions than charge reversal or aggregation prone motif disruption.

In the multi-domain protein adenylate kinase, domain insertion facilitates cooperative folding while accommodating function at domain interfaces

Having multiple domains in proteins can lead to partial folding and increased aggregation. Folding cooperativity, the all or nothing folding of a protein, can reduce this aggregation propensity. In agreement with bulk experiments, a coarse-grained structure-based model of the three-domain protein, E. coli Adenylate kinase (AKE), folds cooperatively. Domain interfaces have previously been implicated in the cooperative folding of multi-domain proteins. To understand their role in AKE folding, we computationally create mutants with deleted inter-domain interfaces and simulate their folding. We find that inter-domain interfaces play a minor role in the folding cooperativity of AKE. On further analysis, we find that unlike other multi-domain proteins whose folding has been studied, the domains of AKE are not singly-linked. Two of its domains have two linkers to the third one, i.e., they are inserted into the third one. We use circular permutation to modify AKE chain-connectivity and convert inserted-domains into singly-linked domains. We find that domain insertion in AKE achieves the following: (1) It facilitates folding cooperativity even when domains have different stabilities. Insertion constrains the N- and C-termini of inserted domains and stabilizes their folded states. Therefore, domains that perform conformational transitions can be smaller with fewer stabilizing interactions. (2) Inter-domain interactions are not needed to promote folding cooperativity and can be tuned for function. In AKE, these interactions help promote conformational dynamics limited catalysis. Finally, using structural bioinformatics, we suggest that domain insertion may also facilitate the cooperative folding of other multi-domain proteins.

Most individual protein domains fold in an all or nothing fashion. This cooperative folding is important because it reduces the existence of partially folded proteins which can stick to each other and create disease causing aggregates. However, numerous proteins have multiple domains, independent units of folding, stability and/or function. Several such proteins also fold cooperatively. It is thought that strong interactions between individual domains allow the folding to propagate from a nucleating domain to neighbouring ones and this enables cooperative folding in multi-domain proteins. Here, we computationally study the folding of the three-domain protein AKE and find instead that the topology of the protein, wherein the two less stable domains are inserted into the more stable one, promotes folding cooperativity. When the more stable domain is folded, the ends of the inserted domains are constrained and this allows them to fold easily. In such a protein topology, strong inter-domain interactions are not needed to promote folding cooperativity. Interface amino acids which would have been involved in ensuring that the domains fit together correctly can now be tuned for binding or catalysis or conformational transitions. Thus, inserted domains may be present in multi-domain proteins to promote both function and folding.

Understanding the mechanism of atovaquone drug resistance in Plasmodium falciparum cytochrome b mutation Y268S using computational methods

The rapid appearance of resistant malarial parasites after introduction of atovaquone (ATQ) drug has prompted the search for new drugs as even single point mutations in the active site of Cytochrome b protein can rapidly render ATQ ineffective. The presence of Y268 mutations in the Cytochrome b (Cyt b) protein is previously suggested to be responsible for the ATQ resistance in Plasmodium falciparum (P. falciparum). In this study, we examined the resistance mechanism against ATQ in P. falciparum through computational methods. Here, we reported a reliable protein model of Cyt bc1 complex containing Cyt b and the Iron-Sulphur Protein (ISP) of P. falciparum using composite modeling method by combining threading, ab initio modeling and atomic-level structure refinement approaches. The molecular dynamics simulations suggest that Y268S mutation causes ATQ resistance by reducing hydrophobic interactions between Cyt bc1 protein complex and ATQ. Moreover, the important histidine contact of ATQ with the ISP chain is also lost due to Y268S mutation. We noticed the induced mutation alters the arrangement of active site residues in a fashion that enforces ATQ to find its new stable binding site far away from the wild-type binding pocket. The MM-PBSA calculations also shows that the binding affinity of ATQ with Cyt bc1 complex is enough to hold it at this new site that ultimately leads to the ATQ resistance.

On the role of aggregation prone regions in protein evolution, stability, and enzymatic catalysis: insights from diverse analyses

The various roles that aggregation prone regions (APRs) are capable of playing in proteins are investigated here via comprehensive analyses of multiple non-redundant datasets containing randomly generated amino acid sequences, monomeric proteins, intrinsically disordered proteins (IDPs) and catalytic residues. Results from this study indicate that the aggregation propensities of monomeric protein sequences have been minimized compared to random sequences with uniform and natural amino acid compositions, as observed by a lower average aggregation propensity and fewer APRs that are shorter in length and more often punctuated by gate-keeper residues. However, evidence for evolutionary selective pressure to disrupt these sequence regions among homologous proteins is inconsistent. APRs are less conserved than average sequence identity among closely related homologues (≥80% sequence identity with a parent) but APRs are more conserved than average sequence identity among homologues that have at least 50% sequence identity with a parent. Structural analyses of APRs indicate that APRs are three times more likely to contain ordered versus disordered residues and that APRs frequently contribute more towards stabilizing proteins than equal length segments from the same protein. Catalytic residues and APRs were also found to be in structural contact significantly more often than expected by random chance. Our findings suggest that proteins have evolved by optimizing their risk of aggregation for cellular environments by both minimizing aggregation prone regions and by conserving those that are important for folding and function. In many cases, these sequence optimizations are insufficient to develop recombinant proteins into commercial products. Rational design strategies aimed at improving protein solubility for biotechnological purposes should carefully evaluate the contributions made by candidate APRs, targeted for disruption, towards protein structure and activity.

Biotechnology requires the large-scale expression, yield, and storage of recombinant proteins. Each step in protein production has the potential to cause aggregation as proteins, not evolved to exist outside the cell, endure the various steps involved in commercial manufacturing processes. Mechanistic studies into protein aggregation have revealed that certain sequence regions contribute more to the aggregation propensity of a protein than other sequence regions do. Efforts to disrupt these regions have thus far indicated that rational sequence engineering is a useful technique to reduce the aggregation of biotechnologically relevant proteins. To improve our ability to rationally engineer proteins with enhanced expression, solubility, and shelf-life we conducted extensive analyses of aggregation prone regions (APRs) within protein sequences to characterize the various roles these regions play in proteins. Findings from this work indicate that protein sequences have evolved by minimizing their aggregation propensities. However, we also found that many APRs are conserved in protein families and are essential to maintain protein stability and function. Therefore, the contributions that APRs, targeted for disruption, make towards protein stability and function should be carefully evaluated when improving protein solubility via rational design.

Structure prediction and binding sites analysis of curcin protein of Jatropha curcas using computational approaches

Ribosome inactivating proteins (RIPs) are defense proteins in a number of higher-plant species that are directly targeted toward herbivores. Jatropha curcas is one of the biodiesel plants having RIPs. The Jatropha seed meal, after extraction of oil, is rich in curcin, a highly toxic RIP similar to ricin, which makes it unsuitable for animal feed. Although the toxicity of curcin is well documented in the literature, the detailed toxic properties and the 3D structure of curcin has not been determined by X-ray crystallography, NMR spectroscopy or any in silico techniques to date. In this pursuit, the structure of curcin was modeled by a composite approach of 3D structure prediction using threading and ab initio modeling. Assessment of model quality was assessed by methods which include Ramachandran plot analysis and Qmean score estimation. Further, we applied the protein-ligand docking approach to identify the r-RNA binding residue of curcin. The present work provides the first structural insight into the binding mode of r-RNA adenine to the curcin protein and forms the basis for designing future inhibitors of curcin. Cloning of a future peptide inhibitor within J. curcas can produce non-toxic varieties of J. curcas, which would make the seed-cake suitable as animal feed without curcin detoxification.

Evidence of an intact N-terminal translocation sequence of human mitochondrial adenylate kinase 4

Adenylate kinases are abundant nucleoside monophosphate kinases, which catalyze the phosphorylation of AMP by using ATP or GTP as phosphate donors. A previously cloned cDNA was named adenylate kinase 4 (AK4) based on its sequence similarity with known AKs but with no confirmed AK enzyme activity. In the present study the AK4 cDNA was expressed in Escherichia coli and the substrate specificity and kinetic properties of the recombinant protein were characterized. The enzyme catalyzed the phosphorylation of AMP, dAMP, CMP and dCMP with ATP or GTP as phosphate donors and AK4 also phosphorylated AMP with UTP as phosphate donor. The kinetic parameters of the enzyme were determined for AMP and dAMP with ATP as phosphate donor and for AMP with GTP as phosphate donor. AK4 showed its highest efficiency when phosphorylating AMP with GTP and a slightly lower efficiency for the phosphorylation of AMP with ATP. Among the three reactions for which kinetics were performed, dAMP was the poorest substrate. The AK4 mitochondrial localization was confirmed by expression of AK4 as a fusion protein with GFP in HeLa cells. The mitochondrial import sequence was shown to be located within the first N-terminal 11 amino acid residues, very close to the ATP-binding region of the enzyme. Import analysis suggested that the mitochondrial import sequence was not cleaved and thus the enzyme retained its activity upon entering the mitochondria. Site directed mutagenesis of amino acids Lys 4 and Arg 7 showed that these two residues were essential for mitochondrial import.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Intra-molecular chaperone: the role of the N-terminal in conformational selection and kinetic control

The vast majority of the proteins in nature are under thermodynamic control, consistent with the universally accepted notion that proteins exist in their thermodynamically most stable state. Yet, recently a number of examples of proteins whose fold is under kinetic control have come to light. Their functions and environments vary. The first among these are some proteases, discovered in the early 1990s. There, an N-terminal proregion is self-cleaved after the protein folded, leaving the remainder of the chain in a kinetically trapped state. A related scenario was observed for microcin J25, an antibacterial peptide. This peptide presents a trapped covalently knotted conformation. The third and the most recently discovered case is the multidrug-resistant transporter protein, P-glycoprotein. There, a synonymous 'silent' mutation leads to ribosome stalling with a consequent altered kinetically trapped state. Here we argue that in all three examples, the N-terminal plays the role of an intra-molecular chaperone, that is, the N-terminal conformation selects among all competing local conformations of a downstream segment. By providing a pattern, the N-terminal chaperone segment assists the protein folding process. If the N-terminal is subsequently cleaved, the protein can be under kinetic control, since it is trapped in a thermodynamically less-stable state.

Molecular dynamics of apo-adenylate kinase: a distance replica exchange method for the free energy of conformational fluctuations

A large domain motion in adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes a large-scale rearrangement of its lid and AMP-binding domains when the open form closes over its substrates, AMP, and Mg2+-ATP, whereby the AMP-binding and lid domains come closer to the core. The third domain, the core, is relatively stable during this motion. A reaction coordinate that monitors the distance between the AMP-binding and core domains is selected to be able to compare with the results of energy transfer experiments. Sampling along this reaction coordinate is carried out by using a distance replica exchange method (DREM), where systems that differ by a restraint potential enforcing different reaction coordinate values are independently simulated with periodic attempts at exchange of these systems. Several methods are used to study the efficiency and convergence properties of the DREM simulation and compared with an analogous non-DREM simulation. The DREM greatly accelerates the rate and extent of configurational sampling and leads to equilibrium sampling as measured by monitoring collective modes obtained from a principal coordinate analysis. The potential of mean force along the reaction coordinate reveals a rather flat region for distances from the open to a relatively closed AKE conformation. The potential of mean force for smaller distances has a distinct minimum that is quite close to that found in the closed form X-ray structure. In concert with a decrease in the reaction coordinate distance (AMP-binding-to-core distance) the lid-to-core distance of AKE also decreases. Therefore, apo AKE can fluctuate from its open form to conformations that are quite similar to its closed form X-ray structure, even in the absence of its substrates.

Molecular Dynamics of Apo-Adenylate Kinase: A Principal Component Analysis

Adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes large-scale motions of its Lid and AMP-binding domains when its open form closes over its substrates, AMP and Mg2+-ATP. The third domain, the Core, is relatively stable during closing. The resulting trajectory is analyzed with a principal component analysis method that decomposes the atom motions into modes ordered by their decreasing contributions to the total protein fluctuation. Simulations at 303 K (normal T) and 500 K (high T) reveal that at both temperatures the first three modes account for 70% of the total fluctuation. The residues that contribute the most to these three modes are concentrated in the Lid and AMP-binding domains. Analysis of the normal T modes indicates that the Lid and AMP-binding domains sample a broad distribution of conformations indicating that AKE is designed to provide its substrates with a large set of conformations. The high T results show that the Lid initially closes toward the Core. Subsequently, the Lid rotates to a new stable conformation that is different from what is observed in the substrate-bound AKE. These results are discussed in the context of experimental data that indicate that adenylate kinases do sample more than one conformational state in solution and that each of these conformational states undergoes substantial fluctuations. A pair of residues is suggested for labeling that would be useful for monitoring distance fluctuations by energy transfer experiments.

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.

A GTP:AMP phosphotransferase, Adk2p, in Saccharomyces cerevisiae. Role of the C terminus in protein folding/stabilization, thermal tolerance, and enzymatic activity

Adenylate kinases participate in maintaining the homeostasis of cellular nucleotides. Depending on the yeast strains, the GTP:AMP phosphotransferase is encoded by the nuclear gene ADK2 with or without a single base pair deletion/insertion near the 3' end of the open reading frame, and the corresponding protein exists as either Adk2p (short) or Adk2p (long) in the mitochondrial matrix. These two forms are identical except that the three C-terminal residues of Adk2p (short) are changed in Adk2p (long), and the latter contains an additional nine amino acids at the C terminus of the protein. The short form of Adk2p has so far been considered to be inactive (Schricker, R., Magdolen, V., Strobel, G., Bogengruber, E., Breitenbach, M., and Bandlow, W. (1995) J. Biol. Chem. 270, 31103-31110). Using purified proteins, we show that at the physiological temperature for yeast growth (30 degrees C), both short and long forms of Adk2p are enzymatically active. However, in contrast to the short form, Adk2p (long) is quite resistant to thermal inactivation, urea denaturation, and degradation by trypsin. Unfolding of the long form by high concentrations of urea greatly stimulated its import into isolated mitochondria. Using an integration-based gene-swapping approach, we found that regardless of the yeast strains used, the steady state levels of endogenous Adk2p (long) in mitochondria were 5-10-fold lower compared with those of Adk2p (short). Together, these results suggest that the modified C-terminal domain in Adk2p (long) is not essential for enzyme activity, but it contributes to and strengthens protein folding and/or stability and is particularly important for maintaining enzyme activity under stress conditions.

Amino acid coupling patterns in thermophilic proteins

Structural analysis is useful in elucidating structural features responsible for enhanced thermal stability of proteins. However, due to the rapid increase of sequenced genomic data, there are far more protein sequences than the corresponding three-dimensional (3D) structures. The usual sequence-based amino acid composition analysis provides useful but simplified clues about the amino acid types related to thermal stability of proteins. In this work, we developed a statistical approach to identify the significant amino acid coupling sequence patterns in thermophilic proteins. The amino acid coupling sequence pattern is defined as any 2 types of amino acids separated by 1 or more amino acids. Using this approach, we construct the rho profiles for the coupling patterns. The rho value gives a measure of the relative occurrence of a coupling pattern in thermophiles compared with mesophiles. We found that thermophiles and mesophiles exhibit significant bias in their amino acid coupling patterns. We showed that such bias is mainly due to temperature adaptation instead of species or GC content variations. Though no single outstanding coupling pattern can adequately account for protein thermostability, we can use a group of amino acid coupling patterns having strong statistical significance (p values < 10(-7)) to distinguish between thermophilic and mesophilic proteins. We found a good correlation between the optimal growth temperatures of the genomes and the occurrences of the coupling patterns (the correlation coefficient is 0.89). Furthermore, we can separate the thermophilic proteins from their mesophilic orthologs using the amino acid coupling patterns. These results may be useful in the study of the enhanced stability of proteins from thermophiles-especially when structural information is scarce. Proteins 2005. (c) 2005 Wiley-Liss, Inc.

Potential folding-function interrelationship in proteins

The possibility is addressed that protein folding and function may be related via regions that are critical for both folding and function. This approach is based on the building blocks folding model that describes protein folding as binding events of conformationally fluctuating building blocks. Within these, we identify building block fragments that are critical for achieving the native fold. A library of such critical building blocks (CBBs) is constructed. Then, it is asked whether the functionally important residues fall in these CBB fragments. We find that for over two-thirds of the proteins in our library with available functional information, the catalytic or binding site residues lie within the CBB regions. From the evolutionary standpoint, a folding-function relationship is advantageous, since the need to guard against mutations is limited to one region. Furthermore, conformationally similar CBBs are found in globally unrelated proteins with different functions. Hence, substituting CBBs may lead to designed proteins with altered functions. We further find that the CBBs in our library are conformationally unstable.

Two structurally distinct and spatially compartmentalized adenylate kinases are expressed from the AK1 gene in mouse brain

Adenylate kinases (AK, EC 2.7.4.3) have been considered important enzymes for energy homeostasis and metabolic signaling. To gain a better understanding of their cell-specific significance we studied the structural and functional aspects of products of one adenylate kinase gene, AK1, in mouse tissues. By combined computer database comparison and Northern analysis of mRNAs, we identified transcripts of 0.7 and 2.0 kilobases with different 5' and 3' non-coding regions which result from alternative use of promoters and polyadenylation sites. These mRNAs specify two distinct proteins, AK1 and a membrane-bound AK1 isoform (AK1beta), which differ in their N-terminal end and are co-expressed in several tissues with high-energy demand, including the brain. Immunohistochemical analysis of brain tissue and primary neurons and astrocytes in culture demonstrated that AK1 isoforms are expressed predominantly in neurons. AK1beta, when tested in transfected COS-1 and N2a neuroblastoma cells, located at the cellular membrane and was able to catalyze phosphorylation of ADP in vitro. In addition, AK1beta mediated AMP-induced activation of recombinant ATP-sensitive potassium channels in the presence of ATP. Thus, two structurally distinct AK1 isoforms co-exist in the mouse brain within distinct cellular locations. These enzymes may function in promoting energy homeostasis in the compartmentalized cytosol and in translating cellular energetic signals to membrane metabolic sensors.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Reducing the computational complexity of protein folding via fragment folding and assembly

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.

Hierarchical protein folding pathways: a computational study of protein fragments

We have previously presented a building block folding model. The model postulates that protein folding is a hierarchical top-down process. The basic unit from which a fold is constructed, referred to as a hydrophobic folding unit, is the outcome of combinatorial assembly of a set of "building blocks." Results obtained by the computational cutting procedure yield fragments that are in agreement with those obtained experimentally by limited proteolysis. Here we show that as expected, proteins from the same family give very similar building blocks. However, different proteins can also give building blocks that are similar in structure. In such cases the building blocks differ in sequence, stability, contacts with other building blocks, and in their 3D locations in the protein structure. This result, which we have repeatedly observed in many cases, leads us to conclude that while a building block is influenced by its environment, nevertheless, it can be viewed as a stand-alone unit. For small-sized building blocks existing in multiple conformations, interactions with sister building blocks in the protein will increase the population time of the native conformer. With this conclusion in hand, it is possible to develop an algorithm that predicts the building block assignment of a protein sequence whose structure is unknown. Toward this goal, we have created sequentially nonredundant databases of building block sequences. A protein sequence can be aligned against these, in order to be matched to a set of potential building blocks.

Comparison of protein fragments identified by limited proteolysis and by computational cutting of proteins

Here we present a comparison between protein fragments produced by limited proteolysis and those identified by computational cutting based on the building block folding model. The principles upon which the two methods are based are different. Limited proteolysis of natively folded proteins occurs at flexible sites and never at the level of chain segments of regular secondary structure such as alpha-helices. Therefore, the targets for limited proteolysis are locally unfolded regions. In contrast, the computational cutting algorithm considers the compactness of the fragments, their nonpolar buried surface area, and their isolatedness, that is, the surface area which was buried prior to the cutting and becomes exposed subsequently. Despite the different criteria, there is an overall correspondence between sites or regions of limited proteolysis with those identified by computational cutting. The computational cutting method has been applied to several model proteins for which detailed limited proteolysis data are available, namely apomyoglobin, cytochrome c, ribonuclease A, alpha-lactalbumin, and thermolysin. As expected, more cuts are obtained computationally than experimentally and the agreement is better when a number of proteolytic enzymes are used. For example, cytochrome c is cleaved by thermolysin at 56-57, 45-46, and at 80-81, and by proteinase K at 48-49 and 50-51. Incubation of the noncovalent and native-like complex of cytochrome c fragments 1-56 and 57-104 with proteinase K yielded the gapped protein species 1-48/57-104 and finally 1-40/57-104. Computational cutting of cytochrome c reproduced the major experimental observations, with cuts at 47, 64-65 or 65-66 and 80-81 and an unstable 32-47 region not assigned to any building block. The next step, not addressed in this work, is to probe the ability of the generated fragments to fold independently. Since both the computational algorithm and limited proteolysis attempt to dissect the protein folding problem, the general agreement between the two procedures is gratifying. This consistency allows us to propose the use of limited proteolysis to produce protein fragments that can adopt an independent folding and, therefore, to study folding intermediates. The results of the present study appear to validate the building block folding model and are in line with the proposal that protein folding is a hierarchical process, where parts constituting local minima of energy fold first, with their subsequent association and mutual stabilization to finally yield the global fold.

Conformational polymorphism of wild-type and mutant prion proteins: Energy landscape analysis

Conformational transitions are thought to be the prime mechanism of prion diseases. In this study, the energy landscapes of a wild-type prion protein (PrP) and the D178N and E200K mutant proteins were mapped, enabling the characterization of the normal isoforms (PrP(C)) and partially unfolded isoforms (PrP(PU)) of the three prion protein analogs. It was found that the three energy landscapes differ in three respects: (i) the relative stability of the PrP(C) and the PrP(PU) states, (ii) the transition pathways from PrP(C) to PrP(PU), and (iii) the relative stability of the three helices in the PrP(C) state. In particular, it was found that although helix 1 (residues 144-156) is the most stable helix in wild-type PrP, its stability is dramatically reduced by both mutations. This destabilization is due to changes in the charge distribution that affects the internal salt bridges responsible for the greater stability of this helix in wild-type PrP. Although both mutations result in similar destabilization of helix 1, they a have different effect on the overall stability of PrP(C) and of PrP(PU) isoforms and on structural properties. The destabilization of helix 1 by mutations provides additional evidences to the role of this helix in the pathogenic transition from the PrP(C) to the pathogenic isoform PrP(SC).

Contributions to the binding free energy of ligands to avidin and streptavidin

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.

The refolding of type II shikimate kinase from Erwinia chrysanthemi after denaturation in urea

Shikimate kinase was chosen as a convenient representative example of the subclass of alpha/beta proteins with which to examine the mechanism of protein folding. In this paper we report on the refolding of the enzyme after denaturation in urea. As shown by the changes in secondary and tertiary structure monitored by far UV circular dichroism (CD) and fluorescence, respectively, the enzyme was fully unfolded in 4 m urea. From an analysis of the unfolding curve in terms of the two-state model, the stability of the folded state could be estimated as 17 kJ.mol-1. Approximately 95% of the enzyme activity could be recovered on dilution of the urea from 4 to 0.36 m. The results of spectroscopic studies indicated that refolding occurred in at least four kinetic phases, the slowest of which (k = 0.009 s-1) corresponded with the regain of shikimate binding and of enzyme activity. The two most rapid phases were associated with a substantial increase in the binding of 8-anilino-1-naphthalenesulfonic acid with only modest changes in the far UV CD, indicating that a collapsed intermediate with only partial native secondary structure was formed rapidly. The relevance of the results to the folding of other alpha/beta domain proteins is discussed.

Building blocks, hinge-bending motions and protein topology

Here we show that the locations of molecular hinges in protein structures fall between building block elements. Building blocks are fragments of the protein chain which constitute local minima. These elements fold first. In the next step they associate through a combinatorial assembly process. While chain-linked building blocks may be expected to trial-associate first, if unstable, alternate more stable associations will take place. Hence, we would expect that molecular hinges will be at such inter-building block locations, or at the less stable, unassigned regions. On the other hand, hinge-bending motions are well known to be critical for protein function. Hence, protein folding and protein function are evolutionarily related. Further, the pathways through which proteins attain their three dimensional folds are determined by protein topology. However, at the same time the locations of the hinges, and hinge-bending motions are also an outcome of protein topology. Thus, protein folding and function appear coupled, and relate to protein topology. Here we provide some results illustrating such a relationship.