BH3-in-groove dimerization initiates and helix 9 dimerization expands Bax pore assembly in membranes.

Pro-apoptotic Bax induces mitochondrial outer membrane permeabilization (MOMP) by forming oligomers through a largely undefined process. Using site-specific disulfide crosslinking, compartment-specific chemical labeling, and mutational analysis, we found that activated integral membrane Bax proteins form a BH3-in-groove dimer interface on the MOM surface similar to that observed in crystals. However, after the α5 helix was released into the MOM, the remaining interface with α2, α3, and α4 helices was rearranged. Another dimer interface was formed inside the MOM by two intersected or parallel α9 helices. Combinations of these interfaces generated oligomers in the MOM. Oligomerization was initiated by BH3-in-groove dimerization, without which neither the other dimerizations nor MOMP occurred. In contrast, α9 dimerization occurred downstream and was required for release of large but not small proteins from mitochondria. Moreover, the release of large proteins was facilitated by α9 insertion into the MOM and localization to the pore rim. Therefore, the BH3-in-groove dimerization on the MOM nucleates the assembly of an oligomeric Bax pore that is enlarged by α9 dimerization at the rim.

A frequent, GxxxG-mediated, transmembrane association motif is optimized for the formation of interhelical Cα-H hydrogen bonds.

Carbon hydrogen bonds between Cα-H donors and carbonyl acceptors are frequently observed between transmembrane helices (Cα-H···O=C). Networks of these interactions occur often at helix-helix interfaces mediated by GxxxG and similar patterns. Cα-H hydrogen bonds have been hypothesized to be important in membrane protein folding and association, but evidence that they are major determinants of helix association is still lacking. Here we present a comprehensive geometric analysis of homodimeric helices that demonstrates the existence of a single region in conformational space with high propensity for Cα-H···O=C hydrogen bond formation. This region corresponds to the most frequent motif for parallel dimers, GASright, whose best-known example is glycophorin A. The finding suggests a causal link between the high frequency of occurrence of GASright and its propensity for carbon hydrogen bond formation. Investigation of the sequence dependency of the motif determined that Gly residues are required at specific positions where only Gly can act as a donor with its "side chain" Hα. Gly also reduces the steric barrier for non-Gly amino acids at other positions to act as Cα donors, promoting the formation of cooperative hydrogen bonding networks. These findings offer a structural rationale for the occurrence of GxxxG patterns at the GASright interface. The analysis identified the conformational space and the sequence requirement of Cα-H···O=C mediated motifs; we took advantage of these results to develop a structural prediction method. The resulting program, CATM, predicts ab initio the known high-resolution structures of homodimeric GASright motifs at near-atomic level.

Backbone dependency further improves side chain prediction efficiency in the Energy-based Conformer Library (bEBL).

Side chain optimization is an integral component of many protein modeling applications. In these applications, the conformational freedom of the side chains is often explored using libraries of discrete, frequently occurring conformations. Because side chain optimization can pose a computationally intensive combinatorial problem, the nature of these conformer libraries is important for ensuring efficiency and accuracy in side chain prediction. We have previously developed an innovative method to create a conformer library with enhanced performance. The Energy-based Library (EBL) was obtained by analyzing the energetic interactions between conformers and a large number of natural protein environments from crystal structures. This process guided the selection of conformers with the highest propensity to fit into spaces that should accommodate a side chain. Because the method requires a large crystallographic data-set, the EBL was created in a backbone-independent fashion. However, it is well established that side chain conformation is strongly dependent on the local backbone geometry, and that backbone-dependent libraries are more efficient in side chain optimization. Here we present the backbone-dependent EBL (bEBL), whose conformers are independently sorted for each populated region of Ramachandran space. The resulting library closely mirrors the local backbone-dependent distribution of side chain conformation. Compared to the EBL, we demonstrate that the bEBL uses fewer conformers to produce similar side chain prediction outcomes, thus further improving performance with respect to the already efficient backbone-independent version of the library.

Structural organization of FtsB, a transmembrane protein of the bacterial divisome.

We report the first structural analysis of an integral membrane protein of the bacterial divisome. FtsB is a single-pass membrane protein with a periplasmic coiled coil. Its heterologous association with its partner FtsL represents an essential event for the recruitment of the late components to the division site. Using a combination of mutagenesis, computational modeling, and X-ray crystallography, we determined that FtsB self-associates, and we investigated its structural organization. We found that the transmembrane domain of FtsB homo-oligomerizes through an evolutionarily conserved interaction interface where a polar residue (Gln 16) plays a critical role through the formation of an interhelical hydrogen bond. The crystal structure of the periplasmic domain, solved as a fusion with Gp7, shows that 30 juxta-membrane amino acids of FtsB form a canonical coiled coil. The presence of conserved Gly residue in the linker region suggests that flexibility between the transmembrane and coiled coil domains is functionally important. We hypothesize that the transmembrane helices of FtsB form a stable dimeric core for its association with FtsL into a higher-order oligomer and that FtsL is required to stabilize the periplasmic domain of FtsB, leading to the formation of a complex that is competent for binding to FtsQ, and to their consequent recruitment to the divisome. The study provides an experimentally validated structural model and identifies point mutations that disrupt association, thereby establishing important groundwork for the functional characterization of FtsB in vivo.

An energy-based conformer library for side chain optimization: improved prediction and adjustable sampling.

Side chain optimization is a fundamental component of protein modeling applications such as docking, structural prediction, and design. In these applications side chain flexibility is often provided by rotamer or conformer libraries, which are collections of representative side chain conformations. Here we demonstrate that the sampling provided by the library can be substantially improved by adding an energetic criterion to its creation. The result of the new procedure is the Energy-Based library, a conformer library selected according to the propensity of its elements to fit energetically into natural protein environments. The new library performs outstandingly well in side chain optimization, producing structures with significantly lower energies and resulting in improved side chain conformation prediction. In addition, because the library was created as an ordered list, its size can be adjusted to any desired level. This feature provides unprecedented versatility in tuning sampling. It allows to precisely balance the number of conformers required by each amino acid type, equalizing their chances to fit into structural environments. It also allows to scale the amount of sampling to the specific requirement of any given side optimization problem. A rotameric version of the library was also produced with the same method to support applications that require a dihedral-only description of side chain conformation. The libraries are available at http://seneslab.org/EBL.

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.