Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

Targeted drug delivery to the skin and deeper tissues: Role of physiology, solute structure and disease

1. Drug delivery through the skin has been used to target the epidermis, dermis and deeper tissues and for systemic delivery, The major barrier for the transport of drugs through the skin is the stratum corneum, with most transport occurring through the intercellular region, The polarity of the intercellular region appears to be similar to butanol, with the diffusion of solutes being hindered by saturable hydrogen bonding to the polar head groups of the ceramides, fatty acids and other intercellular lipids, Accordingly, the permeability of the more lipophilic solutes is greatest from aqueous solutions, whereas polar solute permeability is favoured by hydrocarbon-based vehicles. 2. The skin is capable of metabolizing many substances and, through its microvasculature, limits the transport of most substances into regions below the dermis. 3. Although the flux of solutes through the skin should be identical for different vehicles when the solute exists as a saturated solution, the fluxes vary in accordance with the skin penetration enhancement properties of the vehicle. It is therefore desirable that the regulatory standards required for the bioequivalence of topical products include skin studies. 4. Deep tissue penetration can be related to solute protein binding, solute molecular size and dermal blood flow. 5. Iontophoresis is a promising area of skin drug delivery, especially for ionized solutes and when a rapid effect is required. 6. In general, psoriasis and other skin diseases facilitate drug delivery through the skin. 7. It is concluded that the variability in skin permeability remains an obstacle in optimizing drug delivery by this route.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear H-1-N-15 NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta -strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. N-15 longitudinal and transverse relaxation rates, and {H-1}-N-15 heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85 +/- 0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides. (C) 2001 Academic Press.

The leucine-rich repeat as a protein recognition motif

Leucine-rich repeats (LRRs) are 20-29-residue sequence motifs present in a number of proteins with diverse functions. The primary function of these motifs appears to be to provide a versatile structural framework for the formation of protein-protein interactions. The past two years have seen an explosion of new structural information on proteins with LRRs. The new structures represent different LRR subfamilies and proteins with diverse functions, including GTPase-activating protein rna 1 p from the ribonuclease-inhibitor-like subfamily; spliceosomal protein U2A', Rab geranylgeranyltransferase, internalin B, dynein light chain 1 and nuclear export protein TAP from the SDS22-like subfamily; Skp2 from the cysteine-containing subfamily; and YopM from the bacterial subfamily. The new structural information has increased our understanding of the structural determinants of LRR proteins and our ability to model such proteins with unknown structures, and has shed new light on how these proteins participate in protein-protein interactions.

Analysis of thermodynamic non-ideality in terms of protein solvation

The effects of thermodynamic non-ideality on the forms of sedimentation equilibrium distributions for several isoelectric proteins have been analysed on the statistical-mechanical basis of excluded volume to obtain an estimate of the extent of protein solvation. Values of the effective solvation. parameter delta are reported for ellipsoidal as well as spherical models of the proteins, taken to be rigid, impenetrable macromolecular structures. The dependence of the effective solvated radius upon protein molecular mass exhibits reasonable agreement with the relationship calculated for a model in which the unsolvated protein molecule is surrounded by a 0.52-nm solvation shell. Although the observation that this shell thickness corresponds to a double layer of water molecules may be of questionable relevance to mechanistic interpretation of protein hydration, it augurs well for the assignment of magnitudes to the second virial coefficients of putative complexes in the quantitative characterization of protein-protein interactions under conditions where effects of thermodynamic non-ideality cannot justifiably be neglected. (C) 2001 Elsevier Science B.V. All rights reserved.

Determination of the intramolecular disulfide bond arrangement and biochemical identification of the glycosylation sites of the nonstructural protein NS1 of Murray Valley encephalitis virus

