Development of Microfluidic Devices with the Use of Nanotechnology to Aid in the Analysis of Biological Systems Including Membrane Protein Separation, Single Cell Analysis, and Genetic Markers

Computation technology has dramatically changed the world around us; you can hardly find an area where cell phones have not saturated the market, yet there is a significant lack of breakthroughs in the development to integrate the computer with biological environments. This is largely the result of the incompatibility of the materials used in both environments; biological environments and experiments tend to need aqueous environments. To help aid in these development chemists, engineers, physicists and biologists have begun to develop microfluidics to help bridge this divide. Unfortunately, the microfluidic devices required large external support equipment to run the device. This thesis presents a series of several microfluidic methods that can help integrate engineering and biology by exploiting nanotechnology to help push the field of microfluidics back to its intended purpose, small integrated biological and electrical devices. I demonstrate this goal by developing different methods and devices to (1) separate membrane bound proteins with the use of microfluidics, (2) use optical technology to make fiber optic cables into protein sensors, (3) generate new fluidic devices using semiconductor material to manipulate single cells, and (4) develop a new genetic microfluidic based diagnostic assay that works with current PCR methodology to provide faster and cheaper results. All of these methods and systems can be used as components to build a self-contained biomedical device.

A study of protein targeting reveals insights into mitigating protein aggregation

A unique chloroplast Signal Recognition Particle (SRP) in green plants is primarily dedicated to the post-translational targeting of light harvesting chlorophyll-a/b binding (LHC) proteins. Our study of the thermodynamics and kinetics of the GTPases of the system demonstrates that GTPase complex assembly and activation are highly coupled in the chloroplast GTPases, suggesting they may forego the GTPase activation step as a key regulatory point. This reflects adaptations of the chloroplast SRP to the delivery of their unique substrate protein. Devotion to one highly hydrophobic family of proteins also may have allowed the chloroplast SRP system to evolve an efficient chaperone in the cpSRP43 subunit. To understand the mechanism of disaggregation, we showed that LHC proteins form micellar, disc-shaped aggregates that present a recognition motif (L18) on the aggregate surface. Further molecular genetic and structure-activity analyses reveal that the action of cpSRP43 can be dissected into two steps: (i) initial recognition of L18 on the aggregate surface; and (ii) aggregate remodeling, during which highly adaptable binding interactions of cpSRP43 with hydrophobic transmembrane domains of the substrate protein compete with the packing interactions within the aggregate. We also tested the adaptability of cpSRP43 for alternative substrates, specifically in attempts to improve membrane protein expression and inhibition of amyloid beta fibrillization. These preliminary results attest to cpSRP43’s potential as a molecular chaperone and provides the impetus for further engineering endeavors to address problems that stem from protein aggregation.

The Get3 ATPase drives unidirectional targeting of tail-anchored membrane proteins

Efficient and accurate localization of membrane proteins is essential to all cells and requires a complex cascade of interactions between protein machineries. This is exemplified in the recently discovered Guided Entry of Tail-anchored protein pathway, in which the central targeting factor Get3 must sequentially interact with three distinct binding partners (Get4, Get1 and Get2) to ensure the targeted delivery of Tail-anchored proteins to the endoplasmic reticulum membrane. To understand the molecular and energetic principles that provide the vectorial driving force of these interactions, we used a quantitative fluorescence approach combined with mechanistic enzymology to monitor the effector interactions of Get3 at each stage of Tail-anchored protein targeting. We show that nucleotide and membrane protein substrate generate a gradient of interaction energies that drive the cyclic and ordered transit of Get3 from Get4 to Get2 and lastly to Get1. These data also define how the Get3/Tail-anchored complex is captured, handed over, and disassembled by the Get1/2 receptor at the membrane, and reveal a novel role for Get4/5 in recycling Get3 from the endoplasmic reticulum membrane at the end of the targeting reaction. These results provide general insights into how complex cascades of protein interactions are coordinated and coupled to energy inputs in biological systems.

Understanding Co-translational Protein Targeting and Lithium Dendrite Formation through Free Energy Simulations and Coarse-Grained Models

