Searching for Multiple Folding Pathways of a Nearly Symmetrical Protein: Temperature Dependent Φ-Value Analysis of the B Domain of Protein A

Helical Content Correlations and Hydration Structures of the Folding Ensemble of the B Domain of Protein A

The B domain of protein A (BdpA), a small three-helix bundle, folds on a time scale of a few microseconds with heterogeneous native and unfolded states. It is widely used as a model for understanding protein folding mechanisms. In this work, we use structure-based models (SBMs) and atomistic simulations to comprehensively investigate how BdpA folding is associated with the formation of its secondary structure. The energy landscape visualization method (ELViM) was used to characterize the pathways that connect the folded and unfolded states of BdpA as well as the sets of structures displaying specific ellipticity patterns. We show that the native state conformational diversity is due mainly to the conformational variability of helix I. Helices I, II, and III occur in a weakly correlated manner, with Spearman's rank correlation coefficients of 0.1539 (I and II), 0.1259 (I and III), and 0.2561 (II and III). These results, therefore, suggest the highest cooperativity between helices II and III. Our results allow the clustering of partially folded structures of folding of the B domain of protein A on the basis of its secondary structure, paving the way to an understanding of environmental factors in the relative stability of the basins of the folding ensemble, which are illustrated by the structural dependency of the protein hydration structures, as computed with minimum-distance distribution functions.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔG_u°'(H₂O), of UBA(1) is 2.4 kcal mol^-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, k_f, in the absence of a denaturant of 13,000 s^-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone ¹⁵NH, ¹³CO, and ¹³Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Mutational Analysis of Protein Folding Transition States: Phi Values

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Simultaneous Determination of Two Subdomain Folding Rates Using the "Transfer-Quench" Method

The investigation of the mechanism of protein folding is complicated by the context dependence of the rates of intramolecular contact formation. Methods based on site-specific labeling and ultrafast spectroscopic detection of fluorescence signals were developed for monitoring the rates of individual subdomain folding transitions in situ, in the context of the whole molecule. However, each site-specific labeling modification might affect rates of folding of near-neighbor structural elements, and thus limit the ability to resolve fine differences in rates of folding of these elements. Therefore, it is highly desirable to be able to study the rates of folding of two or more neighboring subdomain structures using a single mutant to facilitate resolution of the order and interdependence of such steps. Here, we report the development of the "Transfer-Quench" method for measuring the rate of formation of two structural elements using a single triple-labeled mutant. This method is based on Förster resonance energy transfer combined with fluorescence quenching. We placed the donor and acceptor at the loop ends, and a quencher at an α-helical element involved in the node forming the loop. The folding of the triple-labeled mutant is monitored by the acceptor emission. The formation of nonlocal contact (loop closure) increases the time-dependent acceptor emission, while the closure of the labeled helix turn reduces this emission. The method was applied in a study of the folding mechanism of the common model protein, the B domain of staphylococcal protein A. Only natural amino acids were used as probes, and thus possible structural perturbations were minimized. Tyr and Trp residues served as donor and acceptor at the ends of a long loop between helices I and II, and a Cys residue as a quencher for the acceptor. We found that the closure of the loop (segment 14-33) occurs with the same rate constant as the nucleation of helix HII (segment 33-29), in line with the nucleation-condensation model.

Packing energetics determine the folding routes of the RNase-H proteins

The RNase-H proteins show a diverse range of folding routes with structurally distinct folding nuclei.

Complexity of the folding transition of the B domain of protein A revealed by the high-speed tracking of single-molecule fluorescence time series

