Engineered protease inhibitors based on sunflower trypsin inhibitor-1 (SFTI-1) provide insights into the role of sequence and conformation in Laskowski mechanism inhibition

Laskowski inhibitors regulate serine proteases by an intriguing mode of action that involves deceiving the protease into synthesizing a peptide bond. Studies exploring naturally occurring Laskowski inhibitors have uncovered several structural features that convey the inhibitor's resistance to hydrolysis and exceptional binding affinity. However, in the context of Laskowski inhibitor engineering, the way that various modifications intended to fine-tune an inhibitor's potency and selectivity impact on its association and dissociation rates remains unclear. This information is important as Laskowski inhibitors are becoming increasingly used as design templates to develop new protease inhibitors for pharmaceutical applications. In this study, we used the cyclic peptide, sunflower trypsin inhibitor-1 (SFTI-1), as a model system to explore how the inhibitor's sequence and structure relate to its binding kinetics and function. Using enzyme assays, MD simulations and NMR spectroscopy to study SFTI variants with diverse sequence and backbone modifications, we show that the geometry of the binding loop mainly influences the inhibitor's potency by modulating the association rate, such that variants lacking a favourable conformation show dramatic losses in activity. Additionally, we show that the inhibitor's sequence (including both the binding loop and its scaffolding) influences its potency and selectivity by modulating both the association and the dissociation rates. These findings provide new insights into protease inhibitor function and design that we apply by engineering novel inhibitors for classical serine proteases, trypsin and chymotrypsin and two kallikrein-related peptidases (KLK5 and KLK14) that are implicated in various cancers and skin diseases.

The application of islanded photovoltaic cells and battery storage on the reliability of distribution networks

Renewable energy resources, in particularly PV and battery storage are increasingly becoming part of residential and agriculture premises to manage their electricity consumption. This thesis addresses the tremendous technical, financial and planning challenges for utilities created by these increases, by offering techniques to examine the significance of various renewable resources in electricity network planning. The outcome of this research should assist utilities and customers for adequate planning that can be financially effective.

Random existence of charge ordered stripes and its influence on the magnetotransport properties of La0.6Sr0.4MnO3 perovskite substituted with diamagnetic ions at Mn sublattice

Phase-singular solid solutions of La0.6Sr0.4Mn1-yMeyO3 (0 <= y <= 0.3) [Me=Li1+, Mg2+, Al3+, Ti4+, Nb5+, Mo6+ or W6+] [LSMey] perovskite of rhombohedral symmetry (space group: R (3) over barc) have been prepared wherein the valence of the diamagnetic substituent at Mn site ranged from 1 to 6. With increasing y-content in LSMey, the metal-insulator (TM-I) transition in resistivity-temperature rho(T) curves shifted to low temperatures. The magnetization studies M(H) as well as the M(T) indicated two groups for LSMey. (1) Group A with Me=Mg, Al, Ti, or Nb which are paramagnetic insulators (PIs) at room temperature with low values of M (< 0.5 mu(B)/Mn); the magnetic transition [ferromagnetic insulator (FMI)-PI] temperature (T-C) shifts to low temperatures and nearly coincides with that of TM-I and the maximum magnetoresistance (MR) of similar to 50% prevails near T-C (approximate to TM-I). (2) Group-B samples with Me=Li, Mo, or W which are FMIs with M-s=3.3-3.58 mu(B)/Mn and marginal reduction in T-C similar to 350 K as compared to the undoped LSMO (T-C similar to 378 K). The latter samples show large temperature differences Delta T=T-c-TM-I, reaching up to similar to 288 K. The maximum MR (similar to 60%) prevails at low temperatures corresponding to the M-I transition TM-I rather than around T-C. High resolution lattice images as well as microscopy analysis revealed the prevalence of inhomogeneous phase mixtures of randomly distributed charge ordered-insulating (COI) bistripes (similar to 3-5 nm width) within FMI charge-disordered regions, yet maintaining crystallographically single phase with no secondary precipitate formation. The averaged ionic radius < r(B)>, valency, or charge/radius ratio < CRR > cannot be correlated with that of large Delta T; hence cannot be used to parametrize the discrepancy between T-C and TM-I. The M-I transition is controlled by the charge conduction within the electronically heterogeneous mixtures (COI bistripes+FMI charge disordered); large MR at TM-I suggests that the spin-ordered FM-insulating regions assist the charge transport, whereas the T-C is associated with the bulk spin ordered regions corresponding to the FMI phase of higher volume fraction of which anchors the T-C to higher temperatures. The present analysis showed that the double-exchange model alone cannot account for the wide bifurcation of the magnetic and electric transitions, contributions from the charge as well as lattice degrees of freedom to be separated from spin/orbital ordering. The heterogeneous phase mixtures (COI+FMI) cannot be treated as of granular composite behavior. (c) 2008 American Institute of Physics.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Seismic Tomography and Earthquake Mechanism Beneath the Central Fennoscandian Shield

