Role of the minor energetic determinants of chicken egg white lysozyme (HEWL) to the stability of the HEWL.antibody scFv-10 complex

Seven of the 13 non-glycine contact amino acids in the hen (chicken) egg white lysozyme (HEWL) epitope for antibody Fab-10 each contribute < or =0.3 kcal/mol to the change in free energy (DeltaDeltaG(D)) from wild type (WT) when replaced by alanine (nullspots), and three others each give (0.7 < DeltaDeltaG(D) < or = 1. 0) kcal/mol (warm spots) (Rajpal et al. Protein Sci 1998;7:1868-1874). The low DeltaDeltaG(D) values introduced by alanine mutations present an opportunity to explore accurately their cumulative effects, as the sum of the combined DeltaDeltaG(D) values is not so large as to destabilize the complex beyond the range of accurate measurement. Substitution of six of the seven null spot residues by alanine leads to a cumulative DeltaDeltaG(D) = 2.25 +/- 0.04 kcal/mol, whereas the sum of the six individual changes is only -0.36 +/- 0.32 kcal/mol. The triple warm spot mutation generates a DeltaDeltaG(D) = 5.11 +/- 0.06 kcal/mol versus DeltaDeltaG(D) = 2.52 +/- 0.22 kcal/mol for the sum of the three individuals. The non-additivity in the individual DeltaDeltaG(D) values for the alanine mutations may indicate that these residues provide a conformationally stabilizing effect on the hot spot residues, each of which exhibits DeltaDeltaG(D) > 4.0 kcal/mol on alanine substitution.     

The association of Duffy binding protein region II polymorphisms and its antigenicity in Plasmodium vivax isolates from Thailand

Plasmodium vivax Duffy binding protein II (DBPII) plays an important role in reticulocyte invasion and is a potential vaccine candidate against vivax malaria. However, polymorphisms in DBPII are a challenge for the successful design of a broadly protective vaccine. In this study, the genetic diversity of DBPII among Thai isolates was analyzed from Plasmodium vivax-infected blood samples and polymorphism characters were defined with the MEGA4 program. Sequence analysis identified 12 variant residues that are common among Thai DBPII haplotypes with variant residues L333F, L424I, W437R and I503K having the highest frequency. Variant residue D384K occurs in combination with either E385K or K386N/Q. Additionally, variant residue L424I occurs in conjunction with W437R in most Thai DBPII alleles and these variants frequently occur in combination with the I503K variant. The polymorphic patterns of Thai isolates were defined into 9 haplotypes (Thai DBL-1, -2, -3, etc.…). Thai DBL-2, -5, -6 haplotypes are the most common DBPII variants in Thai residents. To study the association of these Thai DBPII polymorphisms with antigenic character, the functional inhibition of anti-DBPII monoclonal antibodies against a panel of Thai DBL variants was characterized by an in vitro erythrocyte binding inhibition assay. The functional inhibition of anti-DBPII monoclonal antibodies 3C9, 2D10 and 2C6 against Thai variants was significantly different, suggesting that polymorphisms of Thai DBPII variants alter the antigenic character of the target epitopes. In contrast, anti-DBPII monoclonal antibody 2H2 inhibited all Thai DBPII variants equally well. Our results suggest that the immune efficacy of a DBPII vaccine will depend on the specificity of the anti-DBPII antibodies induced and that it is preferable to optimize responses to conserved epitopes for broadly neutralizing protection against P. vivax.     

Free-energy simulations reveal that both hydrophobic and polar interactions are important for influenza hemagglutinin antibody binding

