From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

Flickering fusion pores comparable with initial exocytotic pores occur in protein-free phospholipid bilayers

For the act of membrane fusion, there are two competing, mutually exclusive molecular models that differ in the structure of the initial pore, the pathway for ionic continuity between formerly separated volumes. Because biological “fusion pores” can be as small as ionic channels or gap junctions, one model posits a proteinaceous initial fusion pore. Because biological fusion pore conductance varies widely, another model proposes a lipidic initial pore. We have found pore opening and flickering during the fusion of protein-free phospholipid vesicles with planar phospholipid bilayers. Fusion pore formation appears to follow the coalescence of contacting monolayers to create a zone of hemifusion where continuity between the two adherent membranes is lipidic, but not aqueous. Hypotonic stress, causing tension in the vesicle membrane, promotes complete fusion. Pores closed soon after opening (flickering), and the distribution of fusion pore conductance appears similar to the distribution of initial fusion pores in biological fusion. Because small flickering pores can form in the absence of protein, the existence of small pores in biological fusion cannot be an argument in support of models based on proteinaceous pores. Rather, these results support the model of a lipidic fusion pore developing within a hemifused contact site.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Efficient Trafficking of TGN38 from the Endosome to the trans-Golgi Network Requires a Free Hydroxyl Group at Position 331 in the Cytosolic Domain

TGN38 is one of the few known resident integral membrane proteins of the trans-Golgi network (TGN). Since it cycles constitutively between the TGN and the plasma membrane, TGN38 is ideally suited as a model protein for the identification of post-Golgi trafficking motifs. Several studies, employing chimeric constructs to detect such motifs within the cytosolic domain of TGN38, have identified the sequence 333YQRL336 as an autonomous signal capable of localizing reporter proteins to the TGN. In addition, one group has found that an upstream serine residue, S331, may also play a role in TGN38 localization. However, the nature and degree of participation of S331 in the localization of TGN38 remain uncertain, and the effect has been studied in chimeric constructs only. Here we investigate the role of S331 in the context of full-length TGN38. Mutations that abolish the hydroxyl moiety at position 331 (A, D, and E) lead to missorting of endocytosed TGN38 to the lysosome. Conversely, mutation of S331 to T has little effect on the endocytic trafficking of TGN38. Together, these findings indicate that the S331 hydroxyl group has a direct or indirect effect on the ability of the cytosolic tail of TGN38 to interact with trafficking and/or sorting machinery at the level of the early endosome. In addition, mutation of S331 to either A or D results in increased levels of TGN38 at the cell surface. The results confirm that S331 plays a critical role in the intracellular trafficking of TGN38 and further reveal that TGN38 undergoes a signal-mediated trafficking step at the level of the endosome.

Biparental Inheritance of γ-Tubulin during Human Fertilization: Molecular Reconstitution of Functional Zygotic Centrosomes in Inseminated Human Oocytes and in Cell-free Extracts Nucleated by Human Sperm

Human sperm centrosome reconstitution and the parental contributions to the zygotic centrosome are examined in mammalian zygotes and after exposure of spermatozoa to Xenopus laevis cell-free extracts. The presence and inheritance of the conserved centrosomal constituents γ-tubulin, centrin, and MPM-2 (which detects phosphorylated epitopes) are traced, as is the sperm microtubule-nucleating capability on reconstituted centrosomes. γ-Tubulin is biparentally inherited in humans (maternal >> than paternal): Western blots detect the presence of paternal γ-tubulin. Recruitment of maternal γ-tubulin to the sperm centrosome occurs after sperm incorporation in vivo or exposure to cell-free extract, especially after sperm “priming” induced by disulfide bond reduction. Centrin is found in the proximal sperm centrosomal region, demonstrates expected calcium sensitivity, but appears absent from the zygotic centrosome after sperm incorporation or exposure to extracts. Sperm centrosome phosphorylation is detected after exposure of primed sperm to egg extracts as well as during the early stages of sperm incorporation after fertilization. Finally, centrosome reconstitution in cell-free extracts permits sperm aster microtubule assembly in vitro. Collectively, these results support a model of a blended zygotic centrosome composed of maternal constituents attracted to an introduced paternal template after insemination.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

Genetically altered AMPA-type glutamate receptor kinetics in interneurons disrupt long-range synchrony of gamma oscillation

Gamma oscillations synchronized between distant neuronal populations may be critical for binding together brain regions devoted to common processing tasks. Network modeling predicts that such synchrony depends in part on the fast time course of excitatory postsynaptic potentials (EPSPs) in interneurons, and that even moderate slowing of this time course will disrupt synchrony. We generated mice with slowed interneuron EPSPs by gene targeting, in which the gene encoding the 67-kDa form of glutamic acid decarboxylase (GAD67) was altered to drive expression of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate receptor subunit GluR-B. GluR-B is a determinant of the relatively slow EPSPs in excitatory neurons and is normally expressed at low levels in γ-aminobutyric acid (GABA)ergic interneurons, but at high levels in the GAD-GluR-B mice. In both wild-type and GAD-GluR-B mice, tetanic stimuli evoked gamma oscillations that were indistinguishable in local field potential recordings. Remarkably, however, oscillation synchrony between spatially separated sites was severely disrupted in the mutant, in association with changes in interneuron firing patterns. The congruence between mouse and model suggests that the rapid time course of AMPA receptor-mediated EPSPs in interneurons might serve to allow gamma oscillations to synchronize over distance.

