Revised mechanism of carboxylation of ribulose-1,5-biphosphate by rubisco from large scale quantum chemical calculations.

Here, we describe a computational approach for studying enzymes that catalyze complex multi-step reactions and apply it to Ribulose 1,5-bisphosphate carboxylase-oxygenase (Rubisco), the enzyme that fixes atmospheric carbon dioxide within photosynthesis. In the 5-step carboxylase reaction, the substrate Ribulose-1,5-bisphosphate (RuBP) first binds Rubisco and undergoes enolization before binding the second substrate, CO2 . Hydration of the RuBP.CO2 complex is followed by CC bond scission and stereospecific protonation. However, details of the roles and protonation states of active-site residues, and sources of protons and water, remain highly speculative. Large-scale computations on active-site models provide a means to better understand this complex chemical mechanism. The computational protocol comprises a combination of hybrid semi-empirical quantum mechanics and molecular mechanics within constrained molecular dynamics simulations, together with constrained gradient minimization calculations using density functional theory. Alternative pathways for hydration of the RuBP.CO2 complex and associated active-site protonation networks and proton and water sources were investigated. The main findings from analysis of the resulting energetics advocate major revision to existing mechanisms such that: hydration takes place anti to the CO2 ; both hydration and CC bond scission require early protonation of CO2 in the RuBP.CO2 complex; CC bond scission and stereospecific protonation reactions are concerted and, effectively, there is only one stable intermediate, the C3-gemdiolate complex. Our main conclusions for interpreting enzyme kinetic results are that the gemdiolate may represent the elusive Michaelis-Menten-like complex corresponding to the empirical Km (=Kc ) with turnover to product via bond scission concerted with stereospecific protonation consistent with the observed catalytic rate. © 2018 Wiley Periodicals, Inc.

Optimization of parameters for molecular dynamics simulation using smooth particle-mesh Ewald in GROMACS 4.5.

Based on our critique of requirements for performing an efficient molecular dynamics simulation with the particle-mesh Ewald (PME) implementation in GROMACS 4.5, we present a computational tool to enable the discovery of parameters that produce a given accuracy in the PME approximation of the full electrostatics. Calculations on two parallel computers with different processor and communication structures showed that a given accuracy can be attained over a range of parameter space, and that the attributes of the hardware and simulation system control which parameter sets are optimal. This information can be used to find the fastest available PME parameter sets that achieve a given accuracy. We hope that this tool will stimulate future work to assess the impact of the quality of the PME approximation on simulation outcomes, particularly with regard to the trade-off between cost and scientific reliability in biomolecular applications.

Kohn-Sham Density Functional Calculations Reveal Proton Wires in the Enolization and Carboxylase Reactions Catalyzed by Rubisco.

Ribulose 1,5-bisphosphate (RuBP) carboxylase-oxygenase (Rubisco) plays a fundamental role in the carbon cycle by fixing the atmospheric CO2 used in photosynthesis. Rubisco is all the more remarkable because it must catalyze some difficult multistep reaction chemistry involving proton transfers within the one active site. In the present study, we have used Kohn-Sham density functional theory at the B3LYP/6-31G* level with basis set superposition error and dispersion corrections (B3LYP-gCP-D3) to examine the possibility that the proton transfers can take place through molecular wires (including active-site water molecules) via the classical Grotthuss proton-shuttle mechanism. The results support an essential role for water molecules found in the crystal structures of Rubisco complexes as facilitators of proton transport in all the rate-limiting (catalytic) reaction steps through a network of short proton wires within the Rubisco active site. We suggest that completion of the initial product turnover (cycle) requires two excess protons produced in the initial carbamylation that is required for Rubisco activation. By use of proton wires, a large number of reaction steps may be accommodated within a single active site without necessitating the input of excessive conformational strain energy arising from the movement of residue side chains into positions where direct protonation of substrates can occur. The involvement of the identified types of proton wires in the kinetic mechanism is capable of providing a unique explanation for various experimental observations, including deuterium isotope effects and the results of site-directed mutagenesis experiments, and may thus provide a realistic solution to the problem of Rubisco's challenging chemistry.

Putative functional role for the invariant aspartate 263 residue of Rhodospirillum rubrum Rubisco.

