Long-Timescale Molecular Dynamics Simulations Elucidate the Dynamics and Kinetics of Exposure of the Hydrophobic Patch in Troponin C

Troponin (Tn) is an important regulatory protein in the thin-filament complex of cardiomyocytes. Calcium binding to the troponin C (TnC) subunit causes a change in its dynamics that leads to the transient opening of a hydrophobic patch on TnC’s surface, to which a helix of another subunit, troponin I (TnI), binds. This process initiates contraction, making it an important target for studies investigating the detailed molecular processes that underlie contraction. Here we use microsecond-timescale Anton molecular dynamics simulations to investigate the dynamics and kinetics of the opening transition of the TnC hydrophobic patch. Free-energy differences for opening are calculated for wild-type Ca²⁺-bound TnC (∼8 kcal/mol), V44Q Ca²⁺-bound TnC (3.2 kcal/mol), E40A Ca²⁺-bound TnC (∼12 kcal/mol), and wild-type apo TnC (∼20 kcal/mol). These results suggest that the mutations have a profound impact on the frequency with which the hydrophobic patch presents to TnI. In addition, these simulations corroborate that cardiac wild-type TnC does not open on timescales relevant to contraction without calcium being bound.     

Refined structure of chicken skeletal muscle troponin C in the two-calcium state at 2-A resolution

The structure of troponin C has been refined at 2A resolution to an R value of 0.172 using a total of 8,100 reflections. Troponin C has an unusual dumbbell shape with only the two C-domain high affinity sites III and IV occupied with metals, while the pair of N-domain low affinity sites I and II are devoid of metals. The coordination of the Ca2+ approaches seven with the last glutamic acid residue in each site forming an asymmetric bidentate ligand. The flanking helices in the metal-bound EF hands are in similar orientation (both 113 degrees) while in the apo sites they are more obtuse (134 and 149 degrees). The EF hands of holo sites III and IV are similar while the apo sites I and II are less similar (rms for backbone atoms, 0.78 and 1.44). The half-loops of the 12-residue holo and apo sites show better agreement than the full loops themselves, suggesting a hinge motion at the midpoint of the loops. The long central helix is stabilized by electrostatic interactions and salt bridges between charged side chains spaced at 3 or 4 residues along the helix. A cluster of water molecules encircle the long helix and hydrogen bond to the backbone carbonyls. At the beginning of the B-helix, a water molecule is interposed at each of two consecutive backbone NH...OC hydrogen bonds. The terminal pair of helices A/D (apo) match with E/H (holo), and the internal pair of helices B/C (apo) match with F/G (holo). Thus, muscle contraction may be triggered by Ca2+ binding to loops I and II which results in a concerted rearrangement of residues in the loops, including the essential Gly at position 6 in each loop. This rearrangement than causes a reorientation of helices B and C along with the BC linker.     

Parallel measurement of Ca2+ binding and fluorescence emission upon Ca2+ titration of recombinant skeletal muscle troponin C. Measurement of sequential calcium binding to the regulatory sites

Calcium binding to chicken recombinant skeletal muscle TnC (TnC) and its mutants containing tryptophan (F29W), 5-hydroxytryptophan (F29HW), or 7-azatryptophan (F29ZW) at position 29 was measured by flow dialysis and by fluorescence. Comparative analysis of the results allowed us to determine the influence of each amino acid on the calcium binding properties of the N-terminal regulatory domain of the protein. Compared with TnC, the Ca(2+) affinity of N-terminal sites was: 1) increased 6-fold in F29W, 2) increased 3-fold in F29ZW, and 3) decreased slightly in F29HW. The Ca(2+) titration of F29ZW monitored by fluorescence displayed a bimodal curve related to sequential Ca(2+) binding to the two N-terminal Ca(2+) binding sites. Single and double mutants of TnC, F29W, F29HW, and F29ZW were constructed by replacing aspartate by alanine at position 30 (site I) or 66 (site II) or both. Ca(2+) binding data showed that the Asp --> Ala mutation at position 30 impairs calcium binding to site I only, whereas the Asp --> Ala mutation at position 66 impairs calcium binding to both sites I and II. Furthermore, the Asp --> Ala mutation at position 30 eliminates the differences in Ca(2+) affinity observed for replacement of Phe at position 29 by Trp, 5-hydroxytryptophan, or 7-azatryptophan. We conclude that position 29 influences the affinity of site I and that Ca(2+) binding to site I is dependent on the previous binding of metal to site II.     