Antibodies binding to conserved epitopes can provide a broad range of neutralization to existing influenza subtypes and may also prevent the propagation of potential pandemic viruses by fighting against emerging strands. Here we propose a computational framework to study structural binding patterns and detailed molecular mechanisms of viral surface glycoprotein hemagglutinin (HA) binding with a broad spectrum of neutralizing monoclonal antibody fragments (Fab). We used rigorous free-energy perturbation (FEP) methods to calculate the antigen-antibody binding affinities, with an aggregate underlying molecular-dynamics simulation time of several microseconds (～2 μs) using all-atom, explicit-solvent models. We achieved a high accuracy in the validation of our FEP protocol against a series of known binding affinities for this complex system, with <0.5 kcal/mol errors on average. We then introduced what to our knowledge are novel mutations into the interfacial region to further study the binding mechanism. We found that the stacking interaction between Trp-21 in HA2 and Phe-55 in the CDR-H2 of Fab is crucial to the antibody-antigen association. A single mutation of either W21A or F55A can cause a binding affinity decrease of ΔΔG > 4.0 kcal/mol (equivalent to an ～1000-fold increase in the dissociation constant K_d). Moreover, for group 1 HA subtypes (which include both the H1N1 swine flu and the H5N1 bird flu), the relative binding affinities change only slightly (< ±1 kcal/mol) when nonpolar residues at the αA helix of HA mutate to conservative amino acids of similar size, which explains the broad neutralization capability of antibodies such as F10 and CR6261. Finally, we found that the hydrogen-bonding network between His-38 (in HA1) and Ser-30/Gln-64 (in Fab) is important for preserving the strong binding of Fab against group 1 HAs, whereas the lack of such hydrogen bonds with Asn-38 in most group 2 HAs may be responsible for the escape of antibody neutralization. These large-scale simulations may provide new insight into the antigen-antibody binding mechanism at the atomic level, which could be essential for designing more-effective vaccines for influenza.     

The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins

We recently found that Asp-His-Ser/Thr-Trp hydrogen-bonded tetrads are widely and uniquely present in the WD40-repeat proteins. WDR5 protein is a seven WD40-repeat propeller with five such tetrads. To explore the effect of the tetrad on the structure and stability of WD40-repeat proteins, the wild-type WDR5 and its seven mutants involving the substitutions of tetrad residues have been isolated. The crystal structures of the wild-type WDR5 and its three WDR5 mutants have been determined by X-ray diffraction method. The mutations of the tetrad residues are found not to change the basic structural features. The denaturing profiles of the wild type and the seven mutants with the use of denaturant guanidine hydrochloride have been studied by circular dichroism spectroscopy to determine the folding free energies of these proteins. The folding free energies of the wild type and the S62A, S146A, S188A, D192E, W330F, W330Y, and D324E mutants are measured to be about -11.6, -2.7, -3.1, -2.9, -3.6, -7.1, -7.0, and -7.5 kcal/mol, respectively. These suggest that (1) the hydrogen bonds in these hydrogen bond networks are unusually strong; (2) each hydrogen-bonded tetrad provides over 12 kcal/mol stability to the protein; thus, the removal of any single tetrad would cause unfolding of the protein; (3) since there are five tetrads, the protein must be in a highly unstable state without the tetrads, which might be related to its biological functions.     

Molecular Docking and Interaction Studies of Identified Abscisic Acid Receptors in Oryza sativa: An In-Silico Perspective on Comprehending Stress Tolerance Mechanisms

BackgroundAbiotic stresses affect plants in several ways and as such, phytohormones such as abscisic acid (ABA) play an important role in conferring tolerance towards these stresses. Hence, to comprehend the role of ABA and its interaction with receptors of the plants, a thorough investigation is essential.AimThe current study aimed to identify the ABA receptors in Oryza sativa, to find the receptor that binds best with ABA and to examine the mutations present to help predict better binding of the receptors with ABA.MethodsProtein sequences of twelve PYL (Pyrabactin resistance 1) and seven PP2C (type 2C protein phosphatase) receptors were retrieved from the Rice Annotation Project database and their 3D structures were predicted using RaptorX. Protein-ligand molecular docking studies between PYL and ABA were performed using AutoDock 1.5.6, followed by 100ns molecular dynamic simulation studies using Desmond to determine the acceptable conformational changes after docking via root mean square deviation RMSD plot analysis. Protein-protein docking was then carried out in three sets: PYL-PP2Cs, PYL-ABA-PP2C and PYL(mut)-ABA-PP2C to scrutinize changes in structural conformations and binding energies between complexes. The amino acids of interest were mapped at their respective genomic coordinates using SNP-seek database to ascertain if there were any naturally occurring single nucleotide polymorphisms (SNPs) responsible for triggering rice PYLs mutations.ResultsInitial protein-ligand docking studies revealed good binding between the complexes, wherein PYL6-ABA complex showed the best energy of -8.15 kcal/mol. The 100ns simulation studies revealed changes in the RMSD values after docking, indicating acceptable conformational changes. Furthermore, mutagenesis study performed at specific PYL-ABA interacting residues followed by downstream PYL(mut)-ABA-PP2C protein-protein docking results after induction of mutations demonstrated binding energy of -8.17 kcal/mol for PP2C79-PYL11-ABA complex. No naturally occurring SNPs that were responsible for triggering rice PYL mutations were identified when specific amino acid coordinates were mapped at respective genomic coordinates.ConclusionThus, the present study provides valuable insights on the interactions of ABA receptors in rice and induced mutations in PYL11 that can enhance the downstream interaction with PP2C.     

