Local Interactions That Contribute Minimal Frustration Determine Foldability.

Earlier experiments suggest that the evolutionary information (conservation and coevolution) encoded in protein sequences is necessary and sufficient to specify the fold of a protein family. However, there is no computational work to quantify the effect of such evolutionary information on the folding process. Here we explore the role of early folding steps for sequences designed using coevolution and conservation through a combination of computational and experimental methods. We simulated a repertoire of native and designed WW domain sequences to analyze early local contact formation and found that the N-terminal ß-hairpin turn would not form correctly due to strong non-native local contacts in unfoldable sequences. Through a maximum likelihood approach, we identified five local contacts that play a critical role in folding, suggesting that a small subset of amino acid pairs can be used to solve the "needle in the haystack" problem to design foldable sequences. Thus, using the contact probability of those five local contacts that form during the early stage of folding, we built a classification model that predicts the foldability of a WW sequence with 81% accuracy. This classification model was used to redesign WW domain sequences that could not fold due to frustration and make them foldable by introducing a few mutations that led to the stabilization of these critical local contacts. The experimental analysis shows that a redesigned sequence folds and binds to polyproline peptides with a similar affinity as those observed for native WW domains. Overall, our analysis shows that evolutionary-designed sequences should not only satisfy the folding stability but also ensure a minimally frustrated folding landscape.

A Designed "Nested" Dimer of Cyanovirin-N Increases Antiviral Activity.

Cyanovirin-N (CV-N) is an antiviral lectin with potent activity against enveloped viruses, including HIV. The mechanism of action involves high affinity binding to mannose-rich glycans that decorate the surface of enveloped viruses. In the case of HIV, antiviral activity of CV-N is postulated to require multivalent interactions with envelope protein gp120, achieved through a pseudo-repeat of sequence that adopts two near-identical glycan-binding sites, and possibly involves a 3D-domain-swapped dimeric form of CV-N. Here, we present a covalent dimer of CV-N that increases the number of active glycan-binding sites, and we characterize its ability to recognize four glycans in solution. A CV-N variant was designed in which two native repeats were separated by the "nested" covalent insertion of two additional repeats of CV-N, resulting in four possible glycan-binding sites. The resulting Nested CV-N folds into a wild-type-like structure as assessed by circular dichroism and NMR spectroscopy, and displays high thermal stability with a Tm of 59 °C, identical to WT. All four glycan-binding domains encompassed by the sequence are functional as demonstrated by isothermal titration calorimetry, which revealed two sets of binding events to dimannose with dissociation constants Kd of 25 µM and 900 µM, assigned to domains B and B' and domains A and A' respectively. Nested CV-N displays a slight increase in activity when compared to WT CV-N in both an anti-HIV cellular assay and a fusion assay. This construct conserves the original binding specifityies of domain A and B, thus indicating correct fold of the two CV-N repeats. Thus, rational design can be used to increase multivalency in antiviral lectins in a controlled manner.

A flexible docking scheme efficiently captures the energetics of glycan-cyanovirin binding.

Cyanovirin-N (CVN), a cyanobacterial lectin, exemplifies a class of antiviral agents that inhibit HIV by binding to the highly glycosylated envelope protein gp120. Here, we investigate the energetics of glycan recognition using a computationally inexpensive flexible docking approach, backbone perturbation docking (BP-Dock). We benchmarked our method using two mutants of CVN: P51G-m4-CVN, which binds dimannose with high affinity through domain B, and CVN((mutDB)), in which binding to domain B has been abolished through mutation of five polar residues to small nonpolar side chains. We investigated the energetic contribution of these polar residues along with the additional position 53 by docking dimannose to single-point CVN mutant models. Analysis of the docking simulations indicated that the E41A/G and T57A mutations led to a significant decrease in binding energy scores due to rearrangements of the hydrogen-bond network that reverberated throughout the binding cavity. N42A decreased the binding score to a level comparable to that of CVN((mutDB)) by affecting the integrity of the local protein structure. In contrast, N53S resulted in a high binding energy score, similar to P51G-m4-CVN. Experimental characterization of the five mutants by NMR spectroscopy confirmed the binding affinity pattern predicted by BP-Dock. Despite their mostly conserved fold and stability, E41A, E41G, and T57A displayed dissociation constants in the millimolar range. N53S showed a binding constant in the low micromolar range, similar to that observed for P51G-m4-CVN. No binding was observed for N42A. Our results show that BP-Dock is a useful tool for rapidly screening the relative binding affinity pattern of in silico-designed mutants compared with wild-type, supporting its use to design novel mutants with enhanced binding properties.

The antiviral lectin cyanovirin-N: probing multivalency and glycan recognition through experimental and computational approaches.

CVN (cyanovirin-N), a small lectin isolated from cyanobacteria, exemplifies a novel class of anti-HIV agents that act by binding to the highly glycosylated envelope protein gp120 (glycoprotein 120), resulting in inhibition of the crucial viral entry step. In the present review, we summarize recent work in our laboratory and others towards determining the crucial role of multivalency in the antiviral activity, and we discuss features that contribute to the high specificity and affinity for the glycan ligand observed in CVN. An integrated approach that encompasses structural determination, mutagenesis analysis and computational work holds particular promise to clarify aspects of the interactions between CVN and glycans.

Modulation of function in a minimalist heme-binding membrane protein.

De novo designed heme-binding proteins have been used successfully to recapitulate features of natural hemoproteins. This approach has now been extended to membrane-soluble model proteins. Our group designed a functional hemoprotein, ME1, by engineering a bishistidine binding site into a natural membrane protein, glycophorin A (Cordova et al. in J Am Chem Soc 129:512-518, 2007). ME1 binds iron(III) protoporphyrin IX with submicromolar affinity, has a redox potential of -128 mV, and displays peroxidase activity. Here, we show the effect of aromatic residues in modulating the redox potential in the context of a membrane-soluble model system. We designed aromatic interactions with the heme through a single-point mutant, G25F, in which a phenylalanine is designed to dock against the porphyrin ring. This mutation results in roughly tenfold tighter binding to iron(III) protoporphyrin IX (K(d,app) = 6.5 × 10(-8) M), and lowers the redox potential of the cofactor to -172 mV. This work demonstrates that specific design features aimed at controlling the properties of bound cofactors can be introduced in a minimalist membrane hemoprotein model. The ability to modulate the redox potential of cofactors embedded in artificial membrane proteins is crucial for the design of electron transfer chains across membranes in functional photosynthetic devices.