Although aspartate residue D263 of Rhodospirillum rubrum Rubisco is close to the active site and invariant in all reported Rubiscos, its possible functional and structural roles in Rubisco activity have not been investigated. We have mutagenised D263 to several selected amino acids (asparagine, alanine, serine, glutamate, and glutamine) to probe possible roles in facilitating proton movements within the active site and maintaining structural positioning of key active-site groups. The mutants have been characterized by kinetic methods and by differential scanning calorimetry (DSC) to examine the effects of the substitutions on the stability of the folded state. We show that D263 is essential for maintaining effective levels of catalysis with the mutations reducing carboxylation variously by up to 100-fold but having less than 10% effect on the carboxylase/oxygenase specificity of the catalytic reaction. Removing the charge of the residue 263 side chain significantly strengthens binding of the activating (carbamylating) CO(2) molecule. In contrast, a charge on the 263 site has only a small influence on binding of the positively charged Mg(2+) ion, suggesting that the local protein structure provides different shielding of the formal charges on the Mg(2+) ion and the epsilon-lysine group of K191. Interestingly, introduction of an internal cavity (D263S and D263A) and insertion of an extra -CH(2)- group (D263E and D263Q) have opposite effects on catalysis, the former relatively small and the latter much larger, suggesting that the extra side-chain group induces a specific structural distortion that inhibits formation of the transition state. As the DSC results show that the mutations only slightly increase the kinetic stability of the folded state, we conclude that the rate-limiting (activated) step of unfolding involves substantial unfolding of the structure but not in the region of site 263. In summary, interaction of D263 with H287 of a largely electrostatic nature appears critical for maintaining correct positioning of catalytic groups in the active site. The conservation of D263 can thus be accounted for by its contribution to the maintenance of a finely tuned structure in this region abutting the active site.

Mechanism of Oxygenase-Pathway Reactions Catalyzed by Rubisco from Large-Scale Kohn-Sham Density Functional Calculations.

Ribulose 1,5-bisphosphate carboxylase-oxygenase (Rubisco) is the primary carbon-fixing enzyme in photosynthesis, fixing CO2 to a 5-carbon sugar, RuBP, in a series of five reactions. However, it also catalyzes an oxygenase reaction by O2 addition to the same enolized RuBP substrate in an analogous reaction series in the same active site, producing a waste product and loss of photosynthetic efficiency. Starting from RuBP, the reactions are enolization to the enediolate form, addition of CO2 or O2 to form the carboxy or peroxo adduct, hydration to form a gemdiolate, scission of the C2-C3 bond of the original RuBP, and stereospecific or nonstereospecific protonation to form two molecules of the 3-carbon PGA product, or one molecule of PGA, one of 2-carbon PG (waste product), and one water molecule. Reducing the loss of efficiency from the oxygenase reaction is an attractive means to increase crop productivity. However, lack of understanding of key aspects of the catalytic mechanisms for both the carboxylase and oxygenase reactions, particularly those involving proton exchanges and roles of water molecules, has stymied efforts at re-engineering Rubisco to reduce losses from the oxygenation reaction. As the stable form of molecular oxygen is the triplet biradical state (3O2), its reaction with near-universal singlet-state molecules is formally spin forbidden. Although in oxygenase enzymes, 3O2 activation is usually achieved by one-electron transfers using transition-metal ions or organic cofactors, recently, cofactor-less oxygenases in which the substrate itself is the source of the electron for 3O2 activation have been identified, but in all such cases an aromatic ring stabilizes the substrate's negative charge. Here we present the first large-scale Kohn-Sham density functional theory study of the reaction mechanism of the Rubisco oxygenase pathway. First, we show that the enediolate substrate complexed to Mg2+ and its ligands extends the region for charge delocalization and stabilization of its negative charge to allow formation of a caged biradical enediolate-O2 complex. Thus, Rubisco is a unique type of oxygenase without precedent in the literature. Second, for the O2 addition to proceed to the singlet peroxo-adduct intermediate, the system must undergo an intersystem crossing. We found that the presence of protonated LYS334 is required to stabilize this intermediate and that both factors (strongly stabilized anion and protonated LYS334) facilitate a barrier-less activation of 3O2. This finding supports our recent proposal that deoxygenation, that is, reversal of gas binding, is possible. Third, as neither CO2 nor O2 binds to the enzyme, our findings support the proposal from our recent carboxylase study that the observed KC or KO (Michaelis-Menten constants) in the steady-state kinetics reflect the respective adducts, carboxy or peroxo. Fourth, after computing hydration pathways with water addition both syn and anti to C3, we found, in contrast to the results of our carboxylation study indicating anti addition, that in the oxygenation reaction only syn-hydration is capable of producing a stable gemdiolate that facilitates the rate-limiting C2-C3 bond scission to final products. Fifth, we propose that an excess proton we previously found was required in the carboxylation reaction for activating the C2-C3 bond scission is utilized in the oxygenation reaction for the required elimination of a water molecule. In summary, despite its oxygenase handicap, Rubisco's success in directing 75% of its substrate through the carboxylation pathway can be considered impressively effective. Although native C3 Rubiscos are in a fix with unwanted activity of 3O2 hampering its primary carboxylase function, mechanistic differences presented here with findings in our recent carboxylase study for both the gas-addition and subsequent reactions provide some clues as to how creative Rubisco re-engineering may offer a solution to reducing the oxygenase activity.

