A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Hantavirus infection: Insights into entry, assembly and pathogenesis

Hantaviruses (family Bunyaviridae, genus Hantavirus) are enveloped viruses incorporating a segmented, negative-sense RNA genome. Each hantavirus is carried by its specific host, either a rodent or an insectivore (shrew), in which the infection is asymptomatic and persistent. In humans, hantaviruses cause Hemorrhagic fever with renal syndrome (HFRS) in Eurasia and Hantavirus cardiopulmonary syndrome (HCPS) in the Americas. In Finland, Puumala virus (genus Hantavirus) is the causative agent of NE, a mild form of HFRS. The HFRS-type diseases are often associated with renal failure and proteinuria that might be mechanistically explained by infected kidney tubular cell degeneration in patients. Previously, it has been shown that non-pathogenic hantavirus, Tula virus (TULV), could cause programmed cell death, apoptosis, in cell cultures. This suggested that the infected kidney tubular degeneration could be caused directly by virus replication. In the first paper of this thesis the molecular mechanisms involved in TULV-induced apoptosis was further elucidated. A virus replication-dependent down-regulation of ERK1/2, concomitantly with the induced apoptosis, was identified. In addition, this phenomenon was not restricted to TULV or to non-pathogenic hantaviruses in general since also a pathogenic hantavirus, Seoul virus, could inhibit ERK1/2 activity. Hantaviruses consist of membrane-spanning glycoproteins Gn and Gc, RNA-dependent RNA polymerase (L protein) and nucleocapsid protein N, which encapsidates the viral genome, and thus forms the ribonucleoprotein (RNP). Interaction between the cytoplasmic tails of viral glycoproteins and RNP is assumed to be the only means how viral genetic material is incorporated into infectious virions. In the second paper of this thesis, it was shown by immunoprecipitation that viral glycoproteins and RNP interact in the purified virions. It was further shown that peptides derived from the cytoplasmic tails (CTs) of both Gn and Gc could bind RNP and recombinant N protein. In the fourth paper the cytoplamic tail of Gn but not Gc was shown to interact with genomic RNA. This interaction was probably rather unspecific since binding of Gn-CT with unrelated RNA and even single-stranded DNA were also observed. However, since the RNP consists of both N protein and N protein-encapsidated genomic RNA, it is possible that the viral genome plays a role in packaging of RNPs into virions. On the other hand, the nucleic acid-binding activity of Gn may have importance in the synthesis of viral RNA. Binding sites of Gn-CT with N protein or nucleic acids were also determined by peptide arrays, and they were largely found to overlap. The Gn-CT of hantaviruses contain a conserved zinc finger (ZF) domain with an unknown function. Some viruses need ZFs in entry or post-entry steps of the viral life cycle. Cysteine residues are required for the folding of ZFs by coordinating zinc-ions, and alkylation of these residues can affect virus infectivity. In the third paper, it was shown that purified hantavirions could be inactivated by treatment with cysteine-alkylating reagents, especially N-ethyl maleimide. However, the effect could not be pin-pointed to the ZF of Gn-CT since also other viral proteins reacted with maleimides, and it was, therefore, impossible to exclude the possibility that other cysteines besides those that were essential in the formation of ZF are required for hantavirus infectivity.

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

Cis-trans peptide variations in structurally similar proteins

The presence of energetically less favourable cis peptides in protein structures has been observed to be strongly associated with its structural integrity and function. Inter-conversion between the cis and trans conformations also has an important role in the folding process. In this study, we analyse the extent of conservation of cis peptides among similar folds. We look at both the amino acid preferences and local structural changes associated with such variations. Nearly 34% of the Xaa-Proline cis bonds are not conserved in structural relatives; Proline also has a high tendency to get replaced by another amino acid in the trans conformer. At both positions bounding the peptide bond, Glycine has a higher tendency to lose the cis conformation. The cis conformation of more than 30% of beta turns of type VIb and IV are not found to be conserved in similar structures. A different view using Protein Block-based description of backbone conformation, suggests that many of the local conformational changes are highly different from the general local structural variations observed among structurally similar proteins. Changes between cis and trans conformations are found to be associated with the evolution of new functions facilitated by local structural changes. This is most frequent in enzymes where new catalytic activity emerges with local changes in the active site. Cis-trans changes are also seen to facilitate inter-domain and inter-protein interactions. As in the case of folding, cis-trans conversions have been used as an important driving factor in evolution.

