Three-dimensional structure of human electron transfer flavoprotein to 2.1-Å resolution

Mammalian electron transfer flavoproteins (ETF) are heterodimers containing a single equivalent of flavin adenine dinucleotide (FAD). They function as electron shuttles between primary flavoprotein dehydrogenases involved in mitochondrial fatty acid and amino acid catabolism and the membrane-bound electron transfer flavoprotein ubiquinone oxidoreductase. The structure of human ETF solved to 2.1-Å resolution reveals that the ETF molecule is comprised of three distinct domains: two domains are contributed by the α subunit and the third domain is made up entirely by the β subunit. The N-terminal portion of the α subunit and the majority of the β subunit have identical polypeptide folds, in the absence of any sequence homology. FAD lies in a cleft between the two subunits, with most of the FAD molecule residing in the C-terminal portion of the α subunit. Alignment of all the known sequences for the ETF α subunits together with the putative FixB gene product shows that the residues directly involved in FAD binding are conserved. A hydrogen bond is formed between the N5 of the FAD isoalloxazine ring and the hydroxyl side chain of αT266, suggesting why the pathogenic mutation, αT266M, affects ETF activity in patients with glutaric acidemia type II. Hydrogen bonds between the 4′-hydroxyl of the ribityl chain of FAD and N1 of the isoalloxazine ring, and between αH286 and the C2-carbonyl oxygen of the isoalloxazine ring, may play a role in the stabilization of the anionic semiquinone. With the known structure of medium chain acyl-CoA dehydrogenase, we hypothesize a possible structure for docking the two proteins.

Design of three-dimensional domain-swapped dimers and fibrous oligomers

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 Å resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

Using three-dimensional microfluidic networks for solving computationally hard problems

This paper describes the design of a parallel algorithm that uses moving fluids in a three-dimensional microfluidic system to solve a nondeterministically polynomial complete problem (the maximal clique problem) in polynomial time. This algorithm relies on (i) parallel fabrication of the microfluidic system, (ii) parallel searching of all potential solutions by using fluid flow, and (iii) parallel optical readout of all solutions. This algorithm was implemented to solve the maximal clique problem for a simple graph with six vertices. The successful implementation of this algorithm to compute solutions for small-size graphs with fluids in microchannels is not useful, per se, but does suggest broader application for microfluidics in computation and control.

Computational adaptive optics for live three-dimensional biological imaging

Light microscopy of thick biological samples, such as tissues, is often limited by aberrations caused by refractive index variations within the sample itself. This problem is particularly severe for live imaging, a field of great current excitement due to the development of inherently fluorescent proteins. We describe a method of removing such aberrations computationally by mapping the refractive index of the sample using differential interference contrast microscopy, modeling the aberrations by ray tracing through this index map, and using space-variant deconvolution to remove aberrations. This approach will open possibilities to study weakly labeled molecules in difficult-to-image live specimens.

Three-dimensional image reconstruction of dephosphorylated smooth muscle heavy meromyosin reveals asymmetry in the interaction between myosin heads and placement of subfragment 2

Regulation of the actin-activated ATPase of smooth muscle myosin II is known to involve an interaction between the two heads that is controlled by phosphorylation of the regulatory light chain. However, the three-dimensional structure of this inactivated form has been unknown. We have used a lipid monolayer to obtain two-dimensional crystalline arrays of the unphosphorylated inactive form of smooth muscle heavy meromyosin suitable for structural studies by electron cryomicroscopy of unstained, frozen-hydrated specimens. The three-dimensional structure reveals an asymmetric interaction between the two myosin heads. The ATPase activity of one head is sterically “blocked” because part of its actin-binding interface is positioned onto the converter domain of the second head. ATPase activity of the second head, which can bind actin, appears to be inhibited through stabilization of converter domain movements needed to release phosphate and achieve strong actin binding. When the subfragment 2 domain of heavy meromyosin is oriented as it would be in an actomyosin filament lattice, the position of the heads is very different from that needed to bind actin, suggesting an additional contribution to ATPase inhibition in situ.

Prediction of protein deamidation rates from primary and three-dimensional structure

A method for the quantitative estimation of instability with respect to deamidation of the asparaginyl (Asn) residues in proteins is described. The procedure involves the observation of several simple aspects of the three-dimensional environment of each Asn residue in the protein and a calculation that includes these observations, the primary amino acid residue sequence, and the previously reported complete set of sequence-dependent rates of deamidation for Asn pentapeptides. This method is demonstrated and evaluated for 23 proteins in which 31 unstable and 167 stable Asn residues have been reported and for 7 unstable and 63 stable Asn residues that have been reported in 61 human hemoglobin variants. The relative importance of primary structure and three-dimensional structure in Asn deamidation is estimated.