Ab Initio Molecular Dynamics Simulation and Energetics of the Ribulose-1,5-biphosphate Carboxylation Reaction Catalyzed by Rubisco: Toward Elucidating the Stereospecific Protonation Mechanism.

In the carboxylation reaction catalyzed by ribulose 1,5-bisphosphate (RuBP) carboxylase-oxygenase (Rubisco), which is fundamental to photosynthesis, scission of a C-C bond in the six-carbon gemdiolate intermediate forms a carbanion that must be protonated stereospecifically to form product. It is thought that a conserved lysine side chain (LYS175 in spinach Rubisco), in the immediate vicinity of the carbanion, provides the necessary proton. Here, we endeavor to determine from the electronic-structure calculations whether protonation via this route is energetically possible. The two-dimensional energy surface was mapped to determine the minimum energy path (MEP) using density functional theory (B3LYP) and incorporating basis set superposition and classical (London) dispersion corrections. The potential of mean force (free energy) was then calculated from ab initio molecular dynamics simulations with umbrella sampling in the vicinity of the MEP on the scission-protonation reaction coordinate. MEP calculations were also carried out to evaluate the possibility of an active-site water near the phosphate (P1) of RuBP, with an excess proton positioned at P1, as an alternative facilitator of stereospecific protonation via a classical Grotthuss mechanism. In both cases, the C-C bond scission in the six-carbon intermediate and proton transfer from the donor was found to be concerted and highly asynchronous, without a stable carbanion intermediate. However, the free energy change was unfavorable for direct protonation by the LYS175 side chain. In contrast, the Grotthuss mechanism yielded stable products and an activation energy in good agreement with experiment. It also provides a plausible mechanism for alternative product formed in enzyme mutations at the LYS175 position and is consistent with the observed deuterium isotope effects.

Redefinition of rubisco carboxylase reaction reveals origin of water for hydration and new roles for active-site residues.

Crystallographic, mutagenesis, kinetic, and computational studies on Rubisco over three decades have revealed much about its catalytic mechanism and the role played by several active-site residues. However, key questions remain unanswered. Specific details of the carboxylase and oxygenase mechanisms, required to underpin the rational re-engineering of Rubisco, are still speculative. Here we address critical gaps in knowledge with a definitive comprehensive computational investigation of the mechanism of carboxylase activity at the Rubisco active site. Density functional theory calculations (B3LYP/6-31G(d,p)) were performed on active-site fragment models of a size up to 77 atoms, not previously possible computationally. All amino acid residues suspected to play roles in the acid-base chemistry in the multistep reaction, and interacting directly with the central Mg (2+) atom and the reactive moiety of substrate and intermediates, were included. The results provide a firm basis for us to propose a novel mechanism for the entire sequence of reactions in the carboxylase catalysis and to define precise roles for the active-site residues, singly and in concert. In this mechanism, the carbamylated LYS201 plays a more limited role than previously proposed but is crucial for initiating the reaction by acting as a base in the enolization. We suggest a wider role for HIS294, with involvement in the carboxylation, hydration, and C2-C3 bond-scission steps, consistent with the suggestion of Harpel et al. (1998) but contrary to the consensus view of Cleland et al. (1998). In contrast to the common assumption that the water molecule for the hydration step comes from within the active site, we propose that the Mg-coordinated water is not dissociated at the start of the gas-addition reaction but rather remains coordinated and is used for the hydration of the C3 carbon atom. New roles are also proposed for LYS175, GLU204, and HIS294. The mechanism suggests roles in the gas-addition step for residues in three spatially distinct regions of the active site, HIS294 and LYS334 in the C-terminal domain of the large subunit (LSU), but also hitherto unsuspected roles for a cluster of three residues (ASN123, GLU60, and TYR20) in the N-terminal domain of the partner LSU of the dimer containing the active site. Our new mechanism is supported by existing experimental data, provides new convincing interpretations of previously puzzling data, and allows new insights into mutational strategies for improving Rubisco activity.