Structural and functional studies on Troponin I and Troponin C interactions

Troponin I (TnI) peptides (TnI inhibitory peptide residues 104-115, Ip; TnI regulatory peptide resides 1-30, TnI1-30), recombinant Troponin C (TnC) and Troponin I mutants were used to study the structural and functional relationship between TnI and TnC. Our results reveal that an intact central D/E helix in TnC is required to maintain the ability of TnC to release the TnI inhibition of the acto-S1-TM ATPase activity. Ca(2+)-titration of the TnC-TnI1-30 complex was monitored by circular dichroism. The results show that binding of TnI1-30 to TnC caused a three-folded increase in Ca(2+) affinity in the high affinity sites (III and IV) of TnC. Gel electrophoresis and high performance liquid chromatography (HPLC) studies demonstrate that the sequences of the N- and C-terminal regions of TnI interact in an anti-parallel fashion with the corresponding N- and C-domain of TnC. Our results also indicate that the N- and C-terminal domains of TnI which flank the TnI inhibitory region (residues 104 to 115) play a vital role in modulating the Ca(2+)- sensitive release of the TnI inhibitory region by TnC within the muscle filament. A modified schematic diagram of the TnC/TnI interaction is proposed.     

Effect of hydrophobic residue substitutions with glutamine on Ca(2+) binding and exchange with the N-domain of troponin C

Troponin C (TnC) is an EF-hand Ca(2+) binding protein that regulates skeletal muscle contraction. The mechanisms that control the Ca(2+) binding properties of TnC and other EF-hand proteins are not completely understood. We individually substituted 27 Phe, Ile, Leu, Val, and Met residues with polar Gln to examine the role of hydrophobic residues in Ca(2+) binding and exchange with the N-domain of a fluorescent TnC(F29W). The global N-terminal Ca(2+) affinities of the TnC(F29W) mutants varied approximately 2340-fold, while Ca(2+) association and dissociation rates varied less than 70-fold and more than 45-fold, respectively. Greater than 2-fold increases in Ca(2+) affinities were obtained primarily by slowing of Ca(2+) dissociation rates, while greater than 2-fold decreases in Ca(2+) affinities were obtained by slowing of Ca(2+) association rates and speeding of Ca(2+) dissociation rates. No correlation was found between the Ca(2+) binding properties of the TnC(F29W) mutants and the solvent accessibility of the hydrophobic amino acids in the apo state, Ca(2+) bound state, or the difference between the two states. However, the effects of these hydrophobic mutations on Ca(2+) binding were contextual possibly because of side chain interactions within the apo and Ca(2+) bound states of the N-domain. These results demonstrate that a single hydrophobic residue, which does not directly ligate Ca(2+), can play a crucial role in controlling Ca(2+) binding and exchange within a coupled and functional EF-hand system.     

HD exchange and PLIMSTEX determine the affinities and order of binding of Ca2+ with troponin C

Troponin C (TnC), present in all striated muscle, is the Ca(2+)-activated trigger that initiates myocyte contraction. The binding of Ca(2+) to TnC initiates a cascade of conformational changes involving the constituent proteins of the thin filament. The functional properties of TnC and its ability to bind Ca(2+) have significant regulatory influence on the contractile reaction of muscle. Changes in TnC may also correlate with cardiac and various other muscle-related diseases. We report here the implementation of the PLIMSTEX strategy (protein ligand interaction by mass spectrometry, titration, and H/D exchange) to elucidate the binding affinity of TnC with Ca(2+) and, more importantly, to determine the order of Ca(2+) binding of the four EF hands of the protein. The four equilibrium constants, K(1) = (5 ± 5) × 10(7) M(-1), K(2) = (1.8 ± 0.8) × 10(7) M(-1), K(3) = (4.2 ± 0.9) × 10(6) M(-1), and K(4) = (1.6 ± 0.6) × 10(6) M(-1), agree well with determinations by other methods and serve to increase our confidence in the PLIMSTEX approach. We determined the order of binding to the four EF hands to be III, IV, II, and I by extracting from the H/DX results the deuterium patterns for each EF hand for each state of the protein (apo through fully Ca(2+) bound). This approach, demonstrated for the first time, may be general for determining binding orders of metal ions and other ligands to proteins.     

