New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

Functional dynamics in the active site of the ribonuclease binase

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-μs time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a β-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum kcat (1,500 s−1) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

Molecular cloning and sequence analysis of the complestatin biosynthetic gene cluster

Streptomyces lavendulae produces complestatin, a cyclic peptide natural product that antagonizes pharmacologically relevant protein–protein interactions including formation of the C4b,2b complex in the complement cascade and gp120-CD4 binding in the HIV life cycle. Complestatin, a member of the vancomycin group of natural products, consists of an α-ketoacyl hexapeptide backbone modified by oxidative phenolic couplings and halogenations. The entire complestatin biosynthetic and regulatory gene cluster spanning ca. 50 kb was cloned and sequenced. It consisted of 16 ORFs, encoding proteins homologous to nonribosomal peptide synthetases, cytochrome P450-related oxidases, ferredoxins, nonheme halogenases, four enzymes involved in 4-hydroxyphenylglycine (Hpg) biosynthesis, transcriptional regulators, and ABC transporters. The nonribosomal peptide synthetase consisted of a priming module, six extending modules, and a terminal thioesterase; their arrangement and domain content was entirely consistent with functions required for the biosynthesis of a heptapeptide or α-ketoacyl hexapeptide backbone. Two oxidase genes were proposed to be responsible for the construction of the unique aryl-ether-aryl-aryl linkage on the linear heptapeptide intermediate. Hpg, 3,5-dichloro-Hpg, and 3,5-dichloro-hydroxybenzoylformate are unusual building blocks that repesent five of the seven requisite monomers in the complestatin peptide. Heterologous expression and biochemical analysis of 4-hydroxyphenylglycine transaminon confirmed its role as an aminotransferase responsible for formation of all three precursors. The close similarity but functional divergence between complestatin and chloroeremomycin biosynthetic genes also presents a unique opportunity for the construction of hybrid vancomycin-type antibiotics.

Studies on the biosynthesis of taxol: the taxane carbon skeleton is not of mevalonoid origin.

A cell culture of Taxus chinensis was established to produce the diterpene 2alpha,5alpha,10beta,14beta-tetra-acetoxy4 ++ +(20),11-taxadiene (taxuyunnanine C) in 2.6% (dry weight) yield. The incorporation of [U-13C6]glucose, [1-13C]glucose, and [1,2-13C2]acetate into this diterpene was analyzed by NMR spectroscopy. Label from [1,2-13C2]acetate was diverted to the four acetyl groups of taxuyunnanine C, but not to the taxane ring system. Label from [1-13C]glucose and [U-13C6]glucose was efficiently incorporated into both the taxane ring system and the acetyl groups. The four isoprenoid moieties of the diterpene showed identical labeling patterns. The analysis of long-range 13C13C couplings in taxuyunnanine C obtained from an experiment with [U-13C6]glucose documents the involvement of an intramolecular rearrangement in the biosynthesis of the isoprenoid precursor. The labeling patterns are inconsistent with the mevalonate pathway. The taxoid data share important features with the alternative pathway of isoprenoid biosynthesis operating in certain eubacteria Rohmer, M., Knani, M., Simonin, P., Sutter, B. & Sahm, H. (1993) Biochem. J. 295, 517-524].

Sample size determination in combinatorial chemistry.

Combinatorial chemistry is gaining wide appeal as a technique for generating molecular diversity. Among the many combinatorial protocols, the split/recombine method is quite popular and particularly efficient at generating large libraries of compounds. In this process, polymer beads are equally divided into a series of pools and each pool is treated with a unique fragment; then the beads are recombined, mixed to uniformity, and redivided equally into a new series of pools for the subsequent couplings. The deviation from the ideal equimolar distribution of the final products is assessed by a special overall relative error, which is shown to be related to the Pearson statistic. Although the split/recombine sampling scheme is quite different from those used in analysis of categorical data, the Pearson statistic is shown to still follow a chi2 distribution. This result allows us to derive the required number of beads such that, with 99% confidence, the overall relative error is controlled to be less than a pregiven tolerable limit L1. In this paper, we also discuss another criterion, which determines the required number of beads so that, with 99% confidence, all individual relative errors are controlled to be less than a pregiven tolerable limit L2 (0 < L2 < 1).

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution.

The measurement of dipolar contributions to the splitting of 15N resonances of 1H-15N amide pairs in multidimensional high-field NMR spectra of field-oriented cyanometmyoglobin is reported. The splittings appear as small field-dependent perturbations of normal scalar couplings. Assignment of more than 90 resonances to specific sequential sites in the protein allows correlation of the dipolar contributions with predictions based on the known susceptibility and known structure of the protein. Implications as an additional source of information for protein structure determination in solution are discussed.