Directions for Optimization of Photosynthetic Carbon Fixation: RuBisCO's Efficiency May Not Be So Constrained After All.

The ubiquitous enzyme Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) fixes atmospheric carbon dioxide within the Calvin-Benson cycle that is utilized by most photosynthetic organisms. Despite this central role, RuBisCO's efficiency surprisingly struggles, with both a very slow turnover rate to products and also impaired substrate specificity, features that have long been an enigma as it would be assumed that its efficiency was under strong evolutionary pressure. RuBisCO's substrate specificity is compromised as it catalyzes a side-fixation reaction with atmospheric oxygen; empirical kinetic results show a trend to tradeoff between relative specificity and low catalytic turnover rate. Although the dominant hypothesis has been that the active-site chemistry constrains the enzyme's evolution, a more recent study on RuBisCO stability and adaptability has implicated competing selection pressures. Elucidating these constraints is crucial for directing future research on improving photosynthesis, as the current literature casts doubt on the potential effectiveness of site-directed mutagenesis to improve RuBisCO's efficiency. Here we use regression analysis to quantify the relationships between kinetic parameters obtained from empirical data sets spanning a wide evolutionary range of RuBisCOs. Most significantly we found that the rate constant for dissociation of CO2 from the enzyme complex was much higher than previous estimates and comparable with the corresponding catalytic rate constant. Observed trends between relative specificity and turnover rate can be expressed as the product of negative and positive correlation factors. This provides an explanation in simple kinetic terms of both the natural variation of relative specificity as well as that obtained by reported site-directed mutagenesis results. We demonstrate that the kinetic behaviour shows a lesser rather than more constrained RuBisCO, consistent with growing empirical evidence of higher variability in relative specificity. In summary our analysis supports an explanation for the origin of the tradeoff between specificity and turnover as due to competition between protein stability and activity, rather than constraints between rate constants imposed by the underlying chemistry. Our analysis suggests that simultaneous improvement in both specificity and turnover rate of RuBisCO is possible.

Integron-sequestered dihydrofolate reductase: a recently redeployed enzyme.

The introduction and wide use of antibacterial drugs has resulted in the emergence of resistant organisms. DfrB dihydrofolate reductase (DHFR) is a bacterial enzyme that is uniquely associated with mobile gene cassettes within integrons, and confers resistance to the drug trimethoprim. This enzyme has intrigued microbiologists since it was discovered more than thirty years ago because of its simple structure, enzymatic inefficiency and its virtual insensitivity to trimethoprim. Here, for the first time, a comprehensive discussion of genetic, evolutionary, structural and functional studies of this enzyme is presented together. This information supports the ideas that DfrB DHFR is a poorly adapted catalyst and has recently been recruited to perform a novel enzymatic activity in response to selective pressure.

The prion protein gene: identifying regulatory signals using marsupial sequence.

The function of the prion protein gene (PRNP) and its normal product PrP(C) is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyT1, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrP(C) in signal transduction and synaptic plasticity.

The C-type lectin-like domain superfamily.

The superfamily of proteins containing C-type lectin-like domains (CTLDs) is a large group of extracellular Metazoan proteins with diverse functions. The CTLD structure has a characteristic double-loop ('loop-in-a-loop') stabilized by two highly conserved disulfide bridges located at the bases of the loops, as well as a set of conserved hydrophobic and polar interactions. The second loop, called the long loop region, is structurally and evolutionarily flexible, and is involved in Ca2+-dependent carbohydrate binding and interaction with other ligands. This loop is completely absent in a subset of CTLDs, which we refer to as compact CTLDs; these include the Link/PTR domain and bacterial CTLDs. CTLD-containing proteins (CTLDcps) were originally classified into seven groups based on their overall domain structure. Analyses of the superfamily representation in several completely sequenced genomes have added 10 new groups to the classification, and shown that it is applicable only to vertebrate CTLDcps; despite the abundance of CTLDcps in the invertebrate genomes studied, the domain architectures of these proteins do not match those of the vertebrate groups. Ca2+-dependent carbohydrate binding is the most common CTLD function in vertebrates, and apparently the ancestral one, as suggested by the many humoral defense CTLDcps characterized in insects and other invertebrates. However, many CTLDs have evolved to specifically recognize protein, lipid and inorganic ligands, including the vertebrate clade-specific snake venoms, and fish antifreeze and bird egg-shell proteins. Recent studies highlight the functional versatility of this protein superfamily and the CTLD scaffold, and suggest further interesting discoveries have yet to be made.