Interaction of troponin I and troponin C: 19F NMR studies of the binding of the inhibitory troponin I peptide to turkey skeletal troponin C.

We have used 19F nuclear magnetic resonance spectroscopy to study the interaction of the inhibitory region of troponin (TnI) with apo- and calcium(II)-saturated turkey skeletal troponin C (TnC), using the synthetic TnI analogue N alpha-acetyl[19FPhe106]TnI(104-115)amide. Dissociation constants of Kd = (3.7 +/- 3.1) x 10(-5) M for the apo interaction and Kd = (4.8 +/- 1.8) x 10(-5) M for the calcium(II)-saturated interaction were obtained using a 1:1 binding model of peptide to protein. The 19F NMR chemical shifts for the F-phenylalanine of the bound peptide are different from the apo- and calcium-saturated protein, indicating a different environment for the bound peptide. The possibility of 2:1 binding of the peptide to Ca(II)-saturated TnC was tested by calculating the fit of the experimental titration data to a series of theoretical binding curves in which the dissociation constants for the two hypothetical binding sites were varied. We obtained the best fit for 0.056 mM less than or equal to Kd1 less than or equal to 0.071 mM and 0.5 mM less than or equal to Kd2 less than or equal to 2.0 mM. These results allow the possibility of a second peptide binding site on calcium(II)-saturated TnC with an affinity 10- to 20-fold weaker than that of the first site.     

Biliverdin Reductase B Dynamics Are Coupled to Coenzyme Binding

Biliverdin reductase B (BLVRB) is a newly identified cellular redox regulator that catalyzes the NADPH-dependent reduction of multiple substrates. Through mass spectrometry analysis, we identified high levels of BLVRB in mature red blood cells, highlighting the importance of BLVRB in redox regulation. The BLVRB conformational changes that occur during conezyme/substrate binding and the role of dynamics in BLVRB function, however, remain unknown. Through a combination of NMR, kinetics, and isothermal titration calorimetry studies, we determined that BLVRB binds its coenzyme 500-fold more tightly than its substrate. While the active site of apo BLVRB is highly dynamic on multiple timescales, active site dynamics are largely quenched within holo BLVRB, in which dynamics are redistributed to other regions of the enzyme. We show that a single point mutation of Arg78?Ala leads to both an increase in active site micro-millisecond motions and an increase in the microscopic rate constants of coenzyme binding. This demonstrates that altering BLVRB active site dynamics can directly cause a change in functional characteristics. Our studies thus address the solution behavior of apo and holo BLVRB and identify a role of enzyme dynamics in coenzyme binding.     

Fourier transform infrared spectroscopic study on the binding of Mg2+ to a mutant Akazara scallop troponin C (E142Q)

Troponin C (TnC) is the Ca(2+)-binding regulatory protein of the troponin complex in muscle tissue. Vertebrate fast skeletal muscle TnCs bind four Ca(2+), while Akazara scallop (Chlamys nipponensis akazara) striated adductor muscle TnC binds only one Ca(2+) at site IV, because all the other EF-hand motifs are short of critical residues for the coordination of Ca(2+). Fourier transform infrared (FTIR) spectroscopy was applied to study coordination structure of Mg(2+) bound in a mutant Akazara scallop TnC (E142Q) in D(2)O solution. The result showed that the side-chain COO(-) groups of Asp 131 and Asp 133 in the Ca(2+)-binding site of E142Q bind to Mg(2+) in the pseudo-bridging mode. Mg(2+) titration experiments for E142Q and the wild-type of Akazara scallop TnC were performed by monitoring the band at about 1600 cm(-1), which is due to the pseudo-bridging Asp COO(-) groups. As a result, the binding constants of them for Mg(2+) were the same value (about 6 mM). Therefore, it was concluded that the side-chain COO(-) group of Glu 142 of the wild type has no relation to the Mg(2+) ligation. The effect of Mg(2+) binding in E142Q was also investigated by CD and fluorescence spectroscopy. The on-off mechanism of the activation of Akazara scallop TnC is discussed on the basis of the coordination structures of Mg(2+) as well as Ca(2+).     