We describe the application of alchemical free energy methods and coarse-grained models to study two key problems: (i) co-translational protein targeting and insertion to direct membrane proteins to the endoplasmic reticulum for proper localization and folding, (ii) lithium dendrite formation during recharging of lithium metal batteries. We show that conformational changes in the signal recognition particle, a central component of the protein targeting machinery, confer additional specificity during the the recognition of signal sequences. We then develop a three-dimensional coarse-grained model to study the long-timescale dynamics of membrane protein integration at the translocon and a framework for the calculation of binding free energies between the ribosome and translocon. Finally, we develop a coarse-grained model to capture the dynamics of lithium deposition and dissolution at the electrode interface with time-dependent voltages to show that pulse plating and reverse pulse plating methods can mitigate dendrite growth.

Regulation of mitochondrial division by the Drp1 receptors

Mitochondria can remodel their membranes by fusing or dividing. These processes are required for the proper development and viability of multicellular organisms. At the cellular level, fusion is important for mitochondrial Ca2+ homeostasis, mitochondrial DNA maintenance, mitochondrial membrane potential, and respiration. Mitochondrial division, which is better known as fission, is important for apoptosis, mitophagy, and for the proper allocation of mitochondria to daughter cells during cellular division.

The functions of proteins involved in fission have been best characterized in the yeast model organism Sarccharomyces cerevisiae. Mitochondrial fission in mammals has some similarities. In both systems, a cytosolic dynamin-like protein, called Dnm1 in yeast and Drp1 in mammals, must be recruited to the mitochondrial surface and polymerized to promote membrane division. Recruitment of yeast Dnm1 requires only one mitochondrial outer membrane protein, named Fis1. Fis1 is conserved in mammals, but its importance for Drp1 recruitment is minor. In mammals, three other receptor proteins—Mff, MiD49, and MiD51—play a major role in recruiting Drp1 to mitochondria. Why mammals require three additional receptors, and whether they function together or separately, are fundamental questions for understanding the mechanism of mitochondrial fission in mammals.

We have determined that Mff, MiD49, or MiD51 can function independently of one another to recruit Drp1 to mitochondria. Fis1 plays a minor role in Drp1 recruitment, suggesting that the emergence of these additional receptors has replaced the system used by yeast. Additionally, we found that Fis1/Mff and the MiDs regulate Drp1 activity differentially. Fis1 and Mff promote constitutive mitochondrial fission, whereas the MiDs activate recruited Drp1 only during loss of respiration.

To better understand the function of the MiDs, we have determined the atomic structure of the cytoplasmic domain of MiD51, and performed a structure-function analysis of MiD49 based on its homology to MiD51. MiD51 adopts a nucleotidyl transferase fold, and binds ADP as a co-factor that is essential for its function. Both MiDs contain a loop segment that is not present in other nucleotidyl transferase proteins, and this loop is used to interact with Drp1 and to recruit it to mitochondria.

Cooperative properties of lipid bilayers : collective director fluctuations and the effects of hydrophobic mismatch in protein/lipid membrane systems

Cooperative director fluctuations in lipid bilayers have been postulated for many years. ^2H-NMR T_1^(-1), T_(1P)^(-1) , and T_2^(-1); measurements have been used identify these motions and to determine the origin of increased slow bilayer motion upon addition of unlike lipids or proteins to a pure lipid bilayer.

The contribution of cooperative director fluctuations to NMR relaxation in lipid bilayers has been expressed mathematically using the approach of Doane et al.^1 and Pace and Chan.^2 The T_2^(-1)’s of pure dimyristoyllecithin (DML) bilayers deuterated at the 2, 9 and 10, and all positions on both lipid hydrocarbon chains have been measured. Several characteristics of these measurements indicate the presence of cooperative director fluctuations. First of all, T_2^(-1) exhibits a linear dependence on S2/CD. Secondly, T_2^(-1) varies across the ^2H-NMR powder pattern as sin^2 (2, β), where , β is the angle between the average bilayer director and the external magnetic field. Furthermore, these fluctuations are restricted near the lecithin head group suggesting that the head group does not participate in these motions but, rather, anchors the hydrocarbon chains in the bilayer.

T_2^(-1)has been measured for selectively deuterated liquid crystalline DML hilayers to which a host of other lipids and proteins have been added. The T_2^(-1) of the DML bilayer is found to increase drastically when chlorophyll a (chl a) and Gramicidin A' (GA') are added to the bilayer. Both these molecules interfere with the lecithin head group spacing in the bilayer. Molecules such as myristic acid, distearoyllecithin (DSL), phytol, and cholesterol, whose hydrocarbon regions are quite different from DML but which have small,neutral polar head groups, leave cooperative fluctuations in the DML bilayer unchanged.