The aim of this thesis was to study the seismic tomography structure of the earth s crust together with earthquake distribution and mechanism beneath the central Fennoscandian Shield, mainly in southern and central Finland. The earthquake foci and some fault plane solutions are correlated with 3-D images of the velocity tomography. The results are discussed in relation to the stress field of the Shield and with other geophysical, e.g. geomagnetic, gravimetric, tectonic, and anisotropy studies of the Shield. The earthquake data of the Fennoscandian Shield has been extracted from the Nordic earthquake parameter data base which was founded at the time of inception of the earthquake catalogue for northern Europe. Eight earlier earthquake source mechanisms are included in a pilot study on creating a novel technique for calculating an earthquake fault plane solution. Altogether, eleven source mechanisms of shallow, weak earthquakes are related in the 3-D tomography model to trace stresses of the crust in southern and central Finland. The earthquakes in the eastern part of the Fennoscandian Shield represent low-active, intraplate seismicity. Earthquake mechanisms with NW-SE oriented horizontal compression confirm that the dominant stress field originates from the ridge-push force in the North Atlantic Ocean. Earthquakes accumulate in coastal areas, in intersections of tectonic lineaments, in main fault zones or are bordered by fault lines. The majority of Fennoscandian earthquakes concentrate on the south-western Shield in southern Norway and Sweden. Onwards, epicentres spread via the ridge of the Shield along the west-coast of the Gulf of Bothnia northwards along the Tornio River - Finnmark fault system to the Barents Sea, and branch out north-eastwards via the Kuusamo region to the White Sea Kola Peninsula faults. The local seismic tomographic method was applied to find the terrane distribution within the central parts of the Shield the Svecofennian Orogen. From 300 local explosions a total of 19765 crustal Pg- and Sg-wave arrival times were inverted to create independent 3-D Vp and Vs tomographic models, from which the Vp/Vs ratio was calculated. The 3-D structure of the crust is presented as a P-wave and for the first time as an S-wave velocity model, and also as a Vp/Vs-ratio model of the SVEKALAPKO area that covers 700x800 km2 in southern and central Finland. Also, some P-wave Moho-reflection data was interpolated to image the relief of the crust-mantle boundary (i.e. Moho). In the tomography model, the seismic velocities vary smoothly. The lateral variations are larger for Vp (dVp =0.7 km/s) than for Vs (dVs =0.4 km/s). The Vp/Vs ratio varies spatially more distinctly than P- and S-wave velocities, usually from 1.70 to 1.74 in the upper crust and from 1.72 to 1.78 in the lower crust. Schist belts and their continuations at depth are associated with lower velocities and lower Vp/Vs ratios than in the granitoid areas. The tomography modelling suggests that the Svecofennian Orogen was accreted from crustal blocks ranging in size from 100x100 km2 to 200x200 km2 in cross-sectional area. The intervening sedimentary belts have ca. 0.2 km/s lower P- and S-wave velocities and ca. 0.04 lower Vp/Vs ratios. Thus, the tomographic model supports the concept that the thick Svecofennian crust was accreted from several crustal terranes, some hidden, and that the crust was later modified by intra- and underplating. In conclusion, as a novel approach the earthquake focal mechanism and focal depth distribution is discussed in relation to the 3-D tomography model. The schist belts and the transformation zones between the high- and low-velocity anomaly blocks are characterized by deeper earthquakes than the granitoid areas where shallow events dominate. Although only a few focal mechanisms were solved for southern Finland, there is a trend towards strike-slip and oblique strike-slip movements inside schist areas. The normal dip-slip type earthquakes are typical in the seismically active Kuusamo district in the NE edge of the SVEKALAPKO area, where the Archean crust is ca. 15-20 km thinner than the Proterozoic Svecofennian crust. Two near vertical dip-slip mechanism earthquakes occurred in the NE-SW junction between the Central Finland Granitoid Complex and the Vyborg rapakivi batholith, where high Vp/Vs-ratio deep-set intrusion splits the southern Finland schist belt into two parts in the tomography model.

