Transmembrane orientation and possible role of the fusogenic peptide from parainfluenza virus 5 (PIV5) in promoting fusion.

Membrane fusion is required for diverse biological functions ranging from viral infection to neurotransmitter release. Fusogenic proteins increase the intrinsically slow rate of fusion by coupling energetically downhill conformational changes of the protein to kinetically unfavorable fusion of the membrane-phospholipid bilayers. Class I viral fusogenic proteins have an N-terminal hydrophobic fusion peptide (FP) domain, important for interaction with the target membrane, plus a C-terminal transmembrane (C-term-TM) helical membrane anchor. The role of the water-soluble regions of fusogenic proteins has been extensively studied, but the contributions of the membrane-interacting FP and C-term-TM peptides are less well characterized. Typically, FPs are thought to bind to membranes at an angle that allows helix penetration but not traversal of the lipid bilayer. Here, we show that the FP from the paramyxovirus parainfluenza virus 5 fusogenic protein, F, forms an N-terminal TM helix, which self-associates into a hexameric bundle. This FP also interacts strongly with the C-term-TM helix. Thus, the fusogenic F protein resembles SNARE proteins involved in vesicle fusion by having water-soluble coiled coils that zipper during fusion and TM helices in both membranes. By analogy to mechanosensitive channels, the force associated with zippering of the water-soluble coiled-coil domain is expected to lead to tilting of the FP helices, promoting interaction with the C-term-TM helices. The energetically unfavorable dehydration of lipid headgroups of opposing bilayers is compensated by thermodynamically favorable interactions between the FP and C-term-TM helices as the coiled coils zipper into the membrane phase, leading to a pore lined by both lipid and protein.

Capture and imaging of a prehairpin fusion intermediate of the paramyxovirus PIV5.

During cell entry, enveloped viruses fuse their viral membrane with a cellular membrane in a process driven by energetically favorable, large-scale conformational rearrangements of their fusion proteins. Structures of the pre- and postfusion states of the fusion proteins including paramyxovirus PIV5 F and influenza virus hemagglutinin suggest that this occurs via two intermediates. Following formation of an initial complex, the proteins structurally elongate, driving a hydrophobic N-terminal "fusion peptide" away from the protein surface into the target membrane. Paradoxically, this first conformation change moves the viral and cellular bilayers further apart. Next, the fusion proteins form a hairpin that drives the two membranes into close opposition. While the pre- and postfusion hairpin forms have been characterized crystallographically, the transiently extended prehairpin intermediate has not been visualized. To provide evidence for this extended intermediate we measured the interbilayer spacing of a paramyxovirus trapped in the process of fusing with solid-supported bilayers. A gold-labeled peptide that binds the prehairpin intermediate was used to stabilize and specifically image F-proteins in the prehairpin intermediate. The interbilayer spacing is precisely that predicted from a computational model of the prehairpin, providing strong evidence for its structure and functional role. Moreover, the F-proteins in the prehairpin conformation preferentially localize to a patch between the target and viral membranes, consistent with the fact that the formation of the prehairpin is triggered by local contacts between F- and neighboring viral receptor-binding proteins (HN) only when HN binds lipids in its target membrane.

Neutralizing epitopes in the membrane-proximal external region of HIV-1 gp41 are influenced by the transmembrane domain and the plasma membrane.

Failure to elicit broadly neutralizing (bNt) antibodies (Abs) against the membrane-proximal external region of HIV-1 gp41 (MPER) reflects the difficulty of mimicking its neutralization-competent structure (NCS). Here, we analyzed MPER antigenicity in the context of the plasma membrane and identified a role for the gp41 transmembrane domain (TM) in exposing the epitopes of three bNt monoclonal Abs (MAbs) (2F5, 4E10, and Z13e1). We transiently expressed DNA constructs encoding gp41 ectodomain fragments fused to either the TM of the platelet-derived growth factor receptor (PDGFR) or the gp41 TM and cytoplasmic tail domain (CT). Constructs encoding the MPER tethered to the gp41 TM followed by a 27-residue CT fragment (MPER-TM1) produced optimal MAb binding. Critical binding residues for the three Nt MAbs were identified using a panel of 24 MPER-TM1 mutants bearing single amino acid substitutions in the MPER; many were previously shown to affect MAb-mediated viral neutralization. Moreover, non-Nt mutants of MAbs 2F5 and 4E10 exhibited a reduction in binding to MPER-TM1 and yet maintained binding to synthetic MPER peptides, indicating that MPER-TM1 better approximates the MPER NCS than peptides. Replacement of the gp41 TM and CT of MPER-TM1 with the PDGFR TM reduced binding by MAb 4E10, but not 2F5, indicating that the gp41 TM plays a pivotal role in orienting the 4E10 epitope, and more globally, in affecting MPER exposure.

