Modeling the noncovalent interactions at the metabolite binding site in purine riboswitches

We present gas phase quantum chemical studies on the metabolite binding interactions in two important purine riboswitches, the adenine and guanine riboswitches, at the B3LYP/6-31G(d,p) level of theory. In order to gain insights into the strucutral basis of their discriminative abilities of regulating gene expression, the structural properties and binding energies for the gas phase optimized geometries of the metabolite bound binding pocket are analyzed and compared with their respective crystal geometries. Kitaura-Morokuma analysis has been carried out to calculate and decompose the interaction energy into various components. NBO and AIM analysis has been carried out to understand the strength and nature of binding of the individual aptamer bases with their respective purine metabolites. The Y74 base, U in case of adenine riboswitch and C in case of guanine riboswitch constitutes the only differentiating element between the two binding pockets. As expected, with W:W cis G:C74 interaction contributing more than 50% of the total binding energy, the interaction energy for metabolite binding as calculated for guanine (-46.43 Kcal/mol) is nearly double compared to the corresponding value for that of adenine (-24.73 Kcal/mol) in the crystal context. Variations in the optimized geometries for different models and comparison of relative contribution to metabolite binding involving four conserved bases reveal the possible role of U47:U51 W:H trans pair in the conformational transition of the riboswitch from the metabolite free to metabolite bound state. Our results are also indicative of significant contributions from stacking and magnesium ion interactions toward cooperativity effects in metabolite recognition.     

Requirements for noncovalent binding of vaccinia topoisomerase I to duplex DNA.

Vaccinia DNA topoisomerase binds duplex DNA and forms a covalent adduct at sites containing a conserved sequence element 5'(C/T)CCTT decreases in the scissile strand. Distinctive aspects of noncovalent versus covalent interaction emerge from analysis of the binding properties of Topo(Phe-274), a mutated protein which is unable to cleave DNA, but which binds DNA noncovalently. Whereas DNA cleavage by wild type enzyme is most efficient with 'suicide' substrates containing fewer than 10 base pairs distal to the scissile bond, optimal noncovalent binding by Topo(Phe-274) requires at least 10-bp of DNA 3' of the cleavage site. Thus, the region of DNA flanking the pentamer motif serves to stabilize the noncovalent topoisomerase-DNA complex. This result is consistent with the downstream dimensions of the DNA binding site deduced from nuclease footprinting. Topo(Phe-274) binds to duplex DNA lacking the consensus pentamer with 7-10-fold lower affinity than to CCCTT-containing DNA.     

Covalent and noncovalent DNA binding by mutants of vaccinia DNA topoisomerase I.

Analysis of vaccinia topoisomerase mutants that are impaired in DNA relaxation has allowed the identification of amino acid residues required for the transesterification step of catalysis. Missense mutations of wild-type residues Gly-132----Asp and Arg-223----Gln rendered the protein inert in formation of the covalent enzyme-DNA complex and hence completely inactive in DNA relaxation. Mutations of Thr-147----Ile and Gly-132----Ser caused severe defects in covalent adduct formation that correlated with the extent of inhibition of relaxation. None of these point mutations had an effect on noncovalent DNA binding sufficient to account for the defect in relaxation. Deletion of amino- or carboxyl-terminal portions of the polypeptide abrogated noncovalent DNA binding. Two distinct topoisomerase-DNA complexes were resolved by native gel electrophoresis. One complex, which was unique to those proteins competent in covalent adduct formation, contained topoisomerase bound to the 5'-portion of the incised DNA strand. The 3'-segment of the cleaved strand had dissociated spontaneously. This complex was isolated and shown to catalyze transfer of the covalently bound DNA to a heterologous acceptor oligonucleotide, thereby proving that the covalent adduct between protein and duplex DNA is a true intermediate in strand breakage and reunion. The role of the active site region of eukaryotic topoisomerase in determining sensitivity or resistance to camptothecin was examined by converting the active site region of the resistant vaccinia enzyme (SKRAY274) to that of the drug-sensitive yeast enzyme (SKINY). The SKINY mutation did not alter the resistance of the vaccinia enzyme to the cleavage-enhancing effects of camptothecin.     

Effects of noncovalent and covalent FAD binding on the redox and catalytic properties of p-cresol methylhydroxylase