The 12 cysteine residues in the flavivirus NS1 protein are strictly conserved, suggesting that they form disulfide bonds that are critical for folding the protein into a functional structure. In this study, we examined the intramolecular disulfide bond arrangement of NS1 of Murray Valley encephalitis virus and elucidated three of the six cysteine-pairing arrangements. Disulfide linkages were identified by separating tryptic-digested NS1 by reverse-phase high pressure liquid chromatography and analysing the resulting peptide peaks by protein sequencing, amino acid analysis and/or electrospray mass spectrometry. The pairing arrangements between the six amino-terminal cysteines were identified as follows: Cys(4)-Cys(15), Cys(55)-Cys(143) and Cys(179)-Cys(223). Although the pairing arrangements between the six carboxyterminal cysteines were not determined, we were able to eliminate several cysteine-pairing combinations. Furthermore, we demonstrated that all three putative N-linked glycosylation sites of NS1 are utilized and that the Asn(207) glycosylation site contains a mannose-rich glycan.

Genomic organization of the CC chemokine MIP-3 alpha/CCL20/LARC/EXODUS/SCYA20, showing gene structure, splice variants, and chromosome localization

We describe the genomic organization of a recently identified CC chemokine, MIP3 alpha /CCL20 (HGMW-approved symbol SCYA20). The MIP-3 alpha /CCL20 gene was cloned and sequenced, revealing a four exon, three intron structure, and was localized by FISK analysis to 2q35-q36. Two distinct cDNAs were identified, encoding two forms of MIP-3 alpha /CCL20, Ala MLP-3 alpha /CCL20 and Ser MIP-3 alpha /CCL20, that differ by one amino acid at the predicted signal peptide cleavage site. Examination of the sequence around the boundary of intron 1 and exon 2 showed that use of alternative splice acceptor sites could give rise to Ata MIP-3 alpha /CCL20 or Ser MIP-3 alpha /CCL20. Both forms of MIP-3cr/CCL20 were chemically synthesized and tested for biological activity. Both flu antigen plus IL-a-activated CD4(+) and CD8(+) T lymphoblasts and cord blood-derived dendritic cells responded to Ser and Ala MIP-3 alpha /CCL20. T lymphocytes exposed only to IL-2 responded inconsistently, while no response was detected in naive T lymphocytes, monocytes, or neutrophils. The biological activity of Ser MIP-3 alpha /CCL20 and Ala MIP-3 alpha /CCL20 and the tissue-specific preference of different splice acceptor sites are not yet known. (C) 2001 Academic Press.

Solution structure of BSTI: A new trypsin inhibitor from skin secretions of Bombina bombina

The three-dimensional solution structure of BSTI, a trypsin inhibitor from the European frog Bombina bombina, has been solved using H-1 NMR spectroscopy. The 60 amino acid protein contains five disulfide bonds, which were unambiguously determined to be Cvs (4-38), Cys (13-34), Cys (17-30), Cys (21-60), and Cys (40-54) by experimental restraints and subsequent structure calculations. The main elements of secondary structure are four beta -strands, arranged as two small antiparallel beta -sheets, The overall fold of BSTI is disk shaped and is characterized by the lack of a hydrophobic core. The presumed active site is located on a loop comprising residues 21-34, which is a relatively disordered region similar to that seen in many other protease inhibitors. However, the overall fold is different to other known protease inhibitors with the exception of a small family of inhibitors isolated from nematodes of the family Ascaris and recently also from the haemolymph of Apis mellifera. BSTI may thus be classified as a new member of this recently discovered family of protease inhibitors.

Cloning of a catalytic subunit of cAMP-dependent protein kinase from the honeybee (Apis mellifera) and its localization in the brain

In the honeybee the cAMP-dependent signal transduction cascade has been implicated in processes underlying learning and memory, The cAMP-dependent protein kinase (PKA) is the major mediator of cAMP action. To characterize the PKA system in the honeybee brain we cloned a homologue of a PKA catalytic subunit from the honeybee,The deduced amino acid sequence shows 80-94% identity with catalytic subunits of PKA from Drosophila melanogaster, Aplysia californica and mammals. The corresponding gene is predominantly expressed in the mushroom bodies, a structure that is involved in learning and memory processes. However, expression can also be found in the antennal and optic lobes,The level of expression varies within all three neuropiles.

Synthesis of an analog of the thyroid hormone-binding protein transthyretin via regioselective chemical ligation

