Dissecting contact potentials for proteins: relative contributions of individual amino acids

Mimotopes mimic the three-dimensional topology of an antigen epitope, and are frequently recognized by antibodies with affinities comparable to those obtained for the original antibody-antigen interaction. Peptides and anti-idiotypic antibodies are two classes of protein mimotopes that mimic the topology (but not necessarily the sequence) of the parental antigen. In this study, we combine these two classes by selecting mimotopes based on single domain IgNAR antibodies, which display exceptionally long CDR3 loop regions (analogous to a constrained peptide library) presented in the context of an immunoglobulin framework with adjacent and supporting CDR1 loops. By screening an in vitro phage-display library of IgNAR variable domains (V(NAR)s) against the target antigen monoclonal antibody MAb5G8, we obtained four potential mimotopes. MAb5G8 targets a linear tripeptide epitope (AYP) in the flexible signal sequence of the Plasmodium falciparum Apical Membrane Antigen-1 (AMA1), and this or similar motifs were detected in the CDR loops of all four V(NAR)s. The V(NAR)s, 1-A-2, -7, -11, and -14, were demonstrated to bind specifically to this paratope by competition studies with an artificial peptide and all showed enhanced affinities (3-46 nM) compared to the parental antigen (175 nM). Crystallographic studies of recombinant proteins 1-A-7 and 1-A-11 showed that the SYP motifs on these V(NAR)s presented at the tip of the exposed CDR3 loops, ideally positioned within bulge-like structures to make contact with the MAb5G8 antibody. These loops, in particular in 1-A-11, were further stabilized by inter- and intra- loop disulphide bridges, hydrogen bonds, electrostatic interactions, and aromatic residue packing. We rationalize the higher affinity of the V(NAR)s compared to the parental antigen by suggesting that adjacent CDR1 and framework residues contribute to binding affinity, through interactions with other CDR regions on the antibody, though of course definitive support of this hypothesis will rely on co-crystallographic studies. Alternatively, the selection of mimotopes from a large (<4 x 10(8)) constrained library may have allowed selection of variants with even more favorable epitope topologies than present in the original antigenic structure, illustrating the power of in vivo selection of mimotopes from phage-displayed molecular libraries.     

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Are Current Molecular Dynamics Force Fields too Helical?

Accurate force fields are essential for the success of molecular dynamics simulations. In apparent contrast to the conformational preferences of most force fields, recent NMR experiments suggest that short polyalanine peptides in water populate the polyproline II structure almost exclusively. To investigate this apparent contradiction, with its ramifications for the assessment of molecular force fields and the structure of unfolded proteins, we performed extensive simulations of Ala₅ in water (∼5 μs total time), using twelve different force fields and three different peptide terminal groups. Using either empirical or density-functional-based Karplus relations for the J-couplings, we find that most current force fields do overpopulate the α-region, with quantitative results depending on the choice of Karplus relation and on the peptide termini. Even after reweighting to match experiment, we find that Ala₅ retains significant α- and β-populations. In fact, several force fields match the experimental data well before reweighting and have a significant helical population. We conclude that radical changes to the best current force fields are not necessary, based on the NMR data. Nevertheless, experiments on short peptides open the way toward the systematic improvement of current simulation models.     

Development of novel statistical potentials for protein fold recognition

The need to perform large-scale studies of protein fold recognition, structure prediction and protein-protein interactions has led to novel developments of residue-level minimal models of proteins. A minimum requirement for useful protein force-fields is that they be successful in the recognition of native conformations. The balance between the level of detail in describing the specific interactions within proteins and the accuracy obtained using minimal protein models is the focus of many current protein studies. Recent results suggest that the introduction of explicit orientation dependence in a coarse-grained, residue-level model improves the ability of inter-residue potentials to recognize the native state. New statistical and optimization computational algorithms can be used to obtain accurate residue-dependent potentials for use in protein fold recognition and, more importantly, structure prediction.     