Each flavoprotein subunit (alpha or PchF) of the alpha(2)beta(2) flavocytochrome p-cresol methylhydroxylase (PCMH) from Pseudomonas putida contains FAD covalently attached to Tyr384. PCMH oxidizes p-cresol to 4-hydroxybenzyl alcohol, which is oxidized subsequently by PCMH to 4-hydroxybenzaldehyde. The Y384F mutant form of PchF (apo-PchF[Y384F]) displayed stoichiometric noncovalent FAD binding. PchF[Y384F]FAD associated with the cytochrome subunit (beta or PchC) (producing PCMH[Y384F]), although not as avidly as with wild-type PchF containing covalently bound FAD (PchF(C)). Dramatic increases in the two-electron E(m,7) (NHE) values for FAD were observed when it bound noncovalently to either apo-PchF or apo-PchF[Y384F], and the two-electron E(m,7) value for FAD was increased further by about 75 mV upon covalent binding to PchF, i.e., PchF(C). The E(m,7) values increased by approximately 20 and 45 mV, respectively, when PchF(C) and PchF[Y384F]FAD associated with PchC. The two-electron E(m,7) for covalently bound FAD in PCMH is 84 mV, the highest measured for a flavoprotein. The values for the one-electron redox potentials (E(m,7), NHE) for FAD were measured also for various forms of PchF. Under anaerobiosis, the reduction of PchF[Y384F]FAD by substrates was similar to that observed previously for PchF containing noncovalently bound FAD. Stopped-flow kinetic studies indicated a rapid substrate reduction of the FAD and heme in PCMH[Y384F] which produced PchF[Y384F]FAD(rad) x PchC, the mutant enzyme containing the flavin radical and reduced heme. These experiments also revealed a slow reduction of unassociated PchC(ox) by PchF[Y384F]FAD(rad) x PchC. Steady-state kinetic studies of the reaction of PCMH[Y384F] with p-cresol indicated that the K(m) for this substrate was unchanged relative to that of PCMH, but that the k(cat) was diminished by an order of magnitude. The data indicate that the covalent attachment of FAD to PchF assists catalysis by raising the E(m,7) of the flavin. Contributions to this effect likely result from conformational changes.     

Small,low nanomolar,noncovalent thrombin inhibitors lacking a group to fill the 'distal binding pocket'

Use of a chlorophenoxyacetamide P1 group with a pyridinone acetamide P2/P3 gave an exceptionally potent thrombin inhibitor (K(i)=40 pM). Truncation of the molecule at the N-terminus gave unique, low nanomolar, non-covalent thrombin inhibitors which do not have a group to fill thrombin's 'distal binding pocket'. A co-crystal structure indicates the importance of an intramolecular hydrogen bond between the P1 side chain and P1/P2 amide link in this series.     

Structural plasticity and noncovalent substrate binding in the GroEL apical domain. A study using electrospay ionization mass spectrometry and fluorescence binding studies

Advances in understanding how GroEL binds to non-native proteins are reported. Conformational flexibility in the GroEL apical domain, which could account for the variety of substrates that GroEL binds, is illustrated by comparison of several independent crystallographic structures of apical domain constructs that show conformational plasticity in helices H and I. Additionally, ESI-MS indicates that apical domain constructs have co-populated conformations at neutral pH. To assess the ability of different apical domain conformers to bind co-chaperone and substrate, model peptides corresponding to the mobile loop of GroES and to helix D from rhodanese were studied. Analysis of apical domain-peptide complexes by ESI-MS indicates that only the folded or partially folded apical domain conformations form complexes that survive gas phase conditions. Fluorescence binding studies show that the apical domain can fully bind both peptides independently. No competition for binding was observed, suggesting the peptides have distinct apical domain-binding sites. Blocking the GroES-apical domain-binding site in GroEL rendered the chaperonin inactive in binding GroES and in assisting the folding of denatured rhodanese, but still capable of binding non-native proteins, supporting the conclusion that GroES and substrate proteins have, at least partially, distinct binding sites even in the intact GroEL tetradecamer.     

Role of noncovalent binding of 11-cis-retinal to opsin in dark adaptation of rod and cone photoreceptors

Regeneration of visual pigments of vertebrate rod and cone photoreceptors occurs by the initial noncovalent binding of 11-cis-retinal to opsin, followed by the formation of a covalent bond between the ligand and the protein. Here, we show that the noncovalent interaction between 11-cis-retinal and opsin affects the rate of dark adaptation. In rods, 11-cis-retinal produces a transient activation of the phototransduction cascade that precedes sensitivity recovery, thus slowing dark adaptation. In cones, 11-cis-retinal immediately deactivates phototransduction. Thus, the initial binding of the same ligand to two very similar G protein receptors, the rod and cone opsins, activates one and deactivates the other, contributing to the remarkable difference in the rates of rod and cone dark adaptation.     

Territrem B, a tremorgenic mycotoxin that inhibits acetylcholinesterase with a noncovalent yet irreversible binding mechanism