Comparative analysis of structural properties of the C-type-lectin-like domain (CTLD).

The superfamily of proteins containing the C-type-lectin-like domain (CTLD) is a group of abundant extracellular metazoan proteins characterized by evolutionary flexibility and functional versatility. Several CTLDs are also found in parasitic prokaryotes and viruses. The 37 distinct currently available CTLD structures demonstrate significant structural conservation despite low or undetectable sequence similarity. Our aim in this study was to perform an extensive comparative analysis of all available CTLD structures to establish the most conserved structural features of the fold, and to test and extend the early analysis of Drickamer. By implication, these features should be those critical for maintenance of integrity of the fold. By analyzing CTLD structures superimposed by several methods, we have established groups of conserved structural positions involved in fold maintenance but not in ligand binding; these are consistent with the fold's known functional flexibility. In addition to the well-recognized disulfide bridges, groups of conserved residues are involved in hydrophobic interactions stabilizing the core of the fold and the long loop region, and in an alpha2-beta1-beta5 polar interaction. Evaluation of the conclusions of the structure comparison study compared with alignments of all available human, mouse and Caenorhabditis elegans CTLD sequences showed that conservation patterns are preserved throughout the whole CTLD sequence space. Our observations provide an improved understanding of CTLD structure, and will help in identification of new CTLDs and the mechanisms that drive and constrain the coevolution of the structure and function of the fold.

C-type lectin-like domains in Fugu rubripes.

BACKGROUND: Members of the C-type lectin domain (CTLD) superfamily are metazoan proteins functionally important in glycoprotein metabolism, mechanisms of multicellular integration and immunity. Three genome-level studies on human, C. elegans and D. melanogaster reported previously demonstrated almost complete divergence among invertebrate and mammalian families of CTLD-containing proteins (CTLDcps). RESULTS: We have performed an analysis of CTLD family composition in Fugu rubripes using the draft genome sequence. The results show that all but two groups of CTLDcps identified in mammals are also found in fish, and that most of the groups have the same members as in mammals. We failed to detect representatives for CTLD groups V (NK cell receptors) and VII (lithostathine), while the DC-SIGN subgroup of group II is overrepresented in Fugu. Several new CTLD-containing genes, highly conserved between Fugu and human, were discovered using the Fugu genome sequence as a reference, including a CSPG family member and an SCP-domain-containing soluble protein. A distinct group of soluble dual-CTLD proteins has been identified, which may be the first reported CTLDcp group shared by invertebrates and vertebrates. We show that CTLDcp-encoding genes are selectively duplicated in Fugu, in a manner that suggests an ancient large-scale duplication event. We have verified 32 gene structures and predicted 63 new ones, and make our annotations available through a distributed annotation system (DAS) server http://anz.anu.edu.au:8080/Fugu_rubripes/ and their sequences as additional files with this paper. CONCLUSIONS: The vertebrate CTLDcp family was essentially formed early in vertebrate evolution and is completely different from the invertebrate families. Comparison of fish and mammalian genomes revealed three groups of CTLDcps and several new members of the known groups, which are highly conserved between fish and mammals, but were not identified in the study using only mammalian genomes. Despite limitations of the draft sequence, the Fugu rubripes genome is a powerful instrument for gene discovery and vertebrate evolutionary analysis. The composition of the CTLDcp superfamily in fish and mammals suggests that large-scale duplication events played an important role in the evolution of vertebrates.

Shadoo, a new protein highly conserved from fish to mammals and with similarity to prion protein.

