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.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

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.

Exploring the active site of chorismate mutase by combinatorial mutagenesis and selection: the importance of electrostatic catalysis.

Chorismate mutase (EC 5.4.99.5) catalyzes the intramolecular rearrangement of chorismate to prephenate. Arg-90 in the active site of the enzyme from Bacillus subtilis is in close proximity to the substrate's ether oxygen and may contribute to efficient catalysis by stabilizing the presumed dipolar transition state that would result upon scission of the C--O bond. To test this idea, we have developed a novel complementation system for chorismate mutase activity in Escherichia coli by reengineering parts of the aromatic amino acid biosynthetic pathway. The codon for Arg-90 was randomized, alone and in combination with that for Cys-88, and active clones were selected. The results show that a positively charged residue either at position 88 (Lys) or 90 (Arg or Lys) is essential. Our data provide strong support for the hypothesis that the positive charge is required for stabilization of the transition state of the enzymatic chorismate rearrangement. The new selection system, in conjunction with combinatorial mutagenesis, renders the mechanism of the natural enzyme(s) accessible to further exploration and opens avenues for the improvement of first generation catalytic antibodies with chorismate mutase activity.

Vpr protein of human immunodeficiency virus type 1 forms cation-selective channels in planar lipid bilayers.

A small (96-aa) protein, virus protein R (Vpr), of human immunodeficiency virus type 1 contains one hydrophobic segment that could form a membrane-spanning helix. Recombinant Vpr, expressed in Escherichia coli and purified by affinity chromatography, formed ion channels in planar lipid bilayers when it was added to the cis chamber and when the trans chamber was held at a negative potential. The channels were more permeable to Na+ than to Cl- ions and were inhibited when the trans potential was made positive. Similar channel activity was caused by Vpr that had a truncated C terminus, but the potential dependence of channel activity was no longer seen. Antibody raised to a peptide mimicking part of the C terminus of Vpr (AbC) inhibited channel activity when added to the trans chamber but had no effect when added to the cis chamber. Antibody to the N terminus of Vpr (AbN) increased channel activity when added to the cis chamber but had no effect when added to the trans chamber. The effects of potential and antibodies on channel activity are consistent with a model in which the positive C-terminal end of dipolar Vpr is induced to traverse the bilayer membrane when the opposite (trans) side of the membrane is at a negative potential. The C terminus of Vpr would then be available for interaction with AbC in the trans chamber, and the N terminus would be available for interaction with AbN in the cis chamber. The ability of Vpr to form ion channels in vitro suggests that channel formation by Vpr in vivo is possible and may be important in the life cycle of human immunodeficiency virus type 1 and/or may cause changes in cells that contribute to AIDS-related pathologies.

Distance determination in proteins using designed metal ion binding sites and site-directed spin labeling: application to the lactose permease of Escherichia coli.

As shown in the accompanying paper, the magnetic dipolar interaction between site-directed metal-nitroxide pairs can be exploited to measure distances in T4 lysozyme, a protein of known structure. To evaluate this potentially powerful method for general use, particularly with membrane proteins that are difficult to crystallize, both a paramagnetic metal ion binding site and a nitroxide side chain were introduced at selected positions in the lactose permease of Escherichia coli, a paradigm for polytopic membrane proteins. Thus, three individual cysteine residues were introduced into putative helix IV of a lactose permease mutant devoid of native cysteine residues containing a high-affinity divalent metal ion binding site in the form of six contiguous histidine residues in the periplasmic loop between helices III and IV. In addition, the construct contained a biotin acceptor domain in the middle cytoplasmic loop to facilitate purification. After purification and spin labeling, electron paramagnetic resonance spectra were obtained with the purified proteins in the absence and presence of Cu(II). The results demonstrate that positions 103, 111, and 121 are 8, 14, and > 23 A from the metal binding site. These data are consistent with an alpha-helical conformation of transmembrane domain IV of the permease. Application of the technique to determine helix packing in lactose permease is 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.

Determination of the distance between two spin labels attached to a macromolecule.

An EPR "spectroscopic ruler" was developed using a series of alpha-helical polypeptides, each modified with two nitroxide spin labels. The EPR line broadening due to electron-electron dipolar interactions in the frozen state was determined using the Fourier deconvolution method. These dipolar spectra were then used to estimate the distances between the two nitroxides separated by 8-25 A. Results agreed well with a simple alpha-helical model. The standard deviation from the model system was 0.9 A in the range of 8-25 A. This technique is applicable to complex systems such as membrane receptors and channels, which are difficult to access with high-resolution NMR or x-ray crystallography, and is expected to be particularly useful for systems for which optical methods are hampered by the presence of light-interfering membranes or chromophores.

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.