Territrem B (TRB) is a fungal metabolite isolated from Aspergillus terreus shown previously to be a potent and irreversible inhibitor of acetylcholinesterase (AChE). In the present study, a number of binding and inhibition assays were carried out to further characterize the inhibitory effect of TRB. The results indicate that the binding of TRB (a) is much more selective than a well characterized selective inhibitor of AChE, BW284C51, (b) adopts a one-to-one stoichiometry with the enzyme, (c) cannot be undone by an AChE-regenerating oxime agent, which contrasts the ability of 8 M urea to release AChE-bound TRB, (d) is enhanced by high concentration NaCl but prevented, unless preincubated, by Triton X-100, and (e) exhibits quasi-first order kinetics with an overall inhibition constant of 0.01 nM(-1) min(-1). Together these results suggest a very different irreversible binding (a noncovalent type) from that of the covalent type, which involves typical irreversible AChE inhibitors such as diisopropylfluorophosphate and neostigmine. According to the prediction of a molecular modeling study, the distinct AChE inhibitory characteristics of TRB may arise from the inhibitor being noncovalently trapped within a unique active-site gorge structure of the enzyme. It was predicted that an optimal TRB. AChE binding would position a narrowing connection of the TRB structure at a constricted area near the entrance of the gorge, thereby providing a structural basis for the observed irreversible binding.     

Kinetics of large-ligand binding to one-dimensional lattices: theory of irreversible binding

Exact solutions are obtained for the time dependence of the extent of irreversible binding of ligands that cover more than one lattice site to a homogeneous one-dimensional lattice. The binding may be cooperative or noncooperative and the lattice either finite or infinite. Although the form of the solution is most convenient when the ligand concentration is buffered, exact numerical or approximate analytical solutions, including upper and lower bounds, can be derived for the case of variable ligand concentration as well. The physical reason behind the relative simplicity of the kinetics of irreversible as opposed to reversible binding in such systems is discussed.     

Separation and evaluation of the covalent and noncovalent interactions which contribute to the binding of pyridoxal 5'-phosphate to D-serine apodehydratase.

The equilibrium constant (KX) for the reaction D-serine dehydratase + pyridoxamine-P in equilibrium KX D-serine apodehydratase: pyridoxamine-P + pyridoxal-P was determined. At 25 degreees, pH 7.80, KX increases from 5.4 times 10-minus 5 to 21 times 10-minus 5 as T/2 is increased from 0.33 to 0.66. A value of 1.3 times 10-minus 4 M at 25 degrees, pH 7.80, T/2 0.33 for the equilibrium constant (KPMP) for dissociation of pyridoxamine-P from D-serine apodehydratase was determined from the ratio of the equilibrium constant for dissociation of pyridoxal-P from holoenzyme to KX. Pyridoxamine-P and the thiazolidine, formed from pyridoxal-P and cysteine, were found to have similar affinities for D-serine apodehydratase. Using the affinities of these derivatives as a measure of the noncovalent interactions between cofactor and protein, it was possible to estimate the contribution of the Schiff base linkage to the stability of the complex formed between pyridoxal-P and protein. The covalent Schiff base linkage in the holoenzyme was found to be no more stable than the Schiff base linkage formed between 6-aminocaproic acid and pyridoxal-P. The contribution of noncovalent interactions to the stability of the cofactor-protein complex was shown to be at least 20 to 40 times greater than the contribution of the covalent Schiff base linkage.     

On the theory of niche overlap

The factors which are likely to limit niche overlap are studied for a class of idealized biological communities in which several species compete on a onedimensional continuum of resources, e.g., food size. Although the shape of the resource spectrum restricts the composition of possible equilibrium communities, and although this restriction is increasingly stringent as niche overlap increases, such considerations do not in principle put a limit to niche overlap. In analysing the stability of these equilibrium communities, two qualitatively different circumstances need be distinguished. In a strictly unvarying (deterministic) environment, stability sets no limit to the degree of overlap, short of complete congruence. However, in a randomly fluctuating (stochastic) environment, community stability requires that the average food sizes for species adjacent on the resource spectrum must differ by an amount roughly equal to the standard deviation in the food size taken by either individual species. This limit to species packing has a weak (logarithmic) dependence on the degree of environmental variance.This mathematical result is very robust, as is shown by considering, inter alia, a wide range of resource spectrum shapes, and a variety of shapes for the functions describing how the species utilize the resource. The effect of including additional resource dimensions is estimated.The general conclusion from the model is that there is an effective limit to niche overlap in the real world, and that this limit is insensitive to the degree of environmental fluctuation, unless it be very severe. This conclusion seems in accord with an increasing body of field data.     