Electrostatic surface plasmon resonance: Direct electric field-induced hybridization and denaturation in monolayer nucleic acid films and label-free discrimination of base mismatches

We demonstrate that in situ optical surface plasmon resonance spectroscopy can be used to monitor hybridization kinetics for unlabeled DNA in tethered monolayer nucleic acid films on gold in the presence of an applied electrostatic field. The dc field can enhance or retard hybridization and can also denature surface-immobilized DNA duplexes. Discrimination between matched and mismatched hybrids is achieved by simple adjustment of the electrode potential. Although the electric field at the interface is extremely large, the tethered single-stranded DNA thiol probes remain bound and can be reused for subsequent hybridization reactions without loss of efficiency. Only capacitive charging currents are drawn; redox reactions are avoided by maintaining the gold electrode potential within the ideally polarizable region. Because of potential-induced changes in the shape of the surface plasmon resonance curve, we account for the full curve rather than simply the shift in the resonance minimum.

Quantitative ER ↔ Golgi Transport Kinetics and Protein Separation upon Golgi Exit Revealed by Vesicular Integral Membrane Protein 36 Dynamics in Live CellsV⃞

To quantitatively investigate the trafficking of the transmembrane lectin VIP36 and its relation to cargo-containing transport carriers (TCs), we analyzed a C-terminal fluorescent-protein (FP) fusion, VIP36-SP-FP. When expressed at moderate levels, VIP36-SP-FP localized to the endoplasmic reticulum, Golgi apparatus, and intermediate transport structures, and colocalized with epitope-tagged VIP36. Temperature shift and pharmacological experiments indicated VIP36-SP-FP recycled in the early secretory pathway, exhibiting trafficking representative of a class of transmembrane cargo receptors, including the closely related lectin ERGIC53. VIP36-SP-FP trafficking structures comprised tubules and globular elements, which translocated in a saltatory manner. Simultaneous visualization of anterograde secretory cargo and VIP36-SP-FP indicated that the globular structures were pre-Golgi carriers, and that VIP36-SP-FP segregated from cargo within the Golgi and was not included in post-Golgi TCs. Organelle-specific bleach experiments directly measured the exchange of VIP36-SP-FP between the Golgi and endoplasmic reticulum (ER). Fitting a two-compartment model to the recovery data predicted first order rate constants of 1.22 ± 0.44%/min for ER → Golgi, and 7.68 ± 1.94%/min for Golgi → ER transport, revealing a half-time of 113 ± 70 min for leaving the ER and 1.67 ± 0.45 min for leaving the Golgi, and accounting for the measured steady-state distribution of VIP36-SP-FP (13% Golgi/87% ER). Perturbing transport with AlF4− treatment altered VIP36-SP-GFP distribution and changed the rate constants. The parameters of the model suggest that relatively small differences in the first order rate constants, perhaps manifested in subtle differences in the tendency to enter distinct TCs, result in large differences in the steady-state localization of secretory components.

Kinetics and mechanism of DNA uptake into the cell nucleus

Gene transfer to eukaryotic cells requires the uptake of exogenous DNA into the cell nucleus. Except during mitosis, molecular access to the nuclear interior is limited to passage through the nuclear pores. Here we demonstrate the nuclear uptake of extended linear DNA molecules by a combination of fluorescence microscopy and single-molecule manipulation techniques, using the latter to follow uptake kinetics of individual molecules in real time. The assays were carried out on nuclei reconstituted in vitro from extracts of Xenopus eggs, which provide both a complete complement of biochemical factors involved in nuclear protein import, and unobstructed access to the nuclear pores. We find that uptake of DNA is independent of ATP or GTP hydrolysis, but is blocked by wheat germ agglutinin. The kinetics are much slower than would be expected from hydrodynamic considerations. A fit of the data to a simple model suggests femto-Newton forces and a large friction relevant to the uptake process.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Intracellular pH in adipocytes: effects of free fatty acid diffusion across the plasma membrane, lipolytic agonists, and insulin.

The main function of white adipose tissue is to store nutrient energy in the form of triglycerides. The mechanism by which free fatty acids (FFA) move into and out of the adipocyte has not been resolved. We show here that changes in intracellular pH (pH1) in adipocytes correlate with the movement of FFA across cellular membranes as predicted by the Kamp and Hamilton model of passive diffusion of FFA. Exposure of fat cells to lipolytic agents or external FFA results is a rapid intracellular acidification that is reversed by metabolism of the FFA or its removal by albumin. In contrast, insulin causes an alkalinization of the cell, consistent with its main function to promote esterification. Inhibition of Na+/H+ exchange in adipocytes does not prevent the changes in pHi caused by FFA, lipolytic agents, or insulin. A fatty acid dimer, which diffuses into the cell but is not metabolized, causes an irreversible acidification. Taken together, the data suggest that changes in pHi occur in adipocytes in response to the passive diffusion of un-ionized FFA (flip-flop) into and out of the cell and in response to their metabolism and production within the cell. These changes in pHi may, in turn, modulate hormonal signaling and metabolism with significant impact on cell function.