Substrate Selectivity and Molecular Adaptation in the Outer Membrane Proteases of Yersinia pestis and Salmonella enterica

Proteolysis is important in bacterial pathogenesis and colonization of animal and plant hosts. In this work I have investigated the functions of the bacterial outer membrane proteases, omptins, of Yersinia pestis and Salmonella enterica. Y. pestis is a zoonotic pathogen that causes plague and has evolved from gastroenteritis-causing Yersinia pseudotuberculosis about 13 000 years ago. S. enterica causes gastroenteritis and typhoid fever in humans. Omptins are transmembrane β-barrels with ten antiparallel β-strands and five surface-exposed loops. The loops are important in substrate recognition, and variation in the loop sequences leads to different substrate selectivities between omptins, which makes omptins an ideal platform to investigate functional adaptation and to alter their polypeptide substrate preferences. The omptins Pla of Y. pestis and PgtE of S. enterica are 75% identical in their amino acid sequences. Pla is a multifunctional protein with proteolytic and non-proteolytic functions, and it increases bacterial penetration and proliferation in the host. Functions of PgtE increase migration of S. enterica in vivo and bacterial survival in mouse macrophages, thus enhancing bacterial spread within the host. Mammalian plasminogen/fibrinolytic system maintains the balance between coagulation and fibrinolysis and participates in several cellular processes, e.g., cell migration and degradation of extracellular matrix proteins. This system consists of activation cascades, which are strictly controlled by several regulators, such as plasminogen activator inhibitor 1 (PAI-1), α2-antiplasmin (α2AP), and thrombin-activatable fibrinolysis inhibitor (TAFI). This work reveals novel interactions of the omptins of Y. pestis and S. enterica with the regulators of the plasminogen/fibrinolytic system: Pla and PgtE inactivate PAI-1 by cleavage at the reactive site peptide bond, and degrade TAFI, preventing its activation to TAFIa. Structure-function relationship studies with Pla showed that threonine 259 of Pla is crucial in plasminogen activation, as it prevents degradation of the plasmin catalytic domain by the omptin and thus maintains plasmin stability. In this work I constructed chimeric proteins between Pla and Epo of Erwinia pyrifoliae that share 78% sequence identity to find out which amino acids and regions in Pla are important for its functions. Epo is neither a plasminogen activator nor an invasin, but it degrades α2AP and PAI-1. Cumulative substitutions towards Pla sequence turned Epo into a Pla-like protein. In addition to threonine 259, loops 3 and 5 are critical in plasminogen activation by Pla. Turning Epo into an invasin required substitution of 31 residues located at the extracellular side of the Epo protein above the lipid bilayer, and also of the β1-strand in the N-terminal transmembrane region of the protein. These studies give an example of how omptins adapt to novel functions that advantage their host bacteria in different ecological niches.

Spectral Tuning and Adaptation to Different Light Environments of Mysid Visual Pigments