On the diffusion theory of phyllotaxis

An inhibitor diffusion theory of phyllotaxis is examined in the steady-state approximation for cylindrical shoot apex models. The model calculations give rise naturally to common patterns of spiral phyllotaxis, as well as to higher whorled patterns. The model also predicts commonly observed subpatterns of axillary organs superimposed on primary phyllotaxis patterns. Application of the model to phyllotaxis patterns in other organisms and in flowers is proposed.     

On the theory of amoeboid movements

It is shown how the application of the principle of maximum energy transformation, discussed in previous papers, leads to expressions which give the probability of occurrence of a pseudopod of a given size and duration in function of other physical constants of the amoeba.     

On the theory of lateral inhibition

Zusammenfassung Das mathematische Modell für das Prinzip der lateralen Inhibition in der Theorie der optischen Perzeption führt auf ein System nichtlinearer Gleichungen für n reelle Variable. Dieses System wird auf Lösbarkeit und eindeutige Lösbarkeit untersucht. Es zeigt sich, daß die Gleichung als Bedingung für die stationären Zustände eines geeigneten zeitabhängigen Systems zu deuten ist. Hier kann man ein diskretes und ein kontinuierliches Modell einführen. In beiden Fällen kann die Frage der Existenz der Lösungen und der Stabilität einigermaßen vollständig geklärt werden. Eine Verallgemeinerung auf kontinuierlich viele Raumvariable ist möglich.     

Study of the fragmentation patterns of the phosphate-arginine noncovalent bond

Our previous work has highlighted the role of certain amino acid residues, mainly two or more adjacent arginine on one peptide and two or more adjacent glutamate, or aspartate, or a phosphorylated residue on the other in the formation of noncovalent complexes (NCX) between peptides. In the present study, we employ ESI-MS to investigate the gas-phase stability and dissociation pathways of the NCX of a basic peptide VLRRRRKRVN, an epitope from the third intracellular loop of the dopamine D(2) receptor, with the phosphopetide SVSTDpTpSAE, an epitope from the cannabinoid CB1 carboxyl terminus. ESI-MS/MS analysis of the NCX between VLRRRRKRVN and SVSTDpTpSAE suggests two dissociation pathways for the NCX. The major pathway is the disruption of the electrostatic interactions between the Arg residues and the phosphate groups, while an alternative pathway is also recorded, in which the complex is dissociated along the covalent bond between the oxygen from either Thr or Ser and HPO(3). To verify the alternative pathway, we have used an ion trap instrument to conduct MS(3) analysis on the product ions of both dissociation pathways.     

C1 repressor of phage P1 is inactivated by noncovalent binding of P1 Coi protein.

The temperate phage P1 encodes two genes whose products antagonize the action of the phage's C1 repressor of lytic functions, namely a distantly linked antirepressor gene, ant, and a closely linked c1 inactivator gene, coi. Starting with an inducible coi-recombinant plasmid, Coi protein was overproduced and purified to near homogeneity. By using a DNA mobility shift assay we demonstrate that Coi protein inhibits the operator binding of the C1 repressors of the closely related P1 and P7 phages. Coi protein (Mr = 7,600) exerts its C1-inactivating function by forming a complex with the C1 repressor (Mr = 32,500) at a molar ratio of about 1:1, as shown by density gradient centrifugation and gel filtration. C1 repressor and Coi protein are recovered in active form from the complex, suggesting that noncovalent interactions are the sole requirements for complex formation. The interplay of repressor and antagonists operating in the life cycle of P1 is discussed.     

On the theory of the cardiovascular system

Most theoretical studies of the circulation have focussed on the transmission line properties of arteries. Only a small number of papers have dealt with the circulation as a closed (lumped) system with two pumps connected by the lesser and greater circulation (Beneken, inCirculatory Analog Computers, No. Holland Publ. Co., Amsterdam, 1963; Defares,et al., inCirculatory Analog Computers, No. Holland Publ. Co., Amsterdam, 1963; Grodins,Quart. Rev. of Biology,34, 93, 1959; Guyton,Cardiac Output and its Regulation, Saunders Publ. Co., New York, 1963). F. W. Cope's recent studies in this journal (Bull. Math. Biophysics,22, 19, 1960;23, 337, 1961;24, 137, 1962) deal with essentially the same questions, although here the circuit is not “closed”. We have attempted to extend the analysis of the areflex (closed) circulation. The complete study is reported elsewhere (Defares,et al., Acta Physiol, et Parmac. Neerl., 1963).