An automated multiplex oligonucleotide synthesizer: development of high-throughput, low-cost DNA synthesis.

An automated oligonucleotide synthesizer has been developed that can simultaneously and rapidly synthesize up to 96 different oligonucleotides in a 96-well microtiter format using phosphoramidite synthesis chemistry. A modified 96-well plate is positioned under reagent valve banks, and appropriate reagents are delivered into individual wells containing the growing oligonucleotide chain, which is bound to a solid support. Each well has a filter bottom that enables the removal of spent reagents while retaining the solid support matrix. A seal design is employed to control synthesis environment and the entire instrument is automated via computer control. Synthesis cycle times for 96 couplings are < 11 min, allowing a plate of 96 20-mers to be synthesized in < 5 hr. Oligonucleotide synthesis quality is comparable to commercial machines, with average coupling efficiencies routinely > 98% across the entire 96-well plate. No significant well-to-well variations in synthesis quality have been observed in > 6000 oligonucleotides synthesized to date. The reduced reagent usage and increased capacity allow the overall synthesis cost to drop by at least a factor of 10. With the development of this instrument, it is now practical and cost-effective to synthesize thousands to tens of thousands of oligonucleotides.

Comparison of the defensive elicitation activity of cerato-populin and cerato-platanin: a structural model

Plants can defend themselves from potential pathogenic microorganisms relying on a complex interplay of signaling pathways: activation of the MAPK cascade, transcription of defense related genes, production of reactive oxygen species, nitric oxide and synthesis of other defensive compounds such as phytoalexins. These events are triggered by the recognition of pathogen’s effectors (effector-triggered immunity) or PAMPs (PAMP-triggered immunity). The Cerato Platanin Family (CPF) members are Cys-rich proteins secreted and localized on fungal cell walls, involved in several aspects of fungal development and pathogen-host interactions. Although more than hundred genes of the CPF have been identified and analyzed, the structural and functional characterization of the expressed proteins has been restricted only to few members of the family. Interestingly, those proteins have been shown to bind chitin with diverse affinity and after foliar treatment they elicit defensive mechanisms in host and non-host plants. This property turns cerato platanins into interesting candidates, worth to be studied to develop new fungal elicitors with applications in sustainable agriculture. This study focus on cerato-platanin (CP), core member of the family and on the orthologous cerato-populin (Pop1). The latter shows an identity of 62% and an overall homology of 73% with respect to CP. Both proteins are able to induce MAPKs phosphorylation, production of reactive oxygen species and nitric oxide, overexpression of defense’s related genes, programmed cell death and synthesis of phytoalexins. CP, however, when compared to Pop1, induces a faster response and, in some cases, a stronger activity on plane leaves. Aim of the present research is to verify if the dissimilarities observed in the defense elicitation activity of these proteins can be associated to their structural and dynamic features. Taking advantage of the available CP NMR structure, Pop1’s 3D one was obtained by homology modeling. Experimental residual dipolar couplings and 1H, 15N, 13C resonance assignments were used to validate the model. Previous works on CPF members, addressed the highly conserved random coil regions (loops b1-b2 and b2-b3) as sufficient and necessary to induce necrosis in plants’ leaves: that region was investigated in both Pop1 and CP. In the two proteins the loops differ, in their primary sequence, for few mutations and an insertion with a consequent diversification of the proteins’ electrostatic surface. A set of 2D and 3D NMR experiments was performed to characterize both the spatial arrangement and the dynamic features of the loops. NOE data revealed a more extended network of interactions between the loops in Pop1 than in CP. In addition, in Pop1 we identified a salt bridge Lys25/Asp52 and a strong hydrophobic interaction between Phe26/Trp53. These structural features were expected not only to affect the loops’ spatial arrangement, but also to reduce the degree of their conformational freedom. Relaxation data and the order parameter S2 indeed highlighted reduced flexibility, in particular for loop b1-b2 of Pop1. In vitro NMR experiments, where Pop1 and CP were titrated with oligosaccharides, supported the hypothesis that the loops structural and dynamic differences may be responsible for the different chitin-binding properties of the two proteins: CP selectively binds tetramers of chitin in a shallow groove on one side of the barrel defined by loops b1-b2, b2-b3 and b4-b5, Pop1, instead, interacts in a non-specific fashion with oligosaccharides. Because the region involved in chitin-binding is also responsible for the defense elicitation activity, possibly being recognized by plant's receptors, it is reasonable to expect that those structural and dynamic modifications may also justify the different extent of defense elicitation. To test that hypothesis, the initial steps of a protocol aimed to the identify a receptor for CP, in silico, are presented.

