Conformations of nicotinamide adenine dinucleotide (NAD(+)) in various environments

Conformational Landscape of NADH and Ion Binding in Water/DMSO Mixtures via 31P NMR Spectroscopy

This study reports on the conformational states of nicotinamide adenine dinucleotide (NADH) in water/DMSO mixtures and examines the influence of ion binding. We observe evidence of conformational changes through a series of NMR techniques, including 31P NMR relaxation (R1 and R2), chemical exchange saturation transfer (CEST), and diffusion-ordered spectroscopy (DOSY) experiments. The observed variation of the conformational states is indicative of the solvent's influence on NADH's structural flexibility. The experimental findings, in combination with viscosity data, are shown to be in line with findings from earlier molecular dynamics studies. The reported observations emphasize the critical role of the solvent environment in determining the conformational states of NADH and similar molecules with relevance for the biophysiological context. The results found herein can help in studying the biophysical behavior of NADH and similar biomolecules and their associated metabolic pathways under various solvent conditions.

The use of NADH anisotropy to investigate mitochondrial cristae alignment

Life may be expressed as the flow of electrons, protons, and other ions, resulting in large potential difference. It is also highly photo-sensitive, as a large proportion of the redox capable molecules it relies on are chromophoric. It is thus suggestive that a key organelle in eukaryotes, the mitochondrion, constantly adapt their morphology as part of the homeostatic process. Studying unstained in vivo nano-scale structure in live cells is technically very challenging. One option is to study a central electron carrier in metabolism, reduced nicotinamide adenine dinucleotide (NADH), which is fluorescent and mostly located within mitochondria. Using one and two-photon absorption (340-360 nm and 730 nm, respectively), fluorescence lifetime imaging and anisotropy spectroscopy of NADH in solution and in live cells, we show that mitochondria do indeed appear to be aligned and exhibit high anisotropy (asymmetric directionality). Aqueous solution of NADH showed an anisotropy of ~ 0.20 compared to fluorescein or coumarin of < 0.1 and 0.04 in water respectively and as expected for small organic molecules. The anisotropy of NADH also increased further to 0.30 in the presence of proteins and 0.42 in glycerol (restricted environment) following two-photon excitation, suggesting more ordered structures. Two-photon NADH fluorescence imaging of Michigan Cancer Foundation-7 (MCF7) also showed strong anisotropy of 0.25 to 0.45. NADH has a quantum yield of fluorescence of 2% compared to more than 40% for photoionisation (electron generation), when exposed to light at 360 nm and below. The consequence of such highly ordered and directional NADH patterns with respect to electron ejection upon ultra-violet (UV) excitation could be very informative-especially in relation to ascertaining the extent of quantum effects in biology, including electron and photonic cascade, communication and modulation of effects such as spin and tunnelling.

Ab Initio Studies of NMNH(2-) Conformers in Water-Methanol Solutions: Comparative Analysis of the Biexponential Fluorescence Signals for NMNH(2-) and NADH

Ab initio studies of the structure of reduced β-nicotinamide d-ribonucleotide (NMNH(2-)) conformations in water and methanol solutions have been carried out for clarifying the role of the phosphate groups in fluorescence parameters of the NMNH(2-) molecule and the reduced β-nicotinamide adenine dinucleotide (NADH) molecule. Relaxed potential energy surfaces as a function of the dihedral rotation angle of the amide group in the NMNH(2-) molecule were calculated in the ground electronic state and the first excited electronic state to better understand the effect of phosphate groups on the nonradiative decay rates in the nicotinamide chromophore groups. The differences in the weighting coefficients in the biexponential fluorescence signals for NMNH(2-) and NADH molecules in solution were explained. A strong hydrogen bonding between the amide hydrogen atom and the nearest oxygen O^- atom from the phosphate group was detected by ab initio calculations for the folded NMNH(2-) conformations in the ground electronic state at trans configurations of the nicotinamide ring. This hydrogen bonding turned out to be much weaker for the first excited electronic state. These calculated data show that, after optical excitation of the NMNH(2-) molecule, a rapid change in the geometry of the molecule is possible. The strong interaction of the phosphate group with the amide group in NMNH(2-) molecules in aqueous solution leads to the predominance of the folded NMNH(2-) conformations and trans configurations of the nicotinamide ring. This explains the reason for the dominance of one fluorescence decay time of NMNH(2-) in the aqueous solution. Based on these data, an important conclusion can be drawn that the contribution of the exponent with the short decay time τ ≈ 0.28 ns to the fluorescence signal of NMNH(2-), NADH, and NADPH molecules is related to the trans configuration of the nicotinamide ring.

