Conformational Motions and Functionally Key Residues for Vitamin B12 Transporter BtuCD–BtuF Revealed by Elastic Network Model with a Function-Related Internal Coordinate

BtuCD–BtuF from Escherichia coli is a binding protein-dependent adenosine triphosphate (ATP)-binding cassette (ABC) transporter system that uses the energy of ATP hydrolysis to transmit vitamin B12 across cellular membranes. Experimental studies have showed that during the transport cycle, the transporter undergoes conformational transitions between the “inward-facing” and “outward-facing” states, which results in the open–closed motions of the cytoplasmic gate of the transport channel. The opening–closing of the channel gate play critical roles for the function of the transporter, which enables the substrate vitamin B12 to be translocated into the cell. In the present work, the extent of opening of the cytoplasmic gate was chosen as a function-related internal coordinate. Then the mean-square fluctuation of the internal coordinate, as well as the cross-correlation between the displacement of the internal coordinate and the movement of each residue in the protein, were calculated based on the normal mode analysis of the elastic network model to analyze the function-related motions encoded in the structure of the system. In addition, the key residues important for the functional motions of the transporter were predicted by using a perturbation method. In order to facilitate the calculations, the internal coordinate was introduced as one of the axes of the coordinate space and the conventional Cartesian coordinate space was transformed into the internal/Cartesian space with linear approximation. All the calculations were carried out in this internal/Cartesian space. Our method can successfully identify the functional motions and key residues for the transporter BtuCD–BtuF, which are well consistent with the experimental observations.     

Analysis of Conformational Motions and Residue Fluctuations for Escherichia coli Ribose-Binding Protein Revealed with Elastic Network Models

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The open (ligand-free) and closed (ligand-bound) forms of RBP have been found. Here we investigate the conformational motions and residue fluctuations of the RBP by analyzing the modes of motion with two coarse-grained elastic network models, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). The calculated B-factors in both the calculated models are in good agreement with the experimentally determined B-factors in X-ray crystal structures. The slowest mode analysis by GNM shows that both forms have the same motion hinge axes around residues Ser103, Gln235, Asp264 and the two domains of both structures have similar fluctuation range. The superposition of the first three dominant modes of ANM, consisting of the rotating, bending and twisting motions of the two forms, accounts for large rearrangement of domains from the ligand-free (open) to ligand-bound (closed) conformation and thus constitutes a critical component of the RBP’s functions. By analyzing cross-correlations between residue fluctuation and the difference-distance plot, it is revealed that the conformational change can be described as a rigid rotation of the two domains with respect to each other, whereas the internal structure of the two domains remains largely intact. The results directly indicate that the dominant dynamic characteristics of protein structures can be captured from their static native state using coarse-grained models.     

Identification and Validation of Aspartic Acid Semialdehyde Dehydrogenase as a New Anti-Mycobacterium Tuberculosis Target

Aspartic acid semialdehyde dehydrogenase (ASADH) lies at the first branch point in the essential aspartic acid biosynthetic pathway that is found in bacteria and plants but is absent from animals. Mutations in the asadh gene encoding ASADH produce an inactive enzyme, which is lethal. Therefore, in this study, we investigated the hypothesis that ASADH represents a new anti-Mycobacterium tuberculosis (MTB) target. An asadh promoter-replacement mutant MTB, designated MTB::asadh, in which asadh gene expression is regulated by pristinamycin, was constructed to investigate the physiological functions of ASADH in the host bacteria. Bacterial growth was evaluated by monitoring OD₆₀₀ and ASADH expression was analyzed by Western blotting. The results showed that the growth and survival of MTB::asadh was completely inhibited in the absence of the inducer pristinamycin. Furthermore, the growth of the mutant was rigorously dependent on the presence of the inducer in the medium. The starved mutant exhibited a marked reduction (approximately 80%) in the cell wall materials compared to the wild-type, in addition to obvious morphological differences that were apparent in scanning electron microscopy studies; however, with the addition of pristinamycin, the cell wall contents and morphology similar to those of the wild-type strain were recovered. The starved mutant also exhibited almost no pathogenicity in an in vitro model of infection using mouse macrophage J774A.1 cells. The mutant showed a concentration-dependent recovery of pathogenicity with the addition of the inducer. These findings implicate ASADH as a promising target for the development of novel anti-MTB drugs.     

