PubMed Health⌕ Search

Biomedical subjects

Richard B Sessions

Publications and source records attributed to Richard B Sessions.

At least 19 recordsLinked to original sources

An efficient, path-independent method for free-energy calculations.

Classical free-energy methods depend on the definition of physical or nonphysical integration paths to calculate free-energy differences between states. This procedure can be problematic and computationally expensive when the states of interest do not overlap and are far apart in phase space. Here we introduce a novel method to calculate free-energy differences that is path-independent by transforming each end state into a reference state in which the vibrational entropy is the sole component of the total entropy, thus allowing direct computation of the relative free energy. We apply the method to calculate side-chain entropies of a beta-hairpin-forming peptide in a variety of backbone conformations, demonstrating its importance in determining structural propensities. We find that low-free-energy conformations achieve their stability through optimal trade off between enthalpic gains due to favorable interatomic interactions and entropic losses incurred by the same.

Chemistry, Physical↗

Mutational analysis of human CEACAM1: the potential of receptor polymorphism in increasing host susceptibility to bacterial infection.

A common overlapping site on the N-terminal IgV-like domain of human carcinoembryonic antigen (CEA)-related cell adhesion molecules (CEACAMs) is targeted by several important human respiratory pathogens. These include Neisseria meningitidis (Nm) and Haemophilus influenzae (Hi) that can cause disseminated or persistent localized infections. To define the precise structural features that determine the binding of distinct pathogens with CEACAMs, we have undertaken molecular modelling and mutation of the receptor molecules at previously implicated key target residues required for bacterial binding. These include Ser-32, Tyr-34, Val-39, Gln-44 and Gln-89, in addition to Ile-91, the primary docking site for the pathogens. Most, but not all, of these residues located adjacent to each other in a previous N-domain model of human CEACAM1, which was based on REI, CD2 and CD4. In the current studies, we have refined this model based on the mouse CEACAM1 crystal structure, and observe that all of the above residues form an exposed continuous binding region on the N-domain. Examination of the model also suggested that substitution of two of these residues 34 and 89 could affect the accessibility of Ile-91 for ligand binding. By introducing selected mutations at the positions 91, 34 and 89, we confirmed the primary importance of Ile-91 in all bacterial binding to CEACAM1 despite the inter- and intraspecies structural differences between the bacterial CEACAM-binding ligands. The studies further indicated that the efficiency of binding was significantly enhanced for specific strains by mutations such as Y34F and Q89N, which also altered the hierarchy of Nm versus Hi strain binding. These studies imply that distinct polymorphisms in human epithelial CEACAMs have the potential to decrease or increase the risk of infection by the receptor-targeting pathogens.

Amino Acid Sequence↗

Antibiotic recognition by binuclear metallo-beta-lactamases revealed by X-ray crystallography.

Metallo-beta-lactamases are zinc-dependent enzymes responsible for resistance to beta-lactam antibiotics in a variety of host bacteria, usually Gram-negative species that act as opportunist pathogens. They hydrolyze all classes of beta-lactam antibiotics, including carbapenems, and escape the action of available beta-lactamase inhibitors. Efforts to develop effective inhibitors have been hampered by the lack of structural information regarding how these enzymes recognize and turn over beta-lactam substrates. We report here the crystal structure of the Stenotrophomonas maltophilia L1 enzyme in complex with the hydrolysis product of the 7alpha-methoxyoxacephem, moxalactam. The on-enzyme complex is a 3'-exo-methylene species generated by elimination of the 1-methyltetrazolyl-5-thiolate anion from the 3'-methyl group. Moxalactam binding to L1 involves direct interaction of the two active site zinc ions with the beta-lactam amide and C4 carboxylate, groups that are common to all beta-lactam substrates. The 7beta-[(4-hydroxyphenyl)malonyl]-amino substituent makes limited hydrophobic and hydrogen bonding contacts with the active site groove. The mode of binding provides strong evidence that a water molecule situated between the two metal ions is the most likely nucleophile in the hydrolytic reaction. These data suggest a reaction mechanism for metallo-beta-lactamases in which both metal ions contribute to catalysis by activating the bridging water/hydroxide nucleophile, polarizing the substrate amide bond for attack and stabilizing anionic nitrogen intermediates. The structure illustrates how a binuclear zinc site confers upon metallo-beta-lactamases the ability both to recognize and efficiently hydrolyze a wide variety of beta-lactam substrates.