Modulating catalytic activity by unnatural amino acid residues in a GSH-binding loop of GST P1-1

The loop following helix α2 in glutathione transferase P1-1 has two conserved residues, Cys48 and Tyr50, important for glutathione (GSH) binding and catalytic activity. Chemical modification of Cys48 thwarts the catalytic activity of the enzyme, and mutation of Tyr50 generally decreases the k_cat value and the affinity for GSH in a differential manner. Cys48 and Tyr50 were targeted by site-specific mutations and chemical modifications in order to investigate how the α2 loop modulates GSH binding and catalysis. Mutation of Cys48 into Ala increased K_M^GSH 24-fold and decreased the binding energy of GSH by 1.5 kcal/mol. Furthermore, the protein stability against thermal inactivation and chemical denaturation decreased. The crystal structure of the Cys-free variant was determined, and its similarity to the wild-type structure suggests that the mutation of Cys48 increases the flexibility of the α2 loop rather than dislocating the GSH-interacting residues. On the other hand, replacement of Tyr50 with Cys, producing mutant Y50C, increased the Gibbs free energy of the catalyzed reaction by 4.8 kcal/mol, lowered the affinity for S-hexyl glutathione by 2.2 kcal/mol, and decreased the thermal stability. The targeted alkylation of Cys50 in Y50C increased the affinity for GSH and protein stability. Characterization of the most active alkylated variants, S-n-butyl-, S-n-pentyl-, and S-cyclobutylmethyl-Y50C, indicated that the affinity for GSH is restored by stabilizing the α2 loop through positioning of the key residue into the lock structure of the neighboring subunit. In addition, k_cat can be further modulated by varying the structure of the key residue side chain, which impinges on the rate-limiting step of catalysis.     

Molecular docking analysis of estrogen receptor binding phytocomponents identified from the ethyl acetate extract of Salicornia herbacea (L)

It is of interest to evaluate the secondary metabolites using high performance thin layer chromatography (HPTLC) finger printing and Gas chromatography-Mass spectroscopy (GC-MS) in S. herbaceaextract. The powdered plant material extracted using different solvents were used for the qualitative analysis of alkaloids, flavonoids, terpenoids and saponins followed by HPTLC finger printing and GC-MS analysis. The components identified in the GC-MS were docked with estrogen receptor (ER) to identify the binding specificity of isolated compounds. The ethyl acetate extract of S. herbaceashowed the presence of high number of secondary metabolites when compared to other solvent system. The qualitative analysis of the plant material also showed the presence of carbohydrates, protein, amino acid, phenol, flavonoids, terpenoids, glycosides, saponins and steroids. The HPTLC finger printing analysis revealed the existence of alkaloid, flavonoid, terpenoid and saponin compounds and GC-MS. GC-MS was performed to identify the phytocomponents constituents in the extract. 8 phytocompounds were identified to analyse binding with ER. The binding affinity score (-6.8 kcal/mol) and interacting ER residues (28) the phyto compound di-n-octyl phthalate showed best docking score with ER α than the standard drugs lasofoxifene, and 4-hydroxytamoxifen. The binding affinity and number of interacting ER residues was -6.9 kcal/mol; 10 and -6.2; 11, respectively. The results identified the presence of ER antagonist in S. herbaceaand warrants further investigation to explore for treating ER regulated diseases.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Effect of proline residues on protein folding

Conformational energy calculations have been used to study the role of the proline residues in the folding of bovine pancreatic trypsin inhibitor. In the calculation, each of the four proline residues of this small protein is forced from the trans to cis peptide isomer while still part of the native folded structure. The cis proline residue can always be accommodated by small changes of the native conformation (< 1 Å root-mean-square deviation). For three of the four proline residues, Pro2, Pro9 and Pro 13, being in the cis form is calculated to destabilize the folded conformation by less than 11 kcal/mol, suggesting that rapid folding to a stable native-like conformation can occur with either isomeric form. For one of these three, Pro13, the destabilization is only 1 kcal/mol, suggesting the existence of an alternative folded native conformation with Pro13 cis. The fourth proline residue, Pro8, is calculated to destabilize the native conformation by so much (33 kcal/mol) that it will block folding in the manner proposed by Brandts et al. (1975).     

