Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Acyclic permutants of naturally occurring cyclic proteins - Characterization of cystine knot and beta-sheet formation in the macrocyclic polypeptide kalata B1

Kalata B1 is a prototypic member of the unique cyclotide family of macrocyclic polypeptides in which the major structural features are a circular peptide backbone, a triple stranded beta-sheet, and a cystine knot arrangement of three disulfide bonds. The cyclotides are the only naturally occurring family of circular proteins and have prompted us to explore the concept of acyclic permutation, i.e. opening the backbone of a cross-linked circular protein in topologically permuted ways. We have synthesized the complete suite of acyclic permutants of kalata B1 and examined the effect of acyclic permutation on structure and activity. Only two of six topologically distinct backbone loops are critical for folding into the native conformation, and these involve disruption of the embedded ring in the cystine knot. Surprisingly, it is possible to disrupt regions of the p-sheet and still allow folding into native-like structure, provided the cystine knot is intact. Kalata B1 has mild hemolytic activity, but despite the overall structure of the native peptide being retained in all but two cases, none of the acyclic permutants displayed hemolytic activity. This loss of activity is not localized to one particular region and suggests that cyclization is critical for hemolytic activity.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 1. Isolation, characterization and chemical synthesis

A novel conotoxin belonging to the 'four-loop' structural class has been isolated from the venom of the piscivorous cone snail Conus tulipa. It was identified using a chemical-directed strategy based largely on mass spectrometric techniques. The new toxin, conotoxin TVIIA, consists of 30 amino-acid residues and contains three disulfide bonds. The amino-acid sequence was determined by Edman analysis as SCSGRDSRCOOVCCMGLMCSRGKCVSIYGE where O = 4-transl-hydroxyproline. Two under-hydroxylated analogues, [Pro10]TVIIA and [Pro10,11]TVIIA, were also identified in the venom of C. tulipa. The sequences of TVIIA and [Pro10]TVIIA were further verified by chemical synthesis and coelution studies with native material. Conotoxin TVIIA has a six cysteine/four-loop structural framework common to many peptides from Conus venoms including the omega-, delta- and kappa-conotoxins. However, TVIIA displays little sequence homology with these well-characterized pharmacological classes of peptides, but displays striking sequence homology with conotoxin GS, a peptide from Conus geographus that blocks skeletal muscle sodium channels. These new toxins and GS share several biochemical features and represent a distinct subgroup of the four-loop conotoxins.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

Integrable Kondo impurities in one-dimensional extended Hubbard models

Three kinds of integrable Kondo problems in one-dimensional extended Hubbard models are studied by means of the boundary graded quantum inverse scattering method. The boundary K matrices depending on the local moments of the impurities are presented as a nontrivial realization of the graded reflection equation algebras acting in a (2s alpha + 1)-dimensional impurity Hilbert space. Furthermore, these models are solved using the algebraic Bethe ansatz method, and the Bethe ansatz equations are obtained.

Solution structures in aqueous SDS micelles of two amyloid beta peptides of A beta(1-28) mutated at the alpha-secretase cleavage site (K16E, K16F)

NMR solution structures are reported for two mutants (K16E, K16F) of the soluble amyloid beta peptide A beta(1-28). The structural effects of these mutations of a positively charged residue to anionic and hydrophobic residues at the alpha-secretase cleavage site (Lys16-Leu17) were examined in the membrane-simulating solvent aqueous SDS micelles. Overall the three-dimensional structures were similar to that for the native A beta(1-28) sequence in that they contained an unstructured N-terminus and a helical C-terminus. These structural elements are similar to those seen in the corresponding regions of full-length A beta peptides A beta(1-40) and A beta(1-42), showing that the shorter peptides are valid model systems. The K16E mutation, which might be expected to stabilize the macrodipole of the helix, slightly increased the helix length (residues 13-24) relative to the K16F mutation, which shortened the helix to between residues 16 and 24. The observed sequence-dependent control over conformation in this region provides an insight into possible conformational switching roles of mutations in the amyloid precursor protein from which A beta peptides are derived. In addition, if conformational transitions from helix to random coil to sheet precede aggregation of A beta peptides in vivo, as they do in vitro, the conformation-inducing effects of mutations at Lys16 may also influence aggregation and fibril formation. (C) 2000 Academic Press.