The equilibrium unfolding transition of the B domain of protein A (BdpA) was investigated by using single-molecule fluorescence spectroscopy based on line-confocal detection of fast-flowing samples. The method achieved the time resolution of 120 μs and the observation time of a few milliseconds in the single-molecule time-series measurements of fluorescence resonance energy transfer (FRET). Two samples of BdpA doubly labeled with donor and acceptor fluorophores, the first possessing fluorophores at residues 22 and 55 (sample 1) and the second at residues 5 and 55 (sample 2), were prepared. The equilibrium unfolding transition induced by guanidium chloride (GdmCl) was monitored by bulk measurements and demonstrated that the both samples obey the apparent two-state unfolding. In the absence of GdmCl, the single-molecule FRET measurements for the both samples showed a single peak assignable to the native state (N). The FRET efficiency for N shifts to lower values as the increase of GdmCl concentration, suggesting the swelling of the native state structure. At the higher concentration of GdmCl, the both samples convert to the unfolded state (U). Near the unfolding midpoint for sample 1, the kinetic exchange between N and U causes the averaging of the two states and the higher values of the relative fluctuation. The time series for different molecules in U showed slightly different FRET efficiencies, suggesting the apparent heterogeneity. Thus, the high-speed tracking of fluorescence signals from single molecules revealed a complexity and heterogeneity hidden in the apparent two-state behavior of protein folding.

Protein folding as a complex reaction: a two-component potential for the driving force of folding and its variation with folding scenario

The Helmholtz decomposition of the vector field of probability fluxes in a two-dimensional space of collective variables makes it possible to introduce a potential for the driving force of protein folding [Chekmarev, J. Chem. Phys. 139 (2013) 145103]. The potential has two components: one component (Φ) is responsible for the source and sink of the folding flow, which represent, respectively, the unfolded and native state of the protein, and the other (Ψ) accounts for the flow vorticity inherently generated at the periphery of the flow field and provides the canalization of the flow between the source and sink. Both components obey Poisson’s equations with the corresponding source/sink terms. In the present paper, we consider how the shape of the potential changes depending on the scenario of protein folding. To mimic protein folding dynamics projected onto a two-dimensional space of collective variables, the two-dimensional Müller and Brown potential is employed. Three characteristic scenarios are considered: a single pathway from the unfolded to the native state without intermediates, two parallel pathways without intermediates, and a single pathway with an off-pathway intermediate. To determine the probability fluxes, the hydrodynamic description of the folding reaction is used, in which the first-passage folding is viewed as a steady flow of the representative points of the protein from the unfolded to the native state. We show that despite the possible complexity of the folding process, the Φ-component is simple and universal in shape. The Ψ-component is more complex and reveals characteristic features of the process of folding. The present approach is potentially applicable to other complex reactions, for which the transition from the reactant to the product can be described in a space of two (collective) variables.

Fast helix formation in the B domain of protein A revealed by site-specific infrared probes

Comparison of experimental and computational protein folding studies can be difficult because of differences in structural resolution. Isotope-edited infrared spectroscopy offers a direct measure of structural changes involved in protein folding at the single-residue level. Here we demonstrate the increased resolution of site-specific infrared probes to the peptide backbone in the B domain of staphylococcal protein A (BdpA). (13)C═(18)O-labeled methionine was incorporated into each of the helices using recombinant protein expression. Laser-induced temperature jumps coupled with infrared spectroscopy were used to probe changes in the peptide backbone on the submillisecond time scale. The relaxation kinetics of the buried helices, solvated helices, and labeled positions were measured independently by probing the corresponding bands assigned in the amide I region. Using these wavelength-dependent measurements, we observe a fast nanosecond phase and slower microsecond phase at each position. We find at least partial formation of helices 1-3 in the fast intermediate state that precedes the transition state. These measurements provide direct, time-resolved experimental evidence of the early formation of partial helical structure in helices 1 and 3, supporting folding models proposed by computer simulations.

The folding of a family of three-helix bundle proteins: spectrin R15 has a robust folding nucleus, unlike its homologous neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.

Engineered solubility tag for solution NMR of proteins

