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.

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.

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.

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.

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.

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.

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.

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.

Ensuring Mixing Efficiency of Replica-Exchange Molecular Dynamics Simulations.

We address the question of constructing a protocol for replica-exchange molecular dynamics (REMD) simulations that make efficient use of the replica space, assess whether published applications are achieving such "mixing" efficiency, and provide a how-to guide to assist users to plan efficient REMD simulations. To address our first question, we introduce and discuss three metrics for assessing the number of replica-exchange attempts required to justify the use of a replica scheme and define a "transit number" as the lower bound for the length of an efficient simulation. Our literature survey of applications of REMD simulations of peptides in explicit solvent indicated that authors are not routinely reporting sufficient details of their simulation protocols to allow readers to make independent assessments of the impact of the method on their results, particularly whether mixing efficiency has been achieved. Necessary details include the expected or observed replica-exchange probability, together with the total number of exchange attempts, the exchange period, and estimates of the autocorrelation time of the potential energy. Our analysis of cases where the necessary information was reported suggests that in many of these simulations there are insufficient exchanges attempted or an insufficiently long period between them to provide confidence that the simulation length justifies the size of the replica scheme. We suggest guidelines for designing REMD simulation protocols to ensure mixing efficiency. Two key recommendations are that the exchange period should in general be larger than 1 ps and the number of exchange attempts should be chosen to significantly exceed the transit number for the replica scheme.

Identification of the RGG box motif in Shadoo: RNA-binding and signaling roles?

Using comparative genomics and in-silico analyses, we previously identified a new member of the prion-protein (PrP) family, the gene SPRN, encoding the protein Shadoo (Sho), and suggested its functions might overlap with those of PrP. Extended bioinformatics and conceptual biology studies to elucidate Sho's functions now reveal Sho has a conserved RGG-box motif, a well-known RNA-binding motif characterized in proteins such as FragileX Mental Retardation Protein. We report a systematic comparative analysis of RGG-box containing proteins which highlights the motif's functional versatility and supports the suggestion that Sho plays a dual role in cell signaling and RNA binding in brain. These findings provide a further link to PrP, which has well-characterized RNA-binding properties.

Calculation of a Complete Enzymic Reaction Surface: Reaction and Activation Free Energies for Hydride-Ion Transfer in Dihydrofolate Reductase.

We present a two-dimensional grid method for the calculation of complete free-energy surfaces for enzyme reactions using a hybrid quantum mechanical/molecular mechanical (QM/MM) potential within the semiempirical (PM3) QM approximation. This implementation is novel in that parallel processing with multiple trajectories (replica-exchange molecular dynamics simulations) is used to sample configuration space. The free energies at each grid point are computed using the thermodynamic integration formalism. From the free-energy surface, the minimum free-energy pathway for the reaction can be defined, and the computed activation and reaction energies can be compared with experimental values. We illustrate its use in a study of the hydride-transfer step in the reduction of dihydrofolate to tetrahydrofolate catalyzed by Escherichia coli dihydrofolate reductase with bound nicotinamide adenine dinucleotide phosphate cofactor. We investigated the effects of changing the QM region, ionization state of the conserved active-site Asp27 residue, and polarization contributions to the activation and reaction free energy. The results clearly show the necessity for including the complete substrate and cofactor molecules in the QM region, and the importance of the overall protein (MM) electrostatic environment in determining the free energy of the transition state (TS) and products relative to reactants. For the model with ionized Asp27, its inclusion in the QM region is essential. We found that the reported [Garcia-Viloca, M.; Truhlar, D. G.; Gao, J. J. Mol. Biol. 2003, 327, 549] stabilization of the TS due to polarization is an artifact that can be attributed to truncation of the electrostatic interactions between the QM and MM atoms. For neutral (protonated) Asp27, our calculated reaction free energy of -4 ± 2 kcal/mol agrees well with the experimental value of -4.4 kcal/mol, although the corresponding activation free-energy estimate is still high at 21 ± 2 kcal/mol compared with the experimental value of 13.4 kcal/mol. The results are less supportive of the ionized Asp27 model, which gives rise to a much higher activation barrier and favors the reverse reaction.

Combining docking and molecular dynamic simulations in drug design.

A rational approach is needed to maximize the chances of finding new drugs, and to exploit the opportunities of potential new drug targets emerging from genomic and proteomic initiatives, and from the large libraries of small compounds now readily available through combinatorial chemistry. Despite a shaky early history, computer-aided drug design techniques can now be effective in reducing costs and speeding up drug discovery. This happy outcome results from development of more accurate and reliable algorithms, use of more thoughtfully planned strategies to apply them, and greatly increased computer power to allow studies with the necessary reliability to be performed. Our review focuses on applications and protocols, with the main emphasis on critical analysis of recent studies where docking calculations and molecular dynamics (MD) simulations were combined to dock small molecules into protein receptors. We highlight successes to demonstrate what is possible now, but also point out drawbacks and future directions. The review is structured to lead the reader from the simpler to more compute-intensive methods. Thus, while inexpensive and fast docking algorithms can be used to scan large compound libraries and reduce their size, more accurate but expensive MD simulations can be applied when a few selected ligand candidates remain. MD simulations can be used: during the preparation of the protein receptor before docking, to optimize its structure and account for protein flexibility; for the refinement of docked complexes, to include solvent effects and account for induced fit; to calculate binding free energies, to provide an accurate ranking of the potential ligands; and in the latest developments, during the docking process itself to find the binding site and correctly dock the ligand a priori.

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 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.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase.

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

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.

Computational methods for the study of enzymic reaction mechanisms III: a perturbation plus QM/MM approach for calculating relative free energies of protonation.

We describe a coupling parameter, that is, perturbation, approach to effectively create and annihilate atoms in the quantum mechanical Hamiltonian within the closed shell restricted Hartree-Fock formalism. This perturbed quantum mechanical atom (PQA) method is combined with molecular mechanics (MM) methods (PQA/MM) within a molecular dynamics simulation, to model the protein environment (MM region) effects that also make a contribution to the overall free energy change. Using the semiempirical PM3 method to model the QM region, the application of this PQA/MM method is illustrated by calculation of the relative protonation free energy of the conserved OD2 (Asp27) and the N5 (dihydrofolate) proton acceptor sites in the active site of Escherichia coli dihydrofolate reductase (DHFR) with the bound nicotinamide adenine dinucleotide phosphate (NADPH) cofactor. For a number of choices for the QM region, the relative protonation free energy was calculated as the sum of contributions from the QM region and the interaction between the QM and MM regions via the thermodynamic integration (TI) method. The results demonstrate the importance of including the whole substrate molecule in the QM region, and the overall protein (MM) environment in determining the relative stabilities of protonation sites in the enzyme active site. The PQA/MM free energies obtained by TI were also compared with those estimated by a less computationally demanding nonperturbative method based on the linear response approximation (LRA). For some choices of QM region, the total free energies calculated using the LRA method were in very close agreement with the PQA/MM values. However, the QM and QM/MM component free energies were found to differ significantly between the two methods.