Anisotropic relaxation in NADH excited states studied by polarization-modulation pump-probe transient spectroscopy

We present the results of experimental and theoretical studies of fast anisotropic relaxation and rotational diffusion in the first electron excited state of biological coenzyme NADH in water-ethanol solutions. The experiments have been carried out by means of a novel polarization-modulation transient method and fluorescence polarization spectroscopy. For interpretation of the experimental results a model of the anisotropic relaxation in terms of scalar and vector properties of transition dipole moments has been developed based on the Born-Oppenheimer approximation. This model allows for the description of fast isotropic and anisotropic excited state relaxation under excitation of molecules by ultrafast laser pulses in transient absorption and upconversion experiments. The results obtained suggest that the dynamics of anisotropic rovibronic relaxation in NADH under excitation with 100 fs pump laser pulses can be characterised by a single vibrational relaxation time τv lying in the range of 2-15 ps and a single rotation diffusion time τr lying in the range of 100-450 ps, both depending on ethanol concentration. The dependence of the times τv and τr on the solution polarity (static permittivity) and viscosity has been determined and analyzed. Limiting values of the term P2(cos θ) describing the rotation of the transition dipole moment in the course of vibrational relaxation have been determined from experiments as a function of ethanol concentration and analyzed.

Two-Photon Excited Fluorescence Dynamics in NADH in Water-Methanol Solutions: The Role of Conformation States

The dynamics of polarized fluorescence in reduced nicotinamide adenine dinucleotide (NADH) at 460 nm under two-photon excitation at 720 nm by femtosecond laser pulses in water-methanol solutions has been studied experimentally and theoretically as a function of methanol concentration. A number of fluorescence parameters have been determined from experiment by means of the global fit procedure and then compared with the results reported by other authors. A comprehensive analysis of experimental errors was made. Ab initio calculations of the structure of NADH in water and methanol and of β-nicotinamide mononucleotide (NMNH) in vacuum have been carried out for clarifying the role of decay time heterogeneity. The main results obtained are as follows. An explanation of the heterogeneity in the measured fluorescence decay times in NADH has been suggested based on the influence of the internal molecular electric field in the nicotinamide ring on nonradiative decay rates. We suggest that different charge distributions in the cis and trans configurations result in different internal electrostatic field distributions that lead to the decay time heterogeneity. A slight but noticeable rise of the fluorescence decay times τ₁ and τ₂ with methanol concentration was observed and treated as a minor effect of a nonradiative relaxation slowing due to the decrease in solution polarity. Relative concentrations of the folded and unfolded NADH conformations in solutions have been determined using a new method of analysis of the rotational diffusion time τ_r as a function of methanol concentration on the basis of the Stokes-Einstein-Debye equation. The analysis of the fluorescence anisotropy parameters obtained under linearly and circularly polarized excitation and the parameter Ω has been carried out and resulted in the determination of the two-photon excitation tensor components and suggested the existence of two excitation channels with comparable intensities. These were the longitudinal excitation channel dominated by the diagonal tensor component S_zz and the mixed excitation channel dominated by the off-diagonal tensor components |S_xz² + S_yz²|^1/2.

Ultrafast Fluorescence Signals from β-Dihydronicotinamide Adenine Dinucleotide: Resonant Energy Transfer in the Folded and Unfolded Forms

β-Dihydronicotinamide adenine dinucleotide (NADH) plays a critical role in biological redox processes. Inside the cell, NADH can be in a folded or an unfolded conformation, depending on the chemical environment that surrounds it. It is known that selective excitation of adenine in NADH can induce energy transfer events from this nucleotide to the reduced nicotinamide chromophore. From the anticipated time scales, this process must compete with adenine's internal conversion channel, which is known to occur in the sub-picosecond time scale. In this work, we studied the dynamics of the excited states of both chromophores in NADH through the time resolution of the spontaneous emission from both nucleotides. Through these experiments, we extend the knowledge about the competition between the different photophysical channels both in the folded and unfolded states. The study involved the folded and unfolded states of NADH by experiments in water and methanol solutions. Our femtosecond fluorescence results were complemented by the first nuclear magnetic resonance through space magnetization transfer measurements on NADH, which establish the solvent-dependent folded versus unfolded states. We determined the dynamics of the excited states by direct excitation of dihydronicotinamide at 380 nm and adenine at 266 nm. From this, we were able to measure for the folded state, a time constant of 90 fs for energy transfer. Additionally, we determined that even in what is referred to as an unfolded state in methanol, non-negligible excitation transfer events do take place. Spontaneous emission anisotropy measurements were used in order to confirm the presence of a minor energy transfer channel in the methanol solutions where the unfolded state dominates.