The effect of chl a on cooperative fluctuations in the DML bilayer has been examined in detail using ^2H-NMR T_1^(-1), T_(1P)^(-1) , and T_2^(-1); measurements. Cooperative fluctuations have been modelled using the continuum theory of the nematic state of liquid crystals. Chl a is found to decrease both the correlation length and the elastic constants in the DML bilayer.

A mismatch between the hydrophobic length of a lipid bilayer and that of an added protein has also been found to change the cooperative properties of the lecithin bilayer. Hydrophobic mismatch has been studied in a series GA' / lecithin bilayers. The dependence of 2H-NMR order parameters and relaxation rates on GA' concentration has been measured in selectively deuterated DML, dipalmitoyllecithin (DPL), and DSL systems. Order parameters, cooperative lengths, and elastic constants of the DML bilayer are most disrupted by GA', while the DSL bilayer is the least perturbed by GA'. Thus, it is concluded that the hydrophobic length of GA' best matches that of the DSL bilayer. Preliminary Raman spectroscopy and Differential Scanning Calorimetry experiments of GA' /lecithin systems support this conclusion. Accommodation of hydrophobic mismatch is used to rationalize the absence of H_(II) phase formation in GA' /DML systems and the observation of H_(II) phase in GA' /DPL and GA' /DSL systems.

1. J. W. Doane and D. L. Johnson, Chem. Phy3. Lett., 6, 291-295 (1970). 2. R. J. Pace and S. I. Chan, J. Chem. Phy3., 16, 4217-4227 (1982).

Chemical studies of viral entry mechanisms: I. Hydrophobic protein-lipid interactions during sendai virus membrane fusion. II. Kinetics of bacteriophage λ DNA injection.

Viruses possess very specific methods of targeting and entering cells. These methods would be extremely useful if they could also be applied to drug delivery, but little is known about the molecular mechanisms of the viral entry process. In order to gain further insight into mechanisms of viral entry, chemical and spectroscopic studies in two systems were conducted, examining hydrophobic protein-lipid interactions during Sendai virus membrane fusion, and the kinetics of bacteriophage λ DNA injection.

Sendai virus glycoprotein interactions with target membranes during the early stages of fusion were examined using time-resolved hydrophobic photoaffinity labeling with the lipid-soluble carbene generator3-(trifluoromethyl)-3-(m-^(125 )I] iodophenyl)diazirine (TID). The probe was incorporated in target membranes prior to virus addition and photolysis. During Sendai virus fusion with liposomes composed of cardiolipin (CL) or phosphatidylserine (PS), the viral fusion (F) protein is preferentially labeled at early time points, supporting the hypothesis that hydrophobic interaction of the fusion peptide at the N-terminus of the F_1 subunit with the target membrane is an initiating event in fusion. Correlation of the hydrophobic interactions with independently monitored fusion kinetics further supports this conclusion. Separation of proteins after labeling shows that the F_1 subunit, containing the putative hydrophobic fusion sequence, is exclusively labeled, and that the F_2 subunit does not participate in fusion. Labeling shows temperature and pH dependence consistent with a need for protein conformational mobility and fusion at neutral pH. Higher amounts of labeling during fusion with CL vesicles than during virus-PS vesicle fusion reflects membrane packing regulation of peptide insertion into target membranes. Labeling of the viral hemagglutinin/neuraminidase (HN) at low pH indicates that HN-mediated fusion is triggered by hydrophobic interactions, after titration of acidic amino acids. HN labeling under nonfusogenic conditions reveals that viral binding may involve hydrophobic as well as electrostatic interactions. Controls for diffusional labeling exclude a major contribution from this source. Labeling during reconstituted Sendai virus envelope-liposome fusion shows that functional reconstitution involves protein retention of the ability to undergo hydrophobic interactions.

Examination of Sendai virus fusion with erythrocyte membranes indicates that hydrophobic interactions also trigger fusion between biological membranes, and that HN binding may involve hydrophobic interactions as well. Labeling of the erythrocyte membranes revealed close membrane association of spectrin, which may play a role in regulating membrane fusion. The data show that hydrophobic fusion protein interaction with both artificial and biological membranes is a triggering event in fusion. Correlation of these results with earlier studies of membrane hydration and fusion kinetics provides a more detailed view of the mechanism of fusion.

