Specificity of RSG-1.2 peptide binding to RRE-IIB RNA element of HIV-1 over Rev peptide is mainly enthalpic in origin

Rev is an essential HIV-1 regulatory protein which binds to the Rev responsive element (RRE) present within the env gene of HIV-1 RNA genome. This binding facilitates the transport of the RNA to the cytoplasm, which in turn triggers the switch between viral latency and active viral replication. Essential components of this complex have been localized to a minimal arginine rich Rev peptide and stem IIB region of RRE. A synthetic peptide known as RSG-1.2 binds with high binding affinity and specificity to the RRE-IIB than the Rev peptide, however the thermodynamic basis of this specificity has not yet been addressed. The present study aims to probe the thermodynamic origin of this specificity of RSG-1.2 over Rev Peptide for RRE-IIB. The temperature dependent melting studies show that RSG-1.2 binding stabilizes the RRE structure significantly (ΔT _m = 4.3°C), in contrast to Rev binding. Interestingly the thermodynamic signatures of the binding have also been found to be different for both the peptides. At pH 7.5, RSG-1.2 binds RRE-IIB with a K_a = 16.2±0.6×10⁷ M⁻¹ where enthalpic change ΔH = −13.9±0.1 kcal/mol is the main driving force with limited unfavorable contribution from entropic change TΔS = −2.8±0.1 kcal/mol. A large part of ΔH may be due to specific stacking between U72 and Arg15. In contrast binding of Rev (K_a = 3.1±0.4×10⁷ M⁻¹) is driven mainly by entropy (ΔH = 0 kcal/mol and TΔS = 10.2±0.2 kcal/mol) which arises from major conformational changes in the RNA upon binding.     

Thermodynamic and kinetic basis for recognition and repair of 8-oxoguanine in DNA by human 8-oxoguanine-DNA glycosylase

We have used a stepwise increase in ligand complexity approach to estimate the relative contributions of the nucleotide units of DNA containing 7,8-dihydro-8-oxoguanine (oxoG) to its total affinity for human 8-oxoguanine DNA glycosylase (OGG1) and construct thermodynamic models of the enzyme interaction with cognate and non-cognate DNA. Non-specific OGG1 interactions with 10–13 nt pairs within its DNA-binding cleft provides approximately 5 orders of magnitude of its affinity for DNA (ΔG° approximately −6.7 kcal/mol). The relative contribution of the oxoG unit of DNA (ΔG° approximately −3.3 kcal/mol) together with other specific interactions (ΔG° approximately −0.7 kcal/mol) provide approximately 3 orders of magnitude of the affinity. Formation of the Michaelis complex of OGG1 with the cognate DNA cannot account for the major part of the enzyme specificity, which lies in the k_cat term instead; the rate increases by 6–7 orders of magnitude for cognate DNA as compared with non-cognate one. The k_cat values for substrates of different sequences correlate with the DNA twist, while the K_M values correlate with ΔG° of the DNA fragments surrounding the lesion (position from −6 to +6). The functions for predicting the K_M and k_cat values for different sequences containing oxoG were found.     

AB‐Bind: Antibody binding mutational database for computational affinity predictions

Antibodies (Abs) are a crucial component of the immune system and are often used as diagnostic and therapeutic agents. The need for high‐affinity and high‐specificity antibodies in research and medicine is driving the development of computational tools for accelerating antibody design and discovery. We report a diverse set of antibody binding data with accompanying structures that can be used to evaluate methods for modeling antibody interactions. Our Antibody‐Bind (AB‐Bind) database includes 1101 mutants with experimentally determined changes in binding free energies (ΔΔG) across 32 complexes. Using the AB‐Bind data set, we evaluated the performance of protein scoring potentials in their ability to predict changes in binding free energies upon mutagenesis. Numerical correlations between computed and observed ΔΔG values were low (r = 0.16–0.45), but the potentials exhibited predictive power for classifying variants as improved vs weakened binders. Performance was evaluated using the area under the curve (AUC) for receiver operator characteristic (ROC) curves; the highest AUC values for 527 mutants with |ΔΔG| > 1.0 kcal/mol were 0.81, 0.87, and 0.88 using STATIUM, FoldX, and Discovery Studio scoring potentials, respectively. Some methods could also enrich for variants with improved binding affinity; FoldX and Discovery Studio were able to correctly rank 42% and 30%, respectively, of the 80 most improved binders (those with ΔΔG < −1.0 kcal/mol) in the top 5% of the database. This modest predictive performance has value but demonstrates the continuing need to develop and improve protein energy functions for affinity prediction.     

