Disease-causing mutations in the cellular retinaldehyde binding protein tighten and abolish ligand interactions

Mutations in the human cellular retinaldehyde binding protein (CRALBP) gene cause retinal pathology. To understand the molecular basis of impaired CRALBP function, we have characterized human recombinant CRALBP containing the disease causing mutations R233W or M225K. Protein structures were verified by amino acid analysis and mass spectrometry, retinoid binding properties were evaluated by UV-visible and fluorescence spectroscopy and substrate carrier functions were assayed for recombinant 11-cis-retinol dehydrogenase (rRDH5). The M225K mutant was less soluble than the R233W mutant and lacked retinoid binding capability and substrate carrier function. In contrast, the R233W mutant exhibited solubility comparable to wild type rCRALBP and bound stoichiometric amounts of 11-cis- and 9-cis-retinal with at least 2-fold higher affinity than wild type rCRALBP. Holo-R233W significantly decreased the apparent affinity of rRDH5 for 11-cis-retinoid relative to wild type rCRALBP. Analyses by heteronuclear single quantum correlation NMR demonstrated that the R233W protein exhibits a different conformation than wild type rCRALBP, including a different retinoid-binding pocket conformation. The R233W mutant also undergoes less extensive structural changes upon photoisomerization of bound ligand, suggesting a more constrained structure than that of the wild type protein. Overall, the results show that the M225K mutation abolishes and the R233W mutation tightens retinoid binding and both impair CRALBP function in the visual cycle as an 11-cis-retinol acceptor and as a substrate carrier.     

Single-chain site-specific mutations of fluorescein-amino acid contact residues in high affinity monoclonal antibody 4-4-20.

Previous crystallographic studies of high affinity anti-fluorescein monoclonal antibody 4-4-20 (Ka = 1.7 x 10(10) M-1) complexed with fluorescyl ligand resolved active site contact residues involved in binding. For better definition of the relative roles of three light chain antigen contact residues (L27dhis, L32tyr and L34arg), four site-specific mutations (L27dhis to L27lys, L32tyr to L32phe, and L34arg to L34lys and L34his) were generated and expressed in single-chain antigen binding derivatives of monoclonal antibody 4-4-20 containing two different polypeptide linkers (SCA 4-4-20/205c, 25 amino acids and SCA 4-4-20/212, 14 amino acids). Results showed that L27dhis and L32tyr were necessary for wild type binding affinities, however, were not required for near-wild type Qmax values (where Qmax is the maximum fluoroscein fluorescence quenching expressed as percent). Tyrosine L32 which hydrogen bonds with ligand was also characterized at the haptenic level through the use of 9-hydroxyphenylfluoron which lacks the carboxyl group to which L32 tyrosine forms a hydrogen bond. Results demonstrated that wild type SCA and mutant L32phe possessed similar HPF binding characteristics. Active site contact residue L34arg was important for fluorescein quenching maxima and binding affinity (L34his mutant), however, substitution of lysine for arginine at L34 did not have a significant effect on observed Qmax value. In addition, substitutions had no effect on structural and topological characteristics, since all mutants retained similar idiotypic and metatypic properties. Finally, two linkers were comparatively examined to determine relative contributions to mutant binding properties and stability. No linker effects were observed. Collectively, these results verified the importance of these light chain fluorescein contact residues in the binding pocket of monoclonal antibody 4-4-20.     

Affinity Maturation of Cry1Aa Toxin to the Bombyx mori Cadherin-Like Receptor by Directed Evolution

Improvement of the activity and insecticidal spectrum of cloned Cry toxins of Bacillus thuringiensis should allow for their wider application as biopesticides and a gene source for gene-modified crops. The insecticidal activity of Cry toxins depends on their binding to the receptor. Therefore, as a model, we aimed to generate improved binding affinity mutant toxins against Bombyx mori cadherin-like receptor (BtR175) using methods of directed evolution with the expectation of insecticidal activity improved mutants. Four serial amino acid residues of ⁴³⁹QAAG⁴⁴² or ⁴⁴³AVYT⁴⁴⁶ of Cry1Aa were replaced with random amino acids and were displayed on the T7 phage for library construction. Through five cycles of panning of the phage libraries using BtR175, 11 mutant phage clones were concentrated, and mutant toxin sequences were confirmed. The binding affinities of the three mutants were 42-, 15-, and 13-fold higher than that of the wild type, indicating that mutants with improved binding affinity to cadherin can be easily selected from randomly replaced loop 3 mutant libraries using directed evolution. We discuss the development of a genetic engineering method based on directed evolution to improve the binding affinity of Cry toxin to receptors.     

