Heterogeneity of Equilibrium Molten Globule State of Cytochrome c Induced by Weak Salt Denaturants under Physiological Condition

While many proteins are recognized to undergo folding via intermediate(s), the heterogeneity of equilibrium folding intermediate(s) along the folding pathway is less understood. In our present study, FTIR spectroscopy, far- and near-UV circular dichroism (CD), ANS and tryptophan fluorescence, near IR absorbance spectroscopy and dynamic light scattering (DLS) were used to study the structural and thermodynamic characteristics of the native (N), denatured (D) and intermediate state (X) of goat cytochorme c (cyt-c) induced by weak salt denaturants (LiBr, LiCl and LiClO₄) at pH 6.0 and 25°C. The LiBr-induced denaturation of cyt-c measured by Soret absorption (Δε ₄₀₀) and CD ([θ]₄₀₉), is a three-step process, N ↔ X ↔ D. It is observed that the X state obtained along the denaturation pathway of cyt-c possesses common structural and thermodynamic characteristics of the molten globule (MG) state. The MG state of cyt-c induced by LiBr is compared for its structural and thermodynamic parameters with those found in other solvent conditions such as LiCl, LiClO₄ and acidic pH. Our observations suggest: (1) that the LiBr-induced MG state of cyt-c retains the native Met80-Fe(III) axial bond and Trp59-propionate interactions; (2) that LiBr-induced MG state of cyt-c is more compact retaining the hydrophobic interactions in comparison to the MG states induced by LiCl, LiClO₄ and 0.5 M NaCl at pH 2.0; and (3) that there exists heterogeneity of equilibrium intermediates along the unfolding pathway of cyt-c as highly ordered (X1), classical (X2) and disordered (X3), i.e., D ↔ X3 ↔ X2 ↔ X1 ↔ N.     

Recombinant, Replication-Defective Adenovirus Gene Transfer Vectors Induce Cell Cycle Dysregulation and Inappropriate Expression of Cyclin Proteins

First-generation adenovirus (Ad) vectors that had been rendered replication defective by removal of the E1 region of the viral genome (ΔE1) or lacking the Ad E3 region in addition to E1 sequences (ΔE1ΔE3) induced G₂ cell cycle arrest and inhibited traverse across G₁/S in primary and immortalized human bronchial epithelial cells. Cell cycle arrest was independent of the cDNA contained in the expression cassette and was associated with the inappropriate expression and increase in cyclin A, cyclin B1, cyclin D, and cyclin-dependent kinase p34^cdc2 protein levels. In some instances, infection with ΔE1 or ΔE1ΔE3 Ad vectors produced aneuploid DNA histogram patterns and induced polyploidization as a result of successive rounds of cell division without mitosis. Cell cycle arrest was absent in cells infected with a second-generation ΔE1Ad vector in which all of the early region E4 except the sixth open reading frame was also deleted. Consequently, E4 viral gene products present in ΔE1 or ΔE1ΔE3 Ad vectors induce G₂ growth arrest, which may pose new and unintended consequences for human gene transfer and gene therapy.     

Study of the Affinity between the Protein Kinase PKA and Peptide Substrates Derived from Kemptide Using Molecular Dynamics Simulations and MM/GBSA

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG_E→TS) determined by the catalytic efficiency (k_cat/K_M) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.     

Structural Characteristic of the Initial Unfolded State on Refolding Determines Catalytic Efficiency of the Folded Protein in Presence of Osmolytes