We report evidence from cDNA isolation and expression analysis as well as analyses of genome, expressed sequence tag (EST), cDNA and expression databases for a new gene named SPRN (shadow of prion protein). SPRN comprises two exons, with the open reading frame (ORF) contained within exon 2, and codes for a protein of 130-150 amino acids named Shadoo (Japanese shadow), predicted to be extracellular and GPI-anchored. The SPRN gene was found in fish (zebrafish, Fugu) and mammals (mouse, rat, human). Conservation of order and transcription orientation of two proximal genes between fishes and mammals strongly indicates gene orthology. Sequence comparison shows: a highly conserved N-terminal signal sequence; Arg-rich basic region containing up to six tetrarepeats of consensus XXRG (where X is G, A or S); a hydrophobic region of 20 residues with strong homology to PrP; a less conserved C-terminal domain containing a conserved glycosylation motif; and a C-terminal peptide predicted to be a signal sequence for glycophosphotidylinositol (GPI)-anchor attachment. Fish Shadoos (Sho) show well conserved sequences (identity 54%) over 106 amino acids (zebrafish length), and conservation among the mammalian sequences is very high (identity 81-96%). The fish and mammalian sequences are also well conserved, particularly for zebrafish, to beyond the end of the hydrophobic sequence (identity 41-53%, 78 amino acids, zebrafish length). The overall structure appears closely related to prion proteins (PrPs), although the C-terminal domains of Shos are quite different from those of PrPs, for which conformational changes in mammals are implicated in disease. The structural similarity is particularly interesting given recent reports of three new genes with similarities to PrPs found in Fugu (PrP-like, PrP-461/stPrP-1 and stPrP-2) and other fish, but for which direct evolution to higher vertebrate PrPs is unlikely and for which no other mammalian homologues have been found. Database information indicates expression of SPRN in embryo, brain and retina of mouse and rat, hippocampus of human, and in embryo and retina of zebrafish, and we directly confirmed a strikingly specific expression of the mammalian (human, mouse, rat) transcripts in whole brain. This result together with some common structural features led to the suggestive hypothesis of a possible functional link between mammalian PrP and Sho proteins.

Copper-induced conformational change in a marsupial prion protein repeat peptide probed using FTIR spectroscopy.

We report the first Fourier transform infrared analysis of prion protein (PrP) repeats and the first study of PrP repeats of marsupial origin. Large changes in the secondary structure and an increase in hydrogen bonding within the peptide groups were evident from a red shift of the amide I band by >7 cm(-1) and an approximately five-fold reduction in amide hydrogen-deuterium exchange for peptide interacting with Cu(2+) ions. Changes in the tertiary structure upon copper binding were also evident from the appearance of a new band at 1564 cm(-1), which arises from the ring vibration of histidine. The copper-induced conformational change is pH dependent, and occurs at pH >7.

An efficient Z-score algorithm for assessing sequence alignments.

We describe an alternative method for scoring of the pairwise alignment of two biological sequences. Designed to overcome the bias due to the composition of the alignment, it measures the distance (in standard deviations) between the given alignment and the mean value of all other alignments that can be obtained by a permutation of either sequence. We demonstrate that the standard deviation can be calculated efficiently. By concentrating upon the ungapped case, the mean and standard deviation can be calculated exactly and in two steps, the first being O(N) time, where N is the length of the sequence, the second in a fixed number of calculations, i.e., in O(1) time. We argue that this statistic is a more consistent measure than a similarity score based upon a standard scoring matrix. Even in the ungapped case, the statistic proves in many cases to be more accurate than the commonly used (FASTA) (Pearson and Lipman, 1988) gapped Z-score in which the sequence is matched against a random sample of the database. We demonstrate the use of the POZ-score as a secondary filter which screens out several well-known types of false positive, reducing the amount of manual screening to be done by the biologist.

Synthesis of quaternised 2-aminopyrimido[4,5-d]pyrimidin-4(3H)-ones and their biological activity with dihydrofolate reductase.