Exploring the cause of aggregation and reduced Zn binding affinity by G85R mutation in SOD1 rendering amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS), a lethal neurodegenerative disorder is characterized by the degeneration of upper and lower motor neuron. ALS occurs due to various notably prominent missense mutations, in gene encoding Cu‐Zn superoxide dismutase (SOD1) thereby leading to aggregation, dysfunction and reduced Zn binding affinity. In this study, one such mutation, G85R was explored in comparison with wild type SOD1, using discrete molecular dynamics (DMD). Accordingly, the conformational changes were significantly observed in mutant SOD1, through various geometrical parameters, which substantiated the difference in conformational deviation, flexibility and compactness, thus stipulating a root cause for aggregation. Followed by, analysis of essential dynamics further authenticated the cause behind the protein dysfunction. In particular, the high content of beta sheet with structural deviations, down to dysfunction was established in mutant as compared to wild type, while passing through secondary structure analysis. Subsequently, the deviation of distance in Zn binding residues was distinctly portrayed in mutant as compared to wild type, thus confirming the cause of reduced Zn binding affinity. In addition, the steered molecular dynamics analysis also authenticated the above results indicating the reduced Zn binding affinity in the mutant as compared to that of the wild type. Hence, this work revealed the theoretical mechanism to unravel the mutational effects of cofactor dependent protein. Proteins 2017; 85:1276–1286. © 2017 Wiley Periodicals, Inc.     

Teaching TetR to recognize a new inducer

Tet Repressor (TetR) recognizes the inducer tetracycline (tc) with high affinity. The tc analog 4-de(dimethylamino)-6-deoxy-6-demethyl-tetracycline (cmt3) is not an inducer for TetR. Induction specificity for cmt3 was generated by employing a directed evolution approach to screen appropriate TetR mutants in four successive steps. The specificity of the best TetR mutant is more than 20,000-fold increased for cmt3 over tc as judged by the ratio of their respective binding constants. Two rounds of directed evolution via DNA shuffling revealed His64 as a key residue for inducer specificity. The best TetR mutant with cmt3 specificity contains the H64K exchange, leading to a 300-fold decreased tc and a 20-fold increased cmt3 affinity. Another round of directed evolution made use of randomized oligonucleotides to mutate selected residues close to the tc-binding pocket of TetR and yielded TetR S135L with a 250-fold increased cmt3 affinity. The double mutant TetR H64K S135L was constructed and again subjected to directed evolution using randomized oligonucleotides to alter residues in the "secondary shell" of the tc-binding pocket. The resulting best mutants TetR H64K E114Q S135L, TetR A61V H64K Q109E Q116E S135L and TetR H64K T112K S135L are fully inducible by cmt3 and not by tc. Thus, their inducer specificity has been redesigned. The molecular mechanism of changed inducer recognition is discussed, based on binding constants with several tc analogs and in light of the TetR crystal structure.     

Computational protein design suggests that human PCNA‐partner interactions are not optimized for affinity

Increasing the affinity of binding proteins is invaluable for basic and applied biological research. Currently, directed protein evolution experiments are the main approach for generating such proteins through the construction and screening of large mutant libraries. Proliferating cell nuclear antigen (PCNA) is an essential hub protein that interacts with many different partners to tightly regulate DNA replication and repair in all eukaryotes. Here, we used computational design to generate human PCNA mutants with enhanced affinity for several different partners. We identified double mutations in PCNA, outside the main partner binding site, that were predicted to increase PCNA‐partner binding affinities compared to the wild‐type protein by forming additional hydrophobic interactions with conserved residues in the PCNA partners. Affinity increases were experimentally validated with four different PCNA partners, demonstrating that computational design can reveal unexpected regions where affinity enhancements in natural systems are possible. The designed PCNA mutants can be used as a valuable tool for further examination of the regulation of PCNA‐partner interactions during DNA replication and repair both in vitro and in vivo. More broadly, the ability to engineer affinity increases toward several PCNA partners suggests that interaction affinity is not an evolutionarily optimized trait of this system. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Engineered pH-dependent recycling antibodies enhance elimination of Staphylococcal enterotoxin B superantigen in mice

