Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

Plasma free amino acid kinetics in rainbow trout (<Emphasis Type="Italic">Oncorhynchus mykiss</Emphasis>) using a bolus injection of <Superscript>15</Superscript>N-labeled amino acids

To gain insight into the metabolic design of the amino acid carrier systems in fish, we injected a bolus of ¹⁵N amino acids into the dorsal aorta in mature rainbow trout (Oncorhynchus mykiss). The plasma kinetic parameters including concentration, pool size, rate of disappearance (R _d), half-life and turnover rate were determined for 15 amino acids. When corrected for metabolic rate, the R _d values obtained for trout for most amino acids were largely comparable to human values, with the exception of glutamine (which was lower) and threonine (which was higher). R _d values ranged from 0.9 μmol 100 g⁻¹ h⁻¹ (lysine) to 22.1 μmol 100 g⁻¹ h⁻¹ (threonine) with most values falling between 2 and 6 μmol 100 g⁻¹ h⁻¹. There was a significant correlation between R _d and the molar proportion of amino acids in rainbow trout whole body protein hydrolysate. Other kinetic parameters did not correlate significantly with whole body amino acid composition. This indicates that an important design feature of the plasma-free amino acids system involves proportional delivery of amino acids to tissues for protein synthesis.     

Studies of the chemical basis of the origin of protein synthesis: Initiation and direction of peptide growth

Summary The data presented in this paper show that the ease of non-enzymatic activation of carboxylic acids by ATP at pH 5 varies directly with the pK_a of the carboxyl group, and is consistent with the idea that it is the protonated form of the carboxyl group which participates in the activation reaction. Consequently, since most N-blocked amino acids have higher pK_a's than do their unblocked forms, they are activated more readily, and we have demonstrated that this principle applies to peptides as well,which are activated more rapidly than single amino acids. We propose that this fact may be partly responsible for the origin of two important features still observed in contemporary protein synthesis: (1) initiation in prokaryotes is accomplished with an N-blocked amino acid, and (2) elongation in all living systems occurs at the carboxyl end of the growing peptide.     

Random-coil chemical shifts of phosphorylated amino acids

The ¹H, ¹³C, ¹⁵N and ³¹ P random-coil chemical shifts and phosphate pK_a values of the phosphorylated amino acids pSer, pThr and pTyr in the protected peptide Ac-Gly-Gly-X-Gly-Gly-NH₂ have been obtained in water at 25°C over the pH range 2 to 9. Analysis of ROESY spectra indicates that the peptides are unstructured. Phosphorylation induces changes in random-coil chemical shifts, some of which are comparable to those caused by secondary structure formation, and are therefore significant in structural analyses based on the chemical shift.     

Graphical analysis of pH-dependent properties of proteins predicted using PROPKA

Background  

Charge states of ionizable residues in proteins determine their pH-dependent properties through their pK_a values. Thus, various theoretical methods to determine ionization constants of residues in biological systems have been developed. One of the more widely used approaches for predicting pK_a values in proteins is the PROPKA program, which provides convenient structural rationalization of the predicted pK_a values without any additional calculations.     

Rheostat Re-Wired: Alternative Hypotheses for the Control of Thioredoxin Reduction Potentials

Thioredoxins are small soluble proteins that contain a redox-active disulfide (CXXC). These disulfides are tuned to oxidizing or reducing potentials depending on the function of the thioredoxin within the cell. The mechanism by which the potential is tuned has been controversial, with two main hypotheses: first, that redox potential (E_m) is specifically governed by a molecular ‘rheostat’—the XX amino acids, which influence the Cys pK_a values, and thereby, E_m; and second, the overall thermodynamics of protein folding stability regulates the potential. Here, we use protein film voltammetry (PFV) to measure the pH dependence of the redox potentials of a series of wild-type and mutant archaeal Trxs, PFV and glutathionine-equilibrium to corroborate the measured potentials, the fluorescence probe BADAN to measure pK_a values, guanidinium-based denaturation to measure protein unfolding, and X-ray crystallography to provide a structural basis for our functional analyses. We find that when these archaeal thioredoxins are probed directly using PFV, both the high and low potential thioredoxins display consistent 2H⁺:2e^- coupling over a physiological pH range, in conflict with the conventional ‘rheostat’ model. Instead, folding measurements reveals an excellent correlation to reduction potentials, supporting the second hypothesis and revealing the molecular mechanism of reduction potential control in the ubiquitous Trx family.     