Symmetrical stabilization of bound Ca2+ ions in a cooperative pair of EF-hands through hydrogen bonding of coordinating water molecules in calbindin D(9k)

Water molecules are found to complete the Ca2+ coordination sphere when a protein fails to provide enough ligating oxygens. Hydrogen bonding of these water molecules to the protein backbone or side chains may contribute favorably to the Ca2+ affinity, as suggested in an earlier study of two calbindin D(9k) mutants [E60D and E60Q; Linse et al. (1994) Biochemistry 33, 12478-12486]. To investigate the generality of this conclusion, another side chain, Gln 22, which hydrogen bonds to a Ca2+-coordinating water molecule in calbindin D(9k), was mutated. Two calbindin D(9k) mutants, (Q22E+P43M) and (Q22N+P43M), were constructed to examine the interaction between Gln 22 and the water molecule in the C-terminal calcium binding site II. Shortening of the side chain, as in (Q22N+P43M), reduces the affinity of binding two calcium ions by a factor of 18 at low ionic strength, whereas introduction of a negative charge, as in (Q22E+P43M), leads to a 12-fold reduction. In 0.15 M KCl, a 7-fold reduction in affinity was observed for both mutants. The cooperativity of Ca2+ binding increases for (Q22E+P43M), while it decreases for (Q22N+P43M). The rates of Ca2+ dissociation are 5.5-fold higher for the double mutants than for P43M at low ionic strength. For both mutants, reduced strength of hydrogen bonding to calcium-coordinating water molecules is a likely explanation for the observed effects on Ca2+ affinity and dissociation. In the apo forms, the (Q22E+P43M) mutant has lower stability toward urea denaturation than (Q22N+P43M) and P43M. 2D (1)H NMR and crystallographic experiments suggest that the structure of (Q22E+P43M) and (Q22N+P43M) is unchanged relative to P43M, except for local perturbations in the loop regions.     

Structural characterization of the apo form of a calcium binding protein from Entamoeba histolytica by hydrogen exchange and its folding to the holo state

One of the calcium binding proteins from Entamoeba histolytica (EhCaBP) is a 134 amino acid residue long (M(r) approximately 14.9 kDa) double domain EF-hand protein containing four Ca(2+) binding sites. CD and NMR studies reveal that the Ca(2+)-free form (apo-EhCaBP) exists in a partially collapsed form compared to the Ca(2+)-bound (holo) form, which has an ordered structure (PDB ID ). Deuterium exchange studies on the partially structured apo-EhCaBP reveal that the C-terminal domain is better structured than the N-terminal domain. The protein can be reversibly folded and unfolded upon addition of Ca(2+) and EGTA, respectively. Titration shows a slow initial folding of the apo form with increasing Ca(2+) concentration, followed by a highly cooperative folding to its final state at a certain threshold of Ca(2+). Ca(2+) and the EGTA titration taken together show that site II in the N-terminal domain has the highest affinity for Ca(2+) contrary to earlier studies. Further, this study has thrown light on the relative Ca(2+) binding affinity and specificity of each site in the intact protein. A structural model for the partially collapsed form of apo-EhCaBP and its equilibrium folding to its completely folded holo state has been suggested. Large conformational changes seen in transforming from the apo to holo form of EhCaBP suggest that this protein should be functioning as a sensor protein and might have a significant role in host-parasite recognition.     

Tyrosine and tyrosinate fluorescence of S-100b. A time-resolved nanosecond fluorescence study. The effect of pH, Ca(II), and Zn(II)

The properties of the tyrosine and tyrosinate emissions from brain S-100b have been studied by nanosecond time-resolved fluorescence at emission wavelengths in the range 305 to 365 nm. The effect of pH on the fluorescence has been studied at pH 6.5, 7.5, and 8.5 for the Ca(II) apo and holo forms of the protein, and for the apo and holo forms in the presence and absence of Zn(II) at pH 7.5. The fluorescence decay is biexponential at pH 8.5 and triexponential at pH 6.5 and 7.5. The three components of the decay have wavelength and metal ion dependent lifetimes in the ranges 0.06 to 1.05 ns, 0.49 to 3.76 ns, and 3.60 to 14.5 ns. The observation of a long lifetime component at wavelengths characteristic of emission from tyrosinate suggests that in class A proteins this may be a useful diagnostic of the environment of tyrosine in their native structures. The time-resolved emission spectra provide evidence for efficient, subnanosecond protolysis of the excited state of the single tyrosine (Tyr17) under all conditions studied except in 6 M guanidium chloride in which the protein shows only emission from tyrosine (lambda em 305 nm), suggesting that the tyrosinate emission is a property of the tertiary structure of the native protein. The Zn(II)-dependence of the fluorescence is fully consistent with the earlier suggestion that Tyr17 is near the Zn(II) binding site and remote from the high affinity Ca(II) binding site.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Structural and functional role of nickel ions in urease by molecular dynamics simulation

