Choroid plexus epithelial expression of MDR1 P glycoprotein and multidrug resistance-associated protein contribute to the blood–cerebrospinal-fluid drug-permeability barrier

The blood–brain barrier and a blood–cerebrospinal-fluid (CSF) barrier function together to isolate the brain from circulating drugs, toxins, and xenobiotics. The blood–CSF drug-permeability barrier is localized to the epithelium of the choroid plexus (CP). However, the molecular mechanisms regulating drug permeability across the CP epithelium are defined poorly. Herein, we describe a drug-permeability barrier in human and rodent CP mediated by epithelial-specific expression of the MDR1 (multidrug resistance) P glycoprotein (Pgp) and the multidrug resistance-associated protein (MRP). Noninvasive single-photon-emission computed tomography with 99mTc-sestamibi, a membrane-permeant radiopharmaceutical whose transport is mediated by both Pgp and MRP, shows a large blood-to-CSF concentration gradient across intact CP epithelium in humans in vivo. In rats, pharmacokinetic analysis with 99mTc-sestamibi determined the concentration gradient to be greater than 100-fold. In membrane fractions of isolated native CP from rat, mouse, and human, the 170-kDa Pgp and 190-kDa MRP are identified readily. Furthermore, the murine proteins are absent in CP isolated from their respective mdr1a/1b(−/−) and mrp(−/−) gene knockout littermates. As determined by immunohistochemical and drug-transport analysis of native CP and polarized epithelial cell cultures derived from neonatal rat CP, Pgp localizes subapically, conferring an apical-to-basal transepithelial permeation barrier to radiolabeled drugs. Conversely, MRP localizes basolaterally, conferring an opposing basal-to-apical drug-permeation barrier. Together, these transporters may coordinate secretion and reabsorption of natural product substrates and therapeutic drugs, including chemotherapeutic agents, antipsychotics, and HIV protease inhibitors, into and out of the central nervous system.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen–deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen–deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is Pstruc, which is more related to the structure of the intermediate, and the other is Papp, which is more related to the stability of the intermediate. Kinetic hydrogen–deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

The helical domain of intestinal fatty acid binding protein is critical for collisional transfer of fatty acids to phospholipid membranes

Fatty acid binding proteins (FABPs) exhibit a β-barrel topology, comprising 10 antiparallel β-sheets capped by two short α-helical segments. Previous studies suggested that fatty acid transfer from several FABPs occurs during interaction between the protein and the acceptor membrane, and that the helical domain of the FABPs plays an important role in this process. In this study, we employed a helix-less variant of intestinal FABP (IFABP-HL) and examined the rate and mechanism of transfer of fluorescent anthroyloxy fatty acids (AOFA) from this protein to model membranes in comparison to the wild type (wIFABP). In marked contrast to wIFABP, IFABP-HL does not show significant modification of the AOFA transfer rate as a function of either the concentration or the composition of the acceptor membranes. These results suggest that the transfer of fatty acids from IFABP-HL occurs by an aqueous diffusion-mediated process, i.e., in the absence of the helical domain, effective collisional transfer of fatty acids to membranes does not occur. Binding of wIFABP and IFABP-HL to membranes was directly analyzed by using a cytochrome c competition assay, and it was shown that IFABP-HL was 80% less efficient in preventing cytochrome c from binding to membranes than the native IFABP. Collectively, these results indicate that the α-helical region of IFABP is involved in membrane interactions and thus plays a critical role in the collisional mechanism of fatty acid transfer from IFABP to phospholipid membranes.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Probing the role of packing specificity in protein design

By using a protein-design algorithm that quantitatively considers side-chain packing, the effect of specific steric constraints on protein design was assessed in the core of the streptococcal protein G β1 domain. The strength of packing constraints used in the design was varied, resulting in core sequences that reflected differing amounts of packing specificity. The structural flexibility and stability of several of the designed proteins were experimentally determined and showed a trend from well-ordered to highly mobile structures as the degree of packing specificity in the design decreased. This trend both demonstrates that the inclusion of specific packing interactions is necessary for the design of native-like proteins and defines a useful range of packing specificity for the design algorithm. In addition, an analysis of the modeled protein structures suggested that penalizing for exposed hydrophobic surface area can improve design performance.

CooA, a CO-sensing transcription factor from Rhodospirillum rubrum, is a CO-binding heme protein