Mutation in the Fe-S scaffold protein Isu bypasses frataxin deletion.

Frataxin is a conserved mitochondrial protein deficient in patients with Friedreich's ataxia. Frataxin has been implicated in control of iron homoeostasis and Fe-S cluster assembly. In yeast or human mitochondria, frataxin interacts with components of the Fe-S cluster synthesis machinery, including the cysteine desulfurase Nfs1, accessory protein Isd11 and scaffold protein Isu. In the present paper, we report that a single amino acid substitution (methionine to isoleucine) at position 107 in the mature form of Isu1 restored many deficient functions in Δyfh1 or frataxin-depleted yeast cells. Iron homoeostasis was improved such that soluble/usable mitochondrial iron was increased and accumulation of insoluble/non-usable iron within mitochondria was largely prevented. Cytochromes were returned to normal and haem synthesis was restored. In mitochondria carrying the mutant Isu1 and no frataxin, Fe-S cluster enzyme activities were improved. The efficiency of new Fe-S cluster synthesis in isolated mitochondria was markedly increased compared with frataxin-negative cells, although the response to added iron was minimal. The M107I substitution in the highly conserved Isu scaffold protein is typically found in bacterial orthologues, suggesting that a unique feature of the bacterial Fe-S cluster machinery may be involved. The mechanism by which the mutant Isu bypasses the absence of frataxin remains to be determined, but could be related to direct effects on Fe-S cluster assembly and/or indirect effects on mitochondrial iron availability.

Structural informatics, modeling, and design with an open-source Molecular Software Library (MSL).

We present the Molecular Software Library (MSL), a C++ library for molecular modeling. MSL is a set of tools that supports a large variety of algorithms for the design, modeling, and analysis of macromolecules. Among the main features supported by the library are methods for applying geometric transformations and alignments, the implementation of a rich set of energy functions, side chain optimization, backbone manipulation, calculation of solvent accessible surface area, and other tools. MSL has a number of unique features, such as the ability of storing alternative atomic coordinates (for modeling) and multiple amino acid identities at the same backbone position (for design). It has a straightforward mechanism for extending its energy functions and can work with any type of molecules. Although the code base is large, MSL was created with ease of developing in mind. It allows the rapid implementation of simple tasks while fully supporting the creation of complex applications. Some of the potentialities of the software are demonstrated here with examples that show how to program complex and essential modeling tasks with few lines of code. MSL is an ongoing and evolving project, with new features and improvements being introduced regularly, but it is mature and suitable for production and has been used in numerous protein modeling and design projects. MSL is open-source software, freely downloadable at http://msl-libraries.org. We propose it as a common platform for the development of new molecular algorithms and to promote the distribution, sharing, and reutilization of computational methods.

Identification of interacting hot spots in the beta3 integrin stalk using comprehensive interface design.

