The Use of the Free Metal – Temperature ‘Phase Diagrams’ for Studies of Single Site Metal Binding Proteins

Typical physico-chemical studies of metal binding proteins are usually aimed at determination of the metal binding constant K for a native protein (K _n), while the significance of the K value for the thermally denatured protein (K _u) is usually underestimated. Meanwhile, metal binding induced shift of thermal denaturation transition of a single site metal binding protein is defined by K _n to K _u ratio, implying that knowledge of both K values is required for full characterization of the system. In the present work, the most universal approach to the studies of single site metal binding proteins, namely construction of a protein “phase diagram” in coordinates of free metal ion concentration – temperature, is considered in detail. The detailed algorithm of construction of the phase diagrams along with underlying mathematic procedures developed here may be of use for studies of other simple protein-target type systems, where target represents low molecular weight ligand. Analysis of the simplest protein-ligand system reveals that thermodynamic properties of apo-protein dictate the maximal possible increase of its affinity to any simple ligand upon thermal denaturation of the protein. Experimental and general problems coupled with the use of the phase diagrams are discussed.     

Amaranth Globulin Polypeptide Heterogeneity

The polypeptides integrating amaranth globulin-p and 11S-globulin were characterized by two-dimensional electrophoresis, ion-exchange chromatography and RP-HPLC. All polypeptides exhibited charge and hydrophobic heterogeneity. Almost all acid (A, pI 5–7) and basic (B, pI 9–10) polypeptides were present in both globulins, and the same happened with the unprocessed M polypeptides with pI in the range of 7–7.5 which fits well with a sequence containing both the A and B polypeptides. There were other polypeptides only present in 11S-globulin, like some of 41 and 16 kDa, which might come from another precursor or be the products of a different processing of the propolypeptide. These results suggested that, although amaranth subunits from different subfamilies are interchangeable in different oligomers, some structural differences between them might affect the assembly of globulin molecules. Structural differences arising from this behavior could account for the different physicochemical properties of globulin molecules.     

HbS-Savaria: The Anti-polymerization Effect of a Single Mutation in Human α-chains

Recombinant α-Savaria globin (α^S49R) was assembled with β^S chains by the alloplex intermediate pathway to generate tetrameric rHbS-Sarvaria (α₂^S49Rβ₂^E6V) that exhibited normal O₂ affinity and co-operatively at pH 7.4. Allosteric effectors, 2,3-DPG, L35, and NaCl increased O₂ affinity by 15%. Bohr effects were similar for rHbS-Savaria and HbS (0.38 ± 0.025 vs. 0.46 ± 0.03, respectively). The C_SAT of HbS increased from 16.7 ± 0.8 to 27.0 ± 1.0 g/dL. Co-polymerization demonstrated inhibition predominantly by the Cis-dimer. Molecular modeling indicated that the positive charge at α-49 generated a strong anion-binding site and reduced flexibility of the CD-region by restricting movement in the E and F helices. The molecular distance between Arg-49 and Asn-78 in the neighboring double strand decreased, and electrostatic repulsion between the inter-double strands increased, resulting in inhibition of polymerization. The Savaria mutation may be useful for the design of super-inhibitory α-chains and gene therapy of sickle cell anemia.     

Equimolar Mixture of 2,2,2-Trifluoroethanol and 4-Chloro-1-butanol is a Stronger Inducer of Molten Globule State: Isothermal Titration Calorimetric and Spectroscopic Studies

A mixture of 4-chloro-1-butanol and 2,2,2-Trifluoroethanol (TFE) has been used to generate Molten globule (MG) state of structurally homologous but functionally different proteins bovine α-lactalbumin and hen egg-white lysozyme. The thermal denaturation was done using UV–Visible spectroscopy. From UV–Visible profile, thermal transition was not observed beyond a particular concentration. There was an indication of molten globule state in case of α-lactalbumin from circular dichroism experiments. By intrinsic tryptophan fluorescence, acrylamide and potassium iodide quenching, 8-anilino-naphthalene sulfonic acid (ANS) binding and energy transfer studies the presence of molten globule state was confirmed. Quantitative characterization of MG state and determining the binding thermodynamics of ANS to the MG state was done using Isothermal Titration Calorimetry (ITC). Results show that α-lactalbumin exists in MG state at a particular concentration but lysozyme does not show features of MG state.     

