Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Molecular cloning and characterization of an invertebrate cellular retinoic acid binding protein

We have cloned a cDNA and gene from the tobacco hornworm, Manduca sexta, which is related to the vertebrate cellular retinoic acid binding proteins (CRABPs). CRABPs are members of the superfamily of lipid binding proteins (LBPs) and are thought to mediate the effects of retinoic acid (RA) on morphogenesis, differentiation, and homeostasis. This discovery of a Manduca sexta CRABP (msCRABP) demonstrates the presence of a CRABP in invertebrates. Compared with bovine/murine CRABP I, the deduced amino acid sequence of msCRABP is 71% homologous overall and 88% homologous for the ligand binding pocket. The genomic organization of msCRABP is conserved with other CRABP family members and the larger LBP superfamily. Importantly, the promoter region contains a motif that resembles an RA response element characteristic of the promoter region of most CRABPs analyzed. Three-dimensional molecular modeling based on postulated structural homology with bovine/murine CRABP I shows msCRABP has a ligand binding pocket that can accommodate RA. The existence of an invertebrate CRABP has significant evolutionary implications, suggesting CRABPs appeared during the evolution of the LBP superfamily well before vertebrate/invertebrate divergence, instead of much later in evolution in selected vertebrates.

Molecular basis for the recognition of two structurally different major histocompatibility complex/peptide complexes by a single T-cell receptor

2C is a typical alloreactive cytotoxic T lymphocyte clone that recognizes two different ligands. These ligands are adducts of the allo-major histocompatibility complex (MHC) molecule H-2Ld and an endogenous octapeptide, and of the self-MHC molecule H-2Kb and another peptide. MHC-binding and T-cell assays with synthetic peptides in combination with molecular modeling studies were employed to analyze the structural basis for this crossreactivity. The molecular surfaces of the two complexes differ greatly in densities and distributions of positive and negative charges. However, modifications of the peptides that increase similarity decrease the capacities of the resulting MHC peptide complexes to induce T-cell responses. Moreover, the roles of the peptides in ligand recognition are different for self- and allo-MHC-restricted T-cell responses. The self-MHC-restricted T-cell responses were finely tuned to recognition of the peptide. The allo-MHC-restricted responses, on the other hand, largely ignore modifications of the peptide. The results strongly suggest that adaptation of the T-cell receptor to the different ligand structures, rather than molecular mimicry by the ligands, is the basis for the crossreactivity of 2C. This conclusion has important implications for T-cell immunology and for the understanding of immunological disorders.

A combined kinetic and modeling study of the catalytic center subsites of human angiogenin.

Kinetic analysis and molecular modeling have been used to map the ribonucleolytic center of angiogenin (Ang). Pyrimidine nucleotides were found to interact very weakly with Ang, consistent with the inaccessible B1 pyrimidine binding site revealed by x-ray crystallography. Ang also lacks an effective phosphate binding site on the 5' side of B1. Although the B2 site that preferentially binds purines on the 3' side of B1 is also weak, its associated phosphate subsites make substantial contributions: both 3',5'-ADP and 5'-ADP have Ki values 6-fold lower than for 5'-AMP, and adding a 3'-phosphate to the substrate CpA increases Kcat/Km by 9-fold. Thus Ang has a functional P2 site on the 3' side of B2 and a site for a second phosphate on the 5' side of B2. Modeling of an Ang-d(ApTpApA) complex suggested that Arg-5 forms part of the P2 site and that a 2'-phosphate might bind more tightly than a 3'-phosphate. Both predictions were confirmed kinetically. The subsite map obtained by this combined approach indicated that 5'-diphosphoadenosine 2'-phosphate might be a more potent inhibitor than any of the nucleotides tested thus far. Indeed, its Ki value of 150 microM is 50-fold lower than that for the best nucleotide previously reported and 400-fold lower than the Km for the best dinucleotide substrate. This compound may serve as a suitable starting point for the eventual design of tight-binding inhibitors of Ang as antiangiogenic agents for human therapy.