Protein-protein interfaces are usually large and complementary surfaces, but specific side chains, representing energetic "hot spots," often contribute disproportionately to binding free energy. We used a computational method, comprehensive interface design, to identify hot spots in the interface between the stalk regions of the ß3 and the complementary αIIb and αv integrin subunits. Using the Rosetta alanine-scanning and design algorithms to predict destabilizing, stabilizing, and neutral mutations in the ß3 region extending from residues Lys(532) through Gly(690), we predicted eight alanine mutations that would destabilize the αIIbß3 interface as well as nine predicted to destabilize the αvß3 interface, by at least 0.3 kcal/mol. The mutations were widely and unevenly distributed, with four between residues 552 and 563 and five between 590 and 610, but none between 565 and 589, and 611 and 655. Further, mutations destabilizing the αvß3 and αIIbß3 interfaces were not identical. The predictions were then tested by introducing selected mutations into the full-length integrins expressed in Chinese hamster ovary cells. Five mutations predicted to destabilize αIIb and ß3 caused fibrinogen binding to αIIbß3, whereas three of four predicted to be neutral or stabilizing did not. Conversely, a mutation predicted to destabilize αvß3, but not αIIbß3 (D552A), caused osteopontin binding to αvß3, but not fibrinogen binding to αIIbß3. These results indicate that stability of the distal stalk interface is involved in constraining integrins in stable, inactive conformations. Further, they demonstrate the ability of comprehensive interface design to identify functionally significant integrin mutations.

Salt bridges: geometrically specific, designable interactions.

Salt bridges occur frequently in proteins, providing conformational specificity and contributing to molecular recognition and catalysis. We present a comprehensive analysis of these interactions in protein structures by surveying a large database of protein structures. Salt bridges between Asp or Glu and His, Arg, or Lys display extremely well-defined geometric preferences. Several previously observed preferences are confirmed, and others that were previously unrecognized are discovered. Salt bridges are explored for their preferences for different separations in sequence and in space, geometric preferences within proteins and at protein-protein interfaces, co-operativity in networked salt bridges, inclusion within metal-binding sites, preference for acidic electrons, apparent conformational side chain entropy reduction on formation, and degree of burial. Salt bridges occur far more frequently between residues at close than distant sequence separations, but, at close distances, there remain strong preferences for salt bridges at specific separations. Specific types of complex salt bridges, involving three or more members, are also discovered. As we observe a strong relationship between the propensity to form a salt bridge and the placement of salt-bridging residues in protein sequences, we discuss the role that salt bridges might play in kinetically influencing protein folding and thermodynamically stabilizing the native conformation. We also develop a quantitative method to select appropriate crystal structure resolution and B-factor cutoffs. Detailed knowledge of these geometric and sequence dependences should aid de novo design and prediction algorithms.

SDR: a database of predicted specificity-determining residues in proteins.

The specificity-determining residue database (SDR database) presents residue positions where mutations are predicted to have changed protein function in large protein families. Because the database pre-calculates predictions on existing protein sequence alignments, users can quickly find the predictions by selecting the appropriate protein family or searching by protein sequence. Predictions can be used to guide mutagenesis or to gain a better understanding of specificity changes in a protein family. The database is available on the web at http://paradox.harvard.edu/sdr.

Energetics of protein-DNA interactions.

Protein-DNA interactions are vital for many processes in living cells, especially transcriptional regulation and DNA modification. To further our understanding of these important processes on the microscopic level, it is necessary that theoretical models describe the macromolecular interaction energetics accurately. While several methods have been proposed, there has not been a careful comparison of how well the different methods are able to predict biologically important quantities such as the correct DNA binding sequence, total binding free energy and free energy changes caused by DNA mutation. In addition to carrying out the comparison, we present two important theoretical models developed initially in protein folding that have not yet been tried on protein-DNA interactions. In the process, we find that the results of these knowledge-based potentials show a strong dependence on the interaction distance and the derivation method. Finally, we present a knowledge-based potential that gives comparable or superior results to the best of the other methods, including the molecular mechanics force field AMBER99.

Predicting specificity-determining residues in two large eukaryotic transcription factor families.

Certain amino acid residues in a protein, when mutated, change the protein's function. We present an improved method of finding these specificity-determining positions that uses all the protein sequence data available for a family of homologous proteins. We study in detail two families of eukaryotic transcription factors, basic leucine zippers and nuclear receptors, because of the large amount of sequences and experimental data available. These protein families also have a clear definition of functional specificity: DNA-binding specificity. We compare our results to three other methods, including the evolutionary trace algorithm and a method that depends on orthology relationships. All of the predictions are compared to the available mutational and crystallographic data. We find that our method provides superior predictions of the known specificity-determining residues and also predicts residue positions within these families that deserve further study for their roles in functional specificity.