INDUCTION OF FOLDED STRUCTURES IN DESIGNED PEPTIDES USING CONFORMATIONALLY CONSTRAINED RESIDUES

Folding into compact globular structures, with well-defined modules of secondary structure, appears to be a characteristic of long polypeptide chains, with a specific patterning of coded amino acid residues along the length of sequence. Cooperative hydrogen bond driven secondary structure formation and solvent forces, which contribute favorably to the entropy of folding, by promoting compaction of the polymeric chain, have long been discussed as major determinants of the folding process. First principles design approaches, which use non-coded amino acids, employ an alternative structure directing strategy, by using amino acid residues which exhibit a strong conformational bias for specific regions of the Ramachandran map. This overview of ongoing studies in the authors' laboratory, attempts to explore the use of conformationally restricted amino acid residues in the design of peptides with well-defined secondary structures. Short peptides composed of 20 genetically coded amino acids usually exist in solution as an ensemble of equilibrating conformations. Apolar peptide sequences, which are readily soluble in organic solvents like chloroform and methanol, facilitate formation of structures which are predominately driven by intramolecular hydrogen bond formation. The choice of sequences containing residues with a limited range of conformational choices strongly favors formation of local turn structures, stabilized by short range intramolecular hydrogen bonds. Two residue beta-turns can nucleate either helical or hairpin folding, depending on the precise conformation of the turn segment Restriction of the conformational space available to amino acid residues is easily achieved by introduction of an additional alkyl group at the C alpha carbon atom or by side chain backbone cyclization, as in proline. Studies of synthetic sequences incorporating two prototype residues alpha-aminoisobutyric acid (Aib) and D-proline (DPro) illustrate the utility of the strategy in construction of helices and hairpins. Extensions to the design of conformationally switchable sequences and structurally defined hybrid peptides containing backbone homologated residues are also surveyed.

Identification of local variations within secondary structures of proteins

Secondary-structure elements (SSEs) play an important role in the folding of proteins. Identification of SSEs in proteins is a common problem in structural biology. A new method, ASSP (Assignment of Secondary Structure in Proteins), using only the path traversed by the C atoms has been developed. The algorithm is based on the premise that the protein structure can be divided into continuous or uniform stretches, which can be defined in terms of helical parameters, and depending on their values the stretches can be classified into different SSEs, namely -helices, 3(10)-helices, -helices, extended -strands and polyproline II (PPII) and other left-handed helices. The methodology was validated using an unbiased clustering of these parameters for a protein data set consisting of 1008 protein chains, which suggested that there are seven well defined clusters associated with different SSEs. Apart from -helices and extended -strands, 3(10)-helices and -helices were also found to occur in substantial numbers. ASSP was able to discriminate non--helical segments from flanking -helices, which were often identified as part of -helices by other algorithms. ASSP can also lead to the identification of novel SSEs. It is believed that ASSP could provide a better understanding of the finer nuances of protein secondary structure and could make an important contribution to the better understanding of comparatively less frequently occurring structural motifs. At the same time, it can contribute to the identification of novel SSEs. A standalone version of the program for the Linux as well as the Windows operating systems is freely downloadable and a web-server version is also available at .

Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

A biochemical and structural characterization of Drosophila neuroglian