Intermolecular and Intramolecular Charge Transfer in Functional Molecular Materials

This thesis is devoted to the investigation of inter and intramolecular charge transfer (CT) in molecular functional materials and specifically organic dyes and CT crystals. An integrated approach encompassing quantum-chemical calculations, semiempirical tools, theoretical models and spectroscopic measurements is applied to understand structure-property relationships governing the low-energy physics of these materials. Four main topics were addressed: 1) Spectral properties of organic dyes. Charge-transfer dyes are constituted by electron donor (D) and electron acceptor (A) units linked through bridge(s) to form molecules with different symmetry and dimensionality. Their low-energy physics is governed by the charge resonance between D and A groups and is effectively described by a family of parametric Hamiltonians known as essential-state models. These models account for few electronic states, corresponding to the main resonance structures of the relevant dye, leading to a simple picture that is completed introducing the coupling of the electronic system to molecular vibrations, treated in a non-adiabatic way, and an effective classical coordinate, describing polar solvation. In this work a specific essential-state model was proposed and parametrized for the dye Brilliant Green. The central issue in this work has been the definition of the diabatic states, a not trivial task for a multi-branched chromophore. In a second effort, we have used essential-state models for the description of the early-stage dynamics of excited states after ultrafast excitation. Crucial to this work is the fully non-adiabatic treatment of the coupled electronic and vibrational motion, allowing for a reliable description of the dynamics of systems showing a multistable, broken-symmetry excited state. 2) Mixed-stack CT salts. Mixed-stack (MS) CT crystals are an interesting class of multifunctional molecular materials, where D and A molecules arrange themselves to form stacks, leading to delocalized electrons in one dimension. The interplay between the intermolecular CT, electrostatic interactions, lattice phonons and molecular vibrations leads to intriguing physical properties that include (photoinduced) phase transitions, multistability, antiferromagnetism, ferroelectricity and potential multiferroicity. The standard microscopic model to describe this family of materials is the Modified Hubbard model accounting for electron-phonon coupling (Peierls coupling), electron-molecular vibrations coupling (Holstein coupling) and electrostatic interactions. We adopt and validate a method, based on DFT calculations on dimeric DA structures, to extract relevant model parameters. The approach offers a powerful tool to shed light on the complex physics of MS-CT salts. 3) Charge transfer in organic radical dipolar dyes. In collaboration with the group of Prof. Jaume Veciana (ICMAB- Barcellona), we have studied spectral properties of a special class of CT dyes with D-bridge-A structure where the acceptor group is a stable radical (of the perchlorotriphenylmethyl, PTM, family), leading to an open-shell CT dyes. These materials are of interest since they associate the electronic and optical properties of CT dyes with magnetic properties from the unpaired electron. The first effort was devoted to the parametrization of the relevant essential-state model. Two strategies were adopted, one based on the calculation of the low-energy spectral properties, the other based on the variation of ground state properties with an applied electric field. 4) The spectral properties of organic nanoparticles based on radical species are investigated in collaboration with Dr. I. Ratera (ICMAB- Barcellona). Intriguing spectroscopic behavior was observed pointing to the presence of excimer states. In an attempt to rationalize these findings, extensive calculations (TD-DFT and ZINDO) were performed. The results for the isolated dyes are validated against experimental spectra in solution. To address intermolecular interactions we studied dimeric structures in the gas phase, but the preliminary results obtained do not support excimer formation.

Excitations in Bose-Einstein condensates: collective modes, quantum turbulence and matter wave statistics

In this thesis, we present the generation and studies of a 87Rb Bose-Einstein condensate (BEC) perturbed by an oscillatory excitation. The atoms are trapped in a harmonic magnetic trap where, after an evaporative cooling process, we produce the BEC. In order to study the effect caused by oscillatory excitations, a quadrupole magnetic field time oscillatory is superimposed to the trapping potential. Through this perturbation, collective modes were observed. The dipole mode is excited even for low excitation amplitudes. However, a minimum excitation energy is needed to excite the condensate quadrupole mode. Observing the excited cloud in TOF expansion, we note that for excitation amplitude in which the quadrupole mode is excited, the cloud expands without invert its aspect ratio. By looking these clouds, after long time-of-flight, it was possible to see vortices and, sometimes, a turbulent state in the condensed cloud. We calculated the momentum distribution of the perturbed BECs and a power law behavior, like the law to Kolmogorov turbulence, was observed. Furthermore, we show that using the method that we have developed to calculate the momentum distribution, the distribution curve (including the power law exponent) exhibits a dependence on the quadrupole mode oscillation of the cloud. The randomness distribution of peaks and depletions in density distribution image of an expanded turbulent BEC, remind us to the intensity profile of a speckle light beam. The analogy between matter-wave speckle and light speckle is justified by showing the similarities in the spatial propagation (or time expansion) of the waves. In addition, the second order correlation function is evaluated and the same dependence with distance was observed for the both waves. This creates the possibility to understand the properties of quantum matter in a disordered state. The propagation of a three-dimensional speckle field (as the matter-wave speckle described here) creates an opportunity to investigate the speckle phenomenon existing in dimensions higher than 2D (the case of light speckle).