The Intrinsic Dynamics and Unfolding Process of an Antibody Fab Fragment Revealed by Elastic Network Model

Antibodies have been increasingly used as pharmaceuticals in clinical treatment. Thermal stability and unfolding process are important properties that must be considered in antibody design. In this paper, the structure-encoded dynamical properties and the unfolding process of the Fab fragment of the phosphocholine-binding antibody McPC603 are investigated by use of the normal mode analysis of Gaussian network model (GNM). Firstly, the temperature factors for the residues of the protein were calculated with GNM and then compared with the experimental measurements. A good result was obtained, which provides the validity for the use of GNM to study the dynamical properties of the protein. Then, with this approach, the mean-square fluctuation (MSF) of the residues, as well as the MSF in the internal distance (MSFID) between all pairwise residues, was calculated to investigate the mobility and flexibility of the protein, respectively. It is found that the mobility and flexibility of the constant regions are higher than those of the variable regions, and the six complementarity-determining regions (CDRs) in the variable regions also exhibit relative large mobility and flexibility. The large amplitude motions of the CDRs are considered to be associated with the immune function of the antibody. In addition, the unfolding process of the protein was simulated by iterative use of the GNM. In our method, only the topology of protein native structure is taken into account, and the protein unfolding process is simulated through breaking the native contacts one by one according to the MSFID values between the residues. It is found that the flexible regions tend to unfold earlier. The sequence of the unfolding events obtained by our method is consistent with the hydrogen-deuterium exchange experimental results. Our studies imply that the unfolding behavior of the Fab fragment of antibody McPc603 is largely determined by the intrinsic dynamics of the protein.     

Cloning and Polymorphisms of Yak Lactate Dehydrogenase b Gene

The main objective of this work was to study the unique polymorphisms of the lactate dehydrogenase-1 (LDH1) gene in yak (Bos grunniens). Native polyacrylamide gel electrophoresis revealed three phenotypes of LDH1 (a tetramer of H subunit) in yak heart and longissimus muscle extracts. The corresponding gene, ldhb, encoding H subunits of three LDH1 phenotypes was obtained by RT-PCR. A total of six nucleotide differences were detected in yak ldhb compared with that of cattle, of which five mutations cause amino acid substitutions. Sequence analysis shows that the G896A and C689A, mutations of ldhb gene, result in alterations of differently charged amino acids, and create the three phenotypes (F, M, and S) of yak LDH1. Molecular modeling of the H subunit of LDH indicates that the substituted amino acids are not located within NAD⁺ or substrate binding sites. PCR-RFLP examination of G896A mutation demonstrated that most LDH1-F samples are actually heterozygote at this site. These results help to elucidate the molecular basis and genetic characteristic of the three unique LDH1 phenotypes in yak.     

Biochemical Characterization and Complete Conversion of Coenzyme Specificity of Isocitrate Dehydrogenase from Bifidobacterium longum

Bifidobacterium longum is a very important gram-positive non-pathogenic bacterium in the human gastrointestinal tract for keeping the digestive and immune system healthy. Isocitrate dehydrogenase (IDH) from B. longum (BlIDH), a novel member in Type II subfamily, was overexpressed, purified and biochemically characterized in detail. The active form of BlIDH was an 83-kDa homodimer. Kinetic analysis showed BlIDH was a NADP⁺-dependent IDH (NADP-IDH), with a 567- and 193-fold preference for NADP⁺ over NAD⁺ in the presence of Mg²⁺ and Mn²⁺, respectively. The maximal activity for BlIDH occurred at 60 °C (with Mn²⁺) and 65 °C (with Mg²⁺), and pH 7.5 (with Mn²⁺) and pH 8.0 (with Mg²⁺). Heat-inactivation profiles revealed that BlIDH retained 50% of maximal activity after incubation at 45 °C for 20 min with either Mn²⁺ or Mg²⁺. Furthermore, the coenzyme specificity of BlIDH can be completely reversed from NADP⁺ to NAD⁺ by a factor of 2387 by replacing six residues. This current work, the first report on the coenzyme specificity conversion of Type II NADP-IDHs, would provide better insight into the evolution of NADP⁺ use by the IDH family.     