Design,synthesis and biological evaluation of 4-aminoquinoline-guanylthiourea derivatives as antimalarial agents

Guanylthiourea (GTU) has been identified as an important antifolate antimalarial pharmacophore unit, whereas, 4-amino quinolones are already known for antimalarial activity. In the present work molecules carrying 4-aminoquinoline and GTU moiety have been designed using molecular docking analysis with PfDHFR enzyme and heme unit. The docking results indicated that the necessary interactions (Asp54 and Ile14) and docking score (−9.63 to −7.36 kcal/mmol) were comparable to WR99210 (−9.89 kcal/mol). From these results nine molecules were selected for synthesis. In vitro analysis of these synthesized compounds reveal that out of the nine molecules, eight show antimalarial activity in the range of 0.61–7.55 μM for PfD6 strain and 0.43–8.04 μM for PfW2 strain. Further, molecular dynamics simulations were performed on the most active molecule to establish comparative binding interactions of these compounds and reference ligand with Plasmodium falciparum dihydrofolate reductase (PfDHFR).     

In silico identification and study of potential anti-mosquito juvenile hormone binding protein (MJHBP) compounds as candidates for dengue virus - Vector insecticides

Dengue has become a huge global health burden. It is currently recognized as the most rapidly spreading mosquito-borne viral disease. Yet, there are currently no licensed vaccines or specific therapeutics to manage the virus, thus, scaling up vector control approaches is important in controlling this viral spread. This study aimed to identify and study in silico, potential anti-mosquito compounds targeting Juvenile hormone (JH) mediated pathways via the Mosquito Juvenile Hormone Binding Protein (MJHBP). The study was implemented using series of computational methods. The query compounds included pyrethroids and those derived from ZINC and ANPDB databases using a simple pharmacophore model in Molecular Operating Environment (MOE). Molecular docking of selected compounds’ library was implemented in MOE. The resultant high-score compounds were further validated by molecular dynamics simulation via Maestro 12.3 module and the respective Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) binding energies computed. The study identified compounds-pyrethroids, natural and synthetic - with high docking energy scores (ranging from 10.91–12.34 kcal/mol). On further analysis of the high-ranking (in terms of docking scores) compounds using MD simulation, the compounds - Ekeberin D4, Maesanin, Silafluofen and ZINC16919139- revealed very low binding energies (?122.99, ?72.91 -104.50 and,-74.94 kcal/mol respectively), fairly stable complex and interesting interaction with JH-binding site amino acid residues on MJHBP. Further studies can explore these compounds in vitro/in vivo in the search for more efficient mosquito vector control.     

Sequence and Structure Based Binding Prediction Study of HLA Class I and cTAP Binding Peptides for Japanese Encephalitis Vaccine Development

Japanese encephalitis is a major threat in developing countries, even the availability of several conventional vaccines, which demand development of more effective vaccines. The present study used propred I and Immune Epitope Database Artificial Neural Network (ANN) algorithm (IEDB-ANN) to identify the conserve and promiscuous T cell epitopes from JEV proteome followed by structure based analysis of potential epitopes. Among all identified 102 epitopes, ten epitope were promiscuous but two epitopes of glycoprotein viz. ⁵⁵LVTVNPFVA⁶³ and ³⁸IPIVSVASL⁴⁶ were found most promiscuous, highly conserved and high population coverage in comparison of known antigenic positive control peptides. The B cell epitopes of glycoprotein also share these two T cell epitopes revealed by BCPred algorithm which can be a basis to confer the protection by neutralizing antibody combined with an effective cell-mediated response. Further, Autodock 4.2 and NAMD–VMD molecular dynamics simulation were used for docking and molecular dynamics simulation respectively, to validate epitope and allele complex binding stability. The 3D structure models were generated for epitopes and corresponding HLA allele by Pepstr and Modeller 9.10 respectively. Epitope LVTVNPFVA–B5101 allele complex showed best energy minimization and stability over the time window during simulation. Here we also present the binding sequel of epitope LVTVNPFVA and its eventual transport through cTAP1 (PDB ID: 1JJ7) revealed by Autodock 4.2, which is an essential path for HLA class I binding epitopes to elicit immune response. The docking experiment of epitope LVTVNPFVA and cTAP1 very well show a 2 H-bond with a binding energy of ?1.88 kcal/mol and other binding state of epitope forming no H-bond with a binding energy of ?1.13 kcal/mol in the lower area of cTAP1 cavity. These results show a smooth pass through of the epitope across the channel of cTAP1. Overall, identified peptides have potential application in the design and development of short peptide based vaccines and diagnostic agents for Japanese encephalitis.     

