A Water Mediated Electrostatic Interaction Gives Thermal Stability to the “Tail” Region of the GrpE Protein from <Emphasis Type="Italic">E. coli</Emphasis>

The GrpE protein from E. coli is a homodimer with an unusual structure of two long paired α-helices from each monomer interacting in a parallel arrangement to form a “tail” at the N-terminal end. Using site-directed mutagenesis, we show that there is a key electrostatic interaction involving R57 (mediated by a water molecule) that provides thermal stability to this “tail” region. The R57A mutant showed a drop in T _m of 8.5°C and a smaller ΔH _u (unfolding) compared to wild-type for the first unfolding transition, but no significant decrease in dimer stability as shown through equilibrium analytical ultracentrifugation studies. Another mutant (E94A) at the dimer interface showed a decrease in ΔH _u but no drop in T _m for the second unfolding transition and a slight increase in dimer stability.     

Analysis of Protein Three-Dimension Structure Using Amino Acids Depths

The issue of amino acid depth in proteins gives important insights to our understanding of protein’s three-dimensional structure. There has already been much research done in mathematical and statistical sciences regarding the general definitions, properties and algorithms describing the particle depth of spatially extended systems. We constructed a method of calculating the amino acids depths and applied it to a set of 527 protein structures. We propose the introduction of amino acid depth tendency factors for three-dimensional structures of proteins. The depth tendency factors relate not only to the hydrophobicity indices but also to the electrostatic charge. We found a relationship between the protein size and the number of residues using the distance between the deepest residue and surface residues. We made a prediction regarding the number of residues on the surface of a protein, the deepest amino acid, and the average depth, all of which are fitted well to a linear functional relationship with the length of the protein. Finally, we have predicted the depths of multiple peptides in protein’s three-dimension structure. Electronic supplementary material The online version of this article () contains supplementary material, which is available to authorized users.     

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.     

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.