Reversible pH-induced fluorescence colour change of gold nanoclusters based on pH-regulated surface interactions

A water-dispersible, biocompatible, fluorescent ratiometric pH nanosensor over a broad pH range is reported.

Insights from ion mobility-mass spectrometry, infrared spectroscopy, and molecular dynamics simulations on nicotinamide adenine dinucleotide structural dynamics: NAD+vs. NADH

Nicotinamide adenine dinucleotide (NAD) is found in all living cells where the oxidized (NAD⁺) and reduced (NADH) forms play important roles in many enzymatic reactions. However, little is known about NAD⁺ and NADH conformational changes and kinetics as a function of the cell environment. In the present work, an analytical workflow is utilized to study NAD⁺ and NADH dynamics as a function of the organic content in solution using fluorescence lifetime spectroscopy and in the gas-phase using trapped ion mobility spectrometry coupled to mass spectrometry (TIMS-MS) and infrared multiple photon dissociation (IRMPD) spectroscopy. NAD solution time decay studies showed a two-component distribution, assigned to changes from a "close" to "open" conformation with the increase of the organic content. NAD gas-phase studies using nESI-TIMS-MS displayed two ion mobility bands for NAD⁺ protonated and sodiated species, while four and two ion mobility bands were observed for NADH protonated and sodiated species, respectively. Changes in the mobility profiles were observed for NADH as a function of the starting solution conditions and the time after desolvation, while NAD⁺ profiles showed no dependence. IRMPD spectroscopy of NAD⁺ and NADH protonated species in the 800-1800 and 3200-3700 cm^-1 spectral regions showed common and signature bands between the NAD forms. Candidate structures were proposed for NAD⁺ and NADH kinetically trapped intermediates of the protonated and sodiated species, based on their collision cross sections and IR profiles. Results showed that NAD⁺ and NADH species exist in open, stack, and closed conformations and that the driving force for conformational dynamics is hydrogen bonding of the N-H-O and O-H-O forms with ribose rings.

Mechanistic Analysis of Fluorescence Quenching of Reduced Nicotinamide Adenine Dinucleotide by Oxamate in Lactate Dehydrogenase Ternary Complexes

Fluorescence of Reduced Nicotinamide Adenine Dinucleotide (NADH) is extensively employed in studies of oxidoreductases. A substantial amount of static and kinetic work has focused on the binding of pyruvate or substrate mimic oxamate to the binary complex of lactate dehydrogenase (LDH)-NADH where substantial fluorescence quenching is typically observed. However, the quenching mechanism is not well understood limiting structural interpretation. Based on time-dependent density functional theory (TDDFT) computations with cam-B3LYP functional in conjunction with the analysis of previous experimental results, we propose that bound oxamate acts as an electron acceptor in the quenching of fluorescence of NADH in the ternary complex, where a charge transfer (CT) state characterized by excitation from the highest occupied molecular orbital (HOMO) of the nicotinamide moiety of NADH to the lowest unoccupied molecular orbital (LUMO) of oxamate exists close to the locally excited (LE) state involving only the nicotinamide moiety. Efficient quenching in the encounter complex like in pig heart LDH requires that oxamate forms a salt bridge with Arg-171 and hydrogen bonds with His-195, Thr-246 and Asn-140. Further structural rearrangement and loop closure, which also brings about another hydrogen bond between oxamate and Arg-109, will increase the rate of fluorescence quenching as well.

Towards understanding the unbound state of drug compounds: Implications for the intramolecular reorganization energy upon binding

