Patrick Barth

Toward high-resolution prediction and design of transmembrane helical protein structures

The prediction and design at the atomic level of membrane protein structures and interactions is a critical but unsolved challenge. To address this problem, we have developed an all-atom physical model that describes intraprotein and protein–solvent interactions in the membrane environment. We evaluated the ability of the model to recapitulate the energetics and structural specificities of polytopic membrane proteins by using a battery of in silico prediction and design tests. First, in side-chain packing and design tests, the model successfully predicts the side-chain conformations at 73% of nonexposed positions and the native amino acid identities at 34% of positions in naturally occurring membrane proteins. Second, the model predicts significant energy gaps between native and nonnative structures of transmembrane helical interfaces and polytopic membrane proteins. Third, distortions in transmembrane helices are successfully recapitulated in docking experiments by using fragments of ideal helices judiciously defined around helical kinks. Finally, de novo structure prediction reaches near-atomic accuracy (<2.5 Å) for several small membrane protein domains (<150 residues). The success of the model highlights the critical role of van der Waals and hydrogen-bonding interactions in the stability and structural specificity of membrane protein structures and sets the stage for the high-resolution prediction and design of complex membrane protein architectures.

Prediction of membrane protein structures with complex topologies using limited constraints

Reliable structure-prediction methods for membrane proteins are important because the experimental determination of high-resolution membrane protein structures remains very difficult, especially for eukaryotic proteins. However, membrane proteins are typically longer than 200 aa and represent a formidable challenge for structure prediction. We have developed a method for predicting the structures of large membrane proteins by constraining helix–helix packing arrangements at particular positions predicted from sequence or identified by experiments. We tested the method on 12 membrane proteins of diverse topologies and functions with lengths ranging between 190 and 300 residues. Enforcing a single constraint during the folding simulations enriched the population of near-native models for 9 proteins. In 4 of the cases in which the constraint was predicted from the sequence, 1 of the 5 lowest energy models was superimposable within 4 Å on the native structure. Near-native structures could also be selected for heme-binding and pore-forming domains from simulations in which pairs of conserved histidine-chelating hemes and one experimentally determined salt bridge were constrained, respectively. These results suggest that models within 4 Å of the native structure can be achieved for complex membrane proteins if even limited information on residue-residue interactions can be obtained from protein structure databases or experiments.

The structure of a receptor with two associating transmembrane domains on the cell surface: integrin αIIbβ3

Structures of intact receptors with single-pass transmembrane domains are essential to understand how extracellular and cytoplasmic domains regulate association and signaling through transmembrane domains. A chemical and computational method to determine structures of the membrane regions of such receptors on the cell surface is developed here and validated with glycophorin A. An integrin heterodimer structure reveals association over most of the lengths of the alpha and beta transmembrane domains and shows that the principles governing association of hetero and homo transmembrane dimers differ. A turn at the Gly of the juxtamembrane GFFKR motif caps the alpha TM helix and brings the two Phe of GFFKR into the alpha/beta interface. A juxtamembrane Lys residue in beta also has an important role in the interface. The structure shows how transmembrane association/dissociation regulates integrin signaling. A joint ectodomain and membrane structure shows that substantial flexibility between the extracellular and TM domains is compatible with TM signaling.

Binding of aldose reductase inhibitors : Correlation of crystallographic and mass spectrometric studies

Aldose reductase is a NADP(H)-dependent enzyme, believed to be strongly implicated in the development of degenerative complications of Diabetes Mellitus. The search for specific inhibitors of this enzyme has thus become a major pharmaceutic challenge. In this study, we applied both X-ray crystallography and mass spectrometry to characterize the interactions between aldose reductase and four representative inhibitors: AminoSNM, Imirestat, LCB3071, and IDD384. If crystallography remains obviously the only way to get an extensive description of the contacts between an inhibitor and the enzymatic site, the duration of the crystallographic analysis makes this technique incompatible with high throughput screenings of inhibitors. On the other hand, dissociation experiments monitored by mass spectrometry permitted us to evaluate rapidly the relative gas-phase stabilities of the aldose reductase-inhibitor noncovalent complexes. In our experiments, dissociation in the gas-phase was provoked by increasing the accelerating voltage of the ions (Vc) in the source-analyzer interface region: the Vc value needed to dissociate 50% of the noncovalent complex initially present (Vc50) was taken as a gas-phase stability parameter of the enzyme-inhibitor complex. Interestingly, the Vc50 were found to correlate with the energy of the electrostatic and H-bond interactions involved in the contact aldose reductase/inhibitor (Eel-H), computed from the crystallographic model. This finding may be specially interesting in a context of drug development. Actually, during a drug design optimization phase, the binding of the drug to the target enzyme is often optimized by modifying its interatomic electrostatic and H-bond contacts, because they usually depend on a single atom change on the drug, and are easier to introduce than the hydrophobic interactions. Therefore, the Vc50 may help to monitor the chemical modifications introduced in new inhibitors. X-ray crystallography is clearly needed to get the details of the contacts and to rationalize the design. Nevertheless, once the cycle of chemical modification is engaged, mass spectrometry can be used to select a priori the drug candidates which are worthy of further crystallographic investigation. We thus propose to use the two techniques in a complementary way, to improve the screening of large collections of inhibitors.

Community-wide assessment of GPCR structure modelling and ligand docking