Three-dimensional domain swapping in p13suc1 occurs in the unfolded state and is controlled by conserved proline residues

p13suc1 has two native states, a monomer and a domain-swapped dimer. We show that their folding pathways are connected by the denatured state, which introduces a kinetic barrier between monomer and dimer under native conditions. The barrier is lowered under conditions that speed up unfolding, thereby allowing, to our knowledge for the first time, a quantitative dissection of the energetics of domain swapping. The monomer–dimer equilibrium is controlled by two conserved prolines in the hinge loop that connects the exchanging domains. These two residues exploit backbone strain to specifically direct dimer formation while preventing higher-order oligomerization. Thus, the loop acts as a loaded molecular spring that releases tension in the monomer by adopting its alternative conformation in the dimer. There is an excellent correlation between domain swapping and aggregation, suggesting they share a common mechanism. These insights have allowed us to redesign the domain-swapping propensity of suc1 from a fully monomeric to a fully dimeric protein.

Three-dimensional reconstruction of the recombinant type 3 ryanodine receptor and localization of its amino terminus

Recombinant type 3 ryanodine receptor (RyR3) has been purified in quantities sufficient for structural characterization by cryoelectron microscopy and three-dimensional (3D) reconstruction. Two cDNAs were prepared and expressed in HEK293 cells, one encoding the wild-type RyR3 and the other encoding RyR3 containing glutathione S-transferase (GST) fused to its amino terminus (GST-RyR3). RyR3 was purified from detergent-solubilized transfected cells by affinity chromatography using 12.6-kDa FK506-binding protein in the form of a GST fusion as the affinity ligand. Purification of GST-RyR3 was achieved by affinity chromatography by using glutathione-Sepharose. Purified recombinant RyR3 and GST-RyR3 proteins exhibited high-affinity [3H]ryanodine binding that was sensitive to activation by Ca2+ and caffeine and to inhibition by Mg2+. 3D reconstructions of both recombinant RyR3 and GST-RyR3 appeared very similar to that of the native RyR3 purified from bovine diaphragm. Comparison of the 3D reconstructions of RyR3 and GST-RyR3 revealed that the GST domains and, hence, the amino termini of the RyR3 subunits are located in the “clamp” structures that form the corners of the square-shaped cytoplasmic region of homotetrameric RyR3. This study describes the 3D reconstruction of a recombinant ryanodine receptor and it demonstrates the potential of this technology for characterizing functional and structural perturbations introduced by site-directed mutagenesis.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

Three-dimensional cryoelectron microscopy of dimeric kinesin and ncd motor domains on microtubules.

Kinesin and ncd motor proteins are homologous in sequence yet move in opposite directions along microtubules. We have previously shown that monomeric kinesin and ncd bind in the same orientation on equivalent sites relative to the ends of tubulin sheets of known polarity. We now report cryoelectron microscope images of 16-protofilament microtubules decorated with both single- and double-headed kinesin and double-headed ncd. Three-dimensional density maps and difference maps show that, in adenosine 5'-[beta,gamma-imido]triphosphate, both dimeric motors bind tightly to microtubules via one head, leaving the other free, though apparently in a fixed position. The attached heads of dimers bind to tubulin in the same way as single kinesin heads. The second heads are connected to the tops of the first but, whereas the second kinesin head is closely associated with the first, pairs of ncd heads are splayed apart. There is also a distinct difference in orientation: the second kinesin head is tilted toward the microtubule plus end, while the second head of ncd points toward the minus end.

Characterizations of diverse residue clusters in protein three-dimensional structures.

We present new methods for identifying and analyzing statistically significant residue clusters that occur in three-dimensional (3D) protein structures. Residue clusters of different kinds occur in many contexts. They often feature the active site (e.g., in substrate binding), the interface between polypeptide units of protein complexes, regions of protein-protein and protein-nucleic acid interactions, or regions of metal ion coordination. The methods are illustrated with 3D clusters centering on four themes. (i) Acidic or histidine-acidic clusters associated with metal ions. (ii) Cysteine clusters including coordination of metals such as zinc or iron-sulfur structures, cysteine knots prominent in growth factors, multiple sets of buried disulfide pairings that putatively nucleate the hydrophobic core, or cysteine clusters of mostly exposed disulfide bridges. (iii) Iron-sulfur proteins and charge clusters. (iv) 3D environments of multiple histidine residues. Study of diverse 3D residue clusters offers a new perspective on protein structure and function. The algorithms can aid in rapid identification of distinctive sites, suggest correlations among protein structures, and serve as a tool in the analysis of new structures.

Clusters of charged residues in protein three-dimensional structures.