The low solubility of many proteins hinders large scale expression and purification as well as biophysical measurements. Here, we devised a general strategy to solubilize a protein by conjugating it at a solvent-exposed position to a 6 kDa protein that was re-engineered to be highly soluble. We applied this method to the CARD domain of Apoptosis-associated speck-like protein containing a CARD (ASC), which represents one member of a class of proteins that are notoriously prone to aggregation. Attachment of the tag to a cysteine residue, introduced by site-directed mutagenesis at its self-association interface, improved the solubility of the ASC CARD over 50-fold under physiological conditions. Although it is not possible to use nuclear magnetic resonance (NMR) to obtain a high quality 2D correlation spectrum of the wild type domain under physiological conditions, we demonstrate that NMR relaxation parameters of the solubilized variant are sufficiently improved to facilitate virtually any demanding measurement. The method shown here represents a straightforward approach for dramatically increasing protein solubility, enabled by ease of labeling as well as flexibility in tag placement with minimal perturbation to the target.

Structure-based histidine substitution for optimizing pH-sensitive Staphylococcus protein A

Optimizing antibody purification is crucial to overcoming a bottleneck in the costly manufacturing process for antibody therapy. To address this issue, we designed a pH-sensitive Staphylococcus aureus protein A variant that retained its innate stability and affinity toward antibody. On the basis of structural information and mutation analysis data, we identified candidate positions for accumulative histidine substitutions to cause electrostatic repulsion under acidic conditions. The histidine substitutions effectively decreased the dissociation rate under acidic conditions by three orders of magnitude. Avoiding deleterious effects of the substitutions, we successfully engineered a protein A variant that exhibited high pH sensitivity and maintained affinity, thermal stability, and alkaline tolerance. The variant was capable of serving as an affinity ligand that made affinity chromatography under milder acidic conditions possible; the elution peak shifted from pH 4.2 to 5.6. Only two substitutions were needed to achieve this pH sensitivity. This structure-based approach is applicable to other protein-based ligands.

Take home lessons from studies of related proteins

The 'Fold Approach' involves a detailed analysis of the folding of several topologically, structurally and/or evolutionarily related proteins. Such studies can reveal determinants of the folding mechanism beyond the gross topology, and can dissect the residues required for folding from those required for stability or function. While this approach has not yet matured to the point where we can predict the native conformation of any polypeptide chain in silico, it has been able to highlight, amongst others, the specific residues that are responsible for nucleation, pathway malleability, kinetic intermediates, chain knotting, internal friction and Paracelsus switches. Some of the most interesting discoveries have resulted from the attempt to explain differences between homologues.

From coarse-grained to atomic-level characterization of protein dynamics: transition state for the folding of B domain of protein A

Atomic-level molecular dynamics simulations are widely used for the characterization of the structural dynamics of proteins; however, they are limited to shorter time scales than the duration of most of the relevant biological processes. Properly designed coarse-grained models that trade atomic resolution for efficient sampling allow access to much longer time-scales. In-depth understanding of the structural dynamics, however, must involve atomic details. In this study, we tested a method for the rapid reconstruction of all-atom models from α carbon atom positions in the application to convert a coarse-grained folding trajectory of a well described model system: the B domain of protein A. The results show that the method and the spatial resolution of the resulting coarse-grained models enable computationally inexpensive reconstruction of realistic all-atom models. Additionally, by means of structural clustering, we determined the most persistent ensembles of the key folding step, the transition state. Importantly, the analysis of the overall structural topologies suggests a dominant folding pathway. This, together with the all-atom characterization of the obtained ensembles, in the form of contact maps, matches the experimental results well.

The folding transition state of protein L is extensive with nonnative interactions (and not small and polarized)