Characterization of the kinetic and thermodynamic landscape of RNA folding using a novel application of isothermal titration calorimetry

A novel isothermal titration calorimetry (ITC) method was applied to investigate RNA helical packing driven by the GAAA tetraloop–receptor interaction in magnesium and potassium solutions. Both the kinetics and thermodynamics were obtained in individual ITC experiments, and analysis of the kinetic data over a range of temperatures provided Arrhenius activation energies (ΔH^‡) and Eyring transition state entropies (ΔS^‡). The resulting rich dataset reveals strongly contrasting kinetic and thermodynamic profiles for this RNA folding system when stabilized by potassium versus magnesium. In potassium, association is highly exothermic (ΔH_25°C = −41.6 ± 1.2 kcal/mol in 150 mM KCl) and the transition state is enthalpically barrierless (ΔH^‡ = −0.6 ± 0.5). These parameters are sigificantly positively shifted in magnesium (ΔH_25°C = −20.5 ± 2.1 kcal/mol, ΔH^‡ = 7.3 ± 2.2 kcal/mol in 0.5 mM MgCl₂). Mixed salt solutions approximating physiological conditions exhibit an intermediate thermodynamic character. The cation-dependent thermodynamic landscape may reflect either a salt-dependent unbound receptor conformation, or alternatively and more generally, it may reflect a small per-cation enthalpic penalty associated with folding-coupled magnesium uptake.     

Thermodynamic examination of the pyrophosphate sensor helix in the thiamine pyrophosphate riboswitch

Riboswitches are functional mRNA that control gene expression. Thiamine pyrophosphate (TPP) binds to thi-box riboswitch RNA and allosterically inhibits genes that code for proteins involved in the biosynthesis and transport of thiamine. Thiamine binding to the pyrimidine sensor helix and pyrophosphate binding to the pyrophosphate sensor helix cause changes in RNA conformation that regulate gene expression. Here we examine the thermodynamic properties of the internal loop of the pyrophosphate binding domain by comparing the wild-type construct (RNA WT) with six modified 2 × 2 bulged RNA and one 2 × 2 bulged DNA. The wild-type construct retains five conserved bases of the pyrophosphate sensor domain, two of which are in the 2 × 2 bulge (C65 and G66). The RNA WT construct was among the most stable (ΔG°₃₇ = −7.7 kcal/mol) in 1 M KCl at pH 7.5. Breaking the A•G mismatch of the bulge decreases the stability of the construct ∼0.5–1 kcal/mol, but does not affect magnesium binding to the RNA WT. Guanine at position 48 is important for RNA–Mg²⁺ interactions of the TPP-binding riboswitch at pH 7.5. In the presence of 9.5 mM magnesium at pH 5.5, the bulged RNA constructs gained an average of 1.1 kcal/mol relative to 1 M salt. Formation of a single A+•C mismatch base pair contributes about 0.5 kcal/mol at pH 5.5, whereas two tandem A+•C mismatch base pairs together contribute about 2 kcal/mol.     

Cold denaturation of monoclonal antibodies

The susceptibility of monoclonal antibodies (mAbs) to undergo cold denaturation remains unexplored. In this study, the phenomenon of cold denaturation was investigated for a mAb, mAb1, through thermodynamic and spectroscopic analyses. tryptophan fluorescence and circular dichroism (CD) spectra were recorded for the guanidine hydrochloride (GuHCl)-induced unfolding of mAb1 at pH 6.3 at temperatures ranging from −5 to 50°C. A three-state unfolding model incorporating the linear extrapolation method was fit to the fluorescence data to obtain an apparent free energy of unfolding, ΔG_u, at each temperature. CD studies revealed that mAb1 exhibited polyproline II helical structure at low temperatures and at high GuHCl concentrations. the Gibbs-Helmholtz expression fit to the ΔG_u versus temperature data from fluorescence gave a ΔC_p of 8.0 kcal mol⁻¹ K⁻¹, a maximum apparent stability of 23.7 kcal mol⁻¹ at 18°C, and an apparent cold denaturation temperature (T_CD) of −23°C. ΔG_u values for another mAb (mAb2) with a similar framework exhibited less stability at low temperatures, suggesting a depressed protein stability curve and a higher relative T_CD. Direct experimental evidence of the susceptibility of mAb1 and mAb2 to undergo cold denaturation in the absence of denaturant was confirmed at pH 2.5. thus, mAbs have a potential to undergo cold denaturation at storage temperatures near −20°C (pH 6.3), and this potential needs to be evaluated independently for individual mAbs.Key words: monoclonal antibodies, thermodynamic stability, cold denaturation, free energy, fluorescence     