Contributions of engineered surface salt bridges to the stability of T4 lysozyme determined by directed mutagenesis.

Six designed mutants of T4 lysozyme were created in an attempt to create putative salt bridges on the surface of the protein. The first three of the mutants, T115E (Thr 115 to Glu), Q123E, and N144E, were designed to introduce a new charged side chain close to one or more existing charged groups of the opposite sign on the surface of the protein. In each of these cases the putative electrostatic interactions introduced by the mutation include possible salt bridges between residues within consecutive turns of an alpha-helix. Effects of the mutations ranged from no change in stability to a 1.5 degrees C (0.5 kcal/mol) increase in melting temperature. In two cases, secondary (double) mutants were constructed as controls in which the charge partner was removed from the primary mutant structure. These controls proteins indicate that the contributions to stability from each of the engineered salt bridges is very small (about 0.1-0.25 kcal/mol in 0.15 M KCl). The structures of the three primary mutants were determined by X-ray crystallography and shown to be essentially the same as the wild-type structure except at the site of the mutation. Although the introduced charges in the T115E and Q123E structures are within 3-5 A of their intended partner, the introduced side chains and their intended partners were observed to be quite mobile. It has been shown that the salt bridge between His 31 and Asp 70 in T4 lysozyme stabilizes the protein by 3-5 kcal/mol [Anderson, D. E., Becktel, W. J., & Dahlquist, F. W. (1990) Biochemistry 29, 2403-2408]. To test the effectiveness of His...Asp interactions in general, three additional double mutants, K60H/L13D, K83H/A112D, and S90H/Q122D, were created in order to introduce histidine-aspartate charge pairs on the surface of the protein. Each of these mutants destabilizes the protein by 1-3 kcal/mol in 0.15 M KCl at pH values from 2 to 6.5. The X-ray crystallographic structure of the mutant K83H/A112D has been determined and shows that there are backbone conformational changes of 0.3-0.6 A extending over several residues. The introduction of the histidine and aspartate presumably introduces strain into the folded protein that destabilizes this variant. It is concluded that pairs of oppositely charged residues that are on the surface of a protein and have freedom to adopt different conformations do not tend to come together to form structurally localized salt bridges. Rather, such residues tend to remain mobile, interact weakly if at all, and do not contribute significantly to protein stability. It is argued that the entropic cost of localizing a pair of solvent-exposed charged groups on the surface of a protein largely offsets the interaction energy expected from the formation of a defined salt bridge. There are examples of strong salt bridges in proteins, but such interactions require that the folding of the protein provides the requisite driving energy to hold the interacting partners in the correct rigid alignment.     

Computational study of bindings of HK20 Fab and D5 Fab to HIV-1 gp41

Antibodies HK20 and D5 have been shown to target HIV-1 gp41, thereby inhibiting membrane fusion that facilitates viral entry. The binding picture is static, based on the X-ray crystal structures of the Fab regions and gp41 mimetic five-helix bundle. In this study, we carried out molecular dynamics simulation to provide the dynamic binding picture. Calculated binding free energies are within reasonable range of and follow the trend of the experimental values: -15.28 kcal/mol for HK20 Fab (expt. -11.60 kcal/mol) and -17.90 kcal/mol for D5 Fab (expt. -11.70 kcal/mol). Alanine scanning at protein-protein interface reveals that the highest contributors to binding for HK20 Fab are F54 and I56, both of V(H) region, as well as R30' of V(L) region; whereas for D5 Fab, F54 of V(H) region, as well as W32' and Y94' of V(L) region. HK20 F54 and I56, as well as D5 I52, F54, and T56, bind to the gp41 hydrophobic binding pocket, an important region targeted by many other fusion inhibitors. Hydrogen bonding analysis also identifies high-occupancy hydrogen bonds at the periphery of gp41 hydrophobic pocket. Considering that almost all interface residues are turn residues, further work may be directed to turn mimics. Pre-orientation by the hydrogen bonds to poise this particular turn towards the binding pocket may also be a point worth pursuing.     