Role of 2′,6′-dimethyl-l-tyrosine (Dmt) in some opioid lead compounds

Here we evaluated how the interchange of the amino acids 2′,6′-dimethyl-l-tyrosine (Dmt), 2′,6′-difluoro-l-tyrosine (Dft), and tyrosine in position 1 can affect the pharmacological characterization of some reference opioid peptides and pseudopeptides. Generally, Dft and Tyr provide analogues with a similar pharmacological profile, despite different pK_a values. Dmt/Tyr(Dft) replacement gives activity changes depending on the reference opioid in which the modification was made. Whereas, H-Dmt-Tic-Asp1-Bid is a potent and selective δ agonist (MVD, IC₅₀ = 0.12 nM); H-Dft-Tic-Asp1-Bid and H-Tyr-Tic-Asp1-Bid are potent and selective δ antagonists (pA₂ = 8.95 and 8.85, respectively). When these amino acids are employed in the synthesis of deltorphin B and its Dmt¹ and Dft¹ analogues, the three compounds maintain a very similar δ agonism (MVD, IC₅₀ 0.32–0.53 nM) with a decrease in selectivity relative to the Dmt¹ analogue. In the less selective H-Dmt-Tic-Gly1-Bid the replacement of Dmt with Dft and Tyr retains the δ agonism but with a decrease in potency. Antagonists containing the Dmt-Tic pharmacophore do not support the exchange of Dmt with Dft or Tyr.     

Hemiacetal stabilization in a chymotrypsin inhibitor complex and the reactivity of the hydroxyl group of the catalytic serine residue of chymotrypsin

The aldehyde inhibitor Z-Ala-Ala-Phe-CHO has been synthesized and shown by ¹³C-NMR to react with the active site serine hydroxyl group of alpha-chymotrypsin to form two diastereomeric hemiacetals. For both hemiacetals oxyanion formation occurs with a pK_a value of ~ 7 showing that chymotrypsin reduces the oxyanion pK_a values by ~ 5.6 pK_a units and stabilizes the oxyanions of both diastereoisomers by ~ 32 kJ mol^− 1. As pH has only a small effect on binding we conclude that oxyanion formation does not have a significant effect on binding the aldehyde inhibitor. By comparing the binding of Z-Ala-Ala-Phe-CHO with that of Z-Ala-Ala-Phe-H we estimate that the aldehyde group increases binding ~ 100 fold. At pH 7.2 the effective molarity of the active site serine hydroxy group is ~ 6000 which is ~ 7 × less effective than with the corresponding glyoxal inhibitor. Using ¹H-NMR we have shown that at both 4 and 25 °C the histidine pK_a is ~ 7.3 in free chymotrypsin and it is raised to ~ 8 when Z-Ala-Ala-Phe-CHO is bound. We conclude that oxyanion formation only has a minor role in raising the histidine pK_a and that the aldehyde hydrogen must be replaced by a larger group to raise the histidine pK_a > 10 and give stereospecific formation of tetrahedral intermediates. The results show that a large increase in the pK_a of the active site histidine is not needed for the active site serine hydroxyl group to have an effective molarity of 6000.     

In Depth Analysis on the Binding Sites of Adamantane Derivatives in HCV (Hepatitis C Virus) p7 Channel Based on the NMR Structure

Background

The recently solved solution structure of HCV (hepatitis C virus) p7 ion channel provides a solid structure basis for drug design against HCV infection. In the p7 channel the ligand amantadine (or rimantadine) was determined in a hydrophobic pocket. However the pharmocophore (−NH₂) of the ligand was not assigned a specific binding site.

Results