Calorimetric and spectroscopic investigation of the interaction between the C‐terminal domain of Enzyme I and its ligands

Enzyme I initiates a series of phosphotransfer reactions during sugar uptake in the bacterial phosphotransferase system. Here, we have isolated a stable recombinant C-terminal domain of Enzyme I (EIC) of Escherichia coli and characterized its interaction with the N-terminal domain of Enzyme I (EIN) and also with various ligands. EIC can phosphorylate EIN, but their binding is transient regardless of the presence of phosphoenolpyruvate (PEP). Circular dichroism and NMR indicate that ligand binding to EIC induces changes near aromatic groups but not in the secondary structure of EIC. Binding of PEP to EIC is an endothermic reaction with the equilibrium dissociation constant (K_D) of 0.28 mM, whereas binding of the inhibitor oxalate is an exothermic reaction with K_D of 0.66 mM from calorimetry. The binding thermodynamics of EIC and PEP compared to that of Enzyme I (EI) and PEP reveals that domain–domain motion in EI can contribute as large as ∼−3.2 kcal/mol toward PEP binding.     

Energetics of echinomycin binding to DNA

Differential scanning calorimetry and UV thermal denaturation have been used to determine a complete thermodynamic profile for the bis-intercalative interaction of the peptide antibiotic echinomycin with DNA. The new calorimetric data are consistent with all previously published binding data, and afford the most rigorous and direct determination of the binding enthalpy possible. For the association of echinomycin with DNA, we found ΔG° = –7.6 kcal mol^–1, ΔH = +3.8 kcal mol^–1 and ΔS = +38.9 cal mol^–1 K^–1 at 20°C. The binding reaction is clearly entropically driven, a hallmark of a process that is predominantly stabilized by hydrophobic interactions, though a deeper analysis of the free energy contributions suggests that direct molecular recognition between echinomycin and DNA, mediated by hydrogen bonding and van der Waals contacts, also plays an important role in stabilizing the complex.     

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ∼16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm_app) by 1–6°C resulting in an increase in free energy of unfolding (ΔG_unf) by 0.04–1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG_unf, (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG_unf = 1.57 kcal⁻¹) justifies the robustness of this relation.     

Protein folding rates and thermodynamic stability are key determinants for interaction with the Hsp70 chaperone system

The Hsp70 family of molecular chaperones participates in vital cellular processes including the heat shock response and protein homeostasis. E. coli''s Hsp70, known as DnaK, works in concert with the DnaJ and GrpE co-chaperones (K/J/E chaperone system), and mediates cotranslational and post-translational protein folding in the cytoplasm. While the role of the K/J/E chaperones is well understood in the presence of large substrates unable to fold independently, it is not known if and how K/J/E modulates the folding of smaller proteins able to fold even in the absence of chaperones. Here, we combine experiments and computation to evaluate the significance of kinetic partitioning as a model to describe the interplay between protein folding and binding to the K/J/E chaperone system. First, we target three nonobligatory substrates, that is, proteins that do not require chaperones to fold. The experimentally observed chaperone association of these client proteins during folding is entirely consistent with predictions from kinetic partitioning. Next, we develop and validate a computational model (CHAMP70) that assumes kinetic partitioning of substrates between folding and interaction with K/J/E. CHAMP70 quantitatively predicts the experimentally measured interaction of RNase H^D as it refolds in the presence of various chaperones. CHAMP70 shows that substrates are posed to interact with K/J/E only if they are slow-folding proteins with a folding rate constant k_f <50 s⁻¹, and/or thermodynamically unstable proteins with a folding free energy ΔG⁰_UN ≥−2 kcal mol⁻¹. Hence, the K/J/E system is tuned to use specific protein folding rates and thermodynamic stabilities as substrate selection criteria.     