A new modality in antibody engineering has emerged in which the antigen affinity is designed to be pH dependent (PHD). In particular, combining high affinity binding at neutral pH with low affinity binding at acidic pH leads to a novel antibody that can more effectively neutralize the target antigen while avoiding antibody-mediated antigen accumulation. Here, we studied how the in vivo pharmacokinetics of the superantigen, Staphylococcal enterotoxin B (SEB), is affected by an engineered antibody with pH-dependent binding. PHD anti-SEB antibodies were engineered by introducing mutations into a high affinity anti-SEB antibody, 3E2, by rational design and directed evolution. Three antibody mutants engineered in the study have an affinity at pH 6.0 that is up to 68-fold weaker than the control antibody. The pH dependency of each mutant, measured as the pH-dependent affinity ratio (PAR – ratio of affinity at pH 7.4 and pH 6.0), ranged from 6.7–11.5 compared to 1.5 for the control antibody. The antibodies were characterized in mice by measuring their effects on the pharmacodynamics and pharmacokinetics (PK) of SEB after co-administration. All antibodies were effective in neutralizing the toxin and reducing the toxin-induced cytokine production. However, engineered PHD antibodies led to significantly faster elimination of the toxin from the circulation than wild type 3E2. The area under the curve computed from the SEB PK profile correlated well with the PAR value of antibody, indicating the importance of fine tuning the pH dependency of binding. These results suggest that a PHD recycling antibody may be useful to treat intoxication from a bacterial toxin by accelerating its clearance.     

Improving the affinity of antigens for mutated antibodies by use of statistical molecular design

We demonstrate the use of statistical molecular design (SMD) in the selection of peptide libraries aimed to systematically investigate antigen-antibody binding spaces. Earlier, we derived two novel antibodies by mutating the complementarity-determining region of the anti-p24 (HIV-1) single chain Fv antibody, CB4-1 that had lost their affinity for a p24 epitope-homologous peptide by 8- and 60-fold. The present study was devoted to explore how peptide libraries can be designed under experimental design criteria for effective screening of peptide antigens. Several small peptide-antigen libraries were selected using SMD principles and their activities were evaluated by their binding to SPOT-synthesized peptide membranes and by fluorescence polarization (FP). The approach was able to reveal the most critical residues required for antigen binding, and finally to increase the binding activity by proper modifications of amino acids in the peptide antigen. A model of the active peptide binding pocket formed by the mutated scFv and the antigen was compatible with the information gained from the experimental data. Our results suggest that SMD approaches can be used to explore peptide antigen features essential for their interactions with antibodies.     

Role of a noncanonical disulfide bond in the stability,affinity, and flexibility of a VHH specific for the Listeria virulence factor InlB

A distinguishing feature of camel (Camelus dromedarius) VHH domains are noncanonical disulfide bonds between CDR1 and CDR3. The disulfide bond may provide an evolutionary advantage, as one of the cysteines in the bond is germline encoded. It has been hypothesized that this additional disulfide bond may play a role in binding affinity by reducing the entropic penalty associated with immobilization of a long CDR3 loop upon antigen binding. To examine the role of a noncanonical disulfide bond on antigen binding and the biophysical properties of a VHH domain, we have used the VHH R303, which binds the Listeria virulence factor InlB as a model. Using site directed mutagenesis, we produced a double mutant of R303 (C33A/C102A) to remove the extra disulfide bond of the VHH R303. Antigen binding was not affected by loss of the disulfide bond, however the mutant VHH displayed reduced thermal stability (T_m = 12°C lower than wild‐type), and a loss of the ability to fold reversibly due to heat induced aggregation. X‐ray structures of the mutant alone and in complex with InlB showed no major changes in the structure. B‐factor analysis of the structures suggested that the loss of the disulfide bond elicited no major change on the flexibility of the CDR loops, and revealed no evidence of loop immobilization upon antigen binding. These results suggest that the noncanonical disulfide bond found in camel VHH may have evolved to stabilize the biophysical properties of the domain, rather than playing a significant role in antigen binding.     