The possible binding sites for amino group of adamantane derivatives is studied based on the NMR structure of p7 channel using QM calculation and molecular modeling. In the hydrophobic cavity and nearby three possible binding sites are proposed: His17, Phe20, and Trp21. The ligand binding energies at the three binding sites are studied using high level QM method CCSD(T)/6–311+G(d,p) and AutoDock calculations, and the interaction details are analyzed. The potential application of the binding sites for rational inhibitor design are discussed.

Conclusions

Some useful viewpoints are concluded as follows. (1) The amino group (−NH₂) of adamantane derivatives is protonated (−NH₃ ⁺), and the positively charged cation may form cation-π interactions with aromatic amino acids. (2) The aromatic amino acids (His17, Phe20, and Trp21) are the possible binding sites for the protonated amino group (−NH₃ ⁺) of adamantane derivatives, and the cation-π bond energies are 3 to 5 times stronger than the energies of common hydrogen bonds. (3) The higher inhibition potent of rimantadine than amantadine probably because of its higher pK_a value (pK_a = 10.40) and the higher positive charge in the amino group. The potential application of p7 channel structure for inhibitor design is discussed.     

Characterization and Implications of the Cell Surface Reactivity of Calothrix sp. Strain KC97

The cell surface reactivity of the cyanobacterium Calothrix sp. strain KC97, an isolate from the Krisuvik hot spring, Iceland, was investigated in terms of its proton binding behavior and charge characteristics by using acid-base titrations, electrophoretic mobility analysis, and transmission electron microscopy. Analysis of titration data with the linear programming optimization method showed that intact filaments were dominated by surface proton binding sites inferred to be carboxyl groups (acid dissociation constants [pK_a] between 5.0 and 6.2) and amine groups (mean pK_a of 8.9). Sheath material isolated by using lysozyme and sodium dodecyl sulfate generated pK_a spectra similarly dominated by carboxyls (pK_a of 4.6 to 6.1) and amines (pK_a of 8.1 to 9.2). In both intact filaments and isolated sheath material, the lower ligand concentrations at mid-pK_a values were ascribed to phosphoryl groups. Whole filaments and isolated sheath material displayed total reactive-site densities of 80.3 × 10⁻⁵ and 12.3 × 10⁻⁵ mol/g (dry mass) of cyanobacteria, respectively, implying that much of the surface reactivity of this microorganism is located on the cell wall and not the sheath. This is corroborated by electrophoretic mobility measurements that showed that the sheath has a net neutral charge at mid-pHs. In contrast, unsheathed cells exhibited a stronger negative-charge characteristic. Additionally, transmission electron microscopy analysis of ultrathin sections stained with heavy metals further demonstrated that most of the reactive binding sites are located upon the cell wall. Thus, the cell surface reactivity of Calothrix sp. strain KC97 can be described as a dual layer composed of a highly reactive cell wall enclosed within a poorly reactive sheath.     

Chemical versus Mechanical Perturbations on the Protonation State of Arginine in Complex Lipid Membranes: Insights from Microscopic pKa Calculations

Charged amino acids such as Arginine play important roles in many membrane-mediated biological processes such as voltage gating of ion channels and membrane translocation of cell penetration peptides. It is well established that local membrane deformation and formation of water defects are crucial to the stabilization of charged species in contact with the membrane, which suggests that mechanical properties of the membrane are relevant although a clear connection has not been established. As a quantitative measure, we study how changes in the composition and therefore mechanical properties of a lipid bilayer influence the pK_a of Arg in the membrane center using free energy simulations. Compared to previous studies in a single-component lipid bilayer containing saturated lipids or lipids with a modest degree of unsaturation, substantially larger pK_a shifts are observed in the presence of highly unsaturated lipid tails and cholesterol. Moreover, the underlying molecular mechanisms for the pK_a perturbation are distinct in different systems, with the unsaturated lipid tails mainly destabilizing the charged state of Arg and the cholesterol stabilizing the neutral state of Arg. The observed behaviors in both cases are at odds with predictions based on mechanical considerations at a mesoscopic level—highlighting that, while mechanical considerations are useful for stimulating hypothesis, their applicability to dissecting phenomena at the molecular-length scale is rather limited.     