Compressing the free energy range of substructure stabilities in iso‐1‐cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 ± 0.4 kcal mol⁻¹ to unfold the least stable substructure compared with 5.05 ± 0.30 kcal mol⁻¹ for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only ∼1.2 kcal mol⁻¹ compared with ∼8 kcal mol⁻¹ for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are ∼3 kcal mol⁻¹M⁻¹ versus ∼4.5 kcal mol⁻¹M⁻¹ for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.     

Temperature Dependence of Nonelectrolyte Permeation across Red Cell Membranes

The temperature dependence of permeation across human red cell membranes has been determined for a series of hydrophilic and lipophilic solutes, including urea and two methyl substituted derivatives, all the straight-chain amides from formamide through valeramide and the two isomers, isobutyramide and isovaleramide. The temperature coefficient for permeation by all the hydrophilic solutes is 12 kcal mol^-1 or less, whereas that for all the lipophilic solutes is 19 kcal mol^-1 or greater. This difference is consonant with the view that hydrophilic molecules cross the membrane by a path different from that taken by the lipophilic ones. The thermodynamic parameters associated with lipophile permeation have been studied in detail. ΔG is negative for adsorption of lipophilic amides onto an oil-water interface, whereas it is positive for transfer of the polar head from the aqueous medium to bulk lipid solvent. Application of absolute reaction rate theory makes it possible to make a clear distinction between diffusion across the water-red cell membrane interface and diffusion within the membrane. Diffusion coefficients and apparent activation enthalpies and entropies have been computed for each process. Transfer of the polar head from the solvent into the interface is characterized by ΔG^‡ = 0 kcal mol^-1 and ΔS^‡ negative, whereas both of these parameters have large positive values for diffusion within the membrane. Diffusion within the membrane is similar to what is expected for diffusion through a highly associated viscous fluid.     

Thermodynamic analysis of F1‐ATPase rotary catalysis using high‐speed imaging

F₁-ATPase (F₁) is a rotary motor protein fueled by ATP hydrolysis. Although the mechanism for coupling rotation and catalysis has been well studied, the molecular details of individual reaction steps remain elusive. In this study, we performed high-speed imaging of F₁ rotation at various temperatures using the total internal reflection dark-field (TIRDF) illumination system, which allows resolution of the F₁ catalytic reaction into elementary reaction steps with a high temporal resolution of 72 µs. At a high concentration of ATP, F₁ rotation comprised distinct 80° and 40° substeps. The 80° substep, which exhibited significant temperature dependence, is triggered by the temperature-sensitive reaction, whereas the 40° substep is triggered by ATP hydrolysis and the release of inorganic phosphate (P_i). Then, we conducted Arrhenius analysis of the reaction rates to obtain the thermodynamic parameters for individual reaction steps, that is, ATP binding, ATP hydrolysis, P_i release, and TS reaction. Although all reaction steps exhibited similar activation free energy values, ΔG^‡ = 53–56 kJ mol⁻¹, the contributions of the enthalpy (ΔH^‡), and entropy (ΔS^‡) terms were significantly different; the reaction steps that induce tight subunit packing, for example, ATP binding and TS reaction, showed high positive values of both ΔH^‡ and ΔS^‡. The results may reflect modulation of the excluded volume as a function of subunit packing tightness at individual reaction steps, leading to a gain or loss in water entropy.     

Distal Histidine Stabilizes Bound O2 and Acts as a Gate for Ligand Entry in Both Subunits of Adult Human Hemoglobin

The role of the distal histidine in regulating ligand binding to adult human hemoglobin (HbA) was re-examined systematically by preparing His(E7) to Gly, Ala, Leu, Gln, Phe, and Trp mutants of both Hb subunits. Rate constants for O₂, CO, and NO binding were measured using rapid mixing and laser photolysis experiments designed to minimize autoxidation of the unstable apolar E7 mutants. Replacing His(E7) with Gly, Ala, Leu, or Phe causes 20–500-fold increases in the rates of O₂ dissociation from either Hb subunit, demonstrating unambiguously that the native His(E7) imidazole side chain forms a strong hydrogen bond with bound O₂ in both the α and β chains (ΔG_{His(E7)H-bond} ≈ −8 kJ/mol). As the size of the E7 amino acid is increased from Gly to Phe, decreases in k_O2′, k_NO′, and calculated bimolecular rates of CO entry (k_entry′) are observed. Replacing His(E7) with Trp causes further decreases in k_O2′, k_NO′, and k_entry′ to 1–2 μm⁻¹ s⁻¹ in β subunits, whereas ligand rebinding to αTrp(E7) subunits after photolysis is markedly biphasic, with fast k_O2′, k_CO′, and k_NO′ values ≈150 μm⁻¹ s⁻¹ and slow rate constants ≈0.1 to 1 μm⁻¹ s⁻¹. Rapid bimolecular rebinding to an open α subunit conformation occurs immediately after photolysis of the αTrp(E7) mutant at high ligand concentrations. However, at equilibrium the closed αTrp(E7) side chain inhibits the rate of ligand binding >200-fold. These data suggest strongly that the E7 side chain functions as a gate for ligand entry in both HbA subunits.     

