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 polyamide ("peptide") nucleic acid binding to duplex DNA.

To elucidate the mechanism of recognition of double-stranded DNA (dsDNA) by homopyrimidine polyamide ("peptide") nucleic acid (PNA) leading to the strand-displacement, the kinetics of the sequence-specific PNA/DNA binding have been studied. The binding was monitored with time by the gel retardation and nuclease S1 cleavage assays. The experimental kinetic curves obey pseudo-first-order kinetics and the dependence of the pseudo-first-order rate constant, kps, on PNA concentration, P, obeys a power law kps approximately P gamma with 2 < gamma < 3. The kps values for binding of decamer PNA to dsDNA target sites with one mismatch are hundreds of times slower than for the correct site. A detailed kinetic scheme for PNA/DNA binding is proposed that includes two major steps of the reaction of strand invasion: (i) a transient partial opening of the PNA binding site on dsDNA and incorporation of one PNA molecule with the formation of an intermediate PNA/DNA duplex and (ii) formation of a very stable PNA2/DNA triplex. A simple theoretical treatment of the proposed kinetic scheme is performed. The interpretation of our experimental data in the framework of the proposed kinetic scheme leads to the following conclusions. The sequence specificity of the recognition is essentially provided at the "search" step of the process, which consists in the highly reversible transient formation of duplex between one PNA molecule and the complementary strand of duplex DNA while the other DNA strand is displaced. This search step is followed by virtually irreversible "locking" step via PNA2/DNA triplex formation. The proposed mechanism explains how the binding of homopyrimidine PNA to dsDNA meets two apparently mutually contradictory features: high sequence specificity of binding and remarkable stability of both correct and mismatched PNA/DNA complexes.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

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.

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Kinetics and thermodynamics of T cell receptor– autoantigen interactions in murine experimental autoimmune encephalomyelitis

In the current study, cellular and molecular approaches have been used to analyze the biophysical nature of T cell receptor (TCR)–peptide MHC (pMHC) interactions for two autoreactive TCRs. These two TCRs recognize the N-terminal epitope of myelin basic protein (MBP1–11) bound to the MHC class II protein, I-Au, and are associated with murine experimental autoimmune encephalomyelitis. Mice transgenic for the TCRs have been generated and characterized in other laboratories. These analyses indicate that the mice either develop encephalomyelitis spontaneously (172.10 TCR) or only if immunized with autoantigen in adjuvant (1934.4 TCR). Here, we show that the 172.10 TCR binds MBP1–11:I-Au with a 4–5-fold higher affinity than the 1934.4 TCR. Consistent with the higher affinity, 172.10 T hybridoma cells are significantly more responsive to autoantigen than 1934.4 cells. The interaction of the 172.10 TCR with cognate ligand is more entropically unfavorable than that of the 1934.4 TCR, indicating that the 172.10 TCR undergoes greater conformational rearrangements upon ligand binding. The studies therefore suggest a correlation between the strength and plasticity of a TCR–pMHC interaction and the frequency of spontaneous disease in the corresponding TCR transgenic mice. The comparative analysis of these two TCRs has implications for understanding autoreactive T cell recognition and activation.

Symmetry, stability, and dynamics of multidomain and multicomponent protein systems

Symmetry is commonly observed in many biological systems. Here we discuss representative examples of the role of symmetry in structural molecular biology. Point group symmetries are observed in many protein oligomers whose three-dimensional atomic structures have been elucidated by x-ray crystallography. Approximate symmetry also occurs in multidomain proteins. Symmetry often confers stability on the molecular system and results in economical usage of basic components to build the macromolecular structure. Symmetry is also associated with cooperativity. Mild perturbation from perfect symmetry may be essential in some systems for dynamic functions.

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.

Effects of Hypoxia on 13NH4+ Fluxes in Rice Roots1 : Kinetics and Compartmental Analysis

Techniques of compartmental (efflux) and kinetic influx analyses with the radiotracer 13NH4+ were used to examine the adaptation to hypoxia (15, 35, and 50% O2 saturation) of root N uptake and metabolism in 3-week-old hydroponically grown rice (Oryza sativa L., cv IR72) seedlings. A time-dependence study of NH4+ influx into rice roots after onset of hypoxia (15% O2) revealed an initial increase in the first 1 to 2.5 h after treatment imposition, followed by a decline to less than 50% of influx in control plants by 4 d. Efflux analyses conducted 0, 1, 3, and 5 d after the treatment confirmed this adaptation pattern of NH4+ uptake. Half-lives for NH4+ exchange with subcellular compartments, cytoplasmic NH4+ concentrations, and efflux (as percentage of influx) were unaffected by hypoxia. However, significant differences were observed in the relative amounts of N allocated to NH4+ assimilation and the vacuole versus translocation to the shoot. Kinetic experiments conducted at 100, 50, 35, and 15% O2 saturation showed no significant change in the Km value for NH4+ uptake with varying O2 supply. However, Vmax was 42% higher than controls at 50% O2 saturation, unchanged at 35%, and 10% lower than controls at 15% O2. The significance of these flux adaptations is discussed.