Osmolytes are low molecular weight organic molecules accumulated by organisms to assist proper protein folding, and to provide protection to the structural integrity of proteins under denaturing stress conditions. It is known that osmolyte-induced protein folding is brought by unfavorable interaction of osmolytes with the denatured/unfolded states. The interaction of osmolyte with the native state does not significantly contribute to the osmolyte-induced protein folding. We have therefore investigated if different denatured states of a protein (generated by different denaturing agents) interact differently with the osmolytes to induce protein folding. We observed that osmolyte-assisted refolding of protein obtained from heat-induced denatured state produces native molecules with higher enzyme activity than those initiated from GdmCl- or urea-induced denatured state indicating that the structural property of the initial denatured state during refolding by osmolytes determines the catalytic efficiency of the folded protein molecule. These conclusions have been reached from the systematic measurements of enzymatic kinetic parameters (K _m and k _cat), thermodynamic stability (T _m and ΔH _m) and secondary and tertiary structures of the folded native proteins obtained from refolding of various denatured states (due to heat-, urea- and GdmCl-induced denaturation) of RNase-A in the presence of various osmolytes.     

Human Coding Synonymous Single Nucleotide Polymorphisms at Ramp Regions of mRNA Translation

According to the ramp model of mRNA translation, the first 50 codons favor rare codons and have slower speed of translation. This study aims to detect translational selection on coding synonymous single nucleotide polymorphisms (sSNP) to support the ramp theory. We investigated fourfold degenerate site (FFDS) sSNPs with A↔G or C↔T substitutions in human genome for distribution bias of synonymous codons (SC), grouped by CpG or non-CpG sites. Distribution bias of sSNPs between the 3^rd ∼50^th codons and the 51^st ∼ remainder codons at non-CpG sites were observed. In the 3^rd ∼50^th codons, G→A sSNPs at non-CpG sites are favored than A→G sSNPs [P = 2.89×10⁻³], and C→T at non-CpG sites are favored than T→C sSNPs [P = 8.50×10⁻³]. The favored direction of SC usage change is from more frequent SCs to less frequent SCs. The distribution bias is more obvious in synonymous substitutions CG(G→A), AC(C→T), and CT(C→T). The distribution bias of sSNPs in human genome, i.e. frequent SCs to less frequent SCs is favored in the 3^rd ∼50^th codons, indicates translational selection on sSNPs in the ramp regions of mRNA templates.     

Energetic signatures of single base bulges: thermodynamic consequences and biological implications

DNA bulges are biologically consequential defects that can arise from template-primer misalignments during replication and pose challenges to the cellular DNA repair machinery. Calorimetric and spectroscopic characterizations of defect-containing duplexes reveal systematic patterns of sequence-context dependent bulge-induced destabilizations. These distinguishing energetic signatures are manifest in three coupled characteristics, namely: the magnitude of the bulge-induced duplex destabilization (ΔΔG_Bulge); the thermodynamic origins of ΔΔG_Bulge (i.e. enthalpic versus entropic); and, the cooperativity of the duplex melting transition (i.e. two-state versus non-two state). We find moderately destabilized duplexes undergo two-state dissociation and exhibit ΔΔG_Bulge values consistent with localized, nearest neighbor perturbations arising from unfavorable entropic contributions. Conversely, strongly destabilized duplexes melt in a non-two-state manner and exhibit ΔΔG_Bulge values consistent with perturbations exceeding nearest-neighbor expectations that are enthalpic in origin. Significantly, our data reveal an intriguing correlation in which the energetic impact of a single bulge base centered in one strand portends the impact of the corresponding complementary bulge base embedded in the opposite strand. We discuss potential correlations between these bulge-specific differential energetic profiles and their overall biological implications in terms of DNA recognition, repair and replication.     

Thermodynamic coupling between neighboring binding sites in homo‐oligomeric ligand sensing proteins from mass resolved ligand‐dependent population distributions

Homo‐oligomeric ligand‐activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand‐binding sites. Understanding the basis for this regulation requires both quantifying the free energy ΔG transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring‐shaped homo‐oligomeric trp RNA‐binding attenuation protein (TRAP). First, we developed a nearest‐neighbor statistical thermodynamic binding model comprising microscopic free energies for ligand binding to isolated sites ΔG ₀, and for coupling between adjacent sites, ΔG _α . Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2 ^N possible liganded states. We then experimentally monitored ligand‐dependent population shifts using conventional spectroscopic and calorimetric methods and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand‐dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2 ^N liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure–thermodynamic linkage that drives cooperativity.     