There has been an explosion of structural information for pharmaceutical compounds bound to biological targets, but the conformations and dynamics of compounds free in solution are poorly characterized, if at all. Yet, knowledge of the unbound state is essential to understand the fundamentals of molecular recognition, including the much debated conformational intramolecular reorganization energy of a compound upon binding (ΔEReorg). Also, dependable observation of the unbound compounds is important for ligand-based drug discovery, e.g. with pharmacophore modelling. Here, these questions are addressed with long (⩾0.5μs) state-of-the-art molecular dynamics (MD) simulations of 26 compounds (including 7 approved drugs) unbound in explicit solvent. These compounds were selected to be chemically diverse, with a range of flexibility, and good quality bioactive X-ray structures. The MD-simulated free compounds are compared to their bioactive structure and conformers generated with ad hoc sampling in vacuo or with implicit generalized Born (GB) aqueous solvation models. The GB conformational models clearly depart from those obtained in explicit solvent, and suffer from conformational collapse almost as severe as in vacuo. Thus, the global energy minima in vacuo or with GB are not suitable representations of the unbound state, which can instead be extensively sampled by MD simulations. Many, but not all, MD-simulated compounds displayed some structural similarity to their bioactive structure, supporting the notion of conformational pre-organization for binding. The ligand-protein complexes were also simulated in explicit solvent, to estimate ΔEReorg as an enthalpic difference ΔHReorg between the intramolecular energies in the bound and unbound states. This fresh approach yielded ΔHReorg values⩽6kcal/mol for 18 out of 26 compounds. For three particularly polar compounds 15⩽ΔHReorg⩽20kcal/mol, supporting the notion that ΔHReorg can be substantial. Those large ΔHReorg values correspond to a redistribution of electrostatic interactions upon binding. Overall, the study illustrates how MD simulations offer a promising avenue to characterize the unbound state of medicinal compounds.

Pharmacophore binding motifs for nicotinamide adenine dinucleotide analogues across multiple protein families: a detailed contact-based analysis of the interaction between proteins and NAD(P) cofactors

We have analyzed the protein-binding pharmacophore of NAD and its close analogues in all protein-ligand structures available in the RCSB database as of February 2012; this analysis has then been used to assess the novelty of structures emerging after that date. We show that proteins have evolved diverse pharmacophore motifs for binding the adenine moiety, fewer, but still diverse, motifs for nicotinamide, and a very limited set of motifs for binding the pyrophosphate linker. Our exhaustive analysis includes a pharmacophore contact analysis for over 1900 protein-ligand structures containing NAD analogues; we have benchmarked this set of contacts against nearly 27 000 protein-ligand structures to demonstrate that the diversity of interactions seen with NAD is very similar to that seen for all other ligands. Hence, variation in binding motifs for NAD is not distinct from that observed for other ligands and they show significant variation across protein families.

The mechanism of discrimination between oxidized and reduced coenzyme in the aldehyde dehydrogenase domain of Aldh1l1

Aldh1l1, also known as 10-formyltetrahydrofolate dehydrogenase (FDH), contains the carboxy-terminal domain (Ct-FDH), which is a structural and functional homolog of aldehyde dehydrogenases (ALDHs). This domain is capable of catalyzing the NADP(+)-dependent oxidation of short chain aldehydes to their corresponding acids, and similar to most ALDHs it has two conserved catalytic residues, Cys707 and Glu673. Previously, we demonstrated that in the Ct-FDH mechanism these residues define the conformation of the bound coenzyme and the affinity of its interaction with the protein. Specifically, the replacement of Cys707 with an alanine resulted in the enzyme lacking the ability to differentiate between the oxidized and reduced coenzyme. We suggested that this was due to the loss of a covalent bond between the cysteine and the C4N atom of nicotinamide ring of NADP(+) formed during Ct-FDH catalysis. To obtain further insight into the functional significance of the covalent bond between Cys707 and the coenzyme, and the overall role of the two catalytic residues in the coenzyme binding and positioning, we have now solved crystal structures of Ct-FDH in the complex with thio-NADP(+) and the complexes of the C707S mutant with NADP(+) and NADPH. This study has allowed us to trap the coenzyme in the contracted conformation, which provided a snapshot of the conformational processing of the coenzyme during the transition from oxidized to reduced form. Overall, the results of this study further support the previously proposed mechanism by which Cys707 helps to differentiate between the oxidized and reduced coenzyme during ALDH catalysis.

Study on nicotinamide adenine dinucleotide adsorbed at nano-boehmite/water and nano-corundum/water interfaces