Molecular interactions involved in proton-dependent gating in KcsA potassium channels

The bacterial potassium channel KcsA is gated open by the binding of protons to amino acids on the intracellular side of the channel. We have identified, via channel mutagenesis and x-ray crystallography, two pH-sensing amino acids and a set of nearby residues involved in molecular interactions that influence gating. We found that the minimal mutation of one histidine (H25) and one glutamate (E118) near the cytoplasmic gate completely abolished pH-dependent gating. Mutation of nearby residues either alone or in pairs altered the channel’s response to pH. In addition, mutations of certain pairs of residues dramatically increased the energy barriers between the closed and open states. We proposed a Monod–Wyman–Changeux model for proton binding and pH-dependent gating in KcsA, where H25 is a “strong” sensor displaying a large shift in pK_a between closed and open states, and E118 is a “weak” pH sensor. Modifying model parameters that are involved in either the intrinsic gating equilibrium or the pK_a values of the pH-sensing residues was sufficient to capture the effects of all mutations.     

A multivariate study of the relationship between the genetic code and the physical-chemical properties of amino acids

Summary The 20 naturally occurring amino acids are characterized by 20 variables: pK_{NH 2}, pK_COOH, pI, molecular weight, substituent van der Waals volume, seven¹H and¹³C nuclear magnetic resonance shift variables, and eight hydrophobicity-hydrophilicity scales. The 20-dimensional data set is reduced to a few new dimensions by principal components analysis. The three first principal components reveal relationships between the properties of the amino acids and the genetic code. Thus the amino acids coded for by adenosine (A), uracil (U), or cytosine (C) in their second codon position (corresponding to U, A, or G in the second anticodon position) are grouped in these components. No grouping was detected for the amino acids coded for by guanine (G) in the second codon position (corresponding to C in the second anticodon position). The results show that a relationship exists between the physical-chemical properties of the amino acids and which of the A (U), U (A), or C (G) nucleotide is used in the second codon (anticodon) position. The amino acids coded for by G (C) in the second codon (anticodon) position do not participate in this relationship.     

Individual heme a and heme a3 contributions to the Soret absorption spectrum of the reduced bovine cytochrome c oxidase

Bovine cytochrome c oxidase (CcO) contains two hemes, a and a₃, chemically identical but differing in coordination and spin state. The Soret absorption band of reduced aa₃-type cytochrome c oxidase consists of overlapping bands of the hemes a²⁺ and a₃²⁺. It shows a peak at ～444 nm and a distinct shoulder at ～425 nm. However, attribution of individual spectral lineshapes to hemes a²⁺ and a₃²⁺ in the Soret is controversial. In the present work, we characterized spectral contributions of hemes a²⁺ and a₃²⁺ using two approaches. First, we reconstructed bovine CcO heme a²⁺ spectrum using a selective Ca²⁺-induced spectral shift of the heme a²⁺. Second, we investigated photobleaching of the reduced Thermus thermophilus ba₃- and bovine aa₃-oxidases in the Soret induced by femtosecond laser pulses in the Q-band. The resolved spectra show splitting of the electronic B_0x-, B_0y-transitions of both reduced hemes. The heme a²⁺ spectrum is shifted to the red relative to heme a₃²⁺ spectrum. The ～425 nm shoulder is mostly attributed to heme a₃²⁺.     

Oncogenic KRAS G12C: Kinetic and redox characterization of covalent inhibition