The kinetics of DNA injection by bacteriophage λ. into liposomes bearing reconstituted receptors were measured using fluorescence spectroscopy. LamB, the bacteriophage receptor, was extracted from bacteria and reconstituted into liposomes by detergent removal dialysis. The DNA binding fluorophore ethidium bromide was encapsulated in the liposomes during dialysis. Enhanced fluorescence of ethidium bromide upon binding to injected DNA was monitored, and showed that injection is a rapid, one-step process. The bimolecular rate law, determined by the method of initial rates, revealed that injection occurs several times faster than indicated by earlier studies employing indirect assays.

It is hoped that these studies will increase the understanding of the mechanisms of virus entry into cells, and to facilitate the development of virus-mimetic drug delivery strategies.

Indentification and characterization of a novel CA^(2+)-binding protein in avian erythrocytes

A novel Ca^(2+)-binding protein with Mr of 23 K (designated p23) has been identified in avian erythrocytes and thrombocytes. p23 localizes to the marginal bands (MBs), centrosomes and discrete sites around the nuclear membrane in mature avian erythrocytes. p23 appears to bind Ca^(2+) directly and its interaction with subcellular organelles seems to be modulated by intracellular [Ca^(2+)]. However, its unique protein sequence lacks any known Ca^(2+)-binding motif. Developmental analysis reveals that p23 association to its target structures occurs only at very late stages of bone marrow definitive erythropoeisis. In primitive erythroid cells, p23 distributes diffusely in the cytoplasm and lacks any distinct localization. It is postulated that p23 association to subcellular structures may be induced in part by decreased intracellular [Ca^(2+)]. In vitro and in vivo experiments indicate that p23 does not appear to act as a classical microtubule-associated protein (MAP) but p23 homologues appear to be expressed in MB-containing cells of a variety of species from different vertebrate classes. It has been hypothesized that p23 may play a regulatory role in MB stabilization in a Ca^(2+)-dependent manner.

Binucleated (bnbn) turkey erythrocytes were found to express a truncated p23 variant (designated p21) with identical subcellular localization as p23 except immunostaining reveals the presence of multi-centrosomes in bnbn cells. The p21 sequence has a 62 amino acid deletion at the C-terminus and must therefore have an additional ~40 amino acids at the N-terminus. In addition, p21 seems to have lost the ability to bind Ca^(2+) and its supramolecular interactions are not modulated by intracellular [Ca^(2+)]. These apparent differences between p23 and p21 raised the possibility that the p23/p21 allelism could be the Bn/bn genotype. However, genetic analysis suggested that p23/p21 allelism had no absolute correlation with the Bn/bn genotype.

Structural and mechanistic motifs in membrane proteins: the three-dimensional modelling of rhodopsin, band 3, and the nicotinic acetylcholine receptor

Because so little is known about the structure of membrane proteins, an attempt has been made in this work to develop techniques by which to model them in three dimensions. The procedures devised rely heavily upon the availability of several sequences of a given protein. The modelling procedure is composed of two parts. The first identifies transmembrane regions within the protein sequence on the basis of hydrophobicity, β-turn potential, and the presence of certain amino acid types, specifically, proline and basic residues. The second part of the procedure arranges these transmembrane helices within the bilayer based upon the evolutionary conservation of their residues. Conserved residues are oriented toward other helices and variable residues are positioned to face the surrounding lipids. Available structural information concerning the protein's helical arrangement, including the lengths of interhelical loops, is also taken into account. Rhodopsin, band 3, and the nicotinic acetylcholine receptor have all been modelled using this methodology, and mechanisms of action could be proposed based upon the resulting structures.

Specific residues in the rhodopsin and iodopsin sequences were identified, which may regulate the proteins' wavelength selectivities. A hinge-like motion of helices M3, M4, and M5 with respect to the rest of the protein was proposed to result in the activation of transducin, the G-protein associated with rhodopsin. A similar mechanism is also proposed for signal transduction by the muscarinic acetylcholine and β-adrenergic receptors.

The nicotinic acetylcholine receptor was modelled with four trans-membrane helices per subunit and with the five homologous M2 helices forming the cation channel. Putative channel-lining residues were identified and a mechanism of channel-opening based upon the concerted, tangential rotation of the M2 helices was proposed.

Band 3, the anion exchange protein found in the erythrocyte membrane, was modelled with 14 transmembrane helices. In general the pathway of anion transport can be viewed as a channel composed of six helices that contains a single hydrophobic restriction. This hydrophobic region will not allow the passage of charged species, unless they are part of an ion-pair. An arginine residue located near this restriction is proposed to be responsible for anion transport. When ion-paired with a transportable anion it rotates across the barrier and releases the anion on the other side of the membrane. A similar process returns it to its original position. This proposed mechanism, based on the three-dimensional model, can account for the passive, electroneutral, anion exchange observed for band 3. Dianions can be transported through a similar mechanism with the additional participation of a histidine residue. Both residues are located on M10.