Three-dimensional structure of RK-1: A novel alpha-defensin peptide

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of RK-1, an antimicrobial peptide from rabbit kidney recently discovered from homology screening based on the distinctive physicochemical properties of the corticostatins/defensins. RK-1 consists of 32 residues, including six cysteines arranged into three disulfide bonds. It exhibits antimicrobial activity against Escherichia coli and activates Ca2+ channels in vitro. Through its physicochemical similarity, identical cysteine spacing, and linkage to the corticostatins/defensins, it was presumed to be a member of this family. However, RK-1 lacks both a large number of arginines in the primary sequence and a high overall positive charge, which are characteristic of this family of peptides. The three-dimensional solution structure, determined by NMR, consists of a triple-stranded antiparallel beta -sheet and a series of turns and is similar to the known structures of other alpha -defensins. This has enabled the definitive classification of RK-1 as a member of this family of antimicrobial peptides. Ultracentrifuge measurements confirmed that like rabbit neutrophil defensins, RK-1 is monomeric in solution, in contrast to human neutrophil defensins, which are dimeric.

Efficient discovery of immune response targets by cyclical refinement of QSAR models of peptide binding

Peptides that induce and recall T-cell responses are called T-cell epitopes. T-cell epitopes may be useful in a subunit vaccine against malaria. Computer models that simulate peptide binding to MHC are useful for selecting candidate T-cell epitopes since they minimize the number of experiments required for their identification. We applied a combination of computational and immunological strategies to select candidate T-cell epitopes. A total of 86 experimental binding assays were performed in three rounds of identification of HLA-All binding peptides from the six preerythrocytic malaria antigens. Thirty-six peptides were experimentally confirmed as binders. We show that the cyclical refinement of the ANN models results in a significant improvement of the efficiency of identifying potential T-cell epitopes. (C) 2001 by Elsevier Science Inc.

Cyclic octapeptides containing thiazole. Effect of stereochemistry and degree of flexibility on calcium binding properties

Solution conformation and calcium binding properties have been investigated for the two cyclic octapeptides cyclo(-D-Thr-D-Val(Thz)-Ile-)(2) (4) and cyclo(-Thr-Gly(Thz)-Ile-Ser-Gly(Thz)-Ile-)(5) and the results are compared to those for the cyclic octapeptides previously studied; ascidiacyclamide (1), patellamide D (2), cyclo(-Thr-D-Val(Thz)-Ile-)(2) (3), and cyclo(-Thr-D-Val-alphaAbu-Ile-)2 (6). Both 4 and 5 contain two heterocyclic thiazole ring constraints but the latter has a larger degree of flexibility as a consequence of the glycine residues within the cyclic framework. The solution conformation of 4 and 5 was determined from H-1 NMR spectra and found to be a twisted figure of eight similar to that for 2. Complexation studies using H-1 NMR and CD spectroscopy yielded 1 : 1 calcium-peptide binding constants (logK) for the two peptides (2.3 (4) and 5.7 (5)). For 5 the magnitude of the binding constant was verified by a competition titration using CD. The different calcium-binding affinities of 3 (logK = 4.0) and 4 is attributed to the stereochemistry of the threonine residue. The magnitude of the binding constant for 5 compared to 3 and 4 (all peptides containing two thiazole ring constrains) demonstrates that the increase in flexibility of the cyclic peptide has a dramatic effect on the Ca2+ binding ability. The affinity for Ca2+ thus decreases in the order (6 similar to 5 > 3 > 2 similar to 1 > 4). The number of carbonyl donors available on each peptide has only a limited effect on calcium binding. The most important factor is the flexibility, which allows for a conformation of the peptide capable of binding calcium efficiently.

Methods for prediction of peptide binding to MHC molecules: A comparative study

