Conformational gating of dimannose binding to the antiviral protein cyanovirin revealed from the crystal structure at 1.35 A resolution

Cyanovirin (CV-N) is a small lectin with potent HIV neutralization activity, which could be exploited for a mucosal defense against HIV infection. The wild-type (wt) protein binds with high affinity to mannose-rich oligosaccharides on the surface of gp120 through two quasi-symmetric sites, located in domains A and B. We recently reported on a mutant of CV-N that contained a single functional mannose-binding site, domain B, showing that multivalent binding to oligomannosides is necessary for antiviral activity. The structure of the complex with dimannose determined at 1.8 A resolution revealed a different conformation of the binding site than previously observed in the NMR structure of wt CV-N. Here, we present the 1.35 A resolution structure of the complex, which traps three different binding conformations of the site and provides experimental support for a locking and gating mechanism in the nanoscale time regime observed by molecular dynamics simulations.     

On Calculation of the Electrostatic Potential of a Phosphatidylinositol Phosphate-Containing Phosphatidylcholine Lipid Membrane Accounting for Membrane Dynamics

Many signaling events require the binding of cytoplasmic proteins to cell membranes by recognition of specific charged lipids, such as phosphoinositol-phosphates. As a model for a protein-membrane binding site, we consider one charged phosphoinositol phosphate (PtdIns(3)P) embedded in a phosphatidylcholine bilayer. As the protein-membrane binding is driven by electrostatic interactions, continuum solvent models require an accurate representation of the electrostatic potential of the phosphoinositol phosphate-containing membrane. We computed and analyzed the electrostatic potentials of snapshots taken at regular intervals from molecular dynamics simulations of the bilayer. We observe considerable variation in the electrostatic potential of the bilayer both along a single simulation and between simulations performed with the GAFF or CHARMM c36 force fields. However, we find that the choice of GAFF or CHARMM c36 parameters has little effect on the electrostatic potential of a given configuration of the bilayer with a PtdIns(3)P embedded in it. From our results, we propose a remedian averaging method for calculating the electrostatic potential of a membrane system that is suitable for simulations of protein-membrane binding with a continuum solvent model.     

NMR solution structure of a cyanovirin homolog from wheat head blight fungus

Members of the cyanovirin-N homolog (CVNH) lectin family are found in bacteria, fungi and plants. As part of our ongoing work on CVNH structure-function studies, we determined the high-resolution NMR solution structure of the homolog from the wheat head blight disease causing ascomycetous fungus Gibberella zeae (or Fusarium graminearum), hereafter called GzCVNH. Like cyanovirin-N (CV-N), GzCVNH comprises two tandem sequence repeats and the protein sequence exhibits 30% identity with CV-N. The overall structure is similar to those of other members of the CVNH family, with the conserved pseudo-symmetric halves of the structure, domains A and B, closely resembling recently determined structures of Tuber borchii, Neurospora crassa, and Ceratopteris richardii CVNH proteins. Although GzCVNH exhibits a similar glycan recognition profile to CV-N and specifically binds to Manα(1-2)Manα, its weak carbohydrate binding affinity to only one binding site is insufficient for conferring anti-HIV activity.     

Multiple drugs and multiple targets: An analysis of the electrostatic determinants of binding between non‐nucleoside HIV‐1 reverse transcriptase inhibitors and variants of HIV‐1 RT