Chemical-scale studies of G protein-coupled receptors and ligand-gated ion channels

This dissertation describes studies of G protein-coupled receptors (GPCRs) and ligand-gated ion channels (LGICs) using unnatural amino acid mutagenesis to gain high precision insights into the function of these important membrane proteins.

Chapter 2 considers the functional role of highly conserved proline residues within the transmembrane helices of the D2 dopamine GPCR. Through mutagenesis employing unnatural α-hydroxy acids, proline analogs, and N-methyl amino acids, we find that lack of backbone hydrogen bond donor ability is important to proline function. At one proline site we additionally find that a substituent on the proline backbone N is important to receptor function.

In Chapter 3, side chain conformation is probed by mutagenesis of GPCRs and the muscle-type nAChR. Specific side chain rearrangements of highly conserved residues have been proposed to accompany activation of these receptors. These rearrangements were probed using conformationally-biased β-substituted analogs of Trp and Phe and unnatural stereoisomers of Thr and Ile. We also modeled the conformational bias of the unnatural Trp and Phe analogs employed.

Chapters 4 and 5 examine details of ligand binding to nAChRs. Chapter 4 describes a study investigating the importance of hydrogen bonds between ligands and the complementary face of muscle-type and α4β4 nAChRs. A hydrogen bond involving the agonist appears to be important for ligand binding in the muscle-type receptor but not the α4β4 receptor.

Chapter 5 describes a study characterizing the binding of varenicline, an actively prescribed smoking cessation therapeutic, to the α7 nAChR. Additionally, binding interactions to the complementary face of the α7 binding site were examined for a small panel of agonists. We identified side chains important for binding large agonists such as varenicline, but dispensable for binding the small agonist ACh.

Chapter 6 describes efforts to image nAChRs site-specifically modified with a fluorophore by unnatural amino acid mutagenesis. While progress was hampered by high levels of fluorescent background, improvements to sample preparation and alternative strategies for fluorophore incorporation are described.

Chapter 7 describes efforts toward a fluorescence assay for G protein association with a GPCR, with the ultimate goal of probing key protein-protein interactions along the G protein/receptor interface. A wide range of fluorescent protein fusions were generated, expressed in Xenopus oocytes, and evaluated for their ability to associate with each other.

Protein-protein recognition: The neonatal Fe receptor and immunoglobulin G

The neonatal Fe receptor (FeRn) binds the Fe portion of immunoglobulin G (IgG) at the acidic pH of endosomes or the gut and releases IgG at the alkaline pH of blood. FeRn is responsible for the maternofetal transfer of IgG and for rescuing endocytosed IgG from a default degradative pathway. We investigated how FeRn interacts with IgG by constructing a heterodimeric form of the Fe (hdFc) that contains one FeRn binding site. This molecule was used to characterize the interaction between one FeRn molecule and one Fe and to determine under what conditions FeRn forms a dimer. The hdFc binds one FeRn molecule at pH 6.0 with a K_d of 80 nM. In solution and with FeRn anchored to solid supports, the heterodimeric Fe does not induce a dimer of FeRn molecules. FcRnhdFc complex crystals were obtained and the complex structure was solved to 2.8 Å resolution. Analysis of this structure refined the understanding of the mechanism of the pH-dependent binding, shed light on the role played by carbohydrates in the Fe binding, and provided insights on how to design therapeutic IgG antibodies with longer serum half-lives. The FcRn-hdFc complex in the crystal did not contain the FeRn dimer. To characterize the tendency of FeRn to form a dimer in a membrane we analyzed the tendency of the hdFc to induce cross-phosphorylation of FeRn-tyrosine kinase chimeras. We also constructed FeRn-cyan and FeRn-yellow fluorescent proteins and have analyzed the tendency of these molecules to exhibit fluorescence resonance energy transfer. As of now, neither of these analyses have lead to conclusive results. In the process of acquiring the context to appreciate the structure of the FcRn-hdFc interface, we developed a study of 171 other nonobligate protein-protein interfaces that includes an original principal component analysis of the quantifiable aspects of these interfaces.