Nickel ions play several roles in the biological processes of microorganisms and plants. Urease has a nickel-containing active site and catalyzes the hydrolysis of urea to yield ammonia and carbamate. In the present study, the role of nickel ions is examined using molecular dynamics simulations of the holo and apo forms. Nonbonded models used for the nickel ions provide good reproduction of the active-site structure as indicated in the crystallized structure. The results confirm that urease has a rigid active site in either its holo or its apo form. A new conformation of the flap is observed in apo urease. The connection between the metal center and His^α323 is proposed to be responsible for maintaining the flap conformation. The binding free energy of acetohydroxamic acid and urease is estimated using the molecular mechanics–generalized Born/surface area method. The binding free energy is primarily driven by electrostatic interactions in the presence of nickel ions. Normal mode analysis is employed to characterize the movements of the flap in apo urease.     

Evidence that both Ca(2+)-specific sites of skeletal muscle TnC are required for full activity

To investigate the role of the Ca2(+)-specific (I and II) sites of fast skeletal muscle troponin C (TnC) in the regulation of contraction, we have produced two TnC mutants which have lost the ability to bind Ca2+ at either site I (VG1) or at site II (VG2). Both mutants were able to partially restore force to TnC-depleted skinned muscle fibers (approximately 25% for VG1 and approximately 50% for VG2). In contrast, bovine cardiac TnC (BCTnC), which like VG1 binds Ca2+ only at site II, could fully reactivate the contraction of TnC-depleted fibers. Higher concentrations of both mutants were required to restore force to the TnC-depleted fibers than with wild type TnC (WTnC) or BCTnC. VG1 and VG2 substituted fibers could not bind additional WTnC, indicating that all of the TnC-binding sites were saturated with the mutant TnC's. The Ca2+ concentration required for force activation was much higher for VG1 and VG2 substituted fibers than for WTnC or BCTnC substituted fibers. Also, the steepness of force activation was much less in VG1 and VG2 versus WTnC and BCTnC substituted fibers. These results suggest cooperative interactions between sites I and II in WTnC. In contrast, BCTnC has essentially the same apparent Ca2+ affinity and steepness of force activation as does WTnC. Thus, cardiac TnC must have structural differences from WTnC which compensate for the lack of site I, while in WTnC, both Ca2(+)-specific sites are probably crucial for full functional activity.     

Kinetic stability of Cu/Zn superoxide dismutase is dependent on its metal ligands: implications for ALS

Over 100 mutants of the enzyme Cu/Zn superoxide dismutase (SOD) have been implicated in the neurodegenerative disease familial amyotrophic lateral sclerosis (FALS). Growing evidence suggests that the aggregation of SOD mutants may play a causative role in FALS and that aberrant copper chemistry, decreased thermodynamic stability, and decreased affinity for metals may contribute independently or synergistically to this process. Since the loss of the copper and zinc ions significantly decreases the thermodynamic stability of SOD, it is expected that this would also decrease its kinetic stability, thereby facilitating partial or global unfolding transitions that may lead to misfolding and aggregation. Here we used wild-type (WT) SOD and five FALS-related mutants (G37R, H46R, G85R, D90A, and L144F) to show that the metals contribute significantly to the kinetic stability of the protein, with demetalated (apo) SOD showing acid-induced unfolding rates about 60-fold greater than the metalated (holo) protein. However, the unfolding rates of SOD WT and mutants were similar to each other in both the holo and apo states, indicating that regardless of the effect of mutation on thermodynamic stability, the kinetic barrier toward SOD unfolding is dependent on the presence of metals. Thus, these results suggest that pathogenic SOD mutations that do not significantly alter the stability of the protein may still lead to SOD aggregation by compromising its ability to bind or retain its metals and thereby decrease its kinetic stability. Furthermore, the mutant-like decrease in the kinetic stability of apo WT SOD raises the possibility that the loss of metals in WT SOD may be involved in nonfamilial forms of ALS.     