Charge density based classification of intermolecular interactions in molecular crystals

Evaluation of intermolecular interactions in terms of both experimental and theoretical charge density analyses has produced a unified picture with which to classify strong and weak hydrogen bonds, along with van der Waals interactions, into three regions.

Synthesis, Structure, Spirocyclization Mechanism, and Glutathione Peroxidase-like Antioxidant Activity of Stable Spirodiazaselenurane and Spirodiazatellurane

The first examples of stable spirodiazaselenurane and spirodiazatellurane were synthesized by oxidative spirocyclization of the corresponding diaryl selenide and telluride and were structurally characterized. X-ray crystal structures of the spirodiazaselenurane and spirodiazatellurane suggest that the structures are distorted trigonal bipyramidal (TBP) with the electronegative nitrogen atoms occupying the apical positions and two carbon atoms and the lone pair of Se/Te occupying the equatorial positions. Interestingly, the spirodiazatellurane underwent spontaneous chiral resolution during crystallization, and the absolute configurations of its enantiomers were confirmed by single-crystal X-ray analyses. A detailed mechanistic study indicates that the cyclization to spirodiazaselenurane and spirodiazatellurane occurs via selenoxide and telluroxide intermediates. The chalcogenoxides cyclize to the corresponding spiro compounds in a stepwise manner via the involvement of hydroxyl chalcogenurane intermediates, and the activation energy for them spirocyclization reaction decreases in the order S > Se > Te. In addition to the synthesis, characterization, and mechanism of cyclization, the glutathione peroxidase (GPx) mimetic activity of the newly synthesized compounds was evaluated. These studies suggest that the tellurium compounds are more effective as GPx mimics than their selenium counterparts due to the fast oxidation of the tellurium center in the presence of peroxide and the involvement of an efficient redox cycle between the telluride and telluroxide intermediate.

Design of mechanism-based inhibitors of transthyretin amyloidosis: Studies with biphenyl ethers and new structural templates

Transthyretin (TTR), a tetrameric thyroxine (T4) carrier protein, is associated with a variety of amyloid diseases. In this study, we explore the potential of biphenyl ethers (BPE), which are shown to interact with a high affinity to its T4 binding site thereby preventing its aggregation and fibrillogenesis. They prevent fibrillogenesis by stabilizing the tetrameric ground state of transthyretin. Additionally, we identify two new structural templates (2-(5-mercapto-[1,3,4]oxadiazol-2-yl)-phenol and 2,3,6-trichloro-N-(4H-[1,2,4]triazol-3-yl) represented as compounds 11 and 12, respectively, throughout the manuscript) exhibiting the ability to arrest TTR amyloidosis. The dissociation constants for the binding of BPEs and compound 11 and 12 to TTR correlate with their efficacies of inhibiting amyloidosis. They also have the ability to inhibit the elongation of intermediate fibrils as well as show nearly complete (> 90%) disruption of the preformed fibrils. The present study thus establishes biphenyl ethers and compounds 11 and 12 as very potent inhibitors of TTR fibrillization and inducible cytotoxicity.

Solution structures of conformationally equilibrium forms of holo-acyl carrier protein (PfACP) from Plasmodium falciparum provides insight into the mechanism of activation of ACPs

Acyl Carrier Protein (ACP) from the malaria parasite, Plasmodium falciparum (PfACP) in its holo form is found to exist in two conformational states in solution. Unique 3D solution structures of holo-PfACP have been determined for both equilibrium conformations, using high-resolution NMR methods. Twenty high-resolution solution structures for each of the two forms of holo-PfACP have been determined on the basis of 1226 and 1218 unambiguously assigned NOEs (including NOEs between 4 '-phosphopantetheine prosthetic group (4 '-PP) and protein), 55 backbone dihedral angles and 26 hydrogen bonds. The atomic rmsd values of the determined structures of two equilibrium forms, about the mean coordinates of the backbone and heavy atoms, are 0.48 +/- 0.09 and 0.92 +/- 0.10 and 0.49 +/- 0.08 and 0.97 +/- 0.11 angstrom, respectively. The interaction of 4 '-PP with the polypeptide backbone is reported here for the first time for any of the ACPs. The structures of holo-PfACP consist of three well-defined helices that are tightly packed. The structured regions of the molecule are stabilized by extensive hydrophobic interactions. The difference between the two forms arises from a reorientation of the 4 '-PP group. The enthalpy difference between the two forms, although small, implies that a conformational switch is essential for the activation of holo-ACP. Sequence and structures of holo-PfACP have been compared with those of the ACPs from type I and type II fatty acid biosynthesis pathways (FAS), in particular with the ACP from rat and the butyryl-ACP from E. coli. The PfACP structure, thus determined has several novel features hitherto not seen in other ACPs.