Bacillus subtilis polynucleotide phosphorylase 3′-to-5′ DNase activity is involved in DNA repair

In the presence of Mn²⁺, an activity in a preparation of purified Bacillus subtilis RecN degrades single-stranded (ss) DNA with a 3′ → 5′ polarity. This activity is not associated with RecN itself, because RecN purified from cells lacking polynucleotide phosphorylase (PNPase) does not show the exonuclease activity. We show here that, in the presence of Mn²⁺ and low-level inorganic phosphate (P_i), PNPase degrades ssDNA. The limited end-processing of DNA is regulated by ATP and is inactive in the presence of Mg²⁺ or high-level P_i. In contrast, the RNase activity of PNPase requires Mg²⁺ and P_i, suggesting that PNPase degradation of RNA and ssDNA occur by mutually exclusive mechanisms. A null pnpA mutation (ΔpnpA) is not epistatic with ΔrecA, but is epistatic with ΔrecN and Δku, which by themselves are non-epistatic. The addA5, ΔrecO, ΔrecQ (ΔrecJ), ΔrecU and ΔrecG mutations (representative of different epistatic groups), in the context of ΔpnpA, demonstrate gain- or loss-of-function by inactivation of repair-by-recombination, depending on acute or chronic exposure to the damaging agent and the nature of the DNA lesion. Our data suggest that PNPase is involved in various nucleic acid metabolic pathways, and its limited ssDNA exonuclease activity plays an important role in RecA-dependent and RecA-independent repair pathways.     

Predicting Thermodynamic Properties of PBXTHs with New Quantum Topological Indexes

Novel group quantitative structure-property relationship (QSPR) models on the thermodynamic properties of PBXTHs were presented, by the multiple linear regression (MLR) analysis method. Four thermodynamic properties were studied: the entropy (S^θ), the standard enthalpy of formation (Δ_fH^θ), the standard Gibbs energy of formation (Δ_fG^θ), and the relative standard Gibbs energy of formation (Δ_RG^θ). The results by the formula indicate that the calculated and predicted data in this study are in good agreement with those in literature and the deviation is within the experimental errors. To validate the estimation reliability for internal samples and the predictive ability for other samples, leave-one-out (LOO) cross validation (CV) and external validation were performed, and the results show that the models are satisfactory.     

Thermodynamically stable aggregation-resistant antibody domains through directed evolution

Protein aggregates are usually formed by interactions between unfolded or partially unfolded species, and often occur when a protein is denatured by, for example, heat or low pH. In earlier work, we used a Darwinian selection strategy to create human antibody variable domains that resisted heat aggregation. The repertoires of domains were displayed on filamentous phage and denatured (at 80 °C in pH 7.4), and folded domains were selected by binding to a generic ligand after cooling. This process appeared to select for domains with denatured states that resisted aggregation, but the domains only had low free energies of folding (ΔG_N-D^o = 15-20 kJ/mol at 25 °C in pH 7.4). Here, using the same phage repertoire, we have extended the method to the selection of domains resistant to acid aggregation. In this case, however, the thermodynamic stabilities of selected domains were higher than those selected by thermal denaturation (under both neutral and acidic conditions; ΔG_N-D^o = 26-47 kJ/mol at 25 °C in pH 7.4, or ΔG_N-D^o = 27-34 kJ/mol in pH 3.2). Furthermore, we identified a key determinant (Arg28) that increased the aggregation resistance of the denatured states of the domains at low pH without compromising their thermodynamic stabilities. Thus, the selection process yielded domains that combined thermodynamic stability and aggregation-resistant unfolded states. We suggest that changes to these properties are controlled by the extent to which the folding equilibrium is displaced during the process of selection.     

Neutralization of the Charge on Asp369 of Na+,K+-ATPase Triggers E1 ? E2 Conformational Changes