Similar antigenic surfaces, rather than sequence homology, dictate T-cell epitope molecular mimicry.

Molecular mimicry, normally defined by the level of primary-sequence similarities between self and foreign antigens, has been considered a key element in the pathogenesis of autoimmunity. Here we describe an example of molecular mimicry between two overlapping peptides within a single self-antigen, both of which are recognized by the same human self-reactive T-cell clone. Two intervening peptides did not stimulate the T-cell clone, even though they share nine amino acids with the stimulatory peptides. Molecular modeling of major histocompatibility complex class II-peptide complexes suggests that both of the recognized peptides generate similar antigenic surfaces, although these are composed of different sets of amino acids. The molecular modeling of a peptide shifted one residue from the stimulatory peptide, which was recognized in the context of the same HLA molecule by another T-cell clone, generated a completely different antigenic surface. Functional studies using truncated peptides confirmed that the anchor residues of the two "mimicking" epitopes in the HLA groove differ. Our results show, for two natural epitopes, how molecular mimicry can occur and suggest that studies of potential antigenic surfaces, rather than sequence similarity, are necessary for analyzing suspected peptide mimicry.

A molecular anchor for stabilizing triple-helical DNA.

Molecular modeling has been used to predict that 2,6-disubstituted amidoanthraquinones, and not the 1,4 series, should preferentially interact with and stabilize triple-stranded DNA structures over duplex DNA. This is due to marked differences in the nature of chromophore-base stacking and groove accessibility for the two series. A DNA foot-printing method that monitors the extent of protection from DNase I cleavage on triplex formation has been used to examine the effects of a number of synthetic isomer compounds in the 1,4 and 2,6 series. The experimental results are in accord with the predicted behavior and confirm that the 1,4 series bind preferentially to double- rather than triple-stranded DNA, whereas the isomeric 2,6 derivatives markedly favor binding to triplex DNA.

Antichagásicos potenciais: síntese e modelagem molecular de híbridos de hidrazonas e liberadores de óxido nítrico

A doença de Chagas é uma parasitose extremamente negligenciada, cujo agente etiológico é o protozoário Trypanosoma cruzi. Atualmente, 21 países da América Latina são considerados regiões endêmicas, onde 75-90 milhões de pessoas estão expostas à infecção, 6-7 milhões estão infectadas e mais de 41 mil novos casos surgem por ano. Entretanto, apenas os fármacos nifurtimox e benznidazol estão disponíveis no mercado. Estes, além da baixa eficácia na fase crônica da parasitose, apresentam diversos efeitos adversos, sendo que no Brasil apenas o benznidazol é utilizado. Este fato mostra a importância de se ampliar o número de fármacos disponíveis e propor quimioterapia mais eficaz para o tratamento da doença de Chagas. Como forma de contribuir para essa busca, este trabalho objetiva a síntese de compostos híbridos bioisostéricos N-acilidrazônicos e sulfonilidrazônicos, contendo grupo liberador de óxido nítrico, com potencial de interação com cisteíno-proteases parasitárias, tais como a cruzaína. Nestes derivados, os grupos liberadores de óxido nítrico utilizados foram os grupos furoxano (contendo substituinte metílico e fenílico) e éster nitrato. Propôs-se a variação de anéis aromáticos substituídos e não-substituídos, com o intuito de avaliar a possível relação estrutura-atividade (REA) desses análogos. Até o momento, somente os compostos da série N-acilidrazônica tiveram avaliação biológica realizada. Os valores de IC50 dos compostos na forma amastigota do parasita variaram entre >100 a 2,88 µM, sendo este último valor comparável ao fármaco de referência. A atividade inibitória frente à cruzaína foi de 25,2 µM a 2,2 µM. Já a liberação de óxido nítrico foi avaliada pelo método indireto de detecção de nitrato e os valores variaram entre 52,0 µM e 4.232,0 µM. Estes são bem inferiores ao composto padrão, além de não se identificar correlação direta entre a atividade biológica e a liberação de NO. Na sequência, os dois compostos mais ativos (6 e 14) foram submetidos a estudos de permeabilidade e de citotoxicidade. O composto 6 foi considerado o de maior permeabilidade segundo o Sistema de Classificação Biofarmacêutica (SCB) e todos os compostos apresentaram a taxa de fluxo menor que 2, indicando a ausência de mecanismo de efluxo. Na avaliação do potencial citotóxico desses compostos em células humanas, o derivado 6 apresentou índice de seletividade superior ao do benznidazol. Em estudos de modelagem molecular usando análise exploratória de dados (HCA e PCA), propriedades estéricas/geométricas e eletrônicas foram consideradas as mais relevantes para a atividade biológica. Além disso, estudos de docking mostraram que a posição do grupo nitro no anel aromático é importante para a interação com a cruzaína. Ademais o composto 6 não provocou mudanças significativas no ciclo celular e na fragmentação de DNA em células humanas, mostrando-se como líder promissor para futuros estudos in vivo. Atividade tripanomicida, citotoxicidade, potencial de liberação de NO e estudos de permeabilidade dos 23 derivados sulfonilidrazônicos e ésteres nitrato estão sendo avaliados.