Kinetics and the Mechanism of Interaction of the Endoplasmic Reticulum Chaperone, Calreticulin, with Monoglucosylated (Glc1Man9GlcNAc2) Substrate

Calreticulin is a lectin-like molecular chaperone of the endoplasmic reticulum in eukaryotes. Its interaction with N-glycosylated polypeptides is mediated by the glycan, Glc(1)Man(9)GlcNAc(2), present on the target glycoproteins. In this work, binding of monoglucosyl IgG (chicken) substrate to calreticulin has been studied using real time association kinetics of the interaction with the biosensor based on surface plasmon resonance (SPR). By SPR, accurate association and dissociation rate constants were determined, and these yielded a micromolar association constant. The nature of reaction was unaffected by immobilization of either of the reactants. The Scatchard analysis values for K-a agreed web crith the one obtained by the ratio k(1)/k(-1). The interaction was completely inhibited by free oligosaccharide, Glc(1)Man(9)GlcNAc(2), whereas Man(9)GlcNAc(2) did not bind to the calreticulin-substrate complex, attesting to the exquisite specificity of this interaction. The binding of calreticulin to IgG was used for the development of immunoassay and the relative affinity of the lectin-substrate association was indirectly measured. The values are in agreement with those obtained with SPR. Although the reactions are several orders of magnitude slower than the diffusion controlled processes, the data are qualitatively and quantitatively consistent with single-step bimolecular association and dissociation reaction. Analyses of the activation parameters indicate that reaction is enthalpically driven and does not involve a highly ordered transition state. Based on these data, the mechanism of its chaperone activity is briefly discussed.

Thermodynamic and kinetic studies on the mechanism of binding of methylumbelliferyl galactosides to the basic agglutinin from winged bean (Psophocarpus tetragonolobus)

The binding of winged bean basic agglutinin (WBA I) to 4-methylumbelliferyl (MeUmb) galactosides was examined by extrinsic fluorescence titration and stopped-flow spectrofluorimetry. Upon binding to WBA I, MeUmb alpha-galactosides show quenching in fluorescence intensity, decrease in UV absorbance with a concomitant blue shift, and decrease in fluorescence excited-state lifetimes. However, their beta-analogues show enhancement in fluorescence intensity, increase in UV absorbance with a red shift, and an increase in fluorescence excited-state lifetimes. This implies that the umbelliferyl groups of alpha- and beta-galactosides experience non-polar and polar microenvironments, respectively, upon binding to WBA I. Replacement of the anomeric hydroxyl group of galactose by 4-methylumbelliferyl moiety increases the affinity of resulting saccharides. Substitution of C-2 hydroxyl of galactose by an acetamido group leads to increased affinity due to a favorable entropy change. This suggests that acetamido group of MeUmb-alpha/beta-GalNAc binds to a relatively non-polar subsite of WBA I. Most interestingly, this substitution also reduces the association rate constants dramatically. Inspection of the activation parameters reveals that the enthalpy of activation is the limiting factor for the differences in the forward rate constants for these saccharides and the entropic contribution to the activation energy is small

Thermodynamic and kinetic studies on the mechanism of binding of methylumbelliferyl glycosides to jacalin