In this study, the adsorption behaviors of nicotinamide adenine dinucleotide (NAD(+)) on nano-boehmite (γ-AlOOH) and nano-corundum (γ-Al(2)O(3)) surfaces were investigated. The results showed that NAD(+) was predominantly adsorbed at the boehmite/water and corundum/water interfaces in outer-sphere fashions by electrostatic interaction between NAD(+) phosphate and surface hydroxyl groups. However, the features of ATR-FTIR spectra suggested that some minor inner-sphere complex should be considered at low pH conditions on corundum surface, which was consistent with the effect of NAD(+) on dissolution rate of corundum. In addition, the adsorption data well fitted with Langmuir and Freundlich isotherms on the boehmite and corundum surfaces, respectively. Also, the Gibbs adsorption energy was negative on the boehmite surface, which indicated that the adsorption behavior was spontaneous.

Linking distinct conformations of nicotinamide adenine dinucleotide with protein fold/function

Nicotinamide adenine dinucleotide (NAD or NADP) are essential cofactor/substrate for enzymes that catalyze redox or nonredox reactions. Because several enzymes involved in NAD(P) metabolism have been implicated in a wide array of diseases, there is great interest in designing inhibitors/activators of these NAD(P)-dependent enzymes based on their structures. Hence, we have elucidated the various distinct enzyme-bound NAD(P) conformations and their correlation with the respective protein fold and function using hierarchical clustering methods. Torsion angles distinguishing enzyme-bound NAD versus NADP conformations and NAD(P) conformations bound to redox versus nonredox enzymes were identified. Although an unusually small χ(N) in diphtheria toxin-bound NAD(+) had been postulated to strain the N-glycosidic bond, thus facilitating catalysis, toxin-bound NAD(+) molecules with χ(N) varying from 0 to 60° were found to exhibit similar C(1D)-N(1N) bond cleavage barriers in water. The findings herein provide useful guidelines in the design of inhibitors/activators of NAD(P)-dependent enzymes that are therapeutic targets.

Getting past the RNA world: the initial Darwinian ancestor

A little-noted result of the confirmation of multiple premises of the RNA-world hypothesis is that we now know something about the dawn organisms that followed the origin of life, perhaps over 4 billion years ago. We are therefore in an improved position to reason about the biota just before RNA times, during the era of the first replicators, the first Darwinian creatures on Earth. An RNA congener still prominent in modern biology is a plausible descendent of these first replicators.

On the diversity of physicochemical environments experienced by identical ligands in binding pockets of unrelated proteins

Most function prediction methods that identify cognate ligands from binding site analyses work on the assumption of molecular complementarity. These approaches build on the conjectured complementarity of geometrical and physicochemical properties between ligands and binding sites so that similar binding sites will bind similar ligands. We found that this assumption does not generally hold for protein-ligand interactions and observed that it is not the chemical composition of ligand molecules that dictates the complementarity between protein and ligand molecules, but that the ligand's share within the functional mechanism of a protein determines the degree of complementarity. Here, we present for a set of cognate ligands a descriptive analysis and comparison of the physicochemical properties that each ligand experiences in various nonhomologous binding pockets. The comparisons in each ligand set reveal large variations in their experienced physicochemical properties, suggesting that the same ligand can bind to distinct physicochemical environments. In some protein ligand complexes, the variation was found to correlate with the electrochemical characteristic of ligand molecules, whereas in others it was disclosed as a prerequisite for the biochemical function of the protein. To achieve binding, proteins were observed to engage in subtle balancing acts between electrostatic and hydrophobic interactions to generate stabilizing free energies of binding. For the presented analysis, a new method for scoring hydrophobicity from molecular environments was developed showing high correlations with experimental determined desolvation energies. The presented results highlight the complexities of molecular recognition and underline the challenges of computational structural biology in developing methods to detect these important subtleties.

Parallel synthesis of a series of non-functional ATP/NAD analogs with activity against trypanosomatid parasites

Non-functional analogs of the cofactors ATP and NAD are putative inhibitors of ATP- or NAD-dependant enzymes. Since pathogenic protozoa rely heavily on the salvage of purine nucleosides from the bloodstream of their host, such compounds are of interest as antiplasmodial and antitrypanosomal agents with a multitude of molecular targets. By replacing the negatively charged phosphate residues with a constrained unsaturated amide spacer and the nicotinamide moiety of NAD with various lipophilic substituents, 15 new ATP/NAD analogs were obtained in screening quantities. In these compounds, a 5'-desoxyadenosine moiety was conserved as key molecular recognition motif. The inhibition of P. falciparum and T. brucei ssp. in a whole parasite in vitro assay is reported.

