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.

Interacting helical faces of subunits a and c in the F1Fo ATP synthase of Escherichia coli defined by disulfide cross-linking

Subunits a and c of Fo are thought to cooperatively catalyze proton translocation during ATP synthesis by the Escherichia coli F1Fo ATP synthase. Optimizing mutations in subunit a at residues A217, I221, and L224 improves the partial function of the cA24D/cD61G double mutant and, on this basis, these three residues were proposed to lie on one face of a transmembrane helix of subunit a, which then interacted with the transmembrane helix of subunit c anchoring the essential aspartyl group. To test this model, in the present work Cys residues were introduced into the second transmembrane helix of subunit c and the predicted fourth transmembrane helix of subunit a. After treating the membrane vesicles of these mutants with Cu(1,10-phenanthroline)2SO4 at 0°, 10°, or 20°C, strong a–c dimer formation was observed at all three temperatures in membranes of 7 of the 65 double mutants constructed, i.e., in the aS207C/cI55C, aN214C/cA62C, aN214C/cM65C, aI221C/cG69C, aI223C/cL72C, aL224C/cY73C, and aI225C/cY73C double mutant proteins. The pattern of cross-linking aligns the helices in a parallel fashion over a span of 19 residues with the aN214C residue lying close to the cA62C and cM65C residues in the middle of the membrane. Lesser a–c dimer formation was observed in nine other double mutants after treatment at 20°C in a pattern generally supporting that indicated by the seven landmark residues cited above. Cross-link formation was not observed between helix-1 of subunit c and helix-4 of subunit a in 19 additional combinations of doubly Cys-substituted proteins. These results provide direct chemical evidence that helix-2 of subunit c and helix-4 of subunit a pack close enough to each other in the membrane to interact during function. The proximity of helices supports the possibility of an interaction between Arg210 in helix-4 of subunit a and Asp61 in helix-2 of subunit c during proton translocation, as has been suggested previously.

Conformational changes couple Na+ and glucose transport

The mechanism by which cotransport proteins couple their substrates across cell membranes is not known. A commonly proposed model is that cotransport results from ligand-induced conformational transitions that change the accessibility of ligand-binding sites from one side of the membrane to the other. To test this model, we have measured the accessibility of covalent probes to a cysteine residue (Q457C) placed in the putative sugar-translocation domain of the Na+/glucose cotransporter (SGLT1). The mutant protein Q457C was able to transport sugar, but transport was abolished after alkylation by methanethiosulfonate reagents. Alkylation blocked sugar translocation but not sugar binding. Accessibility of Q457C to alkylating reagents required external Na+ and was blocked by external sugar and phlorizin. The voltage dependence of accessibility was directly correlated with the presteady–state charge movement of SGLT1. Voltage-jump experiments with rhodamine-6-maleimide-labeled Q457C showed that the time course and level of changes in fluorescence closely followed the presteady–state charge movement. We conclude that conformational changes are responsible for the coupling of Na+ and sugar transport and that Q457 plays a critical role in sugar translocation by SGLT1.

Interactions among enzymes of the Arabidopsis flavonoid biosynthetic pathway

Flavonoids are secondary metabolites derived from phenylalanine and acetate metabolism that perform a variety of essential functions in higher plants. Studies over the past 30 years have supported a model in which flavonoid metabolism is catalyzed by an enzyme complex localized to the endoplasmic reticulum [Hrazdina, G. & Wagner, G. J. (1985) Arch. Biochem. Biophys. 237, 88–100]. To test this model further we assayed for direct interactions between several key flavonoid biosynthetic enzymes in developing Arabidopsis seedlings. Two-hybrid assays indicated that chalcone synthase, chalcone isomerase (CHI), and dihydroflavonol 4-reductase interact in an orientation-dependent manner. Affinity chromatography and immunoprecipitation assays further demonstrated interactions between chalcone synthase, CHI, and flavonol 3-hydroxylase in lysates from Arabidopsis seedlings. These results support the hypothesis that the flavonoid enzymes assemble as a macromolecular complex with contacts between multiple proteins. Evidence was also found for posttranslational modification of CHI. The importance of understanding the subcellular organization of elaborate enzyme systems is discussed in the context of metabolic engineering.