This work investigates the role of charge of the phosphorylated aspartate, Asp³⁶⁹, of Na⁺,K⁺-ATPase on E₁ ↔ E₂ conformational changes. Wild type (porcine α₁/His₁₀-β₁), D369N/D369A/D369E, and T212A mutants were expressed in Pichia pastoris, labeled with fluorescein 5′-isothiocyanate (FITC), and purified. Conformational changes of wild type and mutant proteins were analyzed using fluorescein fluorescence (Karlish, S. J. (1980) J. Bioenerg. Biomembr. 12, 111–136). One central finding is that the D369N/D369A mutants are strongly stabilized in E₂ compared with wild type and D369E or T212A mutants. Stabilization of E₂(Rb) is detected by a reduced K_0.5Rb for the Rb⁺-induced E₁ ↔ E₂(2Rb) transition. The mechanism involves a greatly reduced rate of E₂(2Rb) → E₁Na with no effect on E₁ → E₂(2Rb). Lowering the pH from 7.5 to 5.5 strongly stabilizes wild type in E₂ but affects the D369N mutant only weakly. Thus, this “Bohr” effect of pH on E₁ ↔ E₂ is due largely to protonation of Asp³⁶⁹. Two novel effects of phosphate and vanadate were observed with the D369N/D369A mutants as follows. (a) E₁ → E₂·P is induced by phosphate without Mg²⁺ ions by contrast with wild type, which requires Mg²⁺. (b) Both phosphate and vanadate induce rapid E₁ → E₂ transitions compared with slow rates for the wild type. With reference to crystal structures of Ca²⁺-ATPase and Na⁺,K⁺-ATPase, negatively charged Asp³⁶⁹ favors disengagement of the A domain from N and P domains (E₁), whereas the neutral D369N/D369A mutants favor association of the A domain (TGES sequence) with P and N domains (E₂). Changes in charge interactions of Asp³⁶⁹ may play an important role in triggering E₁P(3Na) ↔ E₂P and E₂(2K) → E₁Na transitions in native Na⁺,K⁺-ATPase.     

Determination of base and backbone contributions to the thermodynamics of premelting and melting transitions in B DNA

In previous papers of this series the temperature-dependent Raman spectra of poly(dA)·poly(dT) and poly(dA–dT)·poly(dA–dT) were used to characterize structurally the melting and premelting transitions in DNAs containing consecutive A·T and alternating A·T/T·A base pairs. Here, we describe procedures for obtaining thermodynamic parameters from the Raman data. The method exploits base-specific and backbone-specific Raman markers to determine separate thermodynamic contributions of A, T and deoxyribosyl-phosphate moieties to premelting and melting transitions. Key findings include the following: (i) Both poly(dA)·poly(dT) and poly(dA–dT)· poly(dA–dT) exhibit robust premelting transitions, due predominantly to backbone conformational changes. (ii) The significant van’t Hoff premelting enthalpies of poly(dA)·poly(dT) [ΔH_vH^pm = 18.0 ± 1.6 kcal·mol^–1 (kilocalories per mole cooperative unit)] and poly(dA–dT)·poly(dA–dT) (ΔH_vH^pm = 13.4 ± 2.5 kcal·mol^–1) differ by an amount (~4.6 kcal·mol^–1) estimated as the contribution from three-centered inter-base hydrogen bonding in (dA)_n·(dT)_n tracts. (iii) The overall stacking free energy of poly(dA)· poly(dT) [–6.88 kcal·mol_bp^–1 (kilocalories per mole base pair)] is greater than that of poly(dA–dT)· poly(dA–dT) (–6.31 kcal·mol_bp^–1). (iv) The difference between stacking free energies of A and T is significant in poly(dA)·poly(dT) (ΔΔG_st = 0.8 ± 0.3 kcal· mol_bp^–1), but marginal in poly(dA–dT)·poly(dA–dT) (ΔΔG_st = 0.3 ± 0.3 kcal·mol_bp^–1). (v) In poly(dA)· poly(dT), the van’t Hoff parameters for melting of A (ΔH_vH^A = 407 ± 23 kcal·mol^–1, ΔS_vH^A = 1166 ± 67 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 60.0 ± 3.2 kcal·mol^–1) are clearly distinguished from those of T (ΔH_vH^T = 185 ± 38 kcal·mol^–1, ΔS_vH^T = 516 ± 109 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 27.1 ± 5.5 kcal·mol^–1). (vi) Similar relative differences are observed in poly(dA–dT)· poly(dA–dT) (ΔH_vH^A = 333 ± 54 kcal·mol^–1, ΔS_vH^A = 961 ± 157 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 45.0 ± 7.6 kcal· mol^–1; ΔH_vH^T = 213 ± 30 kcal·mol^–1, ΔS_vH^T = 617 ± 86 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 29.3 ± 4.9 kcal·mol^–1). The methodology employed here distinguishes thermodynamic contributions of base stacking, base pairing and backbone conformational ordering in the molecular mechanism of double-helical B DNA formation.     