Transthyretin is an essential protein responsible for the transport of thyroid hormones and retinol in human serum and is also implicated in the amyloid diseases familial amyloidotic polyneuropathy and senile systemic amyloidosis. Its folding properties and stabilization by ligands are of current interest due to their importance in understanding and combating these diseases, Here we report the solid phase synthesis of the monomeric unit of a transthyretin analog (equivalent to 127 amino acids) using t-Boc chemistry and peptide ligation and its folding to form a functional 54-kDa tetramer, The monomeric unit of the protein was chemically synthesized in three parts (positions 1-51, 54-99, and 102-127) and ligated using a chemoselective thioether ligation chemistry. The synthetic protein was folded and assembled to a tetrameric structure in the presence of transthyretin's native ligand, thyroxine, as shown by gel filtration chromatography, native gel electrophoresis, transthyretin antibody recognition, and thyroid hormone binding. Other folding products included a high molecular weight aggregate as well as a transient dimeric species. This represents one of the largest macromolecules chemically synthesized to date and demonstrates the potential of protein chemical synthesis for investigations of protein-ligand interactions.

Human pigmentation genes: identification, structure and consequences of polymorphic variation

The synthesis of the visible pigment melanin by the melanocyte cell is the basis of the human pigmentary system, those genes directing the formation, transport and distribution of the specialised melanosome organelle in which melanin accumulates can legitimately be called pigmentation genes. The genes involved in this process have been identified through comparative genomic studies of mouse coat colour mutations and by the molecular characterisation of human hypopigmentary genetic diseases such as OCA1 and OCA2. The melanocyte responds to the peptide hormones a-MSH or ACTH through the MC1R G-protein coupled receptor to stimulate melanin production through induced maturation or switching of melanin type. The pheomelanosome, containing the key enzyme of the pathway tyrosinase, produces light red/yellowish melanin, whereas the eumelanosome produces darker melanins via induction of additional TYRP1, TYRP2, SILV enzymes, and the P-protein. Intramelanosomal pH governed by the P-protein may act as a critical determinant of tyrosinase enzyme activity to control the initial step in melanin synthesis or TYRP complex formation to facilitate melanogenesis and melanosomal maturation. The search for genetic variation in these candidate human pigmentation genes in various human populations has revealed high levels of polymorphism in the MC1R locus, with over 30 variant alleles so far identified. Functional correlation of MC1R alleles with skin and hair colour provides evidence that this receptor molecule is a principle component underlying normal human pigment variation. (C) 2001 Elsevier Science B.V. All rights reserved.

Getting the adrenaline going: Crystal structure of the adrenaline-synthesizing enzyme PNMT

Background: Adrenaline is localized to specific regions of the central nervous system (CNS), but its role therein is unclear because of a lack of suitable pharmacologic agents. Ideally, a chemical is required that crosses the blood-brain barrier, potently inhibits the adrenaline-synthesizing enzyme PNMT, and does not affect other catecholamine processes. Currently available PNMT inhibitors do not meet these criteria. We aim to produce potent, selective, and CNS-active PNMT inhibitors by structure-based design methods. The first step is the structure determination of PNMT. Results: We have solved the crystal structure of human PNMT complexed with a cofactor product and a submicromolar inhibitor at a resolution of 2.4 Angstrom. The structure reveals a highly decorated methyltransferase fold, with an active site protected from solvent by an extensive cover formed from several discrete structural motifs. The structure of PNMT shows that the inhibitor interacts with the enzyme in a different mode from the (modeled) substrate noradrenaline. Specifically, the position and orientation of the amines is not equivalent. Conclusions: An unexpected finding is that the structure of PNMT provides independent evidence of both backward evolution and fold recruitment in the evolution of a complex enzyme from a simple fold. The proposed evolutionary pathway implies that adrenaline, the product of PNMT catalysis, is a relative newcomer in the catecholamine family. The PNMT structure reported here enables the design of potent and selective inhibitors with which to characterize the role of adrenaline in the CNS. Such chemical probes could potentially be useful as novel therapeutics.

Isolation of a novel G-protein gamma-subunit from Arabidopsis thaliana and its interaction with G beta