Identification of the vitamin K-dependent carboxylase active site: Cys-99 and Cys-450 are required for both epoxidation and carboxylation

The vitamin K-dependent carboxylase modifies and renders active vitamin K-dependent proteins involved in hemostasis, cell growth control, and calcium homeostasis. Using a novel mechanism, the carboxylase transduces the free energy of vitamin K hydroquinone (KH2) oxygenation to convert glutamate into a carbanion intermediate, which subsequently attacks CO2, generating the γ-carboxylated glutamate product. How the carboxylase effects this conversion is poorly understood because the active site has not been identified. Dowd and colleagues [Dowd, P., Hershline, R., Ham, S. W. & Naganathan, S. (1995) Science 269, 1684–1691] have proposed that a weak base (cysteine) produces a strong base (oxygenated KH2) capable of generating the carbanion. To define the active site and test this model, we identified the amino acids that participate in these reactions. N-ethyl maleimide inhibited epoxidation and carboxylation, and both activities were equally protected by KH2 preincubation. Amino acid analysis of 14C- N-ethyl maleimide-modified human carboxylase revealed 1.8–2.3 reactive residues and a specific activity of 7 × 108 cpm/hr per mg. Tryptic digestion and liquid chromatography electrospray mass spectrometry identified Cys-99 and Cys-450 as active site residues. Mutation to serine reduced both epoxidation and carboxylation, to 0.2% (Cys-99) or 1% (Cys-450), and increased the Kms for a glutamyl substrate 6- to 8-fold. Retention of some activity indicates a mechanism for enhancing cysteine/serine nucleophilicity, a property shared by many active site thiol enzymes. These studies, which represent a breakthrough in defining the carboxylase active site, suggest a revised model in which the glutamyl substrate indirectly coordinates at least one thiol, forming a catalytic complex that ionizes a thiol to initiate KH2 oxygenation.

Genetic drift and within-host metapopulation dynamics of HIV-1 infection

Most HIV replication occurs in solid lymphoid tissue, which has prominent architecture at the histological level, which separates groups of productively infected CD4+ cells. Nevertheless, current population models of HIV assume panmixis within lymphoid tissue. We present a simple “metapopulation” model of HIV replication, where the population of infected cells is comprised of a large number of small populations, each of which is established by a few founder viruses and undergoes turnover. To test this model, we analyzed viral genetic variation of infected cell subpopulations within the spleen and demonstrated the action of founder effects as well as significant variation in the extent of genetic differentiation between subpopulations among patients. The combination of founder effects and subpopulation turnover can result in an effective population size much lower than the actual population size and may contribute to the importance of genetic drift in HIV evolution despite a large number of infected cells.

Copy-up mutants of the plasmid RK2 replication initiation protein are defective in coupling RK2 replication origins.

The broad host range plasmid RK2 replicates and regulates its copy number in a wide range of Gram-negative bacteria. The plasmid-encoded trans-acting replication protein TrfA and the origin of replication oriV are sufficient for controlled replication of the plasmid in all Gram-negative bacteria tested. The TrfA protein binds specifically to direct repeat sequences (iterons) at the origin of replication. A replication control model, designated handcuffing or coupling, has been proposed whereby the formation of coupled TrfA-oriV complexes between plasmid molecules results in hindrance of origin activity and, consequently, a shut-down of plasmid replication under conditions of higher than normal copy number. Therefore, according to this model, the coupling activity of an initiation protein is essential for copy number control and a copy-up initiation protein mutant should have reduced ability to form coupled complexes. To test this model for plasmid RK2, two previously characterized copy-up TrfA mutations, trfA-254D and trfA-267L, were combined and the resulting copy-up double mutant TFrfA protein TrfA-254D/267L was characterized. Despite initiating runaway (uncontrolled) replication in vivo, the copy-up double-mutant TrfA protein exhibited replication kinetics similar to the wild-type protein in vitro. Purified TrfA-254D, TrfA-267L, and TrfA-254D/267L proteins were then examined for binding to the iterons and for coupling activity using an in vitro ligase-catalyzed multimerization assay. It was found that both single and double TrfA mutant proteins exhibited substantially reduced (single mutants) or barely detectable (double mutant) levels of coupling activity while not being diminished in their capacity to bind to the origin of replication. These observations provide direct evidence in support of the coupling model of replication control.

