Search

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor., ((c) 2007 Wiley-Liss, Inc.)

Structural studies of symmetric homo-oligomers provide mechanistic insights into their roles in essential biological processes, including cell signaling and cellular regulation. This paper presents a novel algorithm for homo-oligomeric structure determination, given the subunit structure, that is both complete, in that it evaluates all possible conformations, and data-driven, in that it evaluates conformations separately for consistency with experimental data and for quality of packing. Completeness ensures that the algorithm does not miss the native conformation, and being data-driven enables it to assess the structural precision possible from data alone. Our algorithm performs a branch-and-bound search in the symmetry configuration space, the space of symmetry axis parameters (positions and orientations) defining all possible C(n) homo-oligomeric complexes for a given subunit structure. It eliminates those symmetry axes inconsistent with intersubunit nuclear Overhauser effect (NOE) distance restraints and then identifies conformations representing any consistent, well-packed structure to within a user-defined similarity level. For the human phospholamban pentamer in dodecylphosphocholine micelles, using the structure of one subunit determined from a subset of the experimental NMR data, our algorithm identifies a diverse set of complex structures consistent with the nine intersubunit NOE restraints. The distribution of determined structures provides an objective characterization of structural uncertainty: backbone RMSD to the previously determined structure ranges from 1.07 to 8.85 A, and variance in backbone atomic coordinates is an average of 12.32 A(2). Incorporating vdW packing reduces structural diversity to a maximum backbone RMSD of 6.24 A and an average backbone variance of 6.80 A(2). By comparing data consistency and packing quality under different assumptions of oligomeric number, our algorithm identifies the pentamer as the most likely oligomeric state of phospholamban, demonstrating that it is possible to determine the oligomeric number directly from NMR data. Additional tests on a number of homo-oligomers, from dimer to heptamer, similarly demonstrate the power of our method to provide unbiased determination and evaluation of homo-oligomeric complex structures., (Proteins 2006. (c) 2006 Wiley-Liss, Inc.)

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures., (2006 Wiley-Liss, Inc.)

Understanding the energetic factors determining the positioning and orientation of single-helical peptides in membranes is of fundamental interest in structural biology. Here, a simple 5-slab continuum dielectric model for the membrane is examined that distinguishes between the solvent, headgroup, and core regions. An analytical solution for the electrostatic solvation of a single dipole and an all-atom model of N-methylacetamide are used to demonstrate the effect of the dielectric boundaries in the system on peptide dipole orientation. The dipole orientation energy is shown to dominate the electrostatic solvation energy of a polyalanine helix in the membrane. With an additional surface-area-dependent term to account for the cavity formation in the aqueous region, the continuum electrostatics description is used to examine several helical peptides, the atoms of which are explicitly represented with a molecular mechanics force field. The experimentally determined tilt angles of a number of peptides of alternating alanine and leucine residues, and of glycophorin and melittin, are accurately reproduced by the model. The factors determining the tilt angles and their fluctuations are analyzed. The tilt angles of the simpler peptides are found to increase approximately linearly with peptide length; this effect is also rationalized. The analysis and model presented here provide a step toward the prediction of helical membrane protein structure., (Copyright 2005 Wiley-Liss, Inc.)

Specific non-covalent interactions between transmembrane (TM) alpha-helices are important in a variety of biological processes. Experimental and computational studies have shown that van der Waals interactions play an important role in the tight packing between TM alpha-helices, although polar interactions can also be important in some instances. Based on the assumption that van der Waals interaction alone is sufficient for a meso-scale (residue-scale) description of the interaction between TM alpha-helices, we have designed a novel residue-scale scoring function for modeling structures of oligomers of TM alpha-helices. We first calculated atomistic van der Waals interaction energies between two amino acids, X and Y, of a pair of parallel alpha-helices, glycine-X-glycine and glycine-Y-glycine and compiled them according to three variables, the distance between the two C(alpha) atoms and the rotational angles of X and Y about their helical axes. Upon averaging over the rotational angles, we obtained one-dimensional interaction energy profiles that are functions of the distance between C(alpha) atoms only. Each of the interaction energy profiles was fitted with a generic fitting function of the distance between C(alpha) atoms, yielding analytical scoring functions for all possible amino acid pairs. For glycophorin A, neu/erbB-2, and phospholamban, lowest-energy conformations obtained through exhaustive scanning of the entire conformational space using the scoring functions were compatible with available experimental data., ((c) 2004 Wiley-Liss, Inc.)

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., (Copyright 2003 Wiley-Liss, Inc.)