In a program to design and develop mechanism-based compounds active as substrates and inhibitors of dihydrofolate reductase (DHFR), we report the synthesis and physical properties of the 6-methyl- (7), 8-methyl- (8a), and 8-ethyl- (8b) derivatives of the parent 2-aminopyrimido[4,5-d]pyrimidin-4-(3H)-one (6). These compounds are the first members of a class of heterocycles related to 8-alkylpterins (N8-alkyl-2-aminopteridin-4(8H)-ones) (2a-2c), which have been shown to be novel substrates for DHFR. Three methods were developed for the synthesis of target compounds 7, 8a and 8b; however, the optimum yields (1-8%) could not be improved because the products decomposed by ring opening (e.g. to 2,4-diamino-5-methyliminomethylpyrimidin-6(1H)-one (9)) under the reaction conditions. The marked pi-electron deficiency of compounds 7, 8a and 8b is the likely cause for the susceptibility of the quaternised pyrimidine ring in the related cations 10, 15a and 15b, respectively, to add nucleophiles, thus promoting the opening of the pyrimidopyrimidine ring system. 1H-NMR spectroscopic studies of compounds 7, 8a and 8b revealed a fast and reversible covalent hydration of the associated cations across the C7z.sbnd;N8 bond for the N6-methyl derivative 7 and across the N6z.sbnd;C7 bond for the N8-methyl derivative 8a. UV spectroscopic studies of methyl derivatives 7 and 8a as well as the parent heterocycle 6 showed that protonation of the latter occurred at N1, while methylation with iodomethane proceeded at N6 and N8. The basicities of the N-methyl derivatives 7 and 8a (pK(a) ca. 5.5) are similar to those of 8-alkylpterins 2; this is an essential element of the design to promote binding to DHFR in their protonated form. Enzyme kinetics of 7, 8a and 8b with chicken DHFR confirmed our predictions that they are substrates, with apparent K(m) values of 3.8, 0.08, and 0.65 mM, and apparent V(max) values of 0.47, 2.27, and 0.30 nmol L(-1) min(-1) (for enzyme concentration 0.122 micro M), respectively. The parent compound 6 was not a substrate.

Methyltetrahydrofolate:corrinoid/iron-sulfur protein methyltransferase (MeTr): protonation state of the ligand and active-site residues.

Methyltetrahydrofolate:corrinoid/iron-sulfur protein methyltransferase (MeTr) catalyzes the transfer of the N5-methyl group from N5-methyltetrahydrofolate (CH(3)THF) to the cobalt center of a corrinoid/iron-sulfur protein, a reaction similar to that of cobalamin-dependent methionine synthase (MetH). For such a reaction to occur, CH(3)THF is expected to be activated by a stereospecific protonation at the N5 position. It has been shown experimentally that binding to MeTr is associated with a pK(a) increase and proton uptake. The enzyme could achieve this by binding the unprotonated form of CH(3)THF, followed by specific protonation at the correct orientation. Here we have used computational approaches to investigate the protonation state of the ligand and active-site residues in MeTr. First, quantum mechanical (QM) methods with the PCM solvation model were used to predict protonation positions and pK(a) values of pterin, folate, and their analogues in an aqueous environment. After a reliable calibration of computational and experimental results was obtained, the effect of the protein environment was then considered. Different protonation states of CH(3)THF and active-site aspartic residues (D75 and D160) were investigated using QM calculations of active-site fragment complexes and the perturbed quantum atom (PQA) approach within QM/MM simulations. The final free energy results indicate that the N5 position of the tetrahydropterin ring is the preferred protonation position of CH(3)THF when bound to the active site of MeTr, followed by Asp160. We also found that the active-site environment is likely to increase the pK(a) of N5 by about 3 units, leading to proton uptake upon CH(3)THF binding, as observed experimentally for MeTr. Some implications of the results are discussed for the MetH mechanism.

Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes Micromonas.

Picoeukaryotes are a taxonomically diverse group of organisms less than 2 micrometers in diameter. Photosynthetic marine picoeukaryotes in the genus Micromonas thrive in ecosystems ranging from tropical to polar and could serve as sentinel organisms for biogeochemical fluxes of modern oceans during climate change. These broadly distributed primary producers belong to an anciently diverged sister clade to land plants. Although Micromonas isolates have high 18S ribosomal RNA gene identity, we found that genomes from two isolates shared only 90% of their predicted genes. Their independent evolutionary paths were emphasized by distinct riboswitch arrangements as well as the discovery of intronic repeat elements in one isolate, and in metagenomic data, but not in other genomes. Divergence appears to have been facilitated by selection and acquisition processes that actively shape the repertoire of genes that are mutually exclusive between the two isolates differently than the core genes. Analyses of the Micromonas genomes offer valuable insights into ecological differentiation and the dynamic nature of early plant evolution.