Toward Classification of BRCA1 Missense Variants Using a Biophysical Approach

Carriers of germ line mutations in breast cancer susceptibility gene BRCA1 have an increased risk of developing breast and ovarian cancers; missense mutations have, however, been difficult to assess for disease association. Here we have used a biophysical approach to classify these variants. We established an assay for measuring the thermodynamic stability of the BRCA1 BRCT domains and investigated the effects of 36 missense mutations. The mutations show a range of effects. Some do not change the stability, whereas others destabilize the protein by as much as 6 kcal mol⁻¹; one-third of the mutants could not be expressed in soluble form in Escherichia coli, and we conclude that these destabilize the protein by an even greater amount. We tested several computer algorithms for their ability to predict the mutant effects and found that by grouping them into two classes (destabilizing by less than or more than 2.2 kcal mol⁻¹), the algorithms could predict the stability changes. Importantly, with the exception of the few mutants located in the binding site, none showed a significant reduction in affinity for phosphorylated substrate. These results indicate that despite very large losses in stability, the integrity of the structure is not compromised by the mutations. Thus, the majority of mutations cause loss of function by reducing the proportion of BRCA1 molecules that are in the folded state and increasing the proportion of molecules that are unfolded. Consequently, small molecule stabilization of the structure could be a generally applicable preventative therapeutic strategy for rescuing many BRCA1 mutations.     

Balanced biochemical reactions: a new approach to unify chemical and biochemical thermodynamics

A novel procedure is presented which, by balancing elements and electric charge of biochemical reactions which occur at constant pH and pMg, allows assessing the thermodynamics properties of reaction Δ_rG ^′0, Δ_rH ^′0, Δ_rS ^′0 and the change in binding of hydrogen and magnesium ions of these reactions. This procedure of general applicability avoids the complex calculations required by the use of the Legendre transformed thermodynamic properties of formation Δ_fG ^′0, Δ_fH ^′0 and Δ_fS ^′0 hitherto considered an obligatory prerequisite to deal with the thermodynamics of biochemical reactions. As a consequence, the term “conditional” is proposed in substitution of “Legendre transformed” to indicate these thermodynamics properties. It is also shown that the thermodynamic potential G is fully adequate to give a criterion of spontaneous chemical change for all biochemical reactions and then that the use of the Legendre transformed G′ is unnecessary. The procedure proposed can be applied to any biochemical reaction, making possible to re-unify the two worlds of chemical and biochemical thermodynamics, which so far have been treated separately.     

Thermodynamics of phospholipid self-assembly

Negatively charged phospholipids are an important component of biological membranes. The thermodynamic parameters governing self-assembly of anionic phospholipids are deduced here from isothermal titration calorimetry. Heats of demicellization were determined for dioctanoyl phosphatidylglycerol (PG) and phosphatidylserine (PS) at different ionic strengths, and for dioctanoyl phosphatidic acid at different pH values. The large heat capacity (ΔC^o_P ∼ −400 J.mol⁻¹ K⁻¹ for PG and PS), and zero enthalpy at a characteristic temperature near the physiological range (T^∗ ∼ 300 K for PG and PS), demonstrate that the driving force for self-assembly is the hydrophobic effect. The pH and ionic-strength dependences indicate that the principal electrostatic contribution to self-assembly comes from the entropy associated with the electrostatic double layer, in agreement with theoretical predictions. These measurements help define the thermodynamic effects of anionic lipids on biomembrane stability.     