3D-QSAR studies of Dipeptidyl peptidase IV inhibitors using a docking based alignment

Dipeptidyl peptidase IV (DPP-IV) deactivates the incretin hormones GLP-1 and GIP by cleaving the penultimate proline or alanine from the N-terminal (P1-position) of the peptide. Inhibition of this enzyme will prevent the degradation of the incretin hormones and maintain glucose homeostasis; this makes it an attractive target for the development of drugs for diabetes. This paper reports 3D-QSAR analysis of several DPP-IV inhibitors, which were aligned by the receptor-based technique. The conformation of the molecules in the active site was obtained through docking methods. The QSAR models were generated on two training sets composed of 74 and 25 molecules which included phenylalanine, thiazolidine, and fluorinated pyrrolidine analogs. The 3D-QSAR models are robust with statistically significant r², q², and values. The CoMFA and CoMSIA models were used to design some new inhibitors with several fold higher binding affinity. Figure The CoMFA contours around molecule D1T155 (a) steric contours - favored (green); disfavored (yellow) (b) electrostatic contours - electropositive (blue); electronegative (red)     

Multipole electrostatic potential derived atomic charges in NDDO-methods with <Emphasis Type="Italic">spd</Emphasis>-basis sets

The recently introduced multipole approach for computing the molecular electrostatic potential (MEP) within the semiempirical neglect of diatomic differential overlap (NDDO) framework [Horn AHC, Lin Jr-H., Clark T (2005) Theor Chem Acc 114:159–168] has been used to obtain atomic charges of nearly ab initio quality by scaling the semiempirical MEP. The parameterization set comprised a total of 797 compounds and included not only the newly parameterized AM1* elements Al, Si, P, S, Cl, Ti, Zr, and Mo but also the standard AM1 elements H, C, N, O and F. For comparison, the ZDO-approximated MEP was also calculated analytically in the spd-basis. For the AM1*-optimized structures, single-point calculations at the B3LYP, HF and MP2 levels with the 6-31G(d) and LanL2DZP basis sets were performed to obtain the MEP. The regression analysis of all 12 combinations of semiempirical and ab initio MEP data yielded correlation coefficients of at least 0.99 in all cases. Scaling the analytical and multipole-derived semiempirical MEP by the regression coefficients yielded mean unsigned errors below 2.6 and 1.9 kcal mol⁻¹, respectively. Subsequently, for 22 drug molecules from the World Drug Index, atomic charges were computed according to the RESP procedure using XX/6-31G(d) (XX=B3LYP, HF, MP2) and scaled AM1* multipole MEP; the correlation coefficients obtained are 0.83, 0.85 and 0.83, respectively. Figure: Schematic representation of the atomic charge generation: The molecular electrostatic potential (MEP) is calculated using the AM1* Hamiltonian; then the semiempirical MEP is scaled to DFT or ab initio level, and atomic charges are generated subsequently by the restraint electrostatic potential (RESP) fit method. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible to authorized users. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

Intrinsic Amyloidogenic Behavior of Terminally Protected Alzheimer’s Aβ<Superscript>17–21</Superscript> Peptide: Self-Aggregation and Amyloid-Like Fibril Formation

The terminally protected peptide Boc-Leu-Val-Phe-Phe-Ala-OMe bearing sequence similarity with the central hydrophobic cluster (CHC) of Alzheimer’s Aβ^17–21 peptide self-assembles to produce amyloid-like straight unbranched fibrils from organic solvents. The fibrils readily bind with a physiological dye Congo red (CR) and exhibits green gold birefringence under polarized light, a characteristic feature of amyloid plaque obtained from many neurodegenerative diseases. FTIR spectroscopy and in silico energy minimization study shed some light on the antiparallel supramolecular β-sheet aggregation of the peptide.     