Recent breakthroughs in the determination of the crystal structures of G protein-coupled receptors (GPCRs) have provided new opportunities for structure-based drug design strategies targeting this protein family. With the aim of evaluating the current status of GPCR structure prediction and ligand docking, a community-wide, blind prediction assessment — GPCR Dock 2008 — was conducted in coordination with the publication of the crystal structure of the human adenosine A2A receptor bound to the ligand ZM241385. Twenty-nine groups submitted 206 structural models before the release of the experimental structure, which were evaluated for the accuracy of the ligand binding mode and the overall receptor model compared with the crystal structure. This analysis highlights important aspects for success and future development, such as accurate modelling of structurally divergent regions and use of additional biochemical insight such as disulphide bridges in the extracellular loops.

Essential Role of a Single Arginine of Photosystem I in Stabilizing the Electron Transfer Complex with Ferredoxin

PsaE is one of the photosystem I subunits involved in ferredoxin binding. The central role of arginine 39 of this 8-kDa peripheral polypeptide has been established by a series of mutations. The neutral substitution R39Q leads to a 250-fold increase of the dissociation constant K d of the photosystem I-ferredoxin complex, as large as the increase induced by PsaE deletion. At pH 8.0, this K d value strongly depends on the charge of the residue substituting Arg-39: 0.22 μmfor wild type, 1.5 μm for R39K, 56 μm for R39Q, and more than 100 μm for R39D. The consequences of arginine 39 substitution for the titratable histidine were analyzed as a function of pH. The K d value of R39H is increased 140 times at pH 8.0 but only 5 times at pH 5.8, which is assigned to the protonation of histidine at low pH. In the mutant R39Q, the association rate of ferredoxin was decreased 3-fold compared with wild type, whereas an 80-fold increase is calculated for the dissociation rate. We propose that a major contribution of PsaE is to provide a prominent positive charge at position 39 for controlling the electrostatic interaction and lifetime of the complex with ferredoxin.

Accurate, conformation-dependent predictions of solvent effects on protein ionization constants

Predicting how aqueous solvent modulates the conformational transitions and influences the pKa values that regulate the biological functions of biomolecules remains an unsolved challenge. To address this problem, we developed FDPB_MF, a rotamer repacking method that exhaustively samples side chain conformational space and rigorously calculates multibody protein–solvent interactions. FDPB_MF predicts the effects on pKa values of various solvent exposures, large ionic strength variations, strong energetic couplings, structural reorganizations and sequence mutations. The method achieves high accuracy, with root mean square deviations within 0.3 pH unit of the experimental values measured for turkey ovomucoid third domain, hen lysozyme, Bacillus circulans xylanase, and human and Escherichia coli thioredoxins. FDPB_MF provides a faithful, quantitative assessment of electrostatic interactions in biological macromolecules.

Naturally evolved G protein-coupled receptors adopt metastable conformations

A wide range of membrane receptors signal through conformational changes, and the resulting protein conformational flexibility often hinders their structural studies. Because the determinants of membrane receptor conformational stability are still poorly understood, identifying a minimal set of perturbations stabilizing a membrane protein in a given conformation remains a major challenge in membrane protein structure determination. We present a novel approach integrating bioinformatics, computational design and experimental techniques that identifies and stabilizes metastable receptor regions. When applied to the beta1-adrenergic receptor, the method generated 13 novel receptor variants stabilized in the intended inactive state among which two exhibit an apparent thermostability higher than WT and M23 (a receptor variant previously stabilized by extensive scanning mutagenesis) by more than 30 °C and 11 °C, respectively. Targeted regions involve nonconserved unsatisfied polar residues or exhibit significant packing defects, features found in all class A G protein-coupled receptor structures. These findings suggest that natural G protein-coupled receptor sequences have evolved to be conformationally metastable through the design of suboptimal polar and van der Waals tertiary interactions. Given sufficiently accurate structural models, our approach should prove useful for designing stabilized variants of many uncharacterized membrane receptors.

Production of crystals of human aldose reductase with very high resolution diffraction

As the action of human aldose reductase (hAR) is thought to be linked to the pathogenesis of diabetic complications, much effort has been directed towards the analysis of the catalytic mechanism and the development of specific inhibitors. Here, the crystallization of recombinant hAR with its cofactor NADP+ at 277 K in the presence of the precipitating agent PEG 6000 is reported. The crystals diffract to high resolution (1.1 A) and belong to the P21 space group with unit-cell parameters a = 49.97, b = 67.14, c = 48. 02 A, beta = 92.2 degrees with one molecule per asymmetric unit. Seleno-substituted hAR crystals were also produced and diffract to 1. 7 A on a conventional X-ray source.

Survivin-specific T cell receptor targets tumor but not T cells

Survivin is a tumor-associated antigen (TAA) that inhibits apoptosis and is widely overexpressed in cancer cells; therefore, survivin has potential as a target for cancer immunotherapy. Application of HLA-A2-restricted survivin-specific T cell receptors (TCRs) isolated from allogeneic HLA-mismatched TCR repertoires has, however, been impeded by the inability of these TCRs to distinguish healthy cells expressing low levels of survivin from cancer cells with high survivin expression levels. Here, we identified an HLA-A2-restricted survivin-specific TCR isolated from autologous TCR repertoires that targets tumor cells in vitro and in vivo but does not cause fratricidal toxicity. Molecular modeling of the TCR-peptide-HLA ternary complexes and alanine scanning revealed that the autologously derived TCRs had tighter interactions with the survivin peptide than did fratricidal TCRs. Similar recognition patterns were observed among 7 additional TAA-specific TCRs isolated from allogeneic versus autologous repertoires. Together, the results from this study indicate that maximal peptide recognition is key for TCR selectivity and likely critical for reducing unwanted off-target toxicities. Moreover, isolating TCRs from autologous repertoires to maximize TCR selectivity has potential as a useful strategy to identify and select other shared tumor- and self-antigen-specific TCRs and ensure selective antitumor activity.