Knowledge-based potential for positioning membrane-associated structures and assessing residue-specific energetic contributions.

The complex hydrophobic and hydrophilic milieus of membrane-associated proteins pose experimental and theoretical challenges to their understanding. Here, we produce a nonredundant database to compute knowledge-based asymmetric cross-membrane potentials from the per-residue distributions of C(ß), C(γ) and functional group atoms. We predict transmembrane and peripherally associated regions from genomic sequence and position peptides and protein structures relative to the bilayer (available at http://www.degradolab.org/ez). The pseudo-energy topological landscapes underscore positional stability and functional mechanisms demonstrated here for antimicrobial peptides, transmembrane proteins, and viral fusion proteins. Moreover, experimental effects of point mutations on the relative ratio changes of dual-topology proteins are quantitatively reproduced. The functional group potential and the membrane-exposed residues display the largest energetic changes enabling to detect native-like structures from decoys. Hence, focusing on the uniqueness of membrane-associated proteins and peptides, we quantitatively parameterize their cross-membrane propensity, thus facilitating structural refinement, characterization, prediction, and design.

Functional analysis of the transmembrane domain in paramyxovirus F protein-mediated membrane fusion.

To enter cells, enveloped viruses use fusion-mediating glycoproteins to facilitate the merger of the viral and host cell membranes. These glycoproteins undergo large-scale irreversible refolding during membrane fusion. The paramyxovirus parainfluenza virus 5 mediates membrane merger through its fusion protein (F). The transmembrane (TM) domains of viral fusion proteins are typically required for fusion. The TM domain of F is particularly interesting in that it is potentially unusually long; multiple calculations suggest a TM helix length between 25 and 48 residues. Oxidative cross-linking of single-cysteine substitutions indicates the F TM trimer forms a helical bundle within the membrane. To assess the functional role of the paramyxovirus parainfluenza virus 5 F protein TM domain, alanine scanning mutagenesis was performed. Two residues located in the outer leaflet of the bilayer are critical for fusion. Multiple amino acid substitutions at these positions indicate the physical properties of the side chain play a critical role in supporting or blocking fusion. Analysis of intermediate steps in F protein refolding indicated that the mutants were not trapped at the open stalk intermediate or the prehairpin intermediate. Incorporation of a known F protein destabilizing mutation that causes a hyperfusogenic phenotype restored fusion activity to the mutants. Further, altering the curvature of the lipid bilayer by addition of oleic acid promoted fusion of the F protein mutants. In aggregate, these data indicate that the TM domain plays a functional role in fusion beyond merely anchoring the protein in the viral envelope and that it can affect the structures and steady-state concentrations of the various conformational intermediates en route to the final postfusion state. We suggest that the unusual length of this TM helix might allow it to serve as a template for formation of or specifically stabilize the lipid stalk intermediate in fusion.

Determining functional specificity from protein sequences.

MOTIVATION: Given a large family of homologous protein sequences, many methods can divide the family into smaller groups that correspond to the different functions carried out by proteins within the family. One important problem, however, has been the absence of a general method for selecting an appropriate level of granularity, or size of the groups. RESULTS: We propose a consistent way of choosing the granularity that is independent of the sequence similarity and sequence clustering method used. We study three large, well-investigated protein families: basic leucine zippers, nuclear receptors and proteins with three consecutive C2H2 zinc fingers. Our method is tested against known functional information, the experimentally determined binding specificities, using a simple scoring method. The significance of the groups is also measured by randomizing the data. Finally, we compare our algorithm against a popular method of grouping proteins, the TRIBE-MCL method. In the end, we determine that dividing the families at the proposed level of granularity creates very significant and useful groups of proteins that correspond to the different DNA-binding motifs. We expect that such groupings will be useful in studying not only DNA binding but also other protein interactions.

CoC: a database of universally conserved residues in protein folds.

UNLABELLED: The conservatism of conservatism (CoC) database presents statistically analyzed information about the conservation of residue positions in folds across protein families. AVAILABILITY: On the web at http://kulibin.mit.edu/coc/