We present a systematic, computational analysis of the electrostatic component of binding of three HIV‐1 RT inhibitors—nevirapine (NVP), efavirenz (EFV), and the recently approved rilpivirine (RPV)—to wild‐type (WT) and mutant variants of RT. Electrostatic charge optimization was applied to determine how suited each molecule's charge distribution is for binding WT and individual mutants of HIV‐1 RT. Although the charge distributions of NVP and EFV are rather far from being optimal for tight binding, RPVs charge distribution is close to the theoretical, optimal charge distribution for binding WT HIV‐1 RT, although slight changes in charge can dramatically impact binding energetics. Moreover, toward the L100I/K103N double mutant, RPVs charge distribution is quite far from optimal. We also determine the contributions of chemical moieties on each molecule toward the electrostatic component of binding and show that different regions of a drug molecule may be used for recognition by different RT variants. The electrostatic contributions of certain RT residues toward drug binding are also computed to highlight critical residues for each interaction. Finally, the charge distribution of RPV is optimized to promiscuously bind to three RT variants rather than to each one in turn, with the resulting charge distribution being a compromise between the optimal charge distributions to each individual variant. Taken together, this work demonstrates that even in a binding site considered quite hydrophobic, electrostatics play a subtle yet varying role that must be considered in designing next‐generation molecules that recognize rapidly mutating targets. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Multiple antiviral activities of cyanovirin-N: blocking of human immunodeficiency virus type 1 gp120 interaction with CD4 and coreceptor and inhibition of diverse enveloped viruses

Cyanovirin-N (CV-N) is a cyanobacterial protein with potent neutralizing activity against human immunodeficiency virus (HIV). CV-N has been shown to bind HIV type 1 (HIV-1) gp120 with high affinity; moreover, it blocks the envelope glycoprotein-mediated membrane fusion reaction associated with HIV-1 entry. However, the inhibitory mechanism(s) remains unclear. In this study, we show that CV-N blocked binding of gp120 to cell-associated CD4. Consistent with this, pretreatment of gp120 with CV-N inhibited soluble CD4 (sCD4)-dependent binding of gp120 to cell-associated CCR5. To investigate possible effects of CV-N at post-CD4 binding steps, we used an assay that measures sCD4 activation of the HIV-1 envelope glycoprotein for fusion with CCR5-expressing cells. CV-N displayed equivalently potent inhibitory effects when added before or after sCD4 activation, suggesting that CV-N also has blocking action at the level of gp120 interaction with coreceptor. This effect was shown not to be due to CV-N-induced coreceptor down-modulation after the CD4 binding step. The multiple activities against the HIV-1 envelope glycoprotein prompted us to examine other enveloped viruses. CV-N potently blocked infection by feline immunodeficiency virus, which utilizes the chemokine receptor CXCR4 as an entry receptor but is CD4 independent. CV-N also inhibited fusion and/or infection by human herpesvirus 6 and measles virus but not by vaccinia virus. Thus, CV-N has broad-spectrum antiviral activity, both for multiple steps in the HIV entry mechanism and for diverse enveloped viruses. This broad specificity has implications for potential clinical utility of CV-N.     

Substrate recognition by ribosome-inactivating protein studied by molecular modeling and molecular electrostatic potentials

A computer model of dianthin 30, a type 1 ribosome-inactivating protein (RIP), is constructed by homology modeling using two known X-ray structures; a type 1 RIP, pokeweed antiviral protein (PAP), and chain A of a type 2 RIP, ricin. The 3D structure is refined by molecular dynamics and its binding site compared with those of PAP and ricin using molecular electrostatic potential mapping. The differences in the maps obtained clearly show how, despite the similarity of the topology of the binding site, differences in electrostatic potential can account for the experimentally observed differences in substrate recognition and binding. This demonstrates the potential of these techniques for guiding further experimental analyses.     

Architecture of the sugar binding sites in carbohydrate binding proteins-a computer modeling study