Energy for Wild-Type Acetylcholine Receptor Channel Gating from Different Choline Derivatives

Agonists, including the neurotransmitter acetylcholine (ACh), bind at two sites in the neuromuscular ACh receptor channel (AChR) to promote a reversible, global change in protein conformation that regulates the flow of ions across the muscle cell membrane. In the synaptic cleft, ACh is hydrolyzed to acetate and choline. Replacement of the transmitter’s ester acetyl group with a hydroxyl (ACh→choline) results in a +1.8 kcal/mol reduction in the energy for gating generated by each agonist molecule from a low- to high-affinity change of the transmitter binding site (ΔG_B). To understand the distinct actions of structurally related agonist molecules, we measured ΔG_B for 10 related choline derivatives. Replacing the hydroxyl group of choline with different substituents, such as hydrogen, chloride, methyl, or amine, increased the energy for gating (i.e., it made ΔG_B more negative relative to choline). Extending the ethyl hydroxide tail of choline to propyl and butyl hydroxide also increased this energy. Our findings reveal the amount of energy that is available for the AChR conformational change provided by different, structurally related agonists. We speculate that a hydrogen bond between the choline hydroxyl and the backbone carbonyl of αW149 positions this agonist’s quaternary ammonium group so as to reduce the cation-π interaction between this moiety and the aromatic groups at the binding site.     

Metal migration and subunit swapping in ALS-linked SOD1: Zn2+ transfer between mutant and wild-type occurs faster than the rate of heterodimerization

The heterodimerization of WT Cu, Zn superoxide dismutase-1 (SOD1), and mutant SOD1 might be a critical step in the pathogenesis of SOD1-linked amyotrophic lateral sclerosis (ALS). Rates and free energies of heterodimerization (ΔG_Het) between WT and ALS-mutant SOD1 in mismatched metalation states—where one subunit is metalated and the other is not—have been difficult to obtain. Consequently, the hypothesis that under-metalated SOD1 might trigger misfolding of metalated SOD1 by “stealing” metal ions remains untested. This study used capillary zone electrophoresis and mass spectrometry to track heterodimerization and metal transfer between WT SOD1, ALS-variant SOD1 (E100K, E100G, D90A), and triply deamidated SOD1 (modeled with N26D/N131D/N139D substitutions). We determined that rates of subunit exchange between apo dimers and metalated dimers—expressed as time to reach 30% heterodimer—ranged from t_30% = 67.75 ± 9.08 to 338.53 ± 26.95 min; free energies of heterodimerization ranged from ΔG_Het = -1.21 ± 0.31 to -3.06 ± 0.12 kJ/mol. Rates and ΔG_Het values of partially metalated heterodimers were more similar to those of fully metalated heterodimers than apo heterodimers, and largely independent of which subunit (mutant or WT) was metal-replete or metal-free. Mass spectrometry and capillary electrophoresis demonstrated that mutant or WT 4Zn-SOD1 could transfer up to two equivalents of Zn²⁺ to mutant or WT apo-SOD1 (at rates faster than the rate of heterodimerization). This result suggests that zinc-replete SOD1 can function as a chaperone to deliver Zn²⁺ to apo-SOD1, and that WT apo-SOD1 might increase the toxicity of mutant SOD1 by stealing its Zn²⁺.     