Background: A variety of methods for prediction of peptide binding to major histocompatibility complex (MHC) have been proposed. These methods are based on binding motifs, binding matrices, hidden Markov models (HMM), or artificial neural networks (ANN). There has been little prior work on the comparative analysis of these methods. Materials and Methods: We performed a comparison of the performance of six methods applied to the prediction of two human MHC class I molecules, including binding matrices and motifs, ANNs, and HMMs. Results: The selection of the optimal prediction method depends on the amount of available data (the number of peptides of known binding affinity to the MHC molecule of interest), the biases in the data set and the intended purpose of the prediction (screening of a single protein versus mass screening). When little or no peptide data are available, binding motifs are the most useful alternative to random guessing or use of a complete overlapping set of peptides for selection of candidate binders. As the number of known peptide binders increases, binding matrices and HMM become more useful predictors. ANN and HMM are the predictive methods of choice for MHC alleles with more than 100 known binding peptides. Conclusion: The ability of bioinformatic methods to reliably predict MHC binding peptides, and thereby potential T-cell epitopes, has major implications for clinical immunology, particularly in the area of vaccine design.

D-amino acid residue in the C-type natriuretic peptide from the venom of the mammal, Ornithorhynchus anatinus, the Australian platypus

The C-type natriuretic peptide from the platypus venom (OvCNP) exists in two forms, OvCNPa and OvCNPb, whose amino acid sequences are identical. Through the use of nuclear magnetic resonance, mass spectrometry, and peptidase digestion studies, we discovered that OvCNPb incorporates a D-amino acid at position 2 in the primary structure. Peptides containing a D-amino acid have been found in lower forms of organism, but this report is the first for a D-amino acid in a biologically active peptide from a mammal. The result implies the existence of a specific isomerase in the platypus that converts an L-amino acid residue in the protein to the D-configuration. (C) 2002 Federation of European Biochemical Societies. Published by Elsevier Science B.V. All rights reserved.

Algebraic Bethe ansatz for integrable one-dimensional extended Hubbard models with open boundary conditions

Nine classes of integrable open boundary conditions, further extending the one-dimensional U-q (gl (212)) extended Hubbard model, have been constructed previously by means of the boundary Z(2)-graded quantum inverse scattering method. The boundary systems are now solved by using the algebraic Bethe ansatz method, and the Bethe ansatz equations are obtained for all nine cases.

Quantification of prolactin-releasing peptide (PrRP) mRNA expression in specific brain regions of the rat during the oestrous cycle and in lactation

Real-time Taqman(TM) RT-PCR was used to make quantitative comparisons of the levels of PrRP mRNA expression in micropunch brain samples from rats at different stages of the oestrous cycle and in lactation. The nucleus of the solitary tract and ventrolateral reticular nuclei of the medulla oblongata contained significantly (P < 0.05) greater levels of PrRP mRNA than any hypothalamic region. Within the hypothalamus, the highest level of PrRP expression was localised to the dorsomedial aspect of the ventromedial hypothalamus. All other hypothalamic regions exhibited significantly (P < 0.05) lower levels of expression, including the rostral and caudal dorsomedial hypothalamus. Very low levels of PrRP expression were observed in the arcuate nucleus, paraventricular nucleus, medial preoptic nucleus and ventrolateral aspect of the ventromedial hypothalamus. No significant changes in PrRP expression were noted in any sampled region between proestrus, oestrus or dioestrus. Similarly, PrRP expression in hypothalamic regions did not differ between lactating and non-lactating (dioestrous) animals. During validation of RT-PCR techniques we cloned and sequenced a novel splice variant of PrRP from the hypothalamus. This variant arises from alternative splicing of the donor site within exon 2, resulting in an insert of 64 base pairs and shift in the-codon:reading frame with the introduction of an early stop codon. In the hypothalamus and brainstem, mRNA expression of the variant was restricted to regions that expressed PrRP. These results suggest that PrRP expression in the hypothalamus may be more Widespread than previously reported. However, the relatively low level of PrRP in the hypothalamus and the lack of significant changes in expression during the oestrous cycle and lactation provides further evidence that PrRP is unlikely to be involved in the regulation of prolactin, secretion. (C) 2003 Elsevier Science B.V. All rights reserved.