Biological sensing of small molecules such as NO, O2, and CO is an important area of research; however, little is know about how CO is sensed biologically. The photosynthetic bacterium Rhodospirillum rubrum responds to CO by activating transcription of two operons that encode a CO-oxidizing system. A protein, CooA, has been identified as necessary for this response. CooA is a member of a family of transcriptional regulators similar to the cAMP receptor protein and fumavate nitrate reduction from Escherichia coli. In this study we report the purification of wild-type CooA from its native organism, R. rubrum, to greater than 95% purity. The purified protein is active in sequence-specific DNA binding in the presence of CO, but not in the absence of CO. Gel filtration experiments reveal the protein to be a dimer in the absence of CO. Purified CooA contains 1.6 mol heme per mol of dimer. Upon interacting with CO, the electronic spectrum of CooA is perturbed, indicating the direct binding of CO to the heme of CooA. A hypothesis for the mechanism of the protein’s response to CO is proposed.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

Enhanced immunogenicity of HIV-1 vaccine construct by modification of the native peptide sequence

Viral proteins are not naturally selected for high affinity major histocompatibility complex (MHC) binding sequences; indeed, if there is any selection, it is likely to be negative in nature. Thus, one should be able to increase viral peptide binding to MHC in the rational design of synthetic peptide vaccines. The T1 helper peptide from the HIV-1 envelope protein was made more immunogenic for inducing T cell proliferation to the native sequence by replacing a residue that exerts an adverse influence on peptide binding to an MHC class II molecule. Mice immunized with vaccine constructs combining the more potent Th helper (Th) epitope with a cytotoxic T lymphocyte (CTL) determinant developed greatly enhanced CTL responses. Use of class II MHC-congenic mice confirmed that the enhancement of CTL response was due to class II-restricted help. Thus, enhanced T cell help is key for optimal induction of CTL, and, by modification of the native immunogen to increase binding to MHC, it is possible to develop second generation vaccine constructs that enhance both Th cell activation and CTL induction.

Preferential exclusion of sucrose from recombinant interleukin-1 receptor antagonist: Role in restricted conformational mobility and compaction of native state

Understanding the mechanism for sucrose-induced protein stabilization is important in many diverse fields, ranging from biochemistry and environmental physiology to pharmaceutical science. Timasheff and Lee [Lee, J. C. & Timasheff, S. N. (1981) J. Biol. Chem. 256, 7193–7201] have established that thermodynamic stabilization of proteins by sucrose is due to preferential exclusion of the sugar from the protein’s surface, which increases protein chemical potential. The current study measures the preferential exclusion of 1 M sucrose from a protein drug, recombinant interleukin 1 receptor antagonist (rhIL-1ra). It is proposed that the degree of preferential exclusion and increase in chemical potential are directly proportional to the protein surface area and that, hence, the system will favor the protein state with the smallest surface area. This mechanism explains the observed sucrose-induced restriction of rhIL-1ra conformational fluctuations, which were studied by hydrogen–deuterium exchange and cysteine reactivity measurements. Furthermore, infrared spectroscopy of rhlL-1ra suggested that a more ordered native conformation is induced by sucrose. Electron paramagnetic resonance spectroscopy demonstrated that in the presence of sucrose, spin-labeled cysteine 116 becomes more buried in the protein’s interior and that the hydrodynamic diameter of the protein is reduced. The preferential exclusion of sucrose from the protein and the resulting shift in the equilibrium between protein states toward the most compact conformation account for sucrose-induced effects on rhIL-1ra.

Conversion of cysteine to formylglycine: A protein modification in the endoplasmic reticulum

In sulfatases a Cα-formylglycine residue is found at a position where their cDNA sequences predict a cysteine residue. In multiple sulfatase deficiency, an inherited lysosomal storage disorder, catalytically inactive sulfatases are synthesized which retain the cysteine residue, indicating that the Cα-formylglycine residue is required for sulfatase activity. Using in vitro translation in the absence or presence of transport competent microsomes we found that newly synthesized sulfatase polypeptides carry a cysteine residue and that the oxidation of its thiol group to an aldehyde is catalyzed in the endoplasmic reticulum. A linear sequence of 16 residues surrounding the Cys-69 in arylsulfatase A is sufficient to direct the oxidation. This novel protein modification occurs after or at a late stage of cotranslational protein translocation into the endoplasmic reticulum when the polypeptide is not yet folded to its native structure.