Thermodynamic criteria for high hit rate antisense oligonucleotide design

Antisense oligonucleotides are used for therapeutic applications and in functional genomic studies. In practice, however, many of the oligonucleotides complementary to an mRNA have little or no antisense activity. Theoretical strategies to improve the ‘hit rate’ in antisense screens will reduce the cost of discovery and may lead to identification of antisense oligonucleotides with increased potency. Statistical analysis performed on data collected from more than 1000 experiments with phosphorothioate-modified oligonucleotides revealed that the oligo-probes, which form stable duplexes with RNA (ΔG^o₃₇ ≤ –30 kcal/mol) and have small self-interaction potential, are more frequently efficient than molecules that form less stable oligonucleotide–RNA hybrids or more stable self-structures. To achieve optimal statistical preference, the values for self-interaction should be (ΔG^o₃₇) ≥ –8 kcal/mol for inter-oligonucleotide pairing and (ΔG^o₃₇) ≥ –1.1 kcal/mol for intra-molecular pairing. Selection of oligonucleotides with these thermodynamic values in the analyzed experiments would have increased the ‘hit rate’ by as much as 6-fold.     

Thermodynamics of Formate-Oxidizing Metabolism and Implications for H2 Production

Formate-dependent proton reduction to H₂ (HCOO⁻ + H₂O → HCO₃⁻ + H₂) has been reported for hyperthermophilic Thermococcus strains. In this study, a hyperthermophilic archaeon, Thermococcus onnurineus strain NA1, yielded H₂ accumulation to a partial pressure of 1 × 10⁵ to 7 × 10⁵ Pa until the values of Gibbs free energy change (ΔG) reached near thermodynamic equilibrium (−1 to −3 kJ mol⁻¹). The bioenergetic requirement for the metabolism to conserve energy was demonstrated by ΔG values as small as −5 kJ mol⁻¹, which are less than the biological minimum energy quantum, −20 kJ mol⁻¹, as calculated by Schink (B. Schink, Microbiol. Mol. Biol. Rev. 61:262-280, 1997). Considering formate as a possible H₂ storage material, the H₂ production potential of the strain was assessed. The volumetric H₂ production rate increased linearly with increasing cell density, leading to 2,820 mmol liter⁻¹ h⁻¹ at an optical density at 600 nm (OD₆₀₀) of 18.6, and resulted in the high specific H₂ production rates of 404 ± 6 mmol g⁻¹ h⁻¹. The H₂ productivity indicates the great potential of T. onnurineus strain NA1 for practical application in comparison with H₂-producing microbes. Our result demonstrates that T. onnurineus strain NA1 has a highly efficient metabolic system to thrive on formate in hydrothermal systems.     

Denatured states of yeast cytochrome c induced by heat and guanidinium chloride are structurally and thermodynamically different

A sequence alignment of mammalian cytochromes c with yeast iso-1-cytochrome c (y-cyt-c) shows that the yeast protein contains five extra N-terminal residues. We have been interested in understanding the question: What is the role of these five extra N-terminal residues in folding and stability of the protein? To answer this question we have prepared five deletants of y-cyt-c by sequentially removing these extra residues. During our studies on the wild type (WT) protein and its deletants, we observed that the amount of secondary structure in the guanidinium chloride (GdmCl)-induced denatured (D) state of each protein is different from that of the heat-induced denatured (H) state. This finding is confirmed by the observation of an additional cooperative transition curve of optical properties between H and D states on the addition of different concentrations of GdmCl to the already heat denatured WT y-cyt-c and its deletants at pH 6.0 and 68°C. For each protein, analysis of transition curves representing processes, native (N) state ? D state, N state ? H state, and H state ? D state, was done to obtain Gibbs free energy changes associated with all the three processes. This analysis showed that, for each protein, thermodynamic cycle accommodates Gibbs free energies associated with transitions between N and D states, N and H states, and H and D states, the characteristics required for a thermodynamic function. All these experimental observations have been supported by our molecular dynamics simulation studies.     