In the present thesis, questions of spectral tuning, the relation of spectral and thermal properties of visual pigments, and evolutionary adaptation to different light environments were addressed using a group of small crustaceans of the genus Mysis as a model. The study was based on microspectrophotometric measurements of visual pigment absorbance spectra, electrophysiological measurements of spectral sensitivities of dark-adapted eyes, and sequencing of the opsin gene retrieved through PCR. The spectral properties were related to the spectral transmission of the respective light environments, as well as to the phylogentic histories of the species. The photoactivation energy (Ea) was estimated from temperature effects on spectral sensitivity in the long-wavelength range, and calculations were made for optimal quantum catch and optimal signal-to-noise ratio in the different light environments. The opsin amino acid sequences of spectrally characterized individuals were compared to find candidate residues for spectral tuning. The general purpose was to clarify to what extent and on what time scale adaptive evolution has driven the functional properties of (mysid) visual pigments towards optimal performance in different light environments. An ultimate goal was to find the molecular mechanisms underlying the spectral tuning and to understand the balance between evolutionary adaptation and molecular constraints. The totally consistent segregation of absorption maxima (λmax) into (shorter-wavelength) marine and (longer-wavelength) freshwater populations suggests that truly adaptive evolution is involved in tuning the visual pigment for optimal performance, driven by selection for high absolute visual sensitivity. On the other hand, the similarity in λmax and opsin sequence between several populations of freshwater M. relicta in spectrally different lakes highlights the limits to adaptation set by evolutionary history and time. A strong inverse correlation between Ea and λmax was found among all visual pigments studied in these respects, including those of M. relicta and 10 species of vertebrate pigments, and this was used to infer thermal noise. The conceptual signal-to-noise ratios thus calculated for pigments with different λmax in the Baltic Sea and Lake Pääjärvi light environments supported the notion that spectral adaptation works towards maximizing the signal-to-noise ratio rather than quantum catch as such. Judged by the shape of absorbance spectra, the visual pigments of all populations of M. relicta and M. salemaai used exclusively the A2 chromophore (3, 4-dehydroretinal). A comparison of amino acid substitutions between M. relicta and M. salemaai indicated that mysid shrimps have a small number of readily available tuning sites to shift between a shorter - and a longer -wavelength opsin. However, phylogenetic history seems to have prevented marine M. relicta from converting back to the (presumably) ancestral opsin form, and thus the more recent reinvention of marine spectral sensitivity has been accomplished by some other novel mechanism, yet to be found

Gabapentin: A Stereochemically Constrained gamma Amino Acid Residue in Hybrid Peptide Design

Nature has used the all-alpha-polypeptide backbone of proteins to create a remarkable diversity of folded structures. Sequential patterns of 20 distinct amino adds, which differ only in their side chains, determine the shape and form of proteins. Our understanding of these specific secondary structures is over half a century old and is based primarily on the fundamental elements: the Pauling alpha-helix and beta-sheet. Researchers can also generate structural diversity through the synthesis of polypeptide chains containing homologated (omega) amino acid residues, which contain a variable number of backbone atoms. However, incorporating amino adds with more atoms within the backbone introduces additional torsional freedom into the structure, which can complicate the structural analysis. Fortunately, gabapentin (Gpn), a readily available bulk drug, is an achiral beta,beta-disubstituted gamma amino add residue that contains a cyclohexyl ring at the C-beta carbon atom, which dramatically limits the range of torsion angles that can be obtained about the flanking C-C bonds. Limiting conformational flexibility also has the desirable effect of increasing peptide crystallinity, which permits unambiguous structural characterization by X-ray diffraction methods. This Account describes studies carried out in our laboratory that establish Gpn as a valuable residue in the design of specifically folded hybrid peptide structures. The insertion of additional atoms into polypeptide backbones facilitates the formation of intramolecular hydrogen bonds whose directionality is opposite to that observed in canonical alpha-peptide helices. If hybrid structures mimic proteins and biologically active peptides, the proteolytic stability conferred by unusual backbones can be a major advantage in the area of medicinal chemistry. We have demonstrated a variety of internally hydrogen-bonded structures in the solid state for Gpn-containing peptides, including the characterization of the C-7 and C-9 hydrogen bonds, which can lead to ribbons in homo-oligomeric sequences. In hybrid alpha gamma sequences, district C-12 hydrogen-bonded turn structures support formation of peptide helices and hairpins in longer sequences. Some peptides that include the Gpn residue have hydrogen-bond directionality that matches alpha-peptide helices, while others have the opposite directionality. We expect that expansion of the polypeptide backbone will lead to new classes of foldamer structures, which are thus far unknown to the world of alpha-polypeptides. The diversity of internally hydrogen-bonded structures observed in hybrid sequences containing Gpn shows promise for the rational design of novel peptide structures incorporating hybrid backbones.

Computer-Aided Drug Design of Pyranopyrazoles and Related Compounds for Checkpoint Kinase-1

Checkpoint-1 kinase plays an important role in the G(2)M cell cycle control, therefore its inhibition by small molecules is of great therapeutic interest in oncology. In this paper, we have reported the virtual screening of an in-house library of 2499 pyranopyrazole derivatives against the ATP-binding site of Chk1 kinase using Glide 5.0 program, which resulted in six hits. All these ligands were docked into the site forming most crucial interactions with Cys87, Glu91 and Leu15 residues. From the observed results these ligands are suggested to be potent inhibitors of Chk1 kinase with sufficient scope for further elaboration.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