A number of cell-cell interactions in the nervous system are mediated by immunoglobulin gene superfamily members. For example, neuroglian, a homophilic neural cell adhesion molecule in Drosophila, has an extracellular portion comprising six C- 2 type immunoglobulin-like domains followed by five fibronectin type III (FnIII) repeats. Neuroglian shares this domain organization and significant sequence identity with Ll, a murine neural adhesion molecule that could be a functional homologue. Here I report the crystal structure of a proteolytic fragment containing the first two FnIII repeats of neuroglian (NgFn 1,2) at 2.0Å. The interpretation of photomicrographs of rotary shadowed Ng, the entire extracellular portion of neuroglian, and NgFnl-5, the five neuroglian Fn III domains, is also discussed.

The structure of NgFn 1,2 consists of two roughly cylindrical β-barrel structural motifs arranged in a head-to-tail fashion with the domains meeting at an angle of ~120, as defined by the cylinder axes. The folding topology of each domain is identical to that previously observed for single FnIII domains from tenascin and fibronectin. The domains of NgFn1,2 are related by an approximate two fold screw axis that is nearly parallel to the longest dimension of the fragment. Assuming this relative orientation is a general property of tandem FnIII repeats, the multiple tandem FnIII domains in neuroglian and other proteins are modeled as thin straight rods with two domain zig-zag repeats. When combined with the dimensions of pairs of tandem immunoglobulin-like domains from CD4 and CD2, this model suggests that neuroglian is a long narrow molecule (20 - 30 Å in diameter) that extends up to 370Å from the cell surface.

In photomicrographs, rotary shadowed Ng and NgFn1-5 appear to be highly flexible rod-like molecules. NgFn 1-5 is observed to bend in at least two positions and has a mean total length consistent with models generated from the NgFn1,2 structure. Ng molecules have up to four bends and a mean total length of 392 Å, consistent with a head-to-tail packing of neuroglian's C2-type domains.

Capturing protein dynamics with time-resolved luminescence spectroscopy

The presented doctoral research utilizes time-resolved spectroscopy to characterize protein dynamics and folding mechanisms. We resolve millisecond-timescale folding by coupling time-resolved fluorescence energy transfer (trFRET) to a continuous flow microfluidic mixer to obtain intramolecular distance distributions throughout the folding process. We have elucidated the folding mechanisms of two cytochromes---one that exhibits two-state folding (cytochrome cb₅₆₂) and one that has both a kinetic refolding intermediate ensemble and a distinct equilibrium unfolding intermediate (cytochrome c₅₅₂). Our data reveal that the distinct structural features of cytochrome c₅₅₂ contribute to its thermostability.

We have also investigated intrachain contact dynamics in unfolded cytochrome cb₅₆₂ by monitoring electron transfer, which occurs as the heme collides with a ruthenium photosensitizer, covalently bound to residues along the polypeptide. Intrachain diffusion for chemically denatured proteins proceeds on the microsecond timescale with an upper limit of 0.1 microseconds. The power-law dependence (slope = -1.5) of the rate constants on the number of peptide bonds between the heme and Ru complex indicate that cytochrome cb₅₆₂ is minimally frustrated.

In addition, we have explored the pathway dependence of electron tunneling rates between metal sites in proteins. Our research group has converted cytochrome b₅₆₂ to a c-type cytochrome with the porphyrin covalently bound to cysteine sidechains. We have investigated the effects of the changes to the protein structure (i.e., increased rigidity and potential new equatorial tunneling pathways) on the electron transfer rates, measured by transient absorption, in a series of ruthenium photosensitizer-modified proteins.

Perversão: um fazer gozar