The binding of Artocarpus integrifolia lectin (jacalin) to 4-methylumbelliferyl (Meumb)-glycosides, Gal alpha Meumb, Gal beta Meumb, GalNAc alpha Meumb, GalNAc beta-Meumb, and Gal beta 3GalNAc beta Meumb was examined by extrinsic fluorescence quenching titration and stopped flow spectrofluorimetry. The binding was characterized by 100% quenching of fluorescence of Meumb-glycosides. Their association constants range from 2.0 x 10(4) to 1.58 x 10(6) M-1 at 15 degrees C. Entropic contribution is the major stabilizing force for avid binding of Meumb-glycosides indicating the existence of a hydrophobic site that is complementary to their methylumbelliferyl group. The second order association rate constants for interaction of these sugars with lectin at 15 degrees C vary from 8.8 x 10(5) to 3.24 x 10(6) M-1 S-1, at pH 7.2. The first order dissociation rate constants range from 2.30 to 43.0 S-1 at 15 degrees C. Despite the differences in their association rate constants, the overall values of association constants for these saccharides are determined by their dissociation rate constants. The second order rate constant for the association of Meumb-glycosides follows a pattern consistent with the magnitude of the activation energies involved therin. Activation parameters for association of all ligands illustrate that the origin of the barrier between binding of jacalin to Meumb-glycosides is entropic, and the enthalpic contribution is small. A correlation between these parameters and the structure of the ligands on the association rates underscores the importance of steric factors in determining protein saccharide recognitions.

Mechanism of lectin-cell interactions: thermodynamic and kinetic analysis of the binding of winged bean acidic lectin (WBA II) to human erythrocytes

In order to identify the forces involved in the binding and to understand the mechanism involved, equilibrium and kinetic studies were performed on the binding of the winged bean acidic lectin to human erythrocytes. The magnitudes of delta S and delta H were positive and negative respectively, an observation differing markedly from the lectin-simple sugar interactions where delta S and delta H are generally negative. Analysis of the sign and magnitudes of these values indicate that ionic and hydrogen bonded interactions prevail over hydrophobic interactions resulting in net -ve delta H (-37.12 kJ.mol-1) and +ve delta S (14.4 J.mole-1 K-1 at 20 degrees C), thereby suggesting that this entropy driven reaction also reflects conformational changes in the lectin and/or the receptor. Presence of two kinds of receptors for WBA II on erythrocytes, as observed by equilibrium studies, is consistent with the biexponential dissociation rate constants (at 20 degrees C K1 = 1.67 x 10(-3) M-1 sec-1 and K2 = 11.1 x 10(-3) M-1 sec-1). These two rate constants differed by an order of magnitude accounting for the difference in the association constants of the two receptors of WBA II. However, the association process remains monoexponential suggesting no observable difference in the association rates of the lectin molecule with both the receptors, under the experimental conditions studied. The thermodynamic parameters calculated from kinetic data correlate well with those observed by equilibrium. A two-step binding mechanism is proposed based on the kinetic parameters for WBA II-receptor interaction

Analysis of dynamics and mechanism of ligand binding to Artocarpus integrifolia agglutinin. A 13C and 19F NMR study

Binding of 13C-labeled N-acetylgalactosamine (13C-GalNAc) and N-trifluoroacetylgalactosamine (19F-GalNAc) to Artocarpus integrifolia agglutinin has been studied using 13C and 19F nuclear magnetic resonance spectroscopy, respectively. Binding of these saccharides resulted in broadening of the resonances, and no change in chemical shift was observed, suggesting that the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc experience a magnetically equivalent environment in the lectin combining site. The alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc were found to be in slow exchange between free and protein bound states. Binding of 13C-GalNAc was studied as a function of temperature. From the temperature dependence of the line broadening, the thermodynamic and kinetic parameters were evaluated. The association rate constants obtained for the alpha-anomers of 13C-GalNAc and 19F-GalNAc (k+1 = 1.01 x 10(5) M-1.s-1 and 0.698 x 10(5) M-1.s-1, respectively) are in close agreement with those obtained for the corresponding beta-anomers (k+1 = 0.95 x 10(5) M-1.s-1 and 0.65 x 10(5) M-1.s-1, respectively), suggesting that the two anomers bind to the lectin by a similar mechanism. In addition these values are several orders of magnitude slower than those obtained for diffusion controlled processes. The dissociation rate constants obtained are 49.9, 56.9, 42, and 43 s-1, respectively, for the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc. A two-step mechanism has been proposed for the interaction of 13C-GalNAc and 19F-GalNAc with A. integrifolia lectin in view of the slow association rates and high activation entropies. The thermodynamic parameters obtained for the association and dissociation reactions suggest that the binding process is entropically favored and that there is a small enthalpic contribution.