Computer aided prediction and identification of potential epitopes in the receptor binding domain (RBD) of spike (S) glycoprotein of MERS-CoV

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) belongs to the coronaviridae family. In spite of several outbreaks in the very recent years, no vaccine against this deadly virus is developed yet. In this study, the receptor binding domain (RBD) of Spike (S) glycoprotein of MERS-CoV was analyzed through Computational Immunology approach to identify the antigenic determinants (epitopes). In order to do so, the sequences of S glycoprotein that belong to different geographical regions were aligned to observe the conservancy of MERS-CoV RBD. The immune parameters of this region were determined using different in silico tools and Immune Epitope Database (IEDB). Molecular docking study was also employed to check the affinity of the potential epitope towards the binding cleft of the specific HLA allele. The N-terminus RBD (S367-S606) of S glycoprotein was found to be conserved among all the available strains of MERS-CoV. Based on the lower IC₅₀ value, a total of eight potential T-cell epitopes and 19 major histocompatibility complex (MHC) class-I alleles were identified for this conserved region. A 9-mer epitope CYSSLILDY displayed interactions with the maximum number of MHC class-I molecules and projected the highest peak in the B-cell antigenicity plot which concludes that it could be a better choice for designing an epitope based peptide vaccine against MERSCoV considering that it must undergo further in vitro and in vivo experiments. Moreover, in molecular docking study, this epitope was found to have a significant binding affinity of -8.5 kcal/mol towards the binding cleft of the HLA-C*12:03 molecule.     

Lysophosphatidic acid and lipopolysaccharide bind to the PIP2-binding domain of gelsolin

The binding of the gelsolin P2 peptide (residues 150-169) with lysophosphatidic acid (LPA) and lipopolysaccharide (LPS) was investigated by isothermal titration calorimetry. P2 binds to LPS with higher affinity than to LPA. For the interaction of 1-oleoyl-LPA with P2 in the absence of salt, K_d and ΔH° were 920 nM and −2.07 kcal/mol, respectively, at pH 7.4 and 25 °C. For the interaction of lipopolysaccharide (LPS) from P. aeruginosa with P2 under the same conditions, K_d was 177 nM and ΔH° was −7.6 kcal/mol.     

Enhanced Amino Acid Selection in Fully Evolved Tryptophanyl-tRNA Synthetase,Relative to Its Urzyme,Requires Domain Motion Sensed by the D1 Switch,a Remote Dynamic Packing Motif

We previously showed (Li, L., and Carter, C. W., Jr. (2013) J. Biol. Chem. 288, 34736–34745) that increased specificity for tryptophan versus tyrosine by contemporary Bacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) over that of TrpRS Urzyme results entirely from coupling between the anticodon-binding domain and an insertion into the Rossmann-fold known as Connecting Peptide 1. We show that this effect is closely related to a long range catalytic effect, in which side chain repacking in a region called the D1 Switch, accounts fully for the entire catalytic contribution of the catalytic Mg²⁺ ion. We report intrinsic and higher order interaction effects on the specificity ratio, (k_cat/K_m)_Trp/(k_cat/K_m)_Tyr, of 15 combinatorial mutants from a previous study (Weinreb, V., Li, L., and Carter, C. W., Jr. (2012) Structure 20, 128–138) of the catalytic role of the D1 Switch. Unexpectedly, the same four-way interaction both activates catalytic assist by Mg²⁺ ion and contributes −4.4 kcal/mol to the free energy of the specificity ratio. A minimum action path computed for the induced-fit and catalytic conformation changes shows that repacking of the four residues precedes a decrease in the volume of the tryptophan-binding pocket. We suggest that previous efforts to alter amino acid specificities of TrpRS and glutaminyl-tRNA synthetase (GlnRS) by mutagenesis without extensive, modular substitution failed because mutations were incompatible with interdomain motions required for catalysis.     

Antigenic structure of human hepatitis A virus defined by analysis of escape mutants selected against murine monoclonal antibodies.