Esse estudo aborda a perversão, tema de extrema relevância política e clínica no âmbito psicanalítico, que se revela como um dos problemas cruciais da psicanálise. Sua importância política imprime-se quando ressalta, com Lacan, que a política da psicanálise é a política da falta-a-ser correlata à ética do desejo e quando propõe que aceitar a diversidade do gozo, com suas múltiplas modalidades, e levar o sujeito a desvelar e a confrontar-se com o seu mais-de-gozar é uma indicação ética que deve orientar a política e a prática do psicanalista. A pesquisa, elaborada a partir dos fundamentos teórico-clínicos de Freud e Lacan, destaca o movimento lógico que delimita a perversão na obra freudiana e verifica uma convergência nas teses de ambos os autores, no que se refere à diferenciação entre a perversidade e a perversão como estrutura clínica. A construção de uma série de casos clínicos e o estudo da vida e da obra de cinco famosos escritores Marquês de Sade, Sacher-Masoch, André Gide, Jean Genet, Yukio Mishima vêm ilustrar que as práticas de gozo perverso não determinam a estrutura perversa. O matema da fantasia sadiana, forjado por Lacan, é tomado para demonstrar a Verleugnung freudiana, o modo que os sujeitos perversos encontram para lidar com castração da mãe/mulher. As fórmulas quânticas da sexuação, formuladas por Lacan, são utilizadas para evidenciar que a diferença entre a neurose e a perversão se explicita na estratégia de gozo que o sujeito utiliza na relação com o seu parceiro. A pesquisa, que se iniciou bibliográfica e se desenvolveu de cunho teórico-clínico, desvela que a clínica da perversão muito pode ensinar aos psicanalistas sobre os quatro conceitos fundamentais da psicanálise, sobre a relação entre o fetiche, a máscara e o semblante, sobre a sublimação, sobre as estratégias de gozo, em particular, o masoquismo. Esse estudo convoca a comunidade analítica ao debate, uma vez que enuncia, através da apresentação de casos de perversão, os impasses da clínica, além de presentificar o real da clínica na condução do tratamento de um caso de perversão.

Viver nas margens: gênero, crime e regulação de conflitos

Esta tese tem como objetivo principal fazer dialogar gênero, cultura popular e crime. O fio condutor do trabalho são as técnicas de conflitos e gestão da violência como interações possuidoras de imagens sexuais (STRATHERN, 2006) analisadas nos encontros de força (MACHADO DA SILVA, 2008) entre moradores e traficantes na favela Nova Holanda, Maré - Rio de Janeiro. Resulta de pesquisa sobre o universo de valentia de galeras Funk e a pacificação do Comando Vermelho analisado a partir de relatos e percursos etnográficos que me foram possíveis através da experiência de moradora cria. Argumento que este processo realiza uma expansão discursiva do crime (FELTRAN, 2006) produzindo perturbações físico-morais (DUARTE, 1987) - expressas pela categoria êmica neurose. As performances proibidas Funk revelam diferentes técnicas de si afirmativas, ofensivas e defensivas parte de um espectro ampliado de representações articuladas por uma gramática patriarcal de guerra e paz responsável pela oposição entre mundos masculinos e femininos. Neste contexto normalizador sãoperformados(BUTLER, 2004) o sujeito homem, a amante e a fiel, enquanto personas populares de valentia que reatualizam a histórica mediação entre o malandro e o bandido (MISSE, 1999). Nas situações vividas, o idioma de gênero permite explorar a produção de fronteiras entre o uso da força e a regulação do convívio, indicando nexos variados de como se processam os conflitos, seja por anulação, mediação e/ou explicitação

mRNA 很可能携带三维遗传信息