Different sugars, Gal, GalNAc and Man were docked at the monosaccharide binding sites of Erythrina corallodenron (EcorL), peanut lectin (PNA), Lathyrus ochrus (LOLI), and pea lectin (PSL). To study the lectin-carbohydrate interactions, in the complexes, the hydroxymethyl group in Man and Gal favors, gg and gt conformations respectively, and is the dominant recognition determination. The monosaccharide binding site in lectins that are specific to Gal/GalNAc is wider due to the additional amino acid residues in loop D as compared to that in lectins specific to Man/Glc, and affects the hydrogen bonds of the sugar involving residues from loop D, but not its orientation in the binding site. The invariant amino acid residues Asp from loop A, and Asn and an aromatic residue (Phe or Tyr) in loop C provides the basic architecture to recognize the common features in C4 epimers. The invariant Gly in loop B together with one or two residues in the variable region of loop D/A holds the sugar tightly at both ends. Loss of any one of these hydrogen bonds leads to weak interaction. While the subtle variations in the sequence and conformation of peptide fragment that resulted due to the size and location of gaps present in amino acid sequence in the neighborhood of the sugar binding site of loop D/A seems to discriminate the binding of sugars which differ at C4 atom (galacto and gluco configurations). The variations at loop B are important in discriminating Gal and GalNAc binding. The present study thus provides a structural basis for the observed specificities of legume lectins which uses the same four invariant residues for binding. These studies also bring out the information that is important for the design/engineering of proteins with the desired carbohydrate specificity.     

Atomic mapping of the sugar interactions in one-site and two-site mutants of cyanovirin-N by NMR spectroscopy

The details of the interaction between two mutants of Cyanovirin-N (CV-N), an HIV inactivating protein, and di- and trimannosides, substructures of Man-9, were investigated by STD NMR spectroscopy. One mutant, CV-N (mutDB), contains only one carbohydrate-binding site on domain A, whereas in CV-N (mutDA), the specificity of domain A for trimannose was changed while the site in domain B was kept intact, allowing for a dissection of the overall binding. Results of the STD NMR experiments revealed close contact between the protein binding site on domain A and H2, H3, and H4 of the nonreducing terminal mannose unit for Manalpha(1-2)Manalpha OMe, Manalpha(1-2)Manalpha(1-3)Manalpha OMe, and Manalpha(1-2)Manalpha(1-6)Manalpha OMe. The Manalpha(1-2)Manalpha(1-2)Manalpha OMe trisaccharide interacted with CV-N with the highest affinity. Further dissection of the interaction was achieved by NMR experiments with synthetic 2'-, 3'-, 4'-, and 6'-deoxy analogues of the disaccharide Manalpha(1-2)Manalpha OMe. STD and (1)H- (15)N HSQC NMR spectroscopy revealed that the 2'- and 6'-deoxy dimannosides were recognized by CV-N, whereas no binding was detected for the 3'- and 4'-deoxy sugars. These results demonstrate that the 3'- and 4'-hydroxyl groups on the terminal residue are engaged in key polar interactions with the protein and are required for high-affinity binding.     

The domain-swapped dimer of cyanovirin-N contains two sets of oligosaccharide binding sites in solution

The binding of high-mannose oligosaccharides to the domain-swapped dimeric form of the potent HIV-inactivating protein cyanovirin-N (CV-N) was investigated in solution by NMR, complementing recent structural studies by X-ray crystallography on similar complexes [J. Biol. Chem. 277 (2002) 34336]. The crystal structures of CV-N dimer complexed with Man-9 and hexamannoside revealed two carbohydrate binding sites on opposite ends of the molecule. No binding was observed at site 1, previously identified on the solution monomer of CV-N [Structure 9 (2001) 931; Shenoy et al., Chem. Biol. 9 (2002) 1109]. Here, we report the presence of four sugar binding sites on the CV-N dimer in solution, identified by chemical shift mapping with hexamannoside and nonamannoside, synthetic substructures of Man-9. Our results demonstrate that in solution the domain-swapped CV-N dimer, like the CV-N monomer, contains two types of sites that are available for carbohydrate binding, suggesting that the occlusion of the primary sites in the crystal is due to specific features of the solid state.     

New insight on beta-lactoglobulin binding sites by 1-anilinonaphthalene-8-sulfonate fluorescence decay