Function of Interfacial Prolines at the Transmitter-binding Sites of the Neuromuscular Acetylcholine Receptor

The neuromuscular acetylcholine (ACh) receptor has two conserved prolines in loop D of the complementary subunit at each of its two transmitter-binding sites (α-ϵ and α-δ). We used single-channel electrophysiology to estimate the energy changes caused by mutations of these prolines with regard to unliganded gating (ΔG₀) and the affinity change for ACh that increases the open channel probability (ΔG_B). The effects of mutations of ProD2 (ϵPro-121/δPro-123) were greater than those of its neighbor (ϵPro-120/δPro-122) and were greater at α-ϵ versus α-δ. The main consequence of the congenital myasthenic syndrome mutation ϵProD2-L was to impair the establishment of a high affinity for ACh and thus make ΔG_B less favorable. At both binding sites, most ProD2 mutations decreased constitutive activity (increased ΔG₀). LRYHQG and RL substitutions reduced substantially the net binding energy (made ΔG_B^ACh less favorable) by ≥2 kcal/mol at α-ϵ and α-δ, respectively. Mutant cycle analyses were used to estimate energy coupling between the two ProD2 residues and between each ProD2 and glycine residues (αGly-147 and αGly-153) on the primary (α subunit) side of each binding pocket. The distant binding site prolines interact weakly. ProD2 interacts strongly with αGly-147 but only at α-ϵ and only when ACh is present. The results suggest that in the low to-high affinity change there is a concerted inter-subunit strain in the backbones at ϵProD2 and αGly-147. It is possible to engineer receptors having a single functional binding site by using a α-ϵ or α-δ ProD2-R knock-out mutation. In adult-type ACh receptors, the energy from the affinity change for ACh is approximately the same at the two binding sites (approximately −5 kcal/mol).     

Thermal Characterization,Crystal Field Analysis and In-Band Pumped Laser Performance of Er Doped NaY(WO4)2 Disordered Laser Crystals

Undoped and Er-doped NaY(WO₄)₂ disordered single crystals have been grown by the Czochralski technique. The specific heat and thermal conductivity (κ) of these crystals have been characterized from T = 4 K to 700 K and 360 K, respectively. It is shown that κ exhibits anisotropy characteristic of single crystals as well as a κ(T) behavior observed in glasses, with a saturation mean free phonon path of 3.6 Å and 4.5 Å for propagation along a and c crystal axes, respectively. The relative energy positions and irreducible representations of Stark Er³⁺ levels up to ⁴G_7/2 multiplet have been determined by the combination of experimental low (<10 K) temperature optical absorption and photoluminescence measurements and simulations with a single-electron Hamiltonian including both free-ion and crystal field interactions. Absorption, emission and gain cross sections of the ⁴I_13/2↔⁴I_15/2 laser related transition have been determined at 77 K. The ⁴I_13/2 Er³⁺ lifetime (τ) was measured in the temperature range of 77–300 K, and was found to change from τ (77K) ≈ 4.5 ms to τ (300K) ≈ 3.5 ms. Laser operation is demonstrated at 77 K and 300 K by resonantly pumping the ⁴I_13/2 multiplet at λ≈1500 nm with a broadband (FWHM≈20 nm) diode laser source perfectly matching the 77 K crystal ⁴I_15/2 → ⁴I_13/2 absorption profile. At 77 K as much as 5.5 W of output power were obtained in π-polarized configuration with a slope efficiency versus absorbed pump power of 57%, the free running laser wavelength in air was λ≈1611 nm with the laser output bandwidth of 3.5 nm. The laser emission was tunable over 30.7 nm, from 1590.7 nm to 1621.4 nm, for the same π-polarized configuration.