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.

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.

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).

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.

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.

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.

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.

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).

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.

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.

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.

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.

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.