Deciphering the Distinct Role for the Metal Coordination Motif in the Catalytic Activity of Mycobacterium smegmatis Topoisomerase I

Mycobacterium smegmatis topoisomerase I (Mstopol) is distinct from typical type IA topoisomerases. The enzyme binds to both single- and double-stranded DNA with high affinity, making specific contacts. The enzyme comprises conserved regions similar to type IA topoisomerases from Escherichia coli and other eubacteria but lacks the typically found zinc fingers in the carboxy-terminal domain. The enzyme can perform DNA cleavage m the absence of Mg2+ but religation needs exogenously added Mg2+. One molecule of Mg2+ tightly bound to the enzyme has no role in DNA cleavage but is needed only for the religation reaction. The toprim. (topoisomerase-primase) domain in MstopoI comprising the Mg2+ binding pocket, conserved in both type IA and type II topoisomerases, was subjected to mutagenesis to understand the role of Mg2+, in different steps of the reaction. The residues D108, D110, and E112 of the enzyme, which form the acidic triad in the DXDXE motif, were changed to alanines. D108A mutation resulted in an enzyme that is Mg2+ dependent for DNA cleavage unlike Mstopol and exhibited enhanced DNA cleavage property and reduced religation activity. The mutant was toxic for cell growth, most likely due to the imbalance in cleavage-religation equilibrium. In contrast, the E112A mutant behaved like wild-type enzyme, cleaving DNA in a Mg2+-independent fashion, albeit to a reduced extent. Intra- and intermolecular religation assays indicated specific roles for D108 and E112 residues during the reaction. Together, these results indicate that the D108 residue has a major role during cleavage and religation, while E112 is important for enhancing the efficiency of cleavage. Thus, although architecturally and mechanistically similar to topoisomerase I from E. coli, the metal coordination pattern of the mycobacterial enzyme is distinct, opening up avenues to exploit the enzyme to develop inhibitors.

Development and characterization of lysine based tripeptide analogues as inhibitors of Sir2 activity

Sirtuins are NAD(+) dependent deacetylases that modulate various essential cellular functions. Development of peptide based inhibitors of Sir2s would prove useful both as pharmaceutical agents and as effectors by which downstream cellular alterations can be monitored. Click chemistry that utilizes Huisgen's 1,3-dipolar cycloaddition permits attachment of novel modifications onto the side chain of lysine. Herein, we report the synthesis of peptide analogues prepared using click reactions on N epsilon-propargyloxycarbonyl protected lysine residues and their characterization as inhibitors of Plasmodium falciparum Sir2 activity. The peptide based inhibitors exhibited parabolic competitive inhibition with respect to acetylated-peptide substrate and parabolic non-competitive inhibition with NAD(+) supporting the formation of EI2 and E.NAD(+).I-2 complexes. Cross-competition inhibition analysis with the non-competitive inhibitor nicotinamide (NAM) ruled out the possibility of the NAM-binding site being the second inhibitor binding site, suggesting the presence of a unique alternate pocket commodating the inhibitor. One of these compounds was also found to be a potent inhibitor of the intraerythrocytic growth of P. falciparum with 50% inhibitory concentration in the micromolar range.

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

Titration calorimetric studies to elucidate the specificity of the interactions of polymyxin B with lipopolysaccharides and lipid A