Progress in understanding protein folding relies heavily upon an interplay between experiment and theory. In particular, readily interpretable experimental data that can be meaningfully compared to simulations are required. According to standard mutational ϕ analysis, the transition state for Protein L contains only a single hairpin. However, we demonstrate here using ψ analysis with engineered metal ion binding sites that the transition state is extensive, containing the entire four-stranded β sheet. Underreporting of the structural content of the transition state by ϕ analysis also occurs for acyl phosphatase [Pandit, A. D., Jha, A., Freed, K. F. & Sosnick, T. R., (2006). Small proteins fold through transition states with native-like topologies. J. Mol. Biol.361, 755-770], ubiquitin [Sosnick, T. R., Dothager, R. S. & Krantz, B. A., (2004). Differences in the folding transition state of ubiquitin indicated by ϕ and ψ analyses. Proc. Natl Acad. Sci. USA 101, 17377-17382] and BdpA [Baxa, M., Freed, K. F. & Sosnick, T. R., (2008). Quantifying the structural requirements of the folding transition state of protein A and other systems. J. Mol. Biol.381, 1362-1381]. The carboxy-terminal hairpin in the transition state of Protein L is found to be nonnative, a significant result that agrees with our Protein Data Bank-based backbone sampling and all-atom simulations. The nonnative character partially explains the failure of accepted experimental and native-centric computational approaches to adequately describe the transition state. Hence, caution is required even when an apparent agreement exists between experiment and theory, thus highlighting the importance of having alternative methods for characterizing transition states.

Combination of Markov state models and kinetic networks for the analysis of molecular dynamics simulations of peptide folding

Atomistic molecular dynamics simulations of the TZ1 beta-hairpin peptide have been carried out using an implicit model for the solvent. The trajectories have been analyzed using a Markov state model defined on the projections along two significant observables and a kinetic network approach. The Markov state model allowed for an unbiased identification of the metastable states of the system, and provided the basis for commitment probability calculations performed on the kinetic network. The kinetic network analysis served to extract the main transition state for folding of the peptide and to validate the results from the Markov state analysis. The combination of the two techniques allowed for a consistent and concise characterization of the dynamics of the peptide. The slowest relaxation process identified is the exchange between variably folded and denatured species, and the second slowest process is the exchange between two different subsets of the denatured state which could not be otherwise identified by simple inspection of the projected trajectory. The third slowest process is the exchange between a fully native and a partially folded intermediate state characterized by a native turn with a proximal backbone H-bond, and frayed side-chain packing and termini. The transition state for the main folding reaction is similar to the intermediate state, although a more native like side-chain packing is observed.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Induced fit or conformational selection for RNA/U1A folding

The hairpin II of U1 snRNA can bind U1A protein with high affinity and specificity. NMR spectra suggest that the loop region of apo-RNA is largely unstructured and undergoes a transition from unstructured to well-folded upon U1Abinding. However, the mechanism that RNA folding coupled protein binding is poorly understood. To get an insight into the mechanism, we have performed explicit-solvent molecular dynamics (MD) to study the folding kinetics of bound RNA and apo-RNA. Room-temperature MD simulations suggest that the conformation of bound RNA has significant adjustment and becomes more stable upon U1A binding. Kinetic analysis of high-temperature MD simulations shows that bound RNA and apo-RNA unfold via a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound RNA folds in the order of RNA contracting, U1A binding, and tertiary folding. The predicted Phi-values suggest that A8, C10, A11, and G16 are key bases for bound RNA folding. Mutant Arg52Gln analysis shows that electrostatic interaction and hydrogen bonds between RNA and U1A (Arg52Gln) decrease. These results are in qualitative agreement with experiments. Furthermore, this method could be used in other studies about biomolecule folding upon receptor binding.

Nearly symmetrical proteins: folding pathways and transition states

The folding pathways of the B domain of protein A have been the subject of many experimental and computational studies. Based on a statistical mechanical model, it has been suggested that the native state symmetry leads to multiple pathways, highly dependent on temperature and denaturant concentration. Experiments, however, have not confirmed this scenario. By considering four nearly symmetrical proteins, one of them being the above molecule, here we show that, if contact energies are properly taken into account, a different picture emerges from kinetic simulations of the above-mentioned model. This is characterized by a dominant folding pathway, which is consistent with the most recent experimental results. Given the simplicity of the model, we also report on a direct sampling of the transition state.