Orientational potentials extracted from protein structures improve native fold recognition

We develop coarse-grained, distance- and orientation-dependent statistical potentials from the growing protein structural databases. For protein structural classes (alpha, beta, and alpha/beta), a substantial number of backbone-backbone and backbone-side-chain contacts stabilize the native folds. By taking into account the importance of backbone interactions with a virtual backbone interaction center as the 21st anisotropic site, we construct a 21 x 21 interaction scheme. The new potentials are studied using spherical harmonics analysis (SHA) and a smooth, continuous version is constructed using spherical harmonic synthesis (SHS). Our approach has the following advantages: (1) The smooth, continuous form of the resulting potentials is more realistic and presents significant advantages for computational simulations, and (2) with SHS, the potential values can be computed efficiently for arbitrary coordinates, requiring only the knowledge of a few spherical harmonic coefficients. The performance of the new orientation-dependent potentials was tested using a standard database of decoy structures. The results show that the ability of the new orientation-dependent potentials to recognize native protein folds from a set of decoy structures is strongly enhanced by the inclusion of anisotropic backbone interaction centers. The anisotropic potentials can be used to develop realistic coarse-grained simulations of proteins, with direct applications to protein design, folding, and aggregation.     

Alzheimer Aβ peptide interactions with lipid membranes: Fibrils,oligomers and polymorphic amyloid channels

Fibrillar aggregates of misfolded amyloid proteins are involved in a variety of diseases such as Alzheimer disease (AD), type 2 diabetes, Parkinson, Huntington and prion-related diseases. In the case of AD amyloid β (Aβ) peptides, the toxicity of amyloid oligomers and larger fibrillar aggregates is related to perturbing the biological function of the adjacent cellular membrane. We used atomistic molecular dynamics (MD) simulations of Aβ_9–40 fibrillar oligomers modeled as protofilament segments, including lipid bilayers and explicit water molecules, to probe the first steps in the mechanism of Aβ-membrane interactions. Our study identified the electrostatic interaction between charged peptide residues and the lipid headgroups as the principal driving force that can modulate the further penetration of the C-termini of amyloid fibrils or fibrillar oligomers into the hydrophobic region of lipid membranes. These findings advance our understanding of the detailed molecular mechanisms and the effects related to Aβ-membrane interactions, and suggest a polymorphic structural character of amyloid ion channels embedded in lipid bilayers. While inter-peptide hydrogen bonds leading to the formation of β-strands may still play a stabilizing role in amyloid channel structures, these may also present a significant helical content in peptide regions (e.g., termini) that are subject to direct interactions with lipids rather than with neighboring Aβ peptides.     

Molecular dynamics simulations of Alzheimer's beta-amyloid protofilaments

Filamentous amyloid aggregates are central to the pathology of Alzheimer's disease. We use all-atom molecular dynamics (MD) simulations with explicit solvent and multiple force fields to probe the structural stability and the conformational dynamics of several models of Alzheimer's beta-amyloid fibril structures, for both wild-type and mutated amino acid sequences. The structural models are based on recent solid state NMR data. In these models, the peptides form in-register parallel beta-sheets along the fibril axis, with dimers of two U-shaped peptides located in layers normal to the fibril axis. Four different topologies are explored for stacking the beta-strand regions against each other to form a hydrophobic core. Our MD results suggest that all four NMR-based models are structurally stable, and we find good agreement with dihedral angles estimated from solid-state NMR experiments. Asp23 and Lys28 form buried salt-bridges, resulting in an alternating arrangement of the negatively and positively charged residues along the fibril axis that is reminiscent of a one-dimensional ionic crystal. Interior water molecules are solvating the buried salt-bridges. Based on data from NMR measurements and MD simulations of short amyloid fibrils, we constructed structural models of long fibrils. Calculated X-ray fiber diffraction patterns show the characteristics of packed beta-sheets seen in experiments, and suggest new experiments that could discriminate between various fibril topologies.