Cancer-related mutations in BRCA1-BRCT cause long-range structural changes in protein-protein binding sites: a molecular dynamics study

Cancer-associated mutations in the BRCT domain of BRCA1 (BRCA1-BRCT) abolish its tumor suppressor function by disrupting interactions with other proteins such as BACH1. Many cancer-related mutations do not cause sufficient destabilization to lead to global unfolding under physiological conditions, and thus abrogation of function probably is due to localized structural changes. To explore the reasons for mutation-induced loss of function, the authors performed molecular dynamics simulations on three cancer-associated mutants, A1708E, M1775R, and Y1853ter, and on the wild type and benign M1652I mutant, and compared the structures and fluctuations. Only the cancer-associated mutants exhibited significant backbone structure differences from the wild-type crystal structure in BACH1-binding regions, some of which are far from the mutation sites. Backbone differences of the A1708E mutant from the liganded wild type structure in these regions are much larger than those of the unliganded wild type X-ray or molecular dynamics structures. These BACH1-binding regions of the cancer-associated mutants also exhibited increases in their fluctuation magnitudes compared with the same regions in the wild type and M1562I mutant, as quantified by quasiharmonic analysis. Several of the regions of increased fluctuation magnitude correspond to correlated motions of residues in contact that provide a continuous path of fluctuating amino acids in contact from the A1708E and Y1853ter mutation sites to the BACH1-binding sites with altered structure and dynamics. The increased fluctuations in the disease-related mutants suggest an increase in vibrational entropy in the unliganded state that could result in a larger entropy loss in the disease-related mutants upon binding BACH1 than in the wild type. To investigate this possibility, vibrational entropies of the A1708E and wild type in the free state and bound to a BACH1-derived phosphopeptide were calculated using quasiharmonic analysis, to determine the binding entropy difference DeltaDeltaS between the A1708E mutant and the wild type. DeltaDeltaS was determined to be -4.0 cal mol(-1) K(-1), with an uncertainty of 2 cal mol(-1) K(-1); that is, the entropy loss upon binding the peptide is 4.0 cal mol(-1) K(-1) greater for the A1708E mutant, corresponding to an entropic contribution to the DeltaDeltaG of binding (-TDeltaDeltaS) 1.1 kcal mol(-1) more positive for the mutant. The observed differences in structure, flexibility, and entropy of binding likely are responsible for abolition of BACH1 binding, and illustrate that many disease- related mutations could have very long-range effects. The methods described here have potential for identifying correlated motions responsible for other long-range effects of deleterious mutations.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Tertiary core rearrangements in a tight binding transfer RNA aptamer

Guided by an in vitro selection experiment designed to obtain tight binding aptamers of Escherichia coli glutamine specific tRNA (tRNAGln) for glutaminyl-tRNA synthetase (GlnRS), we have engineered a tRNA mutant in which the five-nucleotide variable loop sequence 5'-44CAUUC48-3' is replaced by 5'-44AGGU48-3'. This mutant tRNA binds to GlnRS with 30-fold improved affinity compared to the wild type. The 2.7 A cocrystal structure of the RNA aptamer-GlnRS complex reveals major rearrangements in the central tertiary core of the tRNA, while maintaining an RNA-protein interface identical to the wild type. The repacked RNA core features a novel hydrogen bonding arrangement of the trans Levitt pair G15-U48, a new sulfate binding pocket in the major groove, and increased hydrophobic stacking interactions among the bases. These data suggest that enhanced protein binding to a mutant globular RNA can arise from stabilization of RNA tertiary interactions rather than optimization of RNA-protein contacts.     

UPregulated single-stranded DNA-binding protein 1 induces cell chemoresistance to cisplatin in lung cancer cell lines

