Deciphering the folding kinetics of transmembrane helical proteins

Entangled Motifs in Membrane Protein Structures

Entangled motifs are found in one-third of protein domain structures, a reference set that contains mostly globular proteins. Their properties suggest a connection with co-translational folding. Here, we wish to investigate the presence and properties of entangled motifs in membrane protein structures. From existing databases, we build a non-redundant data set of membrane protein domains, annotated with the monotopic/transmembrane and peripheral/integral labels. We evaluate the presence of entangled motifs using the Gaussian entanglement indicator. We find that entangled motifs appear in one-fifth of transmembrane and one-fourth of monotopic proteins. Surprisingly, the main features of the distribution of the values of the entanglement indicator are similar to the reference case of general proteins. The distribution is conserved across different organisms. Differences with respect to the reference set emerge when considering the chirality of entangled motifs. Although the same chirality bias is found for single-winding motifs in both membrane and reference proteins, the bias is reversed, strikingly, for double-winding motifs only in the reference set. We speculate that these observations can be rationalized in terms of the constraints exerted on the nascent chain by the co-translational bio-genesis machinery, which is different for membrane and globular proteins.

STRUCTURE PREDICTION OF HELICAL TRANSMEMBRANE PROTEINS AT TWO LENGTH SCALES

As the first step toward a multi-scale, hierarchical computational approach for membrane protein structure prediction, the packing of transmembrane helices was modeled at the residue and atom levels, respectively. For predictions at the residue level, the helix-helix and helix-membrane interactions were described by a set of knowledge-based energy functions. For predictions at the atom level, CHARMM19 force field was used. To facilitate the system to overcome energy barriers, the Wang–Landau method was employed, where a random walk is performed in the energy space with a uniform probability. Native-like structures were predicted at both levels for two model systems, each of which consists of two transmembrane helices. Interestingly, consistent results were obtained from simulations at the residue and atom levels for the same system, strongly suggesting the feasibility of a hierarchical approach for membrane protein structure predictions.

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.

Probing the Kinetics of Membrane-Mediated Helix Folding

The kinetics of peptide-membrane association have been studied previously using stopped-flow tryptophan fluorescence; however, such experiments do not directly report the coil-to-helix transition process, which is a hallmark of peptide-membrane interaction. Herein, we report a new method for directly assessing the kinetics of the helix formation accompanied by the peptide-membrane association. This method is based on the technique of fluorescence resonance energy transfer (FRET) and an amino acid FRET pair, p-cyano-L-phenylalanine and tryptophan. To demonstrate the utility of this method, we have studied the membrane-mediated helix folding dynamics of a mutant of magainin 2, an antibiotic peptide found in the skin of the African clawed frog, Xenopus laevis. Our results indicate that the coil-to-helix transition occurs during the binding of the peptide to the lipid vesicle (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)], 3:1, wt/wt) but prior to the full insertion of the peptide into the hydrophobic region of the lipid bilayers.

Elucidation of the disulfide-folding pathway of hirudin by a topology-based approach

A theoretical model for the folding of proteins containing disulfide bonds is introduced. The model exploits the knowledge of the native state to favor the progressive establishment of native interactions. At variance with traditional approaches based on native topology, not all native bonds are treated in the same way; in particular, a suitable energy term is introduced to account for the special strength of disulfide bonds, as well as their ability to undergo intramolecular reshuffling. The model thus possesses the minimal ingredients necessary to investigate the much debated issue of whether the refolding process occurs through partially structured intermediates with native or non-native disulfide bonds. This strategy is applied to a context of particular interest, the refolding process of hirudin, a thrombin-specific protease inhibitor, for which conflicting folding pathways have been proposed. We show that the only two parameters in the model (temperature and disulfide strength) can be tuned to reproduce well a set of experimental transitions between species with different number of formed disulfides. This model is then used to provide a characterization of the folding process and a detailed description of the species involved in the rate-limiting step of hirudin refolding.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Experimental and computational studies of determinants of membrane-protein folding

Recent experiments and analysis have demonstrated the important roles of hydrogen bonding and polar-polar interactions in driving transmembrane helix-helix association. Further details of the energetics of helix-helix association and interhelical packing geometry are being mapped out. Many sequence motifs and a few spatial motifs promoting helical associations have also been identified.

Optimal potentials for predicting inter-helical packing in transmembrane proteins

A set of pairwise contact potentials between amino acid residues in transmembrane helices was determined from the known native structure of the transmembrane protein (TMP) bacteriorhodopsin by the method of perceptron learning, using Monte Carlo dynamics to generate suitable "decoy" structures. The procedure of finding these decoys is simpler than for globular proteins, since it is reasonable to assume that helices behave as independent, stable objects and, therefore, the search in the conformational space is greatly reduced. With the learnt potentials, the association of the helices in bacteriorhodopsin was successfully simulated. The folding of a second TMP (the helix-dimer glycophorin A) was then accomplished with only a refinement of the potentials from a small number of decoys.

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.