Statistically significant charge clusters (basic, acidic, or of mixed charge) in tertiary protein structures are identified by new methods from a large representative collection of protein structures. About 10% of protein structures show at least one charge cluster, mostly of mixed type involving about equally anionic and cationic residues. Positive charge clusters are very rare. Negative (or histidine-acidic) charge clusters often coordinate calcium, or magnesium or zinc ions [e.g., thermolysin (PDB code: 3tln), mannose-binding protein (2msb), aminopeptidase (1amp)]. Mixed-charge clusters are prominent at interchain contacts where they stabilize quaternary protein formation [e.g., glutathione S-transferase (2gst), catalase (8act), and fructose-1,6-bisphosphate aldolase (1fba)]. They are also involved in protein-protein interaction and in substrate binding. For example, the mixed-charge cluster of aspartate carbamoyl-transferase (8atc) envelops the aspartate carbonyl substrate in a flexible manner (alternating tense and relaxed states) where charge associations can vary from weak to strong. Other proteins with charge clusters include the P450 cytochrome family (BM-3, Terp, Cam), several flavocytochromes, neuraminidase, hemagglutinin, the photosynthetic reaction center, and annexin. In each case in Table 2 we discuss the possible role of the charge clusters with respect to protein structure and function.

Suppression of apoptosis by basement membrane requires three-dimensional tissue organization and withdrawal from the cell cycle.

The basement membrane (BM) extracellular matrix induces differentiation and suppresses apoptosis in mammary epithelial cells, whereas cells lacking BM lose their differentiated phenotype and undergo apoptosis. Addition of purified BM components, which are known to induce beta-casein expression, did not prevent apoptosis, indicating that a more complex BM was necessary. A comparison of culture conditions where apoptosis would or would not occur allowed us to relate inhibition of apoptosis to a complete withdrawal from the cell cycle, which was observed only when cells acquired a three-dimensional alveolar structure in response to BM. In the absence of this morphology, both the GI cyclin kinase inhibitor p21/WAF-1 and positive proliferative signals including c-myc and cyclin DI were expressed and the retinoblastoma protein (Rb) continued to be hyperphosphorylated. When we overexpressed either c-myc in quiescent cells or p21 when cells were still cycling, apoptosis was induced. In the absence of three-dimensional alveolar structures, mammary epithelial cells secrete a number of factors including transforming growth factor alpha and tenascin, which when added exogenously to quiescent cells induced expression of c-myc and interleukin-beta1-converting enzyme (ICE) mRNA and led to apoptosis. These experiments demonstrate that a correct tissue architecture is crucial for long-range homeostasis, suppression of apoptosis, and maintenance of differentiated phenotype.

Three-dimensional structure of recombinant human osteogenic protein 1: structural paradigm for the transforming growth factor beta superfamily.

We report the three-dimensional structure of osteogenic protein 1 (OP-1, also known as bone morphogenetic protein 7) to 2.8-A resolution. OP-1 is a member of the transforming growth factor beta (TGF-beta) superfamily of proteins and is able to induce new bone formation in vivo. Members of this superfamily share sequence similarity in their C-terminal regions and are implicated in embryonic development and adult tissue repair. Our crystal structure makes possible the structural comparison between two members of the TGF-beta superfamily. We find that although there is limited sequence identity between OP-1 and TGF-beta 2, they share a common polypeptide fold. These results establish a basis for proposing the OP-1/TGF-beta 2 fold as the primary structural motif for the TGF-beta superfamily as a whole. Detailed comparison of the OP-1 and TGF-beta 2 structures has revealed striking differences that provide insights into how these growth factors interact with their receptors.

Analysis of the relation between the sequence and secondary and three-dimensional structures of immunoglobulin molecules.

Methods of structural and statistical analysis of the relation between the sequence and secondary and three-dimensional structures are developed. About 5000 secondary structures of immunoglobulin molecules from the Kabat data base were predicted. Two statistical analyses of amino acids reveal 47 universal positions in strands and loops. Eight universally conservative positions out of the 47 are singled out because they contain the same amino acid in > 90% of all chains. The remaining 39 positions, which we term universally alternative positions, were divided into five groups: hydrophobic, charged and polar, aromatic, hydrophilic, and Gly-Ala, corresponding to the residues that occupied them in almost all chains. The analysis of residue-residue contacts shows that the 47 universal positions can be distinguished by the number and types of contacts. The calculations of contact maps in the 29 antibody structures revealed that residues in 24 of these 47 positions have contacts only with residues of antiparallel beta-strands in the same beta-sheet and residues in the remaining 23 positions always have far-away contacts with residues from other beta-sheets as well. In addition, residues in 6 of the 47 universal positions are also involved in interactions with residues of the other variable or constant domains.