Allozyme diversity in pheasants (<Emphasis Type="Italic">Phasianus</Emphasis> spp.) from breeding stations in Serbia

The pheasant breeds are widely used for restocking of natural populations depleted by hunting. The pheasant population number decline was detected during the 1970s in many hunting areas of Europe. One of its possible reasons might be the loss of adaptability in populations originating from breeding stations, which was caused by inbreeding depression. The aim of this paper was the analysis of genetic variability in pheasant populations from three breeding stations in Vojvodina province (Serbia) by means of allozyme diversity detection. The allozyme variability analysis of pheasants from all three breeding stations revealed polymorphisms at nine loci: Ldh-1, Mor-1, Mor-2, Es-1, Mod-2, Pgd, Gpi-2, Odh, and Sod. The analysis of individuals from three different breeding stations showed mean values of observed heterozygosity of H _o=0.137, polymorphism P _95%=30%, and H/P ratio H/P=0.430, which indicate a normal level of genetic variability for bird populations. Comparative analysis of three pheasant populations showed a high level of interpopulation differentiation.     

Distribution of haplotypes and microsatellite alleles among Asian elephants (<Emphasis Type="Italic">Elephas maximus</Emphasis>) in Thailand

Habitat fragmentation often promotes increased inbreeding depression due to interrupted gene flow between populations. In this study, we demonstrate that Asian elephants most likely also suffer from outbreeding depression due to cryptic speciation. We analysed mitochondrial and nuclear DNA loci from 78 Asian elephants. Haplotype genealogy and analysis of molecular variance revealed two matrilinear clades (α _h, β _h). Microsatellite analyses of individuals grouped according to their haplotype clade (corresponding group of nuclear genotypes called α _nuc and β _nuc) revealed significant genotypic differentiation between α _nuc and β _nuc. Such genealogically differentiated forms in a morphologically uniform species are considered indicative of cryptic speciation. The differentiation was caused by bulls, whereas considering cows only resulted in no differentiation. Such result is best explained by Haldane’s rule whereby hybrid formation between genealogical forms causes lower viability and fertility of heterogametic hybrids. Although the lack of hybrid-specific morphological characteristics renders direct testing of reduced hybrid fitness under natural conditions unfeasible, the effects of Haldane’s rule are demonstrated by reduced male-mediated gene flow between genealogical forms under sympatric conditions, as was indeed suggested by the data found in Asian elephants: male-mediated gene flow between groups α _nuc and β _nuc was much lower than female-mediated gene flow. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Homology modeling and examination of the effect of the D92E mutation on the H5N1 nonstructural protein NS1 effector domain

Virulent H5N1 strains of influenza virus often harbor a D92E point mutation in the nonstructural protein NS1. This crucial mutation has been correlated with increased virulence and/or cytokine resistance, but the structural implications of such a change are still unclear. Furthermore, NS1 protein could also be a potential target for the development of novel antiviral agents against H5N1 strains. Therefore, a reasonable 3D model of H5N1 NS1 is important for the understanding of the molecular basis of increased virulence and the design of novel antiviral agents. Based on the crystal structure of a non-H5N1 NS1 protein, a model of H5N1 NS1 was developed by homology modeling, molecular mechanics and molecular dynamics simulations. It was found that the D92E mutation could result in weakened interactions of the carboxylate side chain with other phosphorylated residues, thereby activating phosphorylation of NS1. Figure Superposition of snapshots picked from the two molecular dynamic (MD) trajectories: a H5N1 NS1 homology model and b non-H5N1 NS1 crystal structure after 0 (green ribbon), 5 (blue ribbon) and 10 ns (pink ribbon) MD simulation