Molecular dynamics simulation of phosphorylated KID post-translational modification

Background

Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood.

Methodology

To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding.

Conclusions

Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix α_B, folding of helix α_A, completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices α_A and α_B are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Φ-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification.

How general is the nucleation-condensation mechanism?

We investigate the structures of the major folding transition states of nine proteins by correlation of published Phi-values with inter-residue contact maps. Combined with previous studies on six proteins, the analysis suggests that at least 10 of the 15 small globular proteins fold via a nucleation-condensation mechanism with a concurrent build-up of secondary and tertiary structure contacts, but a structural consolidation that is clearly nonuniformly distributed over the molecule and most intense in a single structural region suggesting the occurrence of a single folding nucleus. However, on average helix- and sheet-forming residues show somewhat larger Phi-values in the major transition state, suggesting that secondary structure formation is one important driving force in the nucleation-condensation in many proteins and that secondary-structure forming residues tend to be more prominent in folding nuclei. We synthesize the combined information on these 10 of 15 proteins into a unified nucleation-condensation mechanism which also accounts for effects described by the framework, hydrophobic collapse, zipper, and funnel models.

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.

Analysis of 3D structural differences in the IgG-binding domains based on the interresidue average-distance statistics

It is well-known that the IgG-binding domain from staphylococcal protein A folds into a 3 alpha helix bundle structure, while the IgG-binding domain of streptococcal protein G forms an (alpha + beta) structure. Recently, He et al. (Biochemistry 44:14055-14061, 2005) made mutants of these proteins from the wild types of protein A and protein G strains. These mutants are referred to as protein A219 and protein G311, and it was showed that these two mutants have different 3D structures, i.e., the 3 alpha helix bundle structure and the (alpha + beta) structure, respectively, despite the high sequence identity (59%). The purpose of our study was to clarify how such 3D structural differences are coded in the sequences with high homology. To address this problem, we introduce a predicted contact map constructed based on the interresidue average-distance statistics for prediction of folding properties of a protein. We refer to this map as an average distance map (ADM). Furthermore, the statistics of interresidue distances can be converted to an effective interresidue potential. We calculated the contact frequency of each residue of a protein in random conformations with this effective interresidue potential, and then we obtained values similar to phi values. We refer to this contact frequency of each residue as a p(mu) value. The comparison of the p(mu) values to the phi values for a protein suggests that p(mu) values reveal the information on the folding initiation site. Using these techniques, we try to extract the information on the difference in the 3D structures of protein A219 and protein G311 coded in their amino acid sequences in the present work. The results show that the ADM analyses and the p(mu) value analyses predict the information of folding initiation sites, which can be used to detect the 3D difference in both proteins.

Folding domain B of protein A on a dynamically partitioned free energy landscape

The B domain of staphylococcal protein A (BdpA) is a small helical protein that has been studied intensively in kinetics experiments and detailed computer simulations that include explicit water. The simulations indicate that BdpA needs to reorganize in crossing the transition barrier to facilitate folding its C-terminal helix (H3) onto the nucleus formed from helices H1 and H2. This process suggests frustration between two partially ordered forms of the protein, but recent varphi value measurements indicate that the transition structure is relatively constant over a broad range of temperatures. Here we develop a simplistic model to investigate the folding transition in which properties of the free energy landscape can be quantitatively compared with experimental data. The model is a continuation of the Muñoz-Eaton model to include the intermittency of contacts between structured parts of the protein, and the results compare variations in the landscape with denaturant and temperature to varphi value measurements and chevron plots of the kinetic rates. The topography of the model landscape (in particular, the feature of frustration) is consistent with detailed simulations even though variations in the varphi values are close to measured values. The transition barrier is smaller than indicated by the chevron data, but it agrees in order of magnitude with a similar alpha-carbon type of model. Discrepancies with the chevron plots are investigated from the point of view of solvent effects, and an approach is suggested to account for solvent participation in the model.