There is increasing evidence that heterotrimeric G-proteins (G-proteins) are involved in many plant processes including phytohormone response, pathogen defence and stomatal control. In animal systems, each of the three G-protein subunits belong to large multigene families; however, few subunits have been isolated from plants. Here we report the cloning of a second plant G-protein γ-subunit (AGG2) from Arabidopsis thaliana. The predicted AGG2 protein sequence shows 48% identity to the first identified Arabidopsis Gγ-subunit, AGG1. Furthermore, AGG2 contains all of the conserved characteristics of γ-subunits including a small size (100 amino acids, 11.1 kDa), C-terminal CAAX box and a N-terminal α-helix region capable of forming a coiled-coil interaction with the β-subunit. A strong interaction between AGG2 and both the tobacco (TGB1) and Arabidopsis (AGB1) β-subunits was observed in vivo using the yeast two-hybrid system. The strong association between AGG2 and AGB1 was confirmed in vitro. Southern and Northern analyses showed that AGG2 is a single copy gene in Arabidopsis producing two transcripts that are present in all tissues tested. The isolation of a second γ-subunit from A. thaliana indicates that plant G-proteins, like their mammalian counterparts, may form different heterotrimer combinations that presumably regulate multiple signal transduction pathways.

Characterization of a chemsensory protein (ASP3c) from honeybee (Apis mellifera L.) as a brood pheromone carrier

Chemosensory proteins (CSPs) are ubiquitous soluble small proteins isolated from sensory organs of a wide range of insect species, which are believed to be involved in chemical communication. We report the cloning of a honeybee CSP gene called ASP3c, as well as the structural and functional characterization of the encoded protein. The protein was heterologously secreted by the yeast Pichia pastoris using the native signal peptide. ASP3c disulfide bonds were assigned after trypsinolysis followed by chromatography and mass spectrometry combined with microsequencing. The pairing (Cys(I)-Cys(II), Cys(III)-Cys(IV)) was found to be identical to that of Schistocerca gregaria CSPs, suggesting that this pattern occurs commonly throughout the insect CSPs. CD measurements revealed that ASP3c mainly consists of alpha-helices, like other insect CSPs. Gel filtration analysis showed that ASP3c is monomeric at neutral pH. Using ASA, a fluorescent fatty acid anthroyloxy analogue as a probe, ASP3c was shown to bind specifically to large fatty acids and ester derivatives, which are brood pheromone components, in the micromolar range. It was unable to bind tested general odorants and other tested pheromones (sexual and nonsexual). This is the first report on a natural pheromonal ligand bound by a recombinant CSP with a measured affinity constant.

Copper transfer from the Cu(I) chaperone, CopZ, to the repressor, Zn(II)CopY: Metal coordination environments and protein interactions

Extracellular copper regulates the DNA binding activity of the CopY repressor of Enterococcus hirae and thereby controls expression of the copper homeostatic genes encoded by the cop operon. CopY has a CxCxxxxCxC metal binding motif. CopZ, a copper chaperone belonging to a family of metallochaperones characterized by a MxCxxC metal binding motif, transfers copper to CopY. The copper binding stoichiometries of CopZ and CopY were determined by in vitro metal reconstitutions. The stoichiometries were found to be one copper(l) per CopZ and two copper(l) per CopY monomer. X-ray absorption studies suggested a mixture of two- and three-coordinate copper in Cu(1)CopZ, but a purely three-coordinate copper coordination with a Cu-Cu interaction for Cu(1)(2)CopY. The latter coordination is consistent with the formation of a compact binuclear Cu(l)-thiolate core in the CxCxxxxCxC binding motif of CopY. Displacement of zinc, by copper. from CopY was monitored with 2,4-pyridylazoresorcinol. Two copper(l) ions were required to release the single zinc(II) ion bound per CopY monomer. The specificity of copper transfer between CopZ and CopY was dependent on electrostatic interactions. Relative copper binding affinities of the proteins were investigated using the chelator, diethyldithiocarbamic acid (DDC). These data suggest that CopY has a higher affinity for copper than CopZ. However, this affinity difference is not the sole factor in the copper exchange: a charge-based interaction between the two proteins is required for the transfer reaction to proceed. Gain-of-function mutation of a CopZ homologue demonstrated the necessity of four lysine residues on the chaperone for the interaction with CopY. Taken together, these results suggest a mechanism for copper exchange between CopZ and CopY.