Prediction of interaction sites from apo 3D structures when the holo conformation is different

The predictability of catalytic and binding sites from apo structures is addressed for proteins that undergo significant conformational change upon binding. Theoretical microscopic titration curves (THEMATICS), an electrostatics-based method for the prediction of functional sites, is performed on a test set of 24 proteins with both apo and holo structures available. For 23 of these 24 proteins (96%), THEMATICS predicts the correct catalytic or binding site for both the apo and holo forms. For only one of the 24 proteins, THEMATICS makes the correct prediction for the holo structure but fails for the apo structure. The metrics used by THEMATICS to identify functional residues generally are larger in absolute value for the functional residues in the holo forms compared to the corresponding residues in the apo forms. However, even in the apo forms, these identifying metrics are still statistically significantly larger for functional residues than for residues not involved in catalysis or binding. This indicates that some of the unusual electrostatic properties of functional residues are preserved in the apo conformation. Evidence is presented that certain residues immediately surrounding the active catalytic and binding residues impart functionally important chemical and electrostatic properties to the active residues. At least parts of these microenvironments exist in the unbound conformations, such that THEMATICS is able to distinguish the functional residues even in the apo structures.     

Role of residue 147 in the gene regulatory function of the Escherichia coli purine repressor.

The crystal structures of corepressor-bound and free Escherichia coli purine repressor (PurR) have delineated the roles of several residues in corepressor binding and specificity and the intramolecular signal transduction (allosterism) of this LacI/GalR family member. From these structures, residue W147 was implicated as a key component of the allosteric response, but in many members of the LacI/GalR family, position 147 is occupied by an arginine. To understand the role of this tryptophan at position 147, three proteins, substituted by phenylalanine (W147F), alanine (W147A), or arginine (W147R), were constructed and characterized in vivo and in vitro, and their structures were determined. W147F displays a decreased affinity for corepressor and is a poor repressor in vivo. W147A and W147R, on the other hand, are super repressors and bind corepressor 13.6 and 7.9 times more tightly, respectively, than wild-type. Each mutant PurR-hypoxanthine-purF operator holo complex crystallizes isomorphously to wild-type. Whereas the apo corepressor binding domain (CBD) of W147F crystallizes under those conditions used for the wild-type protein, neither the apo CBD of W147R nor W147A crystallizes, although screened extensively for new crystal forms. Structures of the holo repressor mutants have been solved to resolutions between 2.5 and 2.9 A, and the structure of the apo CBD of W147F has been solved to 2.4 A resolution. These structures provide insight into the altered biochemical properties and physiological functions of these mutants, which appear to depend on the sometimes subtle preference for one conformation (apo vs holo) over the other.     

Hydrophobic interactions of transcobalamin II (TC II) from mammalian sera

The hydrophobic properties of mammalian transcobalamin IIs (TC II) were studied by chromatography of radioactive cyanocobalamin (CN[57Co]Cbl)-labeled serum on phenyl-Sepharose CL-4B. Mammalian holo TC IIs (CN[57Co]Cbl-TC II) exhibited species variability in their affinity for the hydrophobic matrix in the order: dog greater than mouse greater than human greater than rat greater than rabbit. Phenyl-Sepharose chromatography of the isolated CN[57Co]Cbl-TC II peaks from gel filtration of dog and rat serum showed no hydrophobic change in dog TC II, but an increase in hydrophobicity of rat TC II. Phenyl-Sepharose chromatography of CN[57Co]Cbl-labeled rabbit serum (holo TC II) and the unlabeled serum (apo TC II) showed apo TC II to be more hydrophobic than holo TC II as has been shown for human TC II (Begley et al., Biochem Biophys Res Commun 103:434-441, 1981). Thus mammalian holo TC IIs differ in their hydrophobic properties and apo TC II, in man and rabbit, is more hydrophobic than holo TC II. In addition, isolation of the TC II in some animal sera by gel filtration may result in a TC II that is more hydrophobic than the native molecule.     

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca2? dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 ?, 5.5 ?, 6 ?, 6.5 ?, 7.0 ? and 7.5 ? were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.