The recent development of mutant-selective inhibitors for the oncogenic KRAS^G12C allele has generated considerable excitement. These inhibitors covalently engage the mutant C12 thiol located within the phosphoryl binding loop of RAS, locking the KRAS^G12C protein in an inactive state. While clinical trials of these inhibitors have been promising, mechanistic questions regarding the reactivity of this thiol remain. Here, we show by NMR and an independent biochemical assay that the pK_a of the C12 thiol is depressed (pK_a ∼7.6), consistent with susceptibility to chemical ligation. Using a validated fluorescent KRAS^Y137W variant amenable to stopped-flow spectroscopy, we characterized the kinetics of KRAS^G12C fluorescence changes upon addition of ARS-853 or AMG 510, noting that at low temperatures, ARS-853 addition elicited both a rapid first phase of fluorescence change (attributed to binding, K_d = 36.0 ± 0.7 μM) and a second, slower pH-dependent phase, taken to represent covalent ligation. Consistent with the lower pK_a of the C12 thiol, we found that reversible and irreversible oxidation of KRAS^G12C occurred readily both in vitro and in the cellular environment, preventing the covalent binding of ARS-853. Moreover, we found that oxidation of the KRAS^G12C Cys12 to a sulfinate altered RAS conformation and dynamics to be more similar to KRAS^G12D in comparison to the unmodified protein, as assessed by molecular dynamics simulations. Taken together, these findings provide insight for future KRAS^G12C drug discovery efforts, and identify the occurrence of G12C oxidation with currently unknown biological ramifications.     

New pentafluorophenyl complexes with phosphine-amide ligands

A series of nickel(II) and palladium(II) complexes containing one or two pentafluorophenyl ligands and the phosphino-amides o-Ph₂PC₆H₄CONHR [R = ⁱPr (a), Ph (b)] displaying different coordination modes have been synthesised. The chelating ability of these ligands and the influence of both coligands and the metal centre in their potential hemilabile behaviour have been explored. The crystal structure of (b) has been determined and reveals N–H⋯O intermolecular hydrogen bonding. Bis-pentafluorophenyl derivatives [M(C₆F₅)₂(o-Ph₂PC₆H₄CO-NHR)] [M = Ni; R = ⁱPr (1a); R = Ph (1b); M = Pd; R = ⁱPr (2a); R = Ph (2b)] in which (a) and (b) act as rigid P, O-chelating ligands were readily prepared from the labile precursors cis-[M(C₆F₅)₂(PhCN)₂]. X-ray structures of (1a), (1b) and (2a) have been established, allowing an interesting comparative structural discussion. Dinuclear [{Pd(C₆F₅)(tht)(μ-Cl)}₂] reacted with (a) and (b) yielding the monopentafluorophenyl complexes [Pd(C₆F₅)Cl{PPh₂(C₆H₄–CONH–R)}] (R = ⁱPr (3a), Ph (3b)) that showed a P, O-chelating behaviour of the ligands, confirmed by the crystal structure determination of (3a). New cationic palladium(II) complexes in which (a) and (b) behave as P-monodentate ligands have been synthesised by reacting them with [{Pd(C₆F₅)(tht)(μ-Cl)}₂], stoichiometric Ag(O₃SCF₃) and external chelating reagents such as cod [Pd(C₆F₅)(cod){PPh₂(C₆H₄-CONH-R)}](O₃SCF₃)(R = ⁱPr (4a), Ph (4b)) and 2,2′-bipy [Pd(C₆F₅)(bipy){PPh₂(C₆H₄-CONH-R)}](O₃SCF₃) (R = ⁱPr (5a), Ph (5b)). When chloride abstraction in [{Pd(C₆F₅)(tht)(μ-Cl)}₂] is promoted by means of a dithioanionic salt as dimethyl dithiophospate in the presence of (a) or (b), the corresponding neutral complexes [Pd(C₆F₅){S(S)P(OMe)₂}{PPh₂(C₆H₄-CONH-R)}] (R = ⁱPr (6a), Ph (6b)) were obtained.     

Effect of nitrogen deficiency on the ion-exchange properties of cell wall polymers from wheat roots