Molecular dynamics simulation studies of the wild-type, I21V, and I16T mutants of isoniazid-resistant Mycobacterium tuberculosis enoyl reductase (InhA) in complex with NADH: toward the understanding of NADH-InhA different affinities

The increasing prevalence of tuberculosis in many areas of the world, associated with the rise in drug-resistant Mycobacterium tuberculosis (MTB) strains, presents a major threat to global health. InhA, the enoyl-ACP reductase from MTB, catalyzes the nicotinamide adenine dinucleotide (NADH)-dependent reduction of long-chain trans-2-enoyl-ACP fatty acids, an intermediate in mycolic acid biosynthesis. Mutations in the structural gene for InhA are associated with isoniazid resistance in vivo due to a reduced affinity for NADH, suggesting that the mechanism of drug resistance may be related to specific interactions between enzyme and cofactor within the NADH binding site. To compare the molecular events underlying ligand affinity in the wild-type, I21V, and I16T mutant enzymes and to identify the molecular aspects related to resistance, molecular dynamics simulations of fully solvated NADH-InhA (wild-type and mutants) were performed. Although very flexible, in the wild-type InhA-NADH complex, the NADH molecule keeps its extended conformation firmly bound to the enzyme's binding site. In the mutant complexes, the NADH pyrophosphate moiety undergoes considerable conformational changes, reducing its interactions with its binding site and probably indicating the initial phase of ligand expulsion from the cavity. This study should contribute to our understanding of specific molecular mechanisms of drug resistance, which is central to the design of more potent antimycobacterial agents for controlling tuberculosis.

Aplysia californica mediated cyclisation of novel 3'-modified NAD+ analogues: a role for hydrogen bonding in the recognition of cyclic adenosine 5'-diphosphate ribose

Cyclic ADP-ribose mobilizes intracellular Ca2+ in a variety of cells. To elucidate the nature of the interaction between the C3' substituent of cADP-ribose and the cADPR receptor, three analogues of NAD+ modified in the adenosine ribase (xyloNAD+ 3'F-xyloNAD+ and 3'F-NAD+ were chemically synthesised from D-xylose and adenine starting materials. 3'F-NAD+ was readily converted to cyclic 3'F-ADP ribose by the action of the cyclase enzyme derived from the mollusc Aplysia californica. XyloNAD+ and 3'F-xyloNAD+ were cyclised only reluctantly and in poor yield to afford unstable cyclic products. Biological evaluation of cyclic 3'F-ADP ribose for calcium release in sea urchin egg homogenate gave an EC(50) of 1.5+/-0.5 microM. This high value suggests that the ability of the C3' substituent to donate a hydrogen bond is crucial for agonism.

Conformational Analysis of Drug-Like Molecules Bound to Proteins: An Extensive Study of Ligand Reorganization upon Binding

This paper describes a large-scale study on the nature and the energetics of the conformational changes drug-like molecules experience upon binding. Ligand strain energies and conformational reorganization were analyzed with different computational methods on 150 crystal structures of pharmaceutically relevant protein-ligand complexes. The common knowledge that ligands rarely bind in their lowest calculated energy conformation was confirmed. Additionally, we found that over 60% of the ligands do not bind in a local minimum conformation. While approximately 60% of the ligands were calculated to bind with strain energies lower than 5 kcal/mol, strain energies over 9 kcal/mol were calculated in at least 10% of the cases regardless of the method used. A clear correlation was found between acceptable strain energy and ligand flexibility, while there was no correlation between strain energy and binding affinity, thus indicating that expensive conformational rearrangements can be tolerated in some cases without overly penalizing the tightness of binding. On the basis of the trends observed, thresholds for the acceptable strain energies of bioactive conformations were defined with consideration of the impact of ligand flexibility. An analysis of the degree of folding of the bound ligands confirmed the general tendency of small molecules to bind in an extended conformation. The results suggest that the unfolding of hydrophobic ligands during binding, which exposes hydrophobic surfaces to contact with protein residues, could be one of the factors accounting for high reorganization energies. Finally, different methods for conformational analysis were evaluated, and guidelines were defined to maximize the prevalence of bioactive conformations in computationally generated ensembles.