Domain structure and dynamics in the helical filaments formed by RecA and Rad51 on DNA

Both the bacterial RecA protein and the eukaryotic Rad51 protein form helical nucleoprotein filaments on DNA that catalyze strand transfer between two homologous DNA molecules. However, only the ATP-binding cores of these proteins have been conserved, and this same core is also found within helicases and the F1-ATPase. The C-terminal domain of the RecA protein forms lobes within the helical RecA filament. However, the Rad51 proteins do not have the C-terminal domain found in RecA, but have an N-terminal extension that is absent in the RecA protein. Both the RecA C-terminal domain and the Rad51 N-terminal domain bind DNA. We have used electron microscopy to show that the lobes of the yeast and human Rad51 filaments appear to be formed by N-terminal domains. These lobes are conformationally flexible in both RecA and Rad51. Within RecA filaments, the change between the “active” and “inactive” states appears to mainly involve a large movement of the C-terminal lobe. The N-terminal domain of Rad51 and the C-terminal domain of RecA may have arisen from convergent evolution to play similar roles in the filaments.

Correlation between kinetics and RNA splicing of alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptors in neocortical neurons.

In the cortex fast excitatory synaptic currents onto excitatory pyramidal neurons and inhibitory nonpyramidal neurons are mediated by alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) receptors exhibiting cell-type-specific differences in their kinetic properties. AMPA receptors consist of four subunits (GluR1-4), each existing as two splice variants, flip and flop, which critically affect the desensitization properties of receptors expressed in heterologous systems. Using single cell reverse transcription PCR to analyze the mRNA of AMPA receptor subunits expressed in layers I-III neocortical neurons, we find that 90% of the GluR1-4 in nonpyramidal neurons are flop variants, whereas 92% of the GluR1-4 in pyramidal neurons are flip variants. We also find that nonpyramidal neurons predominantly express GluR1 mRNA (GluR1/GluR1-4 = 59%), whereas pyramidal neurons contain mainly GluR2 mRNA (GluR2/GluR1-4 = 59%). However, the neuron-type-specific splicing is exhibited by all four AMPA receptor subunits. We suggest that the predominance of the flop variants contributes to the faster and more extensive desensitization in nonpyramidal neurons, compared to pyramidal cells where flip variants are dominant. Alternative splicing of AMPA receptors may play an important role in regulating synaptic function in a cell-type-specific manner, without changing permeation properties.

G protein activation kinetics and spillover of gamma-aminobutyric acid may account for differences between inhibitory responses in the hippocampus and thalamus.

We have developed a model of gamma-aminobutyric acid (GABA)ergic synaptic transmission mediated by GABAA and GABAB receptors, including cooperativity in the guanine nucleotide binding protein (G protein) cascade mediating the activation of K+ channels by GABAB receptors. If the binding of several G proteins is needed to activate the K+ channels, then only a prolonged activation of GABAB receptors evoked detectable currents. This could occur if strong stimuli evoked release in adjacent terminals and the spillover resulted in prolonged activation of the receptors, leading to inhibitory responses similar to those observed in hippocampal slices. The same model also reproduced thalamic GABAB responses to high-frequency bursts of stimuli. In this case, prolonged activation of the receptors was due to high-frequency release conditions. This model provides insights into the function of GABAB receptors in normal and epileptic discharges.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

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.

Proton uptake by bacterial reaction centers: The protein complex responds in a similar manner to the reduction of either quinone acceptor

In bacterial photosynthetic reaction centers, the protonation events associated with the different reduction states of the two quinone molecules constitute intrinsic probes of both the electrostatic interactions and the different kinetic events occurring within the protein in response to the light-generated introduction of a charge. The kinetics and stoichiometries of proton uptake on formation of the primary semiquinone QA− and the secondary acceptor QB− after the first and second flashes have been measured, at pH 7.5, in reaction centers from genetically modified strains and from the wild type. The modified strains are mutated at the L212Glu and/or at the L213Asp sites near QB; some of them carry additional mutations distant from the quinone sites (M231Arg → Leu, M43Asn → Asp, M5Asn → Asp) that compensate for the loss of L213Asp. Our data show that the mutations perturb the response of the protein system to the formation of a semiquinone, how distant compensatory mutations can restore the normal response, and the activity of a tyrosine residue (M247Ala → Tyr) in increasing and accelerating proton uptake. The data demonstrate a direct correlation between the kinetic events of proton uptake that are observed with the formation of either QA− or QB−, suggesting that the same residues respond to the generation of either semiquinone species. Therefore, the efficiency of transferring the first proton to QB is evident from examination of the pattern of H+/QA− proton uptake. This delocalized response of the protein complex to the introduction of a charge is coordinated by an interactive network that links the Q− species, polarizable residues, and numerous water molecules that are located in this region of the reaction center structure. This could be a general property of transmembrane redox proteins that couple electron transfer to proton uptake/release reactions.