The forming mechanism of the three - dimensional structures of proteins，i.e.the mechanism of protein folding，is a basic problem in molecular biology which is still unsolved unitl now. In which a core problem is whether there is the three – dimensional genetic information that decide the three - dimensional structures of proteins. However, the research on this field has mot yet been reported. Recently，we made a comparative study on the folded structures of more than 70 mature messeneger RNAs (mRNAs) and the three - dimensional structures of the proteins encoded by them，it has been found that there exist marked correspondences between their featured structures in the following aspects: 1.The number of the structural units. An RNA molecule can form a secondary structure(stem and loop structure) by the folding and the base pairing of itself. The elementary structural unit of an RNA secondary structure is hairpin(or compound hair pin).The regular structural unit in the secondary structure of a protein is # alpha # - helix or #beta# - sheet . We have found that the hairpin number in the secondary structure of each mature mRNA is equal or approximately equal to the number of the regular secondary structural unis of the encoded protein. 2 .Turning region. Turn is a main structrual element in the secondary structure of a protein, which decides the backbone orientation of a protein molecule to some extent .Our analysis shows that the nucleotide sequence segments in an mRNA which encode the turns of the corresponding protein are overall situated in the turning regions of the mRNA secondary structure such as haipin，bulge loop or multibaranch loops. 3 .The arrangement of structural elements in space. In order to understand the backbone orientation of an RNA molecule and the arangement of its structural elements in space，we have modeled the three一dimensional structure of the mRNA molecule on SGI workstation based on its secondary structure.The result shows that the spatial arrangement of most of the nucleotide sequence segments encoding the structural elements of a protein is consistent with that of these stretural exements in the protein. For instance，the nucleotide sequences corresponding to each pleated sheet of a # beta # - sheet structure are close to each other in the mRNA secondary stucture and in the three - dimensional structure，although some of the nucleotide segments are far apart from each other in the one - dimensional sequence. For another instance，the two triplet codons of cysteines which form a disulphide bridge geneal1y are very close to each other in the mRNA folded structure. In addition，we also analyzed the locations of the codons proline - coding and the distrbution of the nucleotide sequences #alpha# - helix - coding in the folded structures of mRNAs . Some distribution laws have been found. All of these results suggest that the transfer of the genetic information from mRNA to protein not only is one – dimensional but also is three - dime ns ional. That is，there exists the genetic information that decide the three - dimensional structures of proteins. To a certain extent，we could say that the mRNA folding detemines the protein folding. Based on these results，it would be possible to predict the three - dimensional structures of proteins from the primary，secondary and tertiary structures of the m RNAs at a higher accuracy.And more important is that a new clue has been provided to uncover the“spatial coding" of the genetic information.

Secondary structures of mRNAs and splicing of introns in eukaryotes.

After analyzing the secondary structures of 68 exon-intron-exon and the corresponding exon-exon sequence segments, it is found that about 90% of 5' and 3' terminal bases G (splicing sites) of introns are situated in the loops of secondary structures or at the ends of stems near the loops, and most of "G" s in loops are closed to the ends of loops. Approximately 92% of the connecting sites of the adjoining exons also show the similar features. About 82% of the branch point "A" s are situated in loops or at the ends of stems near the loops. Splicing sites and branch points approach each other in space because of the folding.

Stable RNA hairpins in 88 coding regions of human mRNA

RNA hairpins containing UNCG, GNRA, CUUG (N = A, U, C or G, R = G or A) loops are unusually thermodynamic stable and conserved structures. The structural features of these hairpin loops are very special, and they play very important roles in vivo. They are prevalent in rRNA, catalytic RNA and non-coding mRNA. However, the 5' C(UUCG)G 3' hairpin is not found in the folding structure of 88 human mRNA coding regions. It is also different from rRNA in that there is no preference for certain sequences among tetraloops in these 88 mRNA folding structures.

Two-orbit switch-pitch structures

The Hoberman 'switch-pitch ' ball is a transformable structure with a single folding and unfolding path. The underlying cubic structure has a novel mechanism that retains tetrahedral symmetry during folding. Here, we propose a generalized class of structures of a similar type that retain their full symmetry during folding. The key idea is that we require two orbits of nodes for the structure: within each orbit, any node can be copied to any other node by a symmetry operation. Each member is connected to two nodes, which may be in different orbits, by revolute joints. We will describe the symmetry analysis that reveals the symmetry of the internal mechanism modes for a switch-pitch structure. To follow the complete folding path of the structure, a nonlinear iterative predictor-corrector algorithm based on the Newton method is adopted. First, a simple tetrahedral example of the class of two-orbit structures is presented. Typical configurations along the folding path are shown. Larger members of the class of structures are also presented, all with cubic symmetry. These switch-pitch structures could have useful applications as deployable structures.