Cytokine analysis in Atlantic bottlenose dolphins: Molecular characterization of IL-4, IL-6, and IL-10

The health status of wild and captive Atlantic Bottlenose dolphins ( Tersiops truncatis) is difficult to ascertain. Mass strandings of these animals have been attributed to pollutants, as well as bacterial infections. Using human Enzyme Linked Immuno-Assays (ELISA) for immunological cytokines, I measured soluble cytokine levels with respect to their health status. In a retrospective analysis of dolphin sera, there was a trend of higher cytokine levels in “sick” animals. I cultured dolphin lymphocytes in the presence of a mitogen (PHA), a super antigen (Staph-A), Lipopolysaccharide (LPS), and a calcium flux inducer (PMA). Levels of messenger RNA, from these cultured cells, were assayed with Polymerase Chain Reaction (PCR) using primers for the human cytokines IL-2, IL-4, IL-6, IL-10, Tumor Necrosis Factor, and Interferon gamma. Only IL-4, IL-6, and IL-10 messages were obtained, inferring similar nucleotide homology to the human primer sequences. The PCR products were sequenced. Sixteen IL-4 sequences, twelve IL-6 sequences and seven IL-10 sequences were obtained and analyzed. Each cytokine exhibited the same nucleotide sequence in all dolphins examined. There was no difference in the cytokine profile in response to the various stimuli. The derived amino acid composition for each of the dolphin cytokines was used for molecular modeling, which showed that dolphin IL-4, IL-6, and IL-10 were structurally similar to the corresponding proteins of Perissodactyla. ^

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

A New Platinum Complex With Tryptophan: Synthesis, Structural Characterization, Dft Studies And Biological Assays In Vitro Over Human Tumorigenic Cells.

A new platinum(II) complex with the amino acid L-tryptophan (trp), named Pt-trp, was synthesized and characterized. Elemental, thermogravimetric and ESI-QTOF mass spectrometric analyses led to the composition [Pt(C11H11N2O2)2]⋅6H2O. Infrared spectroscopic data indicate the coordination of trp to Pt(II) through the oxygen of the carboxylate group and also through the nitrogen atom of the amino group. The (13)C CP/MAS NMR spectroscopic data confirm coordination through the oxygen atom of the carboxylate group, while the (15)N CP/MAS NMR data confirm coordination of the nitrogen of the NH2 group to the metal. Density functional theory (DFT) studies were applied to evaluate the cis and trans coordination modes of trp to platinum(II). The trans isomer was shown to be energetically more stable than the cis one. The Pt-trp complex was evaluated as a cytotoxic agent against SK-Mel 103 (human melanoma) and Panc-1 (human pancreatic carcinoma) cell lines. The complex was shown to be cytotoxic over the considered cells.

Synthesis, screening for antiacetylcholinesterase activity and binding mode prediction of a new series of [3-(disubstituted-phosphate)-4,4,4-trifluoro-butyl]-carbamic acid ethyl esters