Crystallography, X-Ray↗

Low energy pathways and non-native interactions: the influence of artificial disulfide bridges on the mechanism of folding.

Four versions of a beta-sheet protein (CD2.d1) have been made, each with a single artificial disulfide bond inserted into hairpin structures. Folding kinetics of reduced and oxidized forms shows bridge position strongly influences its effect on the folding reaction. Bridging residues 58 and 62 does not affect the rapidly formed intermediate (I) or rate-limiting transition (t) state, whereas bridging 33 and 38, or 31 and 41, lowers the t-state energy, with the latter having the stronger influence. Bridging residues 79 and 90 stabilizes both I- and t-states. To assess additivity in the energetic effects of these bridges, four double-bridge variants have also been made. All show precise additivity of overall stability, with two showing additivity when ground states and the rate-limiting t-state are assessed, i.e. no measurable change in the folding mechanism occurs. However, combining 31-41 and 79-90 bridges produces a molecule that folds through a different pathway, with a much more stable intermediate than expected and a much higher t-state barrier. This is explained by the artificial introduction of stabilizing, non-native contacts in the I-state. More surprisingly, for another double-bridge version (58-62 and 79-90) both I- and t-states are less stable than expected, showing that conformational constraints introduced by the two bridges prevent formation of non-native contacts that would otherwise stabilize the I- and t-states, thereby lowering the energy of the folding landscape in the wild-type (unbridged) molecule. We conclude that the lowest energy path for folding has I- and t-state structures that are stabilized by non-native interactions.

Animals↗

Mapping the binding site for gossypol-like inhibitors of Plasmodium falciparum lactate dehydrogenase.

Gossypol is a di-sesquiterpene natural-product in the form of a functionalised binaphthyl and is isolated from cotton plants. The compound has long been known to exhibit anti-malarial and other biological activities. Previous studies have indicated that compounds of this type target Plasmodium falciparum lactate dehydrogenase (pfLDH), an essential enzyme for energy generation within the parasite. In this study, we report that simple naphthalene-based compounds, the core of the gossypol structure, exhibit weak inhibition of the parasite lactate dehydrogenase. Crystal structures of the complexes formed by binding of these naphthalene-based compounds to their target enzyme have been used to delineate the molecular features likely to form the gossypol binding site. Two modes of binding are observed: one overlapping the pyruvate but not the co-factor site, the other bridging the binding sites for the co-factor nicontinamide group and pyruvate substrate. This latter site encompasses molecular features unique to Plasmodium forms of LDH and is likely to represent the mode of binding for gossypol derivatives that show selectivity for the parasite enzymes. We also report a substrate analogue that unexpectedly binds within the adenine pocket of the co-factor groove. Although these core pharmacophore-like molecules only exhibit low levels of inhibitory activity, these molecular snapshots provide a rational basis for renewed structure-based development of naphthalene-based compounds as anti-malarial agents.

Animals↗

Density guided importance sampling: application to a reduced model of protein folding.

MOTIVATION: Monte Carlo methods are the most effective means of exploring the energy landscapes of protein folding. The rugged topography of folding energy landscapes causes sampling inefficiencies however, particularly at low, physiological temperatures. RESULTS: A hybrid Monte Carlo method, termed density guided importance sampling (DGIS), is presented that overcomes these sampling inefficiencies. The method is shown to be highly accurate and efficient in determining Boltzmann weighted structural metrics of a discrete off-lattice protein model. In comparison to the Metropolis Monte Carlo method, and the hybrid Monte Carlo methods, jump-walking, smart-walking and replica-exchange, the DGIS method is shown to be more efficient, requiring no parameter optimization. The method guides the simulation towards under-sampled regions of the energy spectrum and recognizes when equilibrium has been reached, avoiding arbitrary and excessively long simulation times. AVAILABILITY: Fortran code available from authors upon request. CONTACT: m.j.parker@leeds.ac.uk.

Algorithms↗

Protein-4.2 association with band 3 (AE1, SLCA4) in Xenopus oocytes: effects of three natural protein-4.2 mutations associated with hemolytic anemia.