Spin glass phase in the four-state three-dimensional Potts model

We perform numerical simulations, including parallel tempering, a four-state Potts glass model with binary random quenched couplings using the JANUS application-oriented computer. We find and characterize a glassy transition, estimating the critical temperature and the value of the critical exponents. Nevertheless, the extrapolation to infinite volume is hampered by strong scaling corrections. We show that there is no ferromagnetic transition in a large temperature range around the glassy critical temperature. We also compare our results with those obtained recently on the “random permutation” Potts glass.

Prediction of hidden multiferroic order in graphene zigzag ribbons

An electronic phase with coexisting magnetic and ferroelectric order is predicted for graphene ribbons with zigzag edges. The electronic structure of the system is described with a mean-field Hubbard model that yields results very similar to those of density functional calculations. Without further approximations, the mean-field theory is recasted in terms of a BCS wave function for electron-hole pairs in the edge bands. The BCS coherence present in each spin channel is related to spin-resolved electric polarization. Although the total electric polarization vanishes, due to an internal phase locking of the BCS state, strong magnetoelectric effects are expected in this system. The formulation naturally accounts for the two gaps in the quasiparticle spectrun, Δ0 and Δ1, and relates them to the intraband and interband self-energies.

Intramolecular surfaces for vicinal proton–proton coupling constants 3JHH

Equations for the intramolecular surfaces of the 3JHH coupling constants in ethane, ethylene, and acetylene are formulated, and the corresponding coefficients are estimated from calculations at the DFT/B3LYP level. The chosen variables are changes in bond lengths, in the torsion angle φ between the coupled protons Ha and Hb, in bond angles, and in dihedral angles. The 3JHH surface of ethane is formulated as an extended Karplus equation with the coefficients of a truncated Fourier series on the torsion angle φ expanded as second-order Taylor series in the chosen variables taking into account the invariance of 3JHH under reflections and rotations of nuclear coordinates. Partial vibrational contributions from linear and square terms corresponding to changes in the geometry of the Ha − Ca − Cb − Hb fragment are important while those from cross terms are small with a few exceptions. The 3JHH surface of ethane is useful to predict contributions to 3JHH from changes in local geometry of derivatives but vibrational contributions are predicted less satisfactorily. The predicted values at the B3LYP/BS2 level of the 3JHH couplings (vibrational contributions at 300 K) from equilibrium geometries are 9.79 (−0.17) for acetylene, and 17.08 (1.93) and 10.73(0.93) for the trans and cis couplings of ethylene.

O direito administrativo na sociedade : a dimensão societal do direito administrativo entre a autopoiese dos sistemas sociais e a função intersistémica do sistema jurídico

Tese de doutoramento, Direito (Ciências Jurídico-Políticas, vertente de Direito Administrativo), Universidade de Lisboa, Faculdade de Direito, 2016

A Comprehensive Discussion of HMBC Pulse Sequences: 4. Establishing Two-Bond Correlations from HMBC and Related Experiments

The utility of the HMBC experiment for structure elucidation is unquestionable, but the nature of the coupling pathways leading to correlations in an HMBC experiment creates the potential for misinterpretation. This misinterpretation potential is intimately linked to the size of the long-range heteronuclear couplings involved, and may become troublesome in those cases of a particularly strong 2JCH correlation that might be mistaken for a 3JCH correlation or a 4JCH correlation of appreciable strength that could be mistaken for a weaker 3JCH correlation. To address these potential avenues of confusion, work from several laboratories has been focused on the development of what might be considered “coupling pathway edited” long-range heteronuclear correlation experiments that are derived from or related to the HMBC experiment. The first example of an effort to address the problems associated with correlation path length was seen in the heteronucleus-detected XCORFE experiment described by Reynolds and co-workers that predated the development of the HMBC experiment. Proton-detected analogs of the HMBC experiment intended to differentiate 2JCH correlations from nJCH correlations where n = 3, 4, include the 2J,3J-HMBC, HMBC-RELAY, H2BC, edited-HMBC, and HAT H2BC experiments. The principles underlying the critical components of each of these experiments are discussed and experimental verification of the results that can be obtained using model compounds are shown. This contribution concludes with a brief discussion of the 1,1-ADEQUATE experiments that provide an alternative means of identifying adjacent protonated and non-protonated carbon correlations by exploiting 1JCC correlations at natural abundance.