The Effect of C-Terminal Helix on the Stability of FF Domain Studied by Molecular Dynamics Simulation

To investigate the effect of C-terminal helix on the stability of the FF domain, we studied the native domain FF3-71 from human HYPA/FBP11 and the truncated version FF3-60 with C-terminal helix being deleted by molecular dynamics simulations with GROMACS package and GROMOS 43A1 force field. The results indicated that the structures of truncated version FF3-60 were evident different from those of native partner FF3-71. Compared with FF3-71, the FF3-60 lost some native contacts and exhibited some similar structural characters to those of intermediate state. The C-terminal helix played a major role in stabilizing the FF3-71 domain. To a certain degree, the FF domain had a tendency to form an intermediate state without the C-terminal helix. In our knowledge, this was the first study to examine the role of C-terminal helix of FF domain in detail by molecular dynamics simulations, which was useful to understand the three-state folding mechanism of the small FF domain.     

Genetic Diversity of the Critically Endangered Thuja sutchuenensis Revealed by ISSR Markers and the Implications for Conservation

Thuja sutchuenensis Franch. is a critically endangered plant endemic to the North-East Chongqing, China. Genetic variation was studied to assess the distribution of genetic diversity within and among seven populations from the single remnant locations, using inter-simple sequence repeat (ISSR) markers. A total of 15 primers generated 310 well defined bands, with an average of 20.7 bands per primer. The seven populations revealed a relatively high level of genetic diversity in the species. The percentage of polymorphic bands, Nei’s gene diversity and Shannon’s information index at the population and species level were 76.1%, 0.155, 0.252 and 100%, 0.165, 0.295, respectively. A low level of genetic differentiation among populations (G_ST = 0.102), in line with the results of Analyses of Molecular Variance (AMOVA), and a high level of gene flow (N_m = 4.407) were observed. Both the Unweighted Pair Group Method with Arithmatic Mean (UPGMA) cluster analysis and Principal Coordinates Analysis (PCoA) supported the grouping of all seven populations into two groups. In addition, Mantel test revealed no significant correlation between genetic and geographical distances (r = 0.329, p = 0.100). The low genetic differentiation among populations implies that the conservation efforts should aim to preserve all the extant populations of this endangered species.     

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.     

Production of (R)-3-Quinuclidinol by E. coli Biocatalysts Possessing NADH-Dependent 3-Quinuclidinone Reductase (QNR or bacC) from Microbacterium luteolum and Leifsonia Alcohol Dehydrogenase (LSADH)

We found two NADH-dependent reductases (QNR and bacC) in Microbacterium luteolum JCM 9174 (M. luteolum JCM 9174) that can reduce 3-quinuclidinone to optically pure (R)-(−)-3-quinuclidinol. Alcohol dehydrogenase from Leifsonia sp. (LSADH) was combined with these reductases to regenerate NAD⁺ to NADH in situ in the presence of 2-propanol as a hydrogen donor. The reductase and LSADH genes were efficiently expressed in E. coli cells. A number of constructed E. coli biocatalysts (intact or immobilized) were applied to the resting cell reaction and optimized. Under the optimized conditions, (R)-(−)-3-quinuclidinol was synthesized from 3-quinuclidinone (15% w/v, 939 mM) giving a conversion yield of 100% for immobilized QNR. The optical purity of the (R)-(−)-3-quinuclidinol produced by the enzymatic reactions was >99.9%. Thus, E. coli biocatalysis should be useful for the practical production of the pharmaceutically important intermediate, (R)-(−)-3-quinuclidinol.     

Differences between Mice and Humans in Regulation and the Molecular Network of Collagen,Type III,Alpha-1 at the Gene Expression Level: Obstacles that Translational Research Must Overcome

