Predicting the complex structure and functional motions of the outer membrane transporter and signal transducer FecA

Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model

The CABS coarse-grained model is a well-established tool for modeling globular proteins (predicting their structure, dynamics, and interactions). Here we introduce an extension of the CABS representation and force field (CABS-membrane) to the modeling of the effect of the biological membrane environment on the structure of membrane proteins. We validate the CABS-membrane model in folding simulations of 10 short helical membrane proteins not using any knowledge about their structure. The simulations start from random protein conformations placed outside the membrane environment and allow for full flexibility of the modeled proteins during their spontaneous insertion into the membrane. In the resulting trajectories, we have found models close to the experimental membrane structures. We also attempted to select the correctly folded models using simple filtering followed by structural clustering combined with reconstruction to the all-atom representation and all-atom scoring. The CABS-membrane model is a promising approach for further development toward modeling of large protein-membrane systems.

Structure prediction of the second extracellular loop in G-protein-coupled receptors

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.

Dynamics of endoglucanase catalytic domains: implications towards thermostability

Thermostable endoglucanases play a crucial role in the production of biofuels to breakdown plant cellulose. Analyzing their structure-dynamics relationship can inform about the origins of their thermostability. Although tertiary structures of many endoglucanase proteins are available, the relationship between thermostability, structure, and dynamics is not explored fully. We have generated elastic network models for thermostable and mesostable endoglucanases with the (αβ)₈ fold in substrate bound and unbound states. The comparative analyses shed light on the relation between protein dynamics, thermostability, and substrate binding. We observed specific differences in the dynamic behavior of catalytic residues in slow modes: while both the nucleophile and the acid/base donor residues show positively correlated motions in the thermophile, their dynamics is uncoupled in the mesophile. Our proof-of-concept comparison study suggests that global dynamics can be harnessed to further our understanding of thermostability.

Gating and conduction of nano-channel forming proteins: a computational approach

Monitoring conformational changes in ion channels is essential to understand their gating mechanism. Here, we explore the structural dynamics of four outer membrane proteins with different structures and functions in the slowest nonzero modes of vibration. Normal mode analysis was performed on the modified elastic network model of channel in the membrane. According to our results, when membrane proteins were analyzed in the dominant mode, the composed pores, TolC and α-hemolysin showed large motions at the intramembrane β-barrel region while, in other porins, OmpA and OmpF, largest motions observed in the region of external flexible loops. A criterion based on equipartition theorem was used to measure the possible amplitude of vibration in channel forming proteins. The current approach complements theoretical and experimental techniques including HOLE, Molecular Dynamics (MD), and voltage clamp used to address the channel's structure and dynamics and provides the means to conduct a theoretical simultaneous study of the structure and function of the channel. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:3.

Structural and dynamical insights into the opening mechanism of P. aeruginosa OprM channel

Originally described in bacteria, drug transporters are now recognized as major determinants in antibiotics resistance. For Gram-negative bacteria, the reversible assembly consisting of an inner membrane protein responsible for the active transport, a periplasmic protein, and an exit outer membrane channel achieves transport. The opening of the outer membrane protein OprM from Pseudomonas aeruginosa was modeled through normal mode analysis starting from a new X-ray structure solved at 2.4 A resolution in P2(1)2(1)2(1) space group. The three monomers are not linked by internal crystallographic symmetries highlighting the possible functional differences. This structure is closed at both ends, but modeling allowed for an opening that is not reduced to the classically proposed "iris-like mechanism."

Asymmetric conformational flexibility in the ATP-binding cassette transporter HI1470/1

Putative metal-chelate-type ABC transporter HI1470/1 is homologous with vitamin B(12) importer BtuCD but exhibits a distinct inward-facing conformation in contrast to the outward-facing conformation of BtuCD. Normal-mode analysis of HI1470/1 reveals the intrinsic asymmetric conformational flexibility in this transporter and demonstrates that the transition from the inward-facing to the outward-facing conformation is realized through the asymmetric motion of individual subunits of the transporter. This analysis suggests that the asymmetric arrangement of the BtuC dimer in the crystal structure of the BtuCD-F complex represents an intermediate state relating HI1470/1 and BtuCD. Furthermore, a twisting motion between transmembrane domains and nucleotide-binding domains encoded in the lowest-frequency normal mode of this type of importer is found to contribute to the conformational transitions during the whole cycle of substrate transportation. A more complete translocation mechanism of the BtuCD type importer is proposed.