We present the results of applying a novel knowledge-based method (FILM) to the prediction of small membrane protein structures. The basis of the method is the addition of a membrane potential to the energy terms (pairwise, solvation, steric, and hydrogen bonding) of a previously developed ab initio technique for the prediction of tertiary structure of globular proteins (FRAGFOLD). The method is based on the assembly of supersecondary structural fragments taken from a library of highly resolved protein structures using a standard simulated annealing algorithm. The membrane potential has been derived by the statistical analysis of a data set made of 640 transmembrane helices with experimentally defined topology and belonging to 133 proteins extracted from the SWISS-PROT database. Results obtained by applying the method to small membrane proteins of known 3D structure show that the method is able to predict, at a reasonable accuracy level, both the helix topology and the conformations of these proteins., (Copyright 2003 Wiley-Liss, Inc.)

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., (Copyright 2002 Wiley-Liss, Inc.)

Biological membrane fusion involves a highly precise and ordered set of protein-protein interactions. Synaptobrevin is a key player in this process. Mutagenesis studies of its single transmembrane segment suggest that it dimerizes in a sequence specific manner. Using the computational methods developed for the successful structure prediction of the glycophorin A transmembrane dimer, we have calculated a structural model for the synaptobrevin dimer. Our computational search yields a well-populated cluster of right-handed structures consistent with the experimentally determined dimerization motif. The three-dimensional structure contains an interface formed primarily by leucine and isoleucine side-chain atoms and has no interhelical hydrogen bonds. The model is the first three-dimensional picture of the synaptobrevin transmembrane dimer and provides a basis for further focused experimentation on its structure and association thermodynamics., (Copyright 2001 Wiley-Liss, Inc.)

The importance of accurately modeling membrane proteins cannot be overstated, in lieu of the difficulties in solving their structures experimentally. Often, however, modeling procedures (e.g., global searching molecular dynamics) generate several possible candidates rather then pointing to a single model. Herein we present a new approach to select among candidate models based on the general hypothesis that silent amino acid substitutions, present in variants identified from evolutionary conservation data or mutagenesis analysis, do not affect the stability of a native structure but may destabilize the non-native structures also found. The proof of this hypothesis has been tested on the alpha-helical transmembrane domains of two homodimers, human glycophorin A and human CD3-zeta, a component of the T-cell receptor. For both proteins, only one structure was identified using all the variants. For glycophorin A, this structure is virtually identical to the structure determined experimentally by NMR. We present a model for the transmembrane domain of CD3-zeta that is consistent with predictions based on mutagenesis, homology modeling, and the presence of a disulfide bond. Our experiments suggest that this method allows the prediction of transmembrane domain structure based only on widely available evolutionary conservation data.

Stone (St(a)) is a variant antigen carried on human erythrocyte MNSs glycophorins (GPSt(a)) that are genetically associated with splicing mutations in GPA genes or with hybrid formation between GPA and GPB genes. Here we identify the first and rare gene conversion event in which GPE, the third member of the family, recombined with GPA, giving rise to a GPA-E-A hybrid gene encoding the St(a) antigen. Western blot detected expression in the proband of both GPA and GPSt(a) on the plasma membrane. Southern blot showed a new restriction fragment from the GPSt(a) gene, indicating an altered exon III-intron 3 junction. Sequencing of RT-PCR products identified one full-length and two shortened glycophorin cDNAs. The shortened forms were derived from GPSt(a) lacking one (exon III) and two exons (exon III and IV), respectively. To define the molecular basis for exon skipping, the genomic region spanning exon III of the GPSt(a) gene was amplified and sequenced. This revealed transfer from GPE to GPA of a DNA segment containing the pseudoexon III and its silent donor splice site. Thus, the inactivation of GPA exon III by conversion of a silent GPE donor splice site portrays a new molecular mechanism for St(a) antigen expression in human erythrocytes.

A more global search method, using fewer assumptions, has been used to predict the structure of the dimeric transmembrane region of the protein glycophorin A. The resulting model significantly differs from that previously determined. In particular, the arrangement between the two transmembrane helices is now more symmetric resulting in improved interaction energies and an increased buried surface area. An increase in the van der Waals interaction energy due to tighter packing compensates for the loss of the interhelical hydrogen bond observed between Thr-87 of each helix in the previous model.