The ion-exchange properties of cell wall polymers isolated from the roots of wheat (Triticum aestivum L.) plants grown on either nitrate-free (N-deficient) or nitrate-containing (+N) hydroponic nutrient medium have been investigated. Irrespective of the nitrogen nutrition regimen, the studied cell walls contained four types of ion-exchange groups: primary amino groups of structural proteins (pK_a < 3), carboxyl groups of polygalacturonic acid in pectin (pK _a ~4.7), carboxyl groups of hydroxycinnamic acids (pK _a ~7.3), and phenolic OH-groups of lignin (pK_a ~10.2). The quantitative ratio between these types of ion-exchange groups, the mass fraction of cell walls in the dry weight of roots, and the swelling coefficient of cell walls depended on the nitrate presence in the growing medium. Compared to the +N variant, the N-deficient variant was characterized by a 2.4 times higher content of phenolic OH-groups in cell walls and 1.24 times higher mass fraction of cell walls; at the same time, the swelling coefficient for this variant was lower by 10%. The obtained data indicate that nitrogen deficiency results in a formation of thicker root cell walls with a higher degree of polymer cross-linking that may be caused by the increased lignin content.     

Incorporation of [15N]Ammonia by the Cellulolytic Ruminal Bacteria Fibrobacter succinogenes BL2, Ruminococcus albus SY3, and Ruminococcus flavefaciens 17

The origin of cell nitrogen and amino acid nitrogen during growth of ruminal cellulolytic bacteria in different growth media was investigated by using ¹⁵NH₃. At high concentrations of peptides (Trypticase, 10 g/liter) and amino acids (15.5 g/liter), significant amounts of cell nitrogen of Fibrobacter succinogenes BL2 (51%), Ruminococcus flavefaciens 17 (43%), and Ruminococcus albus SY3 (46%) were derived from non-NH₃-N. With peptides at 1 g/liter, a mean of 80% of cell nitrogen was from NH₃. More cell nitrogen was formed from NH₃ during growth on cellobiose compared with growth on cellulose in all media. Phenylalanine was essential for F. succinogenes, and its ¹⁵N enrichment declined more than that of other amino acids in all species when amino acids were added to the medium.     

2,3-Diamino acid modifying 3S-tetrahydroisoquinoline-3-carboxylic acids: Leading to a class of novel agents with highly unfolded conformation,selective in vitro anti-platelet aggregation and potent in vivo anti-thrombotic activity

In the preparation of anti-thrombotic agents the 2- and 3-positions of 3S-tetra-hydroisoquinoline-3-carboxylic acid (THIQA) were simultaneously modified with amino acids to form 20 novel N-(3S-N-aminoacyl-1,2,3,4-tetrahydroisoquinoline-3-carbonyl)amino acids (8a–t). On an in vitro platelet aggregation model 8a–t selectively inhibit ADP-induced platelet aggregation and their IC₅₀ values are leas than 3.5 nM. On an extracorporeal circulation of arterioveinos cannula model of rats both orally and intraveously effective doses of 8a–t are less than 30 nmol/kg. Cerius² based stereoview of explores 8a–t having highly unfolded conformation. 3D QSAR analysis gives the importance of the unfolded conformation to high in vitro anti-platelet aggregation and in vivo anti-thrombotic potency rational understanding.     

Importance of Aerodynamic Resistance to Water Use Efficiency in Three Conifers under Field Conditions

The quantitative importance of aerodynamic resistance to H₂O vapor and CO₂ exchange was determined for shoots from saplings of three conifers (Abies lasiocarpa [Hook] Nutt., Pinus contorta Dougl., Juniperus communis L.) under natural conditions in the field. A combination of relatively low stomatal resistances (<300 seconds per centimeter) and low wind speeds (<30 centimeters per second) led to substantial contributions of the aerodynamic resistance (R_wv^a) to water use efficiency (WUE = photosynthesis/transpiration) for all three species. For A. lasiocarpa, transpiration was calculated to be 44% less and photosynthesis 17% less due to the presence of R_wv^a, which led to a predicted increase in WUE of 57% compared to the calculated WUE when R_wv^a was assumed negligible. Similar increases in WUE were computed for P. contorta (48%) with somewhat smaller values for J. communis (34%). These results are discussed in terms of the estimated importance of R_wv^a on water and photosynthetic relations of plants that have relatively low stomatal resistances and grow in microhabitats with low winds.