Lipopolysaccharide (LPS), the major cell wall constituent of Gram-negative bacteria, evokes a multitude of biological effects in mammals including pyrogenicity and toxic shock syndrome. Polymyxin B (PmB), a polycationic cyclic peptide, is known to neutralize most of its activities. The nature of the interaction of PmB with LPS and lipid A was investigated by isothermal titration calorimetry. PmB binds to LPS as well as lipid A stoichiometrically and non-co-operatively with micromolar affinity. These interactions are driven primarily by a favourable change in entropy (delta S) and are endothermic in nature. These positive changes in enthalpies decrease with increasing temperature, yielding a heat capacity change, delta Cp, of -2385 J.mol-1.degree-1 for PmB-LPS interactions while the binding of PmB to lipid A displays a delta Cp of -2259 J.mol-1.degree-1. The negative heat capacity changes provide strong evidence for the role of hydrophobic interactions as the driving force for the association of PmB with LPS and lipid A. A correlation of the energetics of these interactions with analyses of the molecular models of PmB suggests that a cluster of solvent-exposed non-polar amino acid side-chains that line one surface of the molecule, together with a ring of positively charged residues on its other surface, are responsible for its strong and stoichiometric binding to LPS.

Functional characterization of the precursor and spliced forms of RecA protein of Mycobacterium tuberculosis

The recA locus of pathogenic mycobacteria differs from that of nonpathogenic species because it contains large intervening sequences nested in the RecA homology region that are excised by an unusual protein-splicing reaction. In vivo assays indicated that Mycobacterium tuberculosis recA partially complemented Escherichia coli recA mutants for recombination and mutagenesis. Further, splicing of the 85 kDa precursor to 38 kDa MtRecA protein was necessary for the display of its activity, in vivo. To gain insights into the molecular basis for partial and lack of complementation by MtRecA and 85 kDa proteins, respectively, we purified both of them to homogeneity. MtRecA protein, but not the 85 kDa form, bound stoichiometrically to single-stranded DNA in the presence of ATP. MtRecA protein was cross-linked to 8-azidoadenosine 5'-triphosphate with reduced efficiency, and kinetic analysis of ATPase activity suggested that it is due to decreased affinity for ATP. In contrast, the 85 kDa form was unable to bind ATP, in the presence or absence of ssDNA and, consequently, was entirely devoid of ATPase activity. Molecular modeling studies suggested that the decreased affinity of MtRecA protein for ATP and the reduced efficiency of its hydrolysis might be due to the widening of the cleft which alters the hydrogen bonds and the contact area between the enzyme and the substrate and changes in the disposition of the amino acid residues around the magnesium ion and the gamma-phosphate. The formation of joint molecules promoted by MtRecA protein was stimulated by SSB when the former was added first. The probability of an association between the lack and partial levels of biological activity of RecA protein(s) to that of illegitimate recombination in pathogenic mycobacteria is considered.

Impacts of agriculture on amphibians at multiple scales

Agriculture-mediated habitat loss and degradation together with climate change are among the greatest global threats to species, communities, and ecosystem functioning. During the last century, more than 50% of the world's wetlands have been lost and agricultural activities have subjected wetland species to increased isolation and decreased quality of habitats. Likewise, as a part of agricultural intensification, the use of pesticides has increased notably, and pesticide residues occur frequently in wetlands making the exposure of wetland organisms to pesticides highly probable. In this thesis, a set of ecotoxicological and landscape ecological studies were carried out to investigate pesticide-effects on tadpoles, and species-habitat relationships of amphibians in agricultural landscapes. The results show that the fitness of R. temporaria tadpoles can be negatively affected by sublethal pesticide concentrations, and that pesticides may increase the costs of response to natural environmental stressors. However, tadpoles may also be able to compensate for some of the negative effects of pesticides. The results further demonstrate that both historic and current-day agricultural land use can negatively impact amphibians, but that in some cases the costs of living in agricultural habitats may only become apparent when amphibians face other environmental stressors, such as drought. Habitat heterogeneity may, however, increase the persistence of amphibians in agricultural landscapes. Hence, the results suggest that amphibians are likely to be affected by agricultural processes that operate at several spatial and temporal scales, and that it is probable that various processes related to current-day agriculture will affect both larval and adult amphibians. The results imply that maintaining dense wetland patterns could enhance persistence of amphibian populations in agricultural habitats, and indicate that heterogeneous landscapes may lower the risk of regional amphibian population declines under extreme weather perturbations.

Structure and Dynamics of Coil-like Molecules by Residual Dipolar Couplings