The fluorescence time decay parameters of the beta-lactoglobulin-1-anilinonaphthalene-8-sulfonate complex have been investigated under physical and chemical perturbations (2 < pH < 8 and added electrolyte 0 < NaCl < 0.5 M) to obtain new insight on the nature of the protein binding interactions. A double exponential decay of the bound probe lifetime has been confirmed by the presence of a longer component, 11 to 14.5 ns, and a shorter component, 2.5 to 3.5 ns. The two lifetimes are ascribed to different binding modes associated also with different exposure to the solvent; in particular, the longer component is attributed to binding inside the hydrophobic beta barrel, while a "surface" site is suggested for the shorter component. A detailed analysis of the lifetime fractional intensities correlates the binding constants with ionic strength and supports the presence of electrostatic effects at both sites. A Debye-Hückel approach, applied to extrapolate the electrostatic free energy contribution vs. pH at vanishing ionic strength, gives interesting clues on the effective charge felt by the ANS ligands in the proximity of each site. In particular, binding is found to parallel the aspartate and glutamate titrations between pH 3 and pH 4.5; the "surface" site mainly responds to the presence of these local titrating charges while the "internal" site more closely follows the overall protein net charge.     

An application of CIFAP for predicting the binding affinity of Chk1 inhibitors derived from 2‐aminothiazole‐4‐carboxamide

Investigation of protein‐ligand interactions obtained from experiments has a crucial part in the design of newly discovered and effective drugs. Analyzing the data extracted from known interactions could help scientists to predict the binding affinities of promising ligands before conducting experiments. The objective of this study is to advance the CIFAP (compressed images for affinity prediction) method, which is relevant to a protein‐ligand model, identifying 2D electrostatic potential images by separating the binding site of protein‐ligand complexes and using the images for predicting the computational affinity information represented by pIC₅₀ values. The CIFAP method has 2 phases, namely, data modeling and prediction. In data modeling phase, the separated 3D structure of the binding pocket with the ligand inside is fitted into an electrostatic potential grid box, which is then compressed through 3 orthogonal directions into three 2D images for each protein‐ligand complex. Sequential floating forward selection technique is performed for acquiring prediction patterns from the images. In the prediction phase, support vector regression (SVR) and partial least squares regression are used for testing the quality of the CIFAP method for predicting the binding affinity of 45 CHK1 inhibitors derived from 2‐aminothiazole‐4‐carboxamide. The results show that the CIFAP method using both support vector regression and partial least squares regression is very effective for predicting the binding affinities of CHK1‐ligand complexes with low‐error values and high correlation. As a future work, the results could be improved by working on the pose of the ligands inside the grid.     

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of Uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k_cat and reduced K_m relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K_m) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

Molecular docking of superantigens with class II major histocompatibility complex proteins

The molecular recognition of two superantigens with class II major histocompatibility complex molecules was simulated by using protein– protein docking. Superantigens studied were staphylococcal enterotoxin B (SEB) and toxic shock syndrome toxin-1 (TSST-1) in their crystallographic assemblies with HLA-DR1. Rigid-body docking was performed sampling configurational space of the interfacial surfaces by employing a strategy of partitioning the contact regions on HLA-DR1 into separate molecular recognition units. Scoring of docked conformations was based on an electrostatic continuum model evaluated with the finite-difference Poisson– Boltzmann method. Estimates of nonpolar contributions were derived from the buried molecular surface areas. We found for both superantigens that docking the HLA-DR1 surface complementary with the SEB and TSST-1 contact regions containing a homologous hydrophobic surface loop provided sufficient recognition for the reconstitution of native-like conformers exhibiting the highest-scoring free energies. For the SEB complex, the calculations were successful in reproducing the total association free energy. A comparison of the free-energy determinants of the conserved hydrophobic contact residue indicates functional similarity between the two proteins for this interface. Though both superantigens share a common global association mode, differences in binding topology distinguish the conformational specificities underlying recognition. © 1998 John Wiley & Sons, Ltd.     