We have investigated the effects of coexpression of protein 4.2 and three protein-4.2 variants with band 3 in the Xenopus oocyte expression system. Normal protein 4.2 increased band-3-specific chloride transport in the oocytes. Protein 4.2 also coimmunoprecipitated with band 3 and colocalized with band 3 at the oocyte plasma membrane. The increase in band-3-mediated chloride transport and coimmunoprecipitation of protein 4.2 required the presence of the N-terminal cytoplasmic domain of band 3. Protein 4.2 also localized to the oocyte plasma membrane in the absence of band 3. The protein-4.2 variants 4.2 Tozeur (R310Q) and 4.2 Komatsu (D175Y) had impaired ability to bind to band 3 and these variants did not localize to the oocyte plasma membrane when expressed on their own or when coexpressed with band 3. Unexpectedly, 4.2 Nippon (A142T) behaved similarly to normal protein 4.2. In the absence of a crystal structure of protein 4.2, we propose a homology model of protein 4.2 based on the structure of the sequence-related protein transglutaminase. Using our results in oocytes and this homology model we speculate how these mutations affect protein 4.2 and result in hereditary spherocytosis.

Anemia, Hemolytic↗

Design, expression and characterisation of lysine-rich forms of the barley seed protein CI-2.

Chymotrypsin inhibitor CI-2 is a small (84 residue) barley seed protein that has been used extensively to study protein folding. It also contains eight lysine residues, making it an attractive target for expression in transgenic plants to increase their lysine contents. We have designed three lysine-enriched forms of CI-2 and compared their structures and properties with that of the wild type protein. One mutant containing three additional lysine residues in the inhibitory loop shows high stability to denaturation and reduced inhibitory activity, indicating its suitability for use in genetic engineering.

Amino Acid Sequence↗

Water as a conformational editor in protein folding.

As molecules approach one another in aqueous solution, desolvation free energy barriers to association are encountered. Experiments suggest these (de)solvation effects contribute to the free energy barriers separating the folded and unfolded states of protein molecules. To explore their influence on the energy landscapes of protein folding reactions, we have incorporated desolvation barriers into a semi-realistic, off-lattice protein model that uses a simplified physico-chemical force-field determined solely by the sequence of amino acids. Monte Carlo sampling techniques were used to study the effects on the thermodynamics and kinetics of folding of a number of systems, diverse in structure and sequence. In each case, desolvation barriers increase the stability of the native conformation and the cooperativity of the major folding/unfolding transition. The folding times of these systems are reduced significantly upon inclusion of desolvation barriers, demonstrating that the particulate nature of the solvent engenders a more defined route to the native fold.

Kinetics↗

The low-potency, voltage-dependent HERG blocker propafenone--molecular determinants and drug trapping.

The molecular determinants of high-affinity human ether-a-go-go-related gene (HERG) potassium channel blockade by methanesulfonanilides include two aromatic residues (Phe656 and Tyr652) on the inner helices (S6) and residues on the pore helices that face into the inner cavity, but determinants for lower-affinity HERG blockers may be different. In this study, alanine-substituted HERG channel mutants of inner cavity residues were expressed in Xenopus laevis oocytes and were used to characterize the HERG channel binding site of the antiarrhythmic propafenone. Propafenone's blockade of HERG was strongly dependent on residue Phe656 but was insensitive or weakly sensitive to mutation of Tyr652, Thr623, Ser624, Val625, Gly648, or Val659 and did not require functional inactivation. Homology models of HERG based on KcsA and MthK crystal structures, representing the closed and open forms of the channel, respectively, suggest propafenone is trapped in the inner cavity and is unable to interact exclusively with Phe656 in the closed state (whereas exclusive interactions between propafenone and Phe656 are found in the open-channel model). These findings are supported by very slow recovery of wild-type HERG channels from block at -120 mV, but extremely rapid recovery of D540K channels that reopen at this potential. The experiments and modeling suggest that the open-state propafenone binding-site may be formed by the Phe656 residues alone. The binding site for propafenone (which may involve pi-stacking interactions with two or more Phe656 side-chains) is either perturbed or becomes less accessible because of closed-channel gating. This provides further evidence for the existence of gating-induced changes in the spatial location of Phe656 side chains.