In vivo regulation of muscle glycogen synthase and the control of glycogen synthesis.

The activity of glycogen synthase (GSase; EC 2.4.1.11) is regulated by covalent phosphorylation. Because of this regulation, GSase has generally been considered to control the rate of glycogen synthesis. This hypothesis is examined in light of recent in vivo NMR experiments on rat and human muscle and is found to be quantitatively inconsistent with the data under conditions of glycogen synthesis. Our first experiments showed that muscle glycogen synthesis was slower in non-insulin-dependent diabetics compared to normals and that their defect was in the glucose transporter/hexokinase (GT/HK) part of the pathway. From these and other in vivo NMR results a quantitative model is proposed in which the GT/HK steps control the rate of glycogen synthesis in normal humans and rat muscle. The flux through GSase is regulated to match the proximal steps by "feed forward" to glucose 6-phosphate, which is a positive allosteric effector of all forms of GSase. Recent in vivo NMR experiments specifically designed to test the model are analyzed by metabolic control theory and it is shown quantitatively that the GT/HK step controls the rate of glycogen synthesis. Preliminary evidence favors the transporter step. Several conclusions are significant: (i) glucose transport/hexokinase controls the glycogen synthesis flux; (ii) the role of covalent phosphorylation of GSase is to adapt the activity of the enzyme to the flux and to control the metabolite levels not the flux; (iii) the quantitative data needed for inferring and testing the present model of flux control depended upon advances of in vivo NMR methods that accurately measured the concentration of glucose 6-phosphate and the rate of glycogen synthesis.

The recombination hot spot chi activates RecBCD recombination by converting Escherichia coli to a recD mutant phenocopy.

The products of the recB and recC genes are necessary for conjugal recombination and for repair of chromosomal double-chain breaks in Escherichia coli. The recD gene product combines with the RecB and RecC proteins to comprise RecBCD enzyme but is required for neither recombination nor repair. On the contrary, RecBCD enzyme is an exonuclease that inhibits recombination by destroying linear DNA. The RecD ejection model proposes that RecBCD enzyme enters a DNA duplex at a double-chain end and travels destructively until it encounters the recombination hot spot sequence chi. Chi then alters the RecBCD enzyme by weakening the affinity of the RecD subunit for the RecBC heterodimer. With the loss of the RecD subunit, the resulting protein, RecBC(D-), becomes deficient for exonuclease activity and proficient as a recombinagenic helicase. To test the model, genetic crosses between lambda phage were conducted in cells containing chi on a nonhomologous plasmid. Upon delivering a double-chain break to the plasmid, lambda recombined as if the cells had become recD mutants. The ability of chi to alter lambda recombination in trans was reversed by overproducing the RecD subunit. These results indicate that chi can influence a recombination act without directly participating in it.

A test of the "jigsaw puzzle" model for protein folding by multiple methionine substitutions within the core of T4 lysozyme.

To test whether the structure of a protein is determined in a manner akin to the assembly of a jigsaw puzzle, up to 10 adjacent residues within the core of T4 lysozyme were replaced by methionine. Such variants are active and fold cooperatively with progressively reduced stability. The structure of a seven-methionine variant has been shown, crystallographically, to be similar to wild type and to maintain a well ordered core. The interaction between the core residues is, therefore, not strictly comparable with the precise spatial complementarity of the pieces of a jigsaw puzzle. Rather, a certain amount of give and take in forming the core structure is permitted. A simplified hydrophobic core sequence, imposed without genetic selection or computer-based design, is sufficient to retain native properties in a globular protein.