In order to understand the molecular basis of cold adaptation, we have used directed evolution to transform a thermophilic lipase LipR1 into its psychrophilic counterpart. A single round of random mutagenesis followed by screening for improved variants yielded a mutant with single-point mutation LipR1M1 (S130T), with optimum activity at 20?°C. Its activity at 50?°C is only 20% as compared to wild type (100%). It showed catalytic rate constant (k _cat) 3 times higher and a catalytic efficiency (k _cat/K _m) 4 times that of wild type. Circular dichroism and fluorescence studies also supported our observation of mutant structural flexibility. Structure analysis using homology models showed that Threonine 130 is exposed to solvent and has lost H-bond interaction with neighboring amino acid, thereby increasing flexibility of this lipase structure.     

Mutation of Phe91 to Asn in human carbonic anhydrase I unexpectedly enhanced both catalytic activity and affinity for sulfonamide inhibitors

Site-directed mutagenesis has been used to change one amino acid residue considered non essential (Phe91Asn) to catalysis in carbonic anhydrase (CA, EC 4.2.1.1) isozyme I (hCA I), but which is near the substrate binding pocket of the enzyme. This change led to a steady increase of 16% of the catalytic activity of the mutant hCA I over the wild type enzyme, which is a gain of 50% catalytic efficiency if one compares hCA I and hCA II as catalysts for CO₂ hydration. This effect may be due to the bigger hydrophobic pocket in the mutant enzyme compared to the wild type one, which probably leads to the reorganization of the solvent molecules present in the cavity and to a diverse proton transfer pathway in the mutant over the non mutated enzyme. To our surprise, the mutant CA I was not only a better catalyst for the physiologic reaction, but in many cases also showed higher affinity (2.6–15.9 times) for sulfonamide/sulfamate inhibitors compared to the wild type enzyme. As the residue in position 91 is highly variable among the 13 catalytically active CA isoforms, this study may shed a better understanding of catalysis/inhibition by this superfamily of enzymes.     

Molecular mechanism of the enhanced virulence of 2009 pandemic Influenza A (H1N1) virus from D222G mutation in the hemagglutinin: a molecular modeling study

D222G mutation of the hemagglutinin (HA) is of special interest because of its close association with the enhanced virulence of 2009 pandemic influenza A (H1N1) virus through the increased binding affinity to α2,3-linked sialylated glycan receptors. However, there is still a lack of detailed understanding about the molecular mechanism of this enhanced virulence. Here, molecular dynamics simulation and binding free energy calculation were performed to explore the altered glycan receptor binding mechanism of HA upon the D222G mutation by studying the interaction of one α2,3-linked sialylglycan (sequence: SIA-GAL-NAG) with the wild type and D222G mutated HA. The binding free energy calculation based on the molecular mechanics generalized Born surface area (MM-GBSA) method indicates that the D222G mutated HA has a much stronger binding affinity with the studied α2,3-linked glycan than the wild type. This is consistent with the experimental result. The increased binding free energy of D222G mutant mainly comes from the increased energy contribution of Gln223. The structural analysis proves that the altered electrostatic potential of receptor binding domain (RBD) and the increased flexibility of 220-loop are the essential reasons leading to the increased affinity of HA to α2,3-linked sialic acid glycans. The obtained results of this study have allowed a deeper understanding of the receptor recognition mechanism and the pathogenicity of influenza virus, which will be valuable to the structure-based inhibitors design targeting influenza virus entry process.     

Compensatory Evolution of a WW Domain Variant Lacking the Strictly Conserved Trp Residue

Replacement of conserved amino acid residues during evolution of proteins can lead to divergence and the formation of new families with novel functions, but is often deleterious to both protein structure and function. Using the WW domain, we experimentally examined whether and to what degree second-site mutations can compensate for the reduction of function and loss of structure that accompany substitution of a strictly conserved amino acid residue. The W17F mutant of the WW domain, with substitution of the most strictly conserved Trp residue, is known to lack a specific three-dimensional structure and shows reduced binding affinity in comparison to the wild type. To obtain second-site revertants, we performed a selection experiment based on the proline-rich peptide (PY ligand) binding affinity using the W17F mutant as the initial sequence. After selection by ribosome display, we were able to select revertants that exhibited a maximum ninefold higher affinity to the PY ligand than the W17F mutant and showed an even better affinity than the wild type. In addition, we found that the functional restoration resulted in increased binding specificity in selected revertants, and the structures were more compact, with increased amounts of secondary structure, in comparison to the W17F mutant. Our results suggest that the defective structure and function of the proteins caused by mutations in highly conserved residues occurring through divergent evolution not only can be restored but can be further improved by compensatory mutations.     