Amino Acid Sequence↗

Identification and activity of a series of azole-based compounds with lactate dehydrogenase-directed anti-malarial activity.

Plasmodium falciparum, the causative agent of malaria, relies extensively on glycolysis coupled with homolactic fermentation during its blood-borne stages for energy production. Selective inhibitors of the parasite lactate dehydrogenase (LDH), central to NAD(+) regeneration, therefore potentially provide a route to new antimalarial drugs directed against a novel molecular target. A series of heterocyclic, azole-based compounds are described that preferentially inhibit P. falciparum LDH at sub-micromolar concentrations, typically at concentrations about 100-fold lower than required for human lactate dehydrogenase inhibition. Crystal structures show these competitive inhibitors form a network of interactions with amino acids within the active site of the enzyme, stacking alongside the nicotinamide ring of the NAD(+) cofactor. These compounds display modest activity against parasitized erythrocytes, including parasite strains with known resistance to existing anti-malarials and against Plasmodium berghei in BALB/c mice. Initial toxicity data suggest the azole derivatives have generally low cytotoxicity, and preliminary pharmoco-kinetic data show favorable bioavailability and circulation times. These encouraging results suggest that further enhancement of these structures may yield candidates suitable for consideration as new therapeutics for the treatment of malaria. In combination these studies also provide strong support for the validity of targeting the Plasmodium glycolytic pathway and, in particular, LDH in the search for novel anti-malarials.

Animals↗

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding.

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described. More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.

Amides↗

Cloning, sequence and expression of the lactate dehydrogenase gene from the human malaria parasite, Plasmodium vivax.

Increased drug resistance to anti-malarials highlights the need for the development of new therapeutics for the treatment of malaria. To this end, the lactate dehydrogenase (LDH) gene was cloned and sequenced from genomic DNA of Plasmodium vivax ( PvLDH) Belem strain. The 316 amino acid protein-coding region of the PvLDH gene was inserted into the prokaryotic expression vector pKK223-3 and a 34 kDa protein with LDH activity was expressed in E. coli. Structural differences between human LDHs and PfLDH make the latter an attractive target for inhibitors leading to novel anti-malarial drugs. The sequence similarity between PvLDH and PfLDH (90% residue identity and no insertions or deletions) indicate that the same approach could be applied to Plasmodium vivax, the most common human malaria parasite in the world.

Amino Acid Sequence↗

Estimating the energetic contribution of hydrogen bonding to the stability of Candida methylica formate dehydrogenase by using double mutant cycle.

An homology model of Candida methylica formate dehydrogenase (cm FDH) was constructed based on the Pseudomonas sp. 101 formate dehydrogenase (ps FDH) structure. In wild type cm FDH, Thr169 and Thr226 can form hydrogen bonds with each other. We measured the interaction energy between the two threonines independent of other interactions in the proteins by using a so-called double mutant cycle and assessing the protein stability from the concentration of guanidine hydrochloride needed to denature 50% of the molecules. We conclude that the hydrogen bonds stabilize the wild type protein by -4 kcal mol(-1).

Biotechnology↗

Tyrphostin A23 inhibits internalization of the transferrin receptor by perturbing the interaction between tyrosine motifs and the medium chain subunit of the AP-2 adaptor complex.

Several intracellular membrane trafficking events are mediated by tyrosine-containing motifs within the cytosolic domains of integral membrane proteins. Many such motifs conform to the consensus YXXPhi, where Phi represents a bulky hydrophobic residue. This motif interacts with the medium chain (mu) subunits of adaptor complexes that link the cytosolic domains of integral membrane proteins to the clathrin coat involved in vesicle formation. The YXXPhi motif is similar to motifs in which the tyrosine residue is phosphorylated by tyrosine kinases. Tyrphostins (structural analogs of tyrosine) are inhibitors of tyrosine kinases and function by binding to the active sites of the enzymes. We previously showed that, in vitro and in yeast two-hybrid interaction assays, some tyrphostins can inhibit the interaction between YXXPhi motifs and the mu2 subunit of the AP-2 adaptor complex (Crump, C., Williams, J. L., Stephens, D. J., and Banting, G. (1998) J. Biol. Chem. 273, 28073-28077). A23 is such a tyrphostin. We now show that molecular modeling of tyrphostin A23 into the tyrosine-binding pocket in mu2 provides a structural explanation for A23 being able to inhibit the interaction between YXXPhi motifs and mu2. Furthermore, we show that A23 inhibited the internalization of (125)I-transferrin in Heb7a cells without having any discernible effect on the morphology of compartments of the endocytic pathway. Control tyrphostins, active as inhibitors of tyrosine kinase activity, but incapable of inhibiting the YXXPhi motif/mu2 interaction, did not inhibit endocytosis. These data are consistent with A23 inhibition of the YXXPhi motif/mu2 interaction in intact cells and with the possibility that different tyrphostins may be used to inhibit specific membrane trafficking events in eukaryotic cells.