Collagen, type III, alpha-1 (COL3A1) is essential for normal collagen I fibrillogenesis in many organs. There are differences in phenotypes of mutations in the COL3A1 gene in humans and mutations in mice. In order to investigate whether the regulation and gene network of COL3A1 is the same in healthy populations of mice and humans, we compared the quantitative trait loci (QTL) that regulate the expression level of COL3A1 and the gene network of COL3A1 pathways between humans and mice using whole genome expression profiles. Our results showed that, for the regulation of expression of Col3a1 in mice, an eQTL on chromosome (Chr) 12 regulates the expression of Col3a1. However, expression of genes in the syntenic region on human Chr 7 has no association with the expression level of COL3A1. For the gene network comparison, we identified 44 top genes whose expression levels are strongly associated with that of Col3a1 in mice. We next identified 41 genes strongly associated with the expression level of COL3A1 in humans. There are a few but significant differences in the COL3A1 gene network between humans and mice. Several genes showed opposite association with expression of COL3A1. These genes are known to play important roles in development and function of the extracellular matrix of the lung. Difference in the molecular pathway of key genes in the COL3A1 gene network in humans and mice suggest caution should be used in extrapolating results from models of human lung diseases in mice to clinical lung diseases in humans. These differences may influence the efficacy of drugs in humans whose development employed mouse models.     

Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers

Rheum officinale Baill., an important but endangered medicinal herb, is endemic to China. Inter-simple sequence repeat (ISSR) markers were employed to investigate the genetic diversity and differentiation of 12 populations of R. officinale. Thirteen selected primers yielded 189 bright and discernible bands, with an average of 14.54 per primer. The genetic diversity was low at the population level, but pretty high at the species level (H = 0.1008, I = 0.1505, PPB = 28.95% vs. H = 0.3341, I = 0.5000, PPB = 95.24%, respectively) by POPGENE analysis. Analysis of molecular variance (AMOVA) showed that the genetic variation was found mainly among populations (74.38%), in line with the limited gene flow (N(m) = 0.2766) among populations. Mantel test revealed a significant correlation between genetic and geographic distances (r = 0.5381, P = 0.002), indicating the role of geographic isolation in shaping the present population genetic structure. Both Bayesian analysis and UPGMA cluster analysis demonstrated the similar results. Our results imply that the conservation efforts should aim to preserve all the extant populations of this endangered species, and cultivation is proposed in this study.     

Purification and Properties of an Insecticidal Metalloprotease Produced by Photorhabdus luminescens Strain 0805-P5G,the Entomopathogenic Nematode Symbiont

A total of 13 Photorhabdus luminescens strains were screened for proteolytic activity. The P. luminescens strain 0805-P5G had the highest activity on both skim milk and gelatin plates. The protease was purified to electrophoretical homogeneity by using a two-step column chromatographic procedure. It had a molecular weight of 51.8 kDa, as determined by MALDI-TOF mass spectrometry. The optimum pH, temperature, as well as pH and thermal stabilities were 8, 60 °C, 5–10, and 14–60 °C, respectively. It was completely inhibited by EDTA and 1,10-phenanthroline. Bioassay of the purified protease against Galleria mellonella by injection showed high insecticidal activity. The protease also showed high oral toxicity to the diamondback moth (Plutella xylostella) of a Taiwan field-collected strain, but low toxicity to an American strain. To our knowledge, this is the first report to demonstrate that the purified protease of P. luminescens has direct toxicity to P. xylostella and biopesticide potentiality.     

Determination of azinphos-methyl and parathion-methyl in honey by stripping voltammetry