Specificity and Affinity Quantification of Flexible Recognition from Underlying Energy Landscape Topography

Flexibility in biomolecular recognition is essential and critical for many cellular activities. Flexible recognition often leads to moderate affinity but high specificity, in contradiction with the conventional wisdom that high affinity and high specificity are coupled. Furthermore, quantitative understanding of the role of flexibility in biomolecular recognition is still challenging. Here, we meet the challenge by quantifying the intrinsic biomolecular recognition energy landscapes with and without flexibility through the underlying density of states. We quantified the thermodynamic intrinsic specificity by the topography of the intrinsic binding energy landscape and the kinetic specificity by association rate. We found that the thermodynamic and kinetic specificity are strongly correlated. Furthermore, we found that flexibility decreases binding affinity on one hand, but increases binding specificity on the other hand, and the decreasing or increasing proportion of affinity and specificity are strongly correlated with the degree of flexibility. This shows more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Our work provides a theoretical foundation and quantitative explanation of the previous qualitative studies on the relationship among flexibility, affinity and specificity. In addition, we found that the folding energy landscapes are more funneled with binding, indicating that binding helps folding during the recognition. Finally, we demonstrated that the whole binding-folding energy landscapes can be integrated by the rigid binding and isolated folding energy landscapes under weak flexibility. Our results provide a novel way to quantify the affinity and specificity in flexible biomolecular recognition.     

Exploring the role of the phospholipid ligand in endothelial protein C receptor: A molecular dynamics study

Endothelial protein C receptor (EPCR) is a CD1‐like transmembrane glycoprotein with important regulatory roles in protein C (PC) pathway, enhancing PC's anticoagulant, anti‐inflammatory, and antiapoptotic activities. Similarly to homologous CD1d, EPCR binds a phospholipid [phosphatidylethanolamine (PTY)] in a groove corresponding to the antigen‐presenting site, although it is not clear if lipid exchange can occur in EPCR as in CD1d. The presence of PTY seems essential for PC γ‐carboxyglutamic acid (Gla) domain binding. However, the lipid‐free form of the EPCR has not been characterized. We have investigated the structural role of PTY on EPCR, by multiple molecular dynamics (MD) simulations of ligand bound and unbound forms of the protein. Structural changes, subsequent to ligand removal, led to identification of two stable and folded ligand‐free conformations. Compared with the bound form, unbound structures showed a narrowing of the A′ pocket and a high flexibility of the helices around it, in agreement with CD1d simulation. Thus, a lipid exchange with a mechanism similar to CD1d is proposed. In addition, unbound conformations presented a reduced interaction surface for Gla domain, confirming the role of PTY in establishing the proper EPCR conformation for the interaction with its partner protein. Single MD simulations were also obtained for 29 mutant models with predicted structural stability and impaired binding ability. Ligand affinity calculations, based on linear interaction energy method, showed that substitution‐induced conformational changes affecting helices around the A′ pocket were associated to a reduced binding affinity. Mutants responsible for this effect may represent useful reagents for experimental tests. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

X-ray absorption and molecular dynamics study of cation binding sites in the purple membrane

The present work describes the results of a study aimed at identifying candidate cation binding sites on the extracellular region of bacteriorhodopsin, including a site near the retinal pocket. The approach used is a combined effort involving computational chemistry methods (computation of cation affinity maps and molecular dynamics) together with the Extended X-Ray Absorption Fine Structure (EXAFS) technique to obtain relevant information about the local structure of the protein in the neighborhood of Mn(2+) ions in different affinity binding sites. The results permit the identification of a high-affinity binding site where the ion is coordinated simultaneously to Asp212(-) and Asp85(-). Comparison of EXAFS data of the wild type protein with the quadruple mutant E9Q/E74Q/E194Q/E204Q at pH 7.0 and 10.0 demonstrate that extracellular glutamic acid residues are involved in cation binding.