Electrostatic potential of nucleotide-free protein is sufficient for discrimination between adenine and guanine-specific binding sites

Despite sharing many common features, adenine-binding and guanine-binding sites in proteins often show a clear preference for the cognate over the non-cognate ligand. We have analyzed electrostatic potential (ESP) patterns at adenine and guanine-binding sites of a large number of non-redundant proteins where each binding site was first annotated as adenine/guanine-specific or non-specific from a survey of primary literature. We show that more than 90% of ESP variance at the binding sites is accounted for by only two principal component ESP vectors, each aligned to molecular dipoles of adenine and guanine. Projected on these principal component vectors, the adenine/guanine-specific and non-specific binding sites, including adenine-containing dinucleotides, show non-overlapping distributions. Adenine or guanine specificities of the binding sites also show high correlation with the corresponding electrostatic replacement (cognate by non-cognate ligand) energies. High correlation coefficients (0.94 for 35 adenine-binding sites and 1.0 for 20 guanine-binding sites) were obtained when adenine/guanine specificities were predicted using the replacement energies. Our results demonstrate that ligand-free protein ESP is an excellent indicator for discrimination between adenine and guanine-specific binding sites and that ESP of ligand-free protein can be used as a tool to annotate known and putative purine-binding sites in proteins as adenine or guanine-specific.     

The carboxy terminus of yeast Atg13 binds phospholipid membrane via motifs that overlap with the Vac8-interacting domain

ABSTRACT

Macroautophagy/autophagy is a conserved catabolic recycling pathway involving the sequestration of cytoplasmic components within double-membrane vesicles termed autophagosomes. The autophagy-related (Atg) protein Atg13 is a key member of the autophagy initiation complex. The Atg13 C terminus is an intrinsically disordered region (IDR) harboring a binding site for the vacuolar membrane protein Vac8. Recent reports suggest Atg13 acts as a hub to assemble the initiation complex, and also participates in membrane recognition. Here we show that the Atg13 C terminus directly binds to lipid membranes via electrostatic interactions between positively charged residues in Atg13 and negatively charged phospholipids as well as a hydrophobic insertion of a Phe residue. We identified 2 sets of residues in the Atg13 IDR that affect its phospholipid-binding properties; these residues overlap with the Vac8-binding domain of Atg13. Our data indicate that Atg13 binding to phospholipids and Vac8 is mutually exclusive, and both are required for efficient autophagy.     

Identification and characterization of peptides that bind to cyanovirin-N, a potent human immunodeficiency virus-inactivating protein

Cyanovirin-N (CV-N) exerts a potent human immunodeficiency virus (HIV)-inactivating activity against diverse strains of HIV by binding to the viral surface envelope glycoprotein gp120 and blocking its essential interactions with cellular receptors. Based on previous thermodynamic analyses, it has been speculated that discrete protein-protein interactions might play an important ancillary role in the CV-N/gp120 binding event, in addition to the interactions of CV-N with specific oligosaccharides present on gp120. Here, we report the identification and characterization of CV-N-binding peptides, which were isolated by screening of M13 phage-displayed peptide libraries. After performing three rounds of biopanning of the libraries against biotinylated CV-N, a CV-N-binding motif, X3CX6(W/F)(Y/F)CX2(Y/F), was evident. A vector was designed to express CV-N-binding peptides as a fusion with thioredoxin (Trx) containing a penta-His affinity tag. The CV-N-binding peptides fused with His-tagged Trx inhibited binding of the corresponding peptide-bearing phages to CV-N, confirming that the peptides possessed CV-N-binding activity. Optical biosensor binding studies showed that the one of the CV-N-binding peptide, TN10-1, bound to CV-N with a KD value of 1.9 microM. The results of alanine scanning mutagenesis of the peptide showed that aromatic residues at positions 11, 12, and 16, as well as the conformational structure of the peptide secured by a disulfide bond, were important for the binding interactions. A series of competitive binding assays confirmed that gp120 inhibited CV-N binding of the corresponding peptide-bearing phages, and suggested that TN10-1 peptides were mimicking the protein component of gp120 rather than mimicking specific oligosaccharides present on gp120.     