We examined the antigenic structure of human hepatitis A virus (HAV) by characterizing a series of 21 murine monoclonal-antibody-resistant neutralization escape mutants derived from the HM175 virus strain. The escape phenotype of each mutant was associated with reduced antibody binding in radioimmunofocus assays. Neutralization escape mutations were identified at the Asp-70 and Gln-74 residues of the capsid protein VP3, as well as at Ser-102, Val-171, Ala-176, and Lys-221 of VP1. With the exception of the Lys-221 mutants, substantial cross-resistance was evident among escape mutants tested against a panel of 22 neutralizing monoclonal antibodies, suggesting that the involved residues contribute to epitopes composing a single antigenic site. As mutations at one or more of these residues conferred resistance to 20 of 22 murine antibodies, this site appears to be immunodominant in the mouse. However, multiple mutants selected independently against any one monoclonal antibody had mutations at only one or, at the most, two amino acid residues within the capsid proteins, confirming that there are multiple epitopes within this antigenic site and suggesting that single-amino-acid residues contributing to these epitopes may play key roles in the binding of individual antibodies. A second, potentially independent antigenic site was identified by three escape mutants with different substitutions at Lys-221 of VP1. These mutants were resistant only to antibody H7C27, while H7C27 effectively neutralized all other escape mutants. These data support the existence of an immunodominant neutralization site in the antigenic structure of hepatitis A virus which involves residues of VP3 and VP1 and a second, potentially independent site involving residue 221 of VP1.     

Effects of five-tryptophan mutations on structure, stability and function of Escherichia coli dihydrofolate reductase

To elucidate the roles of tryptophan residues in the structure, stability, and function of Escherichia coli dihydrofolate reductase (DHFR), its five tryptophan residues were replaced by site-directed mutagenesis with leucine, phenylalanine or valine (W22F, W22L, W30L, W47L, W74F, W74L, W133F, and W133V). Far-ultraviolet circular dichroism (CD) spectra of these mutants reveal that exciton coupling between Trp47 and Trp74 strongly affects the peptide CD of wild-type DHFR, and that Trp133 also contributes appreciably. No additivity was observed in the contributions of individual tryptophan residues to the fluorescence spectrum of wild-type DHFR, Trp74 having a dominant effect. These single-tryptophan mutations induce large changes in the free energy of urea unfolding, which showed values of 1.79-7.14 kcal/mol, compared with the value for wild-type DHFR of 6.08 kcal/mol. Analysis of CD and fluorescence spectra suggests that thermal unfolding involves an intermediate with the native-like secondary structure, the disrupted Trp47-Trp74 exciton coupling, and the solvent-exposed Trp30 and Trp47 side chains. All the mutants except W22L (13%) retain more than 50% of the enzyme activity of wild-type DHFR. These results demonstrate that the five tryptophan residues of DHFR play important roles in its structure and stability but do not crucially affect its enzymatic function.     

Investigating the effects of point mutations on the affinity between the cyanobacterial lectin microvirin and high mannose-type glycans present on the HIV envelope glycoprotein

Human immunodeficiency virus (HIV) infections continue to exert an enormous impact on global human health. This led experts to emphasize the importance of new measures for preventing HIV infections, including the development of vaccines and novel drugs. In this context, a promising approach involves the use of lectins that can bind the surface envelope glycoprotein gp120 of HIV with high affinity, preventing viral entry. The cyanobacterial lectin microvirin (MVN) has been proposed as a candidate for development as a topical microbicide because of its ability to bind to high mannose-type glycans, potently inhibiting HIV-1 entry. Thus, the aim of this computational study was to investigate the effects of four point mutations (D53Q, D53E, D53K, and D53W) on the structure and affinity of MVN with di-mannose (MAN). Molecular dynamics simulations followed by binding free energy calculations using MM-GBSA were employed. The calculated binding free energy of ligand-receptor complexation of MVN with MAN was ?26.02 kcal mol^-1. We identified in the wild-type protein that residues I45, T59, and Q81 have a major contribution to the binding free energy of di-mannose. Among the investigated mutants, the most promising one was the D53W mutation, with a theoretical binding free energy value of ?29.16 kcal mol^-1. We suggest that this increased stability is due to the introduction of extra rigidity on the hinge region connecting two key structural elements of the MVN binding site.