An electrochemical method for the determination of azinphos-methyl and parathion-methyl in honey is presented. The determination is established by adsorptive stripping differential pulse voltammetry at hanging mercury working electrode.In contrast to the chromatographic methods for the determination of pesticide residues, the sample preparation of the proposed method is minimal; analytes were extracted from honey samples with a mixture of (acetone):(Britton-Robinson buffer) and then were analyzed without any additional pretreatment.The response of the analytes either individually or as a mixture was studied for a series of deposition time and molar ratio. Two quantitation protocols were compared, using either the external calibration or the standard addition method. Accuracy was tested with spiked honey samples obtaining good recovery values. The limit of detection for the honey sample (for deposition time of 10 s) was calculated 51.71 μg kg⁻¹ for MeP and 65.87 μg kg⁻¹ for AzMet.     

Dual Bioactivities of Essential Oil Extracted from the Leaves of Artemisia argyi as an Antimelanogenic versus Antioxidant Agent and Chemical Composition Analysis by GC/MS

The study was aimed at investigating the antimelanogenic and antioxidant properties of essential oil when extracted from the leaves of Artemisia argyi, then analyzing the chemical composition of the essential oil. The inhibitory effect of the essential oil on melanogenesis was evaluated by a mushroom tyrosinase activity assay and B16F10 melanoma cell model. The antioxidant capacity of the essential oil was assayed by spectrophotometric analysis, and the volatile chemical composition of the essential oil was analyzed with gas chromatography-mass spectrometry (GC/MS). The results revealed that the essential oil significantly inhibits mushroom tyrosinase activity (IC₅₀ = 19.16 mg/mL), down-regulates B16F10 intracellular tyrosinase activity and decreases the amount of melanin content in a dose-dependent pattern. Furthermore, the essential oil significantly scavenged 2,2-diphenyl-1-picryl-hydrazyl (DPPH) and 2,2′-azino-bis (3-ethylbenzthiazoline- 6-sulphonic acid) ABTS radicals, showed an apparent reduction power as compared with metal-ion chelating activities. The chemicals constituents in the essential oil are ether (23.66%), alcohols (16.72%), sesquiterpenes (15.21%), esters (11.78%), monoterpenes (11.63%), ketones (6.09%), aromatic compounds (5.01%), and account for a 90.10% analysis of its chemical composition. It is predicted that eucalyptol and the other constituents, except for alcohols, in the essential oil may contribute to its antioxidant activities. The results indicated that essential oil extracted from A. argyi leaves decreased melanin production in B16F10 cells and showed potent antioxidant activity. The essential oil can thereby be applied as an inhibitor of melanogenesis and could also act as a natural antioxidant in skin care products.     

P5C as an Interface of Proline Interconvertible Amino Acids and Its Role in Regulation of Cell Survival and Apoptosis

Studies of cancer metabolism have focused on the production of energy and the interconversion of carbons between cell cycles. More recently, amino acid metabolism, especially non-essential amino acids (NEAAs), has been investigated, underlining their regulatory role. One of the important mediators in energy production and interconversion of carbons in the cell is Δ¹-pyrroline-5-carboxylate (P5C)—the physiological intracellular intermediate of the interconversion of proline, ornithine, and glutamate. As a central component of these conversions, it links the tricarboxylic acid cycle (TCA), urea cycle (UC), and proline cycle (PC). P5C has a cyclic structure containing a tertiary nitrogen atom (N) and is in tautomeric equilibrium with the open-chain form of L-glutamate-γ-semialdehyde (GSAL). P5C is produced by P5C synthase (P5CS) from glutamate, and ornithine via ornithine δ-amino acid transferase (δOAT). It can also be converted to glutamate by P5C dehydrogenase (P5CDH). P5C is both a direct precursor of proline and a product of its degradation. The conversion of P5C to proline is catalyzed by P5C reductase (PYCR), while proline to P5C by proline dehydrogenase/oxidase (PRODH/POX). P5C-proline-P5C interconversion forms a functional redox couple. Their transformations are accompanied by the transfer of a reducing-oxidizing potential, that affect the NADP+/NADPH ratio and a wide variety of processes, e.g., the synthesis of phosphoribosyl pyrophosphate (PRPP), and purine ribonucleotides, which are crucial for DNA synthesis. This review focuses on the metabolism of P5C in the cell as an interconversion mediator of proline, glutamate, and ornithine and its role in the regulation of survival and death with particular emphasis on the metabolic context.