A delicate interplay of structure, dynamics, and thermodynamics for function: a high pressure NMR study of outer surface protein A

Outer surface protein A (OspA) is a crucial protein in the infection of Borrelia burgdorferi causing Lyme disease. We studied conformational fluctuations of OspA with high-pressure ¹⁵N/¹H two-dimensional NMR along with high-pressure fluorescence spectroscopy. We found evidence within folded, native OspA for rapid local fluctuations of the polypeptide backbone in the nonglobular single layer β-sheet connecting the N- and C-terminal domains with τ << ms, which may give the two domains certain independence in mobility and thermodynamic stability. Furthermore, we found that folded, native OspA is in equilibrium (τ >> ms) with a minor conformer I, which is almost fully disordered and hydrated for the entire C-terminal part of the polypeptide chain from β8 to the C-terminus. Conformer I is characterized with ΔG⁰ = 32 ± 9 kJ/mol and ΔV⁰ = −140 ± 40 mL/mol, populating only ∼0.001% at 40°C at 0.1 MPa, pH 5.9. Because in the folded conformer the receptor binding epitope of OspA is buried in the C-terminal domain, its transition into conformer I under in vivo conditions may be critical for the infection of B. burgdorferi. The formation and stability of the peculiar conformer I are apparently supported by a large packing defect or cavity located in the C-terminal domain.     

Substrate binding modifies the hinge bending characteristics of human 3‐phosphoglycerate kinase: A molecular dynamics study

3‐Phosphogycerate kinase (PGK) is a two domain enzyme, with a binding site of the 1,3‐bisphosphoglycerate on the N‐domain and of the ADP on the C‐domain. To transfer a phosphate group the enzyme has to undergo a hinge bending motion from open to closed conformation to bring the substrates to close proximity. Molecular dynamics simulation was used to elucidate the effect of ligand binding onto the domain motions of this enzyme. The simulation results of the apo form indicate a hinge bending motion in the ns timescale while the time period of the hinge bending motion of the complex form is clearly over the 20 ns simulation time. The apo form exhibits several hinge points that contribute to the hinge bending motion while upon binding the ligands, the hinge bending becomes strictly restrained with one dominant hinge point in the vicinity of the substrates. At the same time, ligand binding results in an enhanced correlation of internal domain motions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Contribution of Residue B5 to the Folding and Function of Insulin and IGF-I: CONSTRAINTS AND FINE-TUNING IN THE EVOLUTION OF A PROTEIN FAMILY*

Proinsulin exhibits a single structure, whereas insulin-like growth factors refold as two disulfide isomers in equilibrium. Native insulin-related growth factor (IGF)-I has canonical cystines (A6—A11, A7–B7, and A20—B19) maintained by IGF-binding proteins; IGF-swap has alternative pairing (A7–A11, A6—B7, and A20—B19) and impaired activity. Studies of mini-domain models suggest that residue B5 (His in insulin and Thr in IGFs) governs the ambiguity or uniqueness of disulfide pairing. Residue B5, a site of mutation in proinsulin causing neonatal diabetes, is thus of broad biophysical interest. Here, we characterize reciprocal B5 substitutions in the two proteins. In insulin, His^B5 → Thr markedly destabilizes the hormone (ΔΔG_u 2.0 ± 0.2 kcal/mol), impairs chain combination, and blocks cellular secretion of proinsulin. The reciprocal IGF-I substitution Thr^B5 → His (residue 4) specifies a unique structure with native ¹H NMR signature. Chemical shifts and nuclear Overhauser effects are similar to those of native IGF-I. Whereas wild-type IGF-I undergoes thiol-catalyzed disulfide exchange to yield IGF-swap, His^B5-IGF-I retains canonical pairing. Chemical denaturation studies indicate that His^B5 does not significantly enhance thermodynamic stability (ΔΔG_u 0.2 ± 0.2 kcal/mol), implying that the substitution favors canonical pairing by destabilizing competing folds. Whereas the activity of Thr^B5-insulin is decreased 5-fold, His^B5-IGF-I exhibits 2-fold increased affinity for the IGF receptor and augmented post-receptor signaling. We propose that conservation of Thr^B5 in IGF-I, rescued from structural ambiguity by IGF-binding proteins, reflects fine-tuning of signal transduction. In contrast, the conservation of His^B5 in insulin highlights its critical role in insulin biosynthesis.