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     

Computational study of enantioselective interaction between C<Subscript>60</Subscript> fullerene and its derivatives with L-histidine

The mechanism of the enantioselective binding of L-histidine with C₆₀ fullerene and its derivatives, (1,2-methanofullerene C₆₀)-61-carboxylic acid, diethyl (1,2-methanofullerene C₆₀)-61-61-dicarboxylate and tert-butyl (1,2-methanofullerene C₆₀)-61-carboxylate based chiral selectors was studied by quantum chemical calculations. All the molecules were fully optimized at RHF/6-31G* basis set. Relative energies between the different complexes were subsequently estimated with single-point electronic energies computed using Møller-Plesset perturbation theory (MP2). Stability and feasibility of all the generated structures were supported by their respective energy minima and fundamental frequencies. It was observed that interaction of fullerene derivatives with L-histidine is due to the existence of hydrogen bonding forces during the complex formation. The intermolecular forces, flow of atomic charges, binding energy, hardness, dipole moment and localization of electrostatic potential are in agreement with enantioselective interaction of L-histidine with C₆₀ fullerene and its derivatives. It is found that theoretical evaluation to be consistent with the experimental data.     

On the relative stabilities of the alkali cations 222 cryptates in the gas phase and in water-methanol solution

The relative stabilities of the alkali [M ⊂ 222]⁺ cryptates (M = Na, K, Rb and Cs) in the gas phase and in solution (80:20 v/v methanol:water mixture) at 298 K, are computed using a combination of ab initio quantum-chemical calculations (HF/6-31G and MP2/6-31+G*//HF/6-31+G*) and explicit-solvent Monte Carlo free-energy simulations. The results suggest that the relative stabilities of the cryptates in solution are due to a combination of steric effects (compression of large ions within the cryptand cavity), electronic effects (delocalization of the ionic charge onto the cryptand atoms) and solvent effects (dominantly the ionic dessolvation penalty). Thus, the relative stabilities in solution cannot be rationalized solely on the basis of a simple match or mismatch between the ionic radius and the cryptand cavity size as has been suggested previously. For example, although the [K ⊂ 222]⁺ cryptate is found to be the most stable in solution, in agreement with experimental data, it is the [Na ⊂ 222]⁺ cryptate that is the most stable in the gas phase. The present results provide further support to the notion that the solvent in which supramolecules are dissolved plays a key role in modulating molecular recognition processes. Figure Alkali cryptates [M ⊂ 222]+ (M = Na, K, Rb and Cs) relative stabilities in gas and methanol:water solution: solvent effects and molecular recognition

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.