Interchange of DNA-binding modes in the deformed and ultrabithorax homeodomains: a structural role for the N-terminal arm

The deformed (Dfd) and ultrabithorax (Ubx) homeoproteins regulate developmental gene expression in Drosophila melanogaster by binding to specific DNA sequences within its genome. DNA binding is largely accomplished via a highly conserved helix-turn-helix DNA-binding domain that is known as a homeodomain (HD). Despite nearly identical DNA recognition helices and similar target DNA sequence preferences, the in vivo functions of the two proteins are quite different. We have previously revealed differences between the two HDs in their interactions with DNA. In an effort to define the individual roles of the HD N-terminal arm and recognition helix in sequence-specific binding, we have characterized the structural details of two Dfd/Ubx chimeric HDs in complex with both the Dfd and Ubx-optimal-binding site sequences. We utilized hydroxyl radical cleavage of DNA to assess the positioning of the proteins on the binding sites. The effects of missing nucleosides and purine methylation on HD binding were also analyzed. Our results show that both the Dfd and Ubx HDs have similar DNA-binding modes when in complex with the Ubx-optimal site. There are subtle but reproducible differences in these modes that are completely interchanged when the Dfd N-terminal arm is replaced with the corresponding region of the Ubx HD. In contrast, we showed previously that the Dfd-optimal site sequence elicits a very different binding mode for the Ubx HD, while the Dfd HD maintains a mode similar to that elicited by the Ubx-optimal site. Our current methylation interference studies suggest that this alternate binding mode involves interaction of the Ubx N-terminal arm with the minor groove on the opposite face of DNA relative to the major groove that is occupied by the recognition helix. As judged by hydroxyl radical footprinting and the missing nucleoside experiment, it appears that interaction of the Ubx recognition helix with the DNA major groove is reduced. Replacing the Dfd N-terminal arm with that of Ubx does not elicit a complete interchange of the DNA-binding mode. Although the position of the chimera relative to DNA, as judged by hydroxyl radical footprinting, is similar to that of the Dfd HD, the missing nucleoside and methylation interference patterns resemble those of the Ubx HD. Repositioning of amino acid side-chains without wholesale structural alteration in the polypeptide appears to occur as a function of N-terminal arm identity and DNA-binding site sequence. Complete interchange of binding modes was achieved only by replacement of the Dfd N-terminal arm and the recognition helix plus 13 carboxyl-terminal residues with the corresponding residues of Ubx. The position of the N-terminal arm in the DNA minor groove appears to differ in a manner that depends on the two base-pair differences between the Dfd and Ubx-optimal-binding sites. Thus, N-terminal arm position dictates the binding mode and the interaction of the recognition helix with nucleosides in the major groove.     

Focusing of the electrostatic potential at EF-hands of calbindin D(9k): titration of acidic residues

Biological functions for a large class of calmodulin-related proteins, such as target protein activation and Ca(2+) buffering, are based on fine-tuned binding and release of Ca(2+) ions by pairs of coupled EF-hand metal binding sites. These are abundantly filled with acidic residues of so far unknown ionization characteristics, but assumed to be essential for protein function in their ionized forms. Here we describe the measurement and modeling of pK(a) values for all aspartic and glutamic acid residues in apo calbindin D(9k), a representative of calmodulin-related proteins. We point out that while all the acidic residues are ionized predominantly at neutral pH, the onset of proton uptake by Ca(2+) ligands with high pK(a) under these conditions may have functional implications. We also show that the negative electrostatic potential is focused at the bidental Ca(2+) ligand of each site, and that the potential is significantly more negative at the N-terminal binding site.