A series of nine new [3-(disubstituted-phosphate)-4,4,4-trifluoro-butyl]-carbamic acid ethyl esters (phosphate-carbamate compounds) was obtained through the reaction of (4,4,4-trifluoro-3-hydroxybut-1-yl)-carbamic acid ethyl esters with phosphorus oxychloride followed by the addition of alcohols. The products were characterized by ¹H, 13C, 31P, and 19F NMR spectroscopy, GC-MS, and elemental analysis. All the synthesized compounds were screened for acetylcholinesterase (AChE) inhibitory activity using the Ellman method. All compounds containing phosphate and carbamate pharmacophores in their structures showed enzyme inhibition, being the compound bearing the diethoxy phosphate group (2b) the most active compound. Molecular modeling studies were performed to investigate the detailed interactions between AChE active site and small-molecule inhibitor candidates, providing valuable structural insights into AChE inhibition.

DFT and electrochemical studies on nortriptyline oxidation sites

A study on the possible sites of oxidation and epoxidation of nortriptyline was performed using electrochemical and quantum chemical methods; these sites are involved in the biological responses (for example, hepatotoxicity) of nortriptyline and other similar antidepressants. Quantum chemical studies and electrochemical experiments demonstrated that the oxidation and epoxidation sites are located on the apolar region of nortriptyline, which will useful for understanding the molecule`s activity. Also, for the determination of the compound in biological fluids or in pharmaceutical formulations, we propose a useful analytical methodology using a graphite-polyurethane composite electrode, which exhibited the best performance when compared with boron-doped diamond or glassy carbon surfaces.

Cruzain inhibition by hydroxymethylnitrofurazone and nitrofurazone: investigation of a new target in Trypanosoma cruzi

Nitrofurazone (NF) and its derivative, hydroxymethylnitrofurazone (NFOH), have presented antichagasic activity. NFOH has higher activity and lower mutagenicity. The aim of this work was to assess whether NF and its derivative NFOH would also be inhibitors of cruzain, besides their trypanothione reductase inhibitory activity. In vitro cruzain inhibition tests were performed for both compounds, and the 50% inhibitory concentration (IC(50)) for NF and NFOH presented values of 22.83 +/- 1.2 mu M and 10.55 +/- 0.81 mu M, respectively. AM1 semi-empirical molecular modeling studies were performed to understand the activity of the compounds, corroborating the observed cruzain inhibitory activity.

Probing the electronic delocalization in a cyclic pyrazine ruthenium cluster hexamer

[Ru(3)O(CH(3)COO)(6)(pz)(CO)](6) is a cyclic hexamer species encompassing six triangular ruthenium cluster centers bridged by pyrazine ligands. The electronic communication among the cluster units strongly depends on their oxidation states, and has been successfully probed by means of cyclic voltammetry and UV-vis spectroelectrochemistry. (C) 2010 Elsevier B.V. All rights reserved.

Enzymatic and structural characterization of a basic phospholipase A(2) from the sea anemone Condylactis gigantea

This work aimed at the isolation and structural/functional characterization of a phospholipase A(2) (CgPLA(2)) from the extract of the anemone Condylactis gigantea. CgPLA2 was isolated with a high purity level through three chromatographic steps, showing pT8.6 and molecular weights of 14,500 and 29,000 for the monomer and dimer, respectively. CgPLA2 showed a high catalytic activity upon fluorescent phospholipids inducing no direct hemolytic activity. This enzyme, which is Ca2+-dependent, showed a lower stability against temperature and pH variations when compared with snake venom enzymes. The enzymatic activity was significantly reduced or completely abolished after chemical modification of CgPLA2 with BPB. Its cDNA was then obtained, with 357 base pairs which codified for a mature protein of 119 amino acid residues. A comparative analysis of the primary structure of CgPLA2 revealed 84%, 61%, 43% and 42% similarity to the PLA2s from Adamsia carciniopados, Nematostella vectensis, Vipera russelli russelli and Both raps jararacussu, respectively. (C) 2010 Elsevier Masson SAS. All rights reserved.