Structural studies on flavin reductase PheA2 reveal binding of NAD in an unusual folded conformation and support novel mechanism of action

The catabolism of toxic phenols in the thermophilic organism Bacillus thermoglucosidasius A7 is initiated by a two-component enzyme system. The smaller flavin reductase PheA2 component catalyzes the NADH-dependent reduction of free FAD according to a ping-pong bisubstrate-biproduct mechanism. The reduced FAD is then used by the larger oxygenase component PheA1 to hydroxylate phenols to the corresponding catechols. We have determined the x-ray structure of PheA2 containing a bound FAD cofactor (2.2 A), which is the first structure of a member of this flavin reductase family. We have also determined the x-ray structure of reduced holo-PheA2 in complex with oxidized NAD (2.1 A). PheA2 is a single domain homodimeric protein with each FAD-containing subunit being organized around a six-stranded beta-sheet and a capping alpha-helix. The tightly bound FAD prosthetic group (K(d) = 10 nm) binds near the dimer interface, and the re face of the FAD isoalloxazine ring is fully exposed to solvent. The addition of NADH to crystalline PheA2 reduced the flavin cofactor, and the NAD product was bound in a wide solvent-accessible groove adopting an unusual folded conformation with ring stacking. This is the first observation of an enzyme that is very likely to react with a folded compact pyridine nucleotide. The PheA2 crystallographic models strongly suggest that reactive exogenous FAD substrate binds in the NADH cleft after release of NAD product. Nanoflow electrospray mass spectrometry data indeed showed that PheA2 is able to bind one FAD cofactor and one FAD substrate. In conclusion, the structural data provide evidence that PheA2 contains a dual binding cleft for NADH and FAD substrate, which alternate during catalysis.

A structurally conserved water molecule in Rossmann dinucleotide-binding domains

A computational comparison of 102 high-resolution (</=1.90 A) enzyme-dinucleotide (NAD, NADP, FAD) complexes was performed to investigate the role of solvent in dinucleotide recognition by Rossmann fold domains. The typical binding site contains about 9-12 water molecules, and about 30% of the hydrogen bonds between the protein and the dinucleotide are water mediated. Detailed inspection of the structures reveals a structurally conserved water molecule bridging dinucleotides with the well-known glycine-rich phosphate-binding loop. This water molecule displays a conserved hydrogen-bonding pattern. It forms hydrogen bonds to the dinucleotide pyrophosphate, two of the three conserved glycine residues of the phosphate-binding loop, and a residue at the C-terminus of strand four of the Rossmann fold. The conserved water molecule is also present in high-resolution structures of apo enzymes. However, the water molecule is not present in structures displaying significant deviations from the classic Rossmann fold motif, such as having nonstandard topology, containing a very short phosphate-binding loop, or having alpha-helix "A" oriented perpendicular to the beta-sheet. Thus, the conserved water molecule appears to be an inherent structural feature of the classic Rossmann dinucleotide-binding domain.

The solute-solvent system: solvent constraints on the conformational dynamics of acetylcholine

The objective of this study was to determine if and how a solvent influences internal motions in a solute molecule. Acetylcholine was chosen as the object of study given its interesting molecular structure and major biological significance. Molecular dynamics simulations were carried out in the vacuum (10 ns), water (5 ns), methanol (5 ns), and octanol (1.5 ns). Seven clusters of conformers were identified, namely, +g+g, -g-g, +gt, -gt, t+g, t-g, and tt, where the gauche and trans labels refer to the dihedral angles tau(2) and tau(3), respectively. As expected, the relative proportion of these conformational clusters was highly solvent-dependent and corresponded to a progressive loss of conformational freedom with increasing molecular weight of the solvent. More importantly, the conformational clusters were used to calculate instantaneous and median angular velocity (omega and omega(M), respectively) and instantaneous and median angular acceleration (alpha and alpha(M), respectively). Angular velocity and angular acceleration were both found to decrease markedly with increasing molecular weight of the solvent, i.e., vacuum (epsilon = 1) > water > methanol > octanol. The decrease from the vacuum to octanol was approximately 40% for tau(2) and approximately 60% for tau(3). Such solvent-dependent constraints on a solute's internal motions may be biologically and pharmacologically relevant.