The effect of the geomagnetic field on cosmic ray energy estimates and large scale anisotropy searches on data from the Pierre Auger Observatory

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than 60 degrees, detected at the Pierre Auger Observatory. the geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the similar to 2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

Fabricación de nanocontactos para células solares: invisibilidad y resonancias plasmónicas

Los dispositivos opto-electrónicos, tales como las células solares, las pantallas planas y los diodos LED (del inglés light emitting diodes), necesitan contactos eléctricos en la cara frontal por la que entra o sale la luz del dispositivo. Estos contactos causan pérdidas por reflexión y absorción de luz (sombra) y por resistencia eléctrica. En una primera aproximación estas pérdidas son contrapuestas, lo que mejora la sombra empeora la resistencia y viceversa. Hasta ahora esto se ha entendido como un compromiso inevitable que limita la eficiencia de conversión energética de los dispositivos opto-electrónicos: disminuir las pérdidas por resistencia eléctrica implica necesariamente aumentar las pérdidas ópticas por sombra. Esta tesis se ha encaminado a tratar de superar esta dificultad a través de la nanoestructuración de la malla de contacto frontal, con especial énfasis en el caso de las células solares de concentración. El objetivo es poder reducir simultáneamente las pérdidas por sombreado y resistencia en serie de la malla. Hemos encontrado, en base a experimentos, teoría y simulaciones, que para tamaños de linea pequeños, en el umbral del régimen de Rayleigh, pero no lo suficientemente pequeños como para que se den las resonancias plasmónicas más intensas (de tipo dipolar), los contactos hacen menos sombra de la que corresponde a su área geométrica. Se puede decir que los contactos se vuelven parcialmente invisibles. En una primera parte de introducción se ha presentado la influencia de la malla en las pérdidas por resistencia en serie producidas en la célula. Se ha analizado el peso de las distintas variables y se ha escogido la reducción del espaciado entre líneas como alternativa a desarrollar. Para no afectar a otras variables, se ha reducido acordemente la anchura de línea manteniendo el factor de sombra geométrico de las células estado del arte. Se ha calculado que para un caso ideal la ganancia puede ser de un 4% absoluto para mallas con líneas de anchura 400-600 nm distribuidas en periodos de 10-20 μm. Se ha visto como otros efectos eléctricos apuntan también a ese rango como óptimo...

Chiral gold(I) vs chiral silver complexes as catalysts for the enantioselective synthesis of the second generation GSK-hepatitis C virus inhibitor

The synthesis of a GSK 2nd generation inhibitor of the hepatitis C virus, by enantioselective 1,3-dipolar cycloaddition between a leucine derived iminoester and tert-butyl acrylate, was studied. The comparison between silver(I) and gold(I) catalysts in this reaction was established by working with chiral phosphoramidites or with chiral BINAP. The best reaction conditions were used for the total synthesis of the hepatitis C virus inhibitor by a four step procedure affording this product in 99% ee and in 63% overall yield. The origin of the enantioselectivity of the chiral gold(I) catalyst was justified according to DFT calculations, the stabilizing coulombic interaction between the nitrogen atom of the thiazole moiety and one of the gold atoms being crucial.

Magnetization reversal of ferromagnetic nanowires studied by magnetic force microscopy

The magnetization reversal of two-dimensional arrays of parallel ferromagnetic Fe nanowires embedded in nanoporous alumina templates has been studied. By combining bulk magnetization measurements (superconducting quantum interference device magnetometry) with field-dependent magnetic force microscopy (MFM), we have been able to decompose the macroscopic hysteresis loop in terms of the irreversible magnetic responses of individual nanowires. The latter are found to behave as monodomain ferromagnetic needles, with hysteresis loops displaced (asymmetric) as a consequence of the strong dipolar interactions between them. The application of field-dependent MFM provides a microscopic method to obtain the hysteresis curve of the array, by simply registering the fraction of up and down magnetized wires as a function of applied field. The observed deviations from the rectangular shape of the macroscopic hysteresis loop of the array can be ascribed to the spatial variation of the dipolar field through the inhomogeneously filled membrane. The system studied proves to be an excellent example of the two-dimensional classical Preisach model, well known from the field of hysteresis modeling and micromagnetism.