Adaptor Protein Complex 2↗

The influence of intramolecular bridges on the dynamics of a protein folding reaction.

Thirteen versions of a beta-sheet protein have been constructed, each with a single, surface-exposed disulfide bridge. A comparison of folding kinetics, in oxidizing and reducing conditions, is used to elucidate the order in which beta-strands become associated during the folding process and, hence, the relationship between topology and folding dynamics. In common with the wild-type molecule, all the proteins fold through a two-step (three state) mechanism with a rapidly formed intermediate which slowly converts to the native state. In a majority of cases, the bridge is seen to stabilize the folded state, and for five of the modified proteins, the additional stability is greater than 3 kcal/mol. Surprisingly, cross-links which connect beta-strands which are distant in sequence predominantly stabilize the rapidly formed intermediate state, suggesting that these strand-strand interactions occur in the initial stages of folding. Cross-links which stabilize local hairpins have their major influence on the second, rate-determining step leading to significant enhancements in the folding rate. We find that enhancement of the folding rate in the second, rate-limiting step is correlated with a reduction in contact order in the same way as in naturally occurring proteins of different folds. The large increases in native-state stability resulting from the insertion of disulfide bridges on the surface of beta-sheet structures have implications for enhancing the robustness of proteins by molecular engineering.

Animals↗

A discrete domain of the human TrkB receptor defines the binding sites for BDNF and NT-4.

TrkB is a member of the Trk family of tyrosine kinase receptors. In vivo, the extracellular region of TrkB is known to bind, with high affinity, the neurotrophin protein brain-derived neurotrophic factor (BDNF) and neurotrophin-4 (NT-4). We describe the expression and purification of the second Ig-like domain of human TrkB (TrkBIg(2)) and show, using surface plasmon resonance, that this domain is sufficient to bind BDNF and NT-4 with subnanomolar affinity. BDNF and NT-4 may have therapeutic implications for a variety of neurodegenerative diseases. The specificity of binding of the neurotrophins to their receptor TrkB is therefore of interest. We examine the specificity of TrkBIg(2) for all the neurotrophins, and use our molecular model of the BDNF-TrkBIg(2) complex to examine the residues involved in binding. It is hoped that the understanding of specific interactions will allow design of small molecule neurotrophin mimetics.

Amino Acid Sequence↗

Characterization and modelling of the hydrophobic domain of a sunflower oleosin.

The oleosins are a group of hydrophobic proteins present on the surface of oil bodies in seeds, where they are thought to prevent coalescence. They contain a central hydrophobic domain of 68-74 residues that is thought to form a loop into the triacylglycerol matrix of the oil body, but the conformation adopted by this sequence is uncertain. We have therefore expressed an oleosin cDNA from sunflower (Helianthus annuus L.) in Escherichia coli as a fusion with maltose-binding protein (MBP) and isolated a peptide corresponding to the hydrophobic domain by sequential digestion with factor Xa (to remove the MBP) followed by trypsin and Staphylococcus V8 protease to remove the N- and C-terminal domains of the oleosin. Circular dichroism spectroscopy of the peptide in two solvent systems chosen to mimic the environment within the oil body (trifluoroethanol and SDS) demonstrated high proportions of alpha-helical structure, with no beta-sheet. A model was therefore developed in which the domain forms an alpha-helical hairpin structure, the two helices being separated by a turn region. We consider that this model is consistent with our current knowledge of oleosin structure and properties.

Amino Acid Sequence↗