Nuclear magnetic resonance (NMR) spectroscopy provides us with many means to study biological macromolecules in solution. Proteins in particular are the most intriguing targets for NMR studies. Protein functions are usually ascribed to specific three-dimensional structures but more recently tails, long loops and non-structural polypeptides have also been shown to be biologically active. Examples include prions, -synuclein, amylin and the NEF HIV-protein. However, conformational preferences in coil-like molecules are difficult to study by traditional methods. Residual dipolar couplings (RDCs) have opened up new opportunities; however their analysis is not trivial. Here we show how to interpret RDCs from these weakly structured molecules. The most notable residual dipolar couplings arise from steric obstruction effects. In dilute liquid crystalline media as well as in anisotropic gels polypeptides encounter nematogens. The shape of a polypeptide conformation limits the encounter with the nematogen. The most elongated conformations may come closest whereas the most compact remain furthest away. As a result there is slightly more room in the solution for the extended than for the compact conformations. This conformation-dependent concentration effect leads to a bias in the measured data. The measured values are not arithmetic averages but essentially weighted averages over conformations. The overall effect can be calculated for random flight chains and simulated for more realistic molecular models. Earlier there was an implicit thought that weakly structured or non-structural molecules would not yield to any observable residual dipolar couplings. However, in the pioneering study by Shortle and Ackerman RDCs were clearly observed. We repeated the study for urea-denatured protein at high temperature and also observed indisputably RDCs. This was very convincing to us but we could not possibly accept the proposed reason for the non-zero RDCs, namely that there would be some residual structure left in the protein that to our understanding was fully denatured. We proceeded to gain understanding via simulations and elementary experiments. In measurements we used simple homopolymers with only two labelled residues and we simulated the data to learn more about the origin of RDCs. We realized that RDCs depend on the position of the residue as well as on the length of the polypeptide. Investigations resulted in a theoretical model for RDCs from coil-like molecules. Later we extended the studies by molecular dynamics. Somewhat surprisingly the effects are small for non-structured molecules whereas the bias may be large for a small compact protein. All in all the work gave clear and unambiguous results on how to interpret RDCs as structural and dynamic parameters of weakly structured proteins.

NUVIEW:software for display and interactive manipulation of nucleic acid models

The NUVIEW software package allows skeletal models of any double helical nucleic acid molecule to be displayed out a graphics monitor and to apply various rotations, translations and scaling transformations interactively, through the keyboard. The skeletal model is generated by connecting any pair of representative points, one from each of the bases in the basepair. In addition to the above mentioned manipulations, the base residues can be identified by using a locator and the distance between any pair of residues can be obtained. A sequence based color coded display allows easy identification of sequence repeats, such as runs of Adenines. The real time interactive manipulation of such skeletal models for large DNA/RNA double helices, can be used to trace the path of the nucleic acid chain in three dimensions and hence get a better idea of its topology, location of linear or curved regions, distances between far off regions in the sequence etc. A physical picture of these features will assist in understanding the relationship between base sequence, structure and biological function in nucleic acids.

Identification of alpha- and beta-Hydroxy Acid Containing Cyclodepsipeptides in Natural Peptide Mixtures Using Negative Ion Mass Spectrometry

Natural peptide libraries often contain cyclodepsipeptides containing alpha or beta hydroxy residues. Extracts of fungal hyphae of Isaria yield a microheterogenous cyclodepsipeptide mixture in which two classes of molecules can be identified by mass spectral fragmentation of negative ions. In the case of isaridins, which contain an alpha-hydroxy residue and a beta-amino acid residue, a characteristic product ion corresponding to a neutral loss of 72 Da is obtained. hi addition, neutral loss of water followed by a 72 Da loss is also observed. Two distinct modes of fragmentation rationalize the observed product ion distribution. The neutral loss of 72 Da has also been obtained for a roseotoxin component, which is also an alpha-hydroxy residue containing cyclodepsipeptide. In the case of isariins, which contain a beta-hydroxy acid residue, ring opening and subsequent loss of the terminal residue as an unsaturated ketene fragment, rationalizes the observed product ion formation. Fragmentation of negative ions provide characteristic neutral losses, which are diagnostic of the presence of alpha-hydroxy or beta-hydroxy residues.