Phospholipids and a protein-conducting channel regulate cotranslational protein targeting

Signal recognition particle (SRP) and signal recognition particle receptor (SR) are evolutionarily conserved GTPases that deliver secretory and membrane proteins to the protein-conducting channel Sec61 complex in the lipid bilayer of the endoplasmic reticulum in eukaryotes or the SecYEG complex in the inner membrane of bacteria. Unlike the canonical Ras-type GTPases, SRP and SR are activated via nucleotide-dependent heterodimerization. Upon formation of the SR•SRP targeting complex, SRP and SR undergo a series of discrete conformational changes that culminate in their reciprocal activation and hydrolysis of GTP. How the SR•SRP GTPase cycle is regulated and coupled to the delivery of the cargo protein to the protein-conducting channel at the target membrane is not well-understood. Here we examine the role of the lipid bilayer and SecYEG in regulation of the SRP-mediated protein targeting pathway and show that they serve as important biological cues that spatially control the targeting reaction.

In the first chapter, we show that anionic phospholipids of the inner membrane activate the bacterial SR, FtsY, and favor the late conformational states of the targeting complex conducive to efficient unloading of the cargo. The results of our studies suggest that the lipid bilayer acts as a spatial cue that weakens the interaction of the cargo protein with SRP and primes the complex for unloading its cargo onto SecYEG.

In the second chapter, we focus on the effect of SecYEG on the conformational states and activity of the targeting complex. While phospholipids prime the complex for unloading its cargo, they are insufficient to trigger hydrolysis of GTP and the release of the cargo from the complex. SecYEG modulates the conformation of the targeting complex and triggers the GTP hydrolysis from the complex, thus driving the targeting reaction to completion. The results of this study suggest that SecYEG is not a passive recipient of the cargo protein; rather, it actively releases the cargo from the targeting complex. Together, anionic phospholipids and SecYEG serve distinct yet complementary roles. They spatially control the targeting reaction in a sequential manner, ensuring efficient delivery and unloading of the cargo protein.

In the third chapter, we reconstitute the transfer reaction in vitro and visualize it in real time. We show that the ribosome-nascent chain complex is transferred to SecYEG via a stepwise mechanism with gradual dissolution and formation of the contacts with SRP and SecYEG, respectively, explaining how the cargo is kept tethered to the membrane during the transfer and how its loss to the cytosol is avoided.

In the fourth chapter, we examine interaction of SecYEG with secretory and membrane proteins and attempt to address the role of a novel insertase YidC in this interaction. We show that detergent-solubilized SecYEG is capable of discriminating between the nascent chains of various lengths and engages a signal sequence in a well-defined conformation in the absence of accessory factors. Further, YidC alters the conformation of the signal peptide bound to SecYEG. The results described in this chapter show that YidC affects the SecYEG-nascent chain interaction at early stages of translocation/insertion and suggest a YidC-facilitated mechanism for lateral exit of transmembrane domains from SecYEG into the lipid bilayer.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Investigating the Role of the GET3-GET4/GET5 Interaction During Tail-Anchor Protein Targeting

The proper targeting of membrane proteins is essential to the viability of all cells. Tail-anchored (TA) proteins, defined as having a single transmembrane helix at their C-terminus, are post-translationally targeted to the endoplasmic reticulum (ER) membrane by the GET pathway (Guided Entry of TA proteins). In the yeast pathway, the handover of TA substrates is mediated by the heterotetrameric Get4/Get5 (Get4/5) complex, which tethers the co-chaperone Sgt2 to the central targeting factor, the Get3 ATPase. Although binding of Get4/5 to Get3 is critical for efficient TA targeting, the mechanisms by which Get4 regulates Get3 are unknown. To understand the molecular basis of Get4 function, we used a combination of structural biology, biochemistry, and cell biology. Get4/5 binds across the Get3 dimer interface, in an orientation only compatible with a closed Get3, providing insight into the role of nucleotide in complex formation. Additionally, this structure reveals two functionally distinct binding interfaces for anchoring and ATPase regulation, and loss of the regulatory interface leads to strong defects in vitro and in vivo. Additional crystal structures of the Get3-Get4/5 complex give rise to an alternate conformation, which represents an initial binding interaction mediated by electrostatics that facilitates the rate of subsequent inhibited complex formation. This interface is supported by an in-depth kinetic analysis of the Get3-Get4/5 interaction confirming the two-step complex formation. These results allow us to generate a refined model for Get4/5 function in TA targeting.