Protein-protein recognition and interaction hot spots in an antigen-antibody complex: free energy decomposition identifies "efficient amino acids"

In silico degradation of fluoroquinolones by a microalgae-based constructed wetland system

Fluoroquinolone antibiotics (FQs) have been used worldwide due to their extended antimicrobial spectrum. However, the overuse of FQs leads to frequent detection in the environment and cannot be efficiently removed. Microalgae-based constructed wetland systems have been proven to be a relatively proper method to treat FQs, mainly by microalgae, plants, microorganisms, and sediments. To improve the removal efficiency of microalgae-constructed wetland, a systematic molecular design, screening, functional, and risk evaluation method was developed using three-dimensional quantitative structure-activity relationship models, molecular dynamics simulation, molecular docking, and TOPKAT approaches. Five designed ciprofloxacin alternatives with improved bactericidal effects and lower human health risks were found to be more easily degraded by microalgae (16.11-167.88 %), plants (6.72-58.86 %), microorganisms (9.10-15.02 %), and sediments (435.83 %-1763.51 %) compared with ciprofloxacin. According to the mechanism analysis, the removal effect of the FQs can be affected via changes in the number, bond energy, and molecular descriptors of favorable and unfavorable amino acids. To the best of our knowledge, this is the first comprehensive study of improving the microalgae, plants, microorganisms, and sediment removal efficiency of FQs in constructed wetlands, which provides theoretical support for the treatment of FQ pollution.

Biosynthesis of a clickable pyoverdine via in vivo enzyme engineering of an adenylation domain

The engineering of non ribosomal peptide synthetases (NRPS) for new substrate specificity is a potent strategy to incorporate non-canonical amino acids into peptide sequences, thereby creating peptide diversity and broadening applications. The non-ribosomal peptide pyoverdine is the primary siderophore produced by Pseudomonas aeruginosa and holds biomedical promise in diagnosis, bio-imaging and antibiotic vectorization. We engineered the adenylation domain of PvdD, the terminal NRPS in pyoverdine biosynthesis, to accept a functionalized amino acid. Guided by molecular modeling, we rationally designed mutants of P. aeruginosa with mutations at two positions in the active site. A single amino acid change results in the successful incorporation of an azido-L-homoalanine leading to the synthesis of a new pyoverdine analog, functionalized with an azide function. We further demonstrated that copper free click chemistry is efficient on the functionalized pyoverdine and that the conjugated siderophore retains the iron chelation properties and its capacity to be recognized and transported by P. aeruginosa. The production of clickable pyoverdine holds substantial biotechnological significance, paving the way for numerous downstream applications.

Biomimetic Design of Peptide Inhibitor to Block CD47/SIRPα Interactions

CD47 on the surface of tumor cells has become a research hot spot in immunotherapy and anticancer therapy, as it can bind to SIRPα protein on the surface of macrophages, which ultimately leads to immune escape of tumor cells. In the present study, molecular interactions between CD47 and human SIRPα proteins (including variant 1, V1 and variant 2, V2) were analyzed through molecular dynamics (MD) simulation and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. Hydrophobic interactions were found as the main driving force for the binding of CD47 on SIRPα. The residues including pyroglutamate acid (Z)1, L2, E35, Y37, E97, L101, and T102 of CD47 were identified with a significant favorable contribution to the binding of CD47 on SIRPα (both V1 and V2). Based on this, a peptide inhibitor library with the sequence ZLXRTLXEXY was designed (X represents the arbitrary residue of 20 standard amino acids) and then screened using molecular docking, MD simulations, and experimental validation. Finally, a peptide ZLIRTLHEWY was determined with high affinity with SIRPα from 8000 candidates, containing 6/10 residues favorable for the binding on SIRPα V1 and 8/10 residues favorable for the binding on SIRPα V2, which was thus considered to have potential anticancer function.

Site-directed modification of multifunctional lignocellulose-degrading enzymes of straw based on homologous modeling

Studying the straw lignocellulose strengthening mechanism during simultaneous degradation has important practical significance for improving resource utilization and reducing environmental pollution. In this paper, the degradation ability of four straw lignocellulose-degrading enzymes was evaluated by molecular docking and molecular dynamics. Using the significantly binds to straw lignocellulose-degrading enzyme as a template, a multifunctional lignocellulose-degrading enzyme 3CBH-1KS5-4XQD-1B85 was constructed based on amino acid recombination and homologous modeling. Five efficient degrading enzymes (3CBH-1, 3CBH-2, 3CBH-3, 3CBH-4, and 3CBH-5) were designed by site-directed mutagenesis of 3CBH-1KS5-4XQD-1B85 amino acid at position 346. Molecular dynamics showed that the degradation ability of 3CBH-1 was significant and it was 1.45 times higher than 3CBH-1KS5-4XQD-1B85. Moreover, the mechanism of enhanced degradability and the stability of the enzymes were explored. With the aid of Taguchi experiments, the suitable external environment for degrading straw was determined. In the presence of inhibitors (organic acids and phenolic compounds), the binding energy of 3CBH-1 (238.46 ± 30.96 kJ/mol) is 36.42% higher than that of 3CBH-1KS5-4XQD-1B85 (174.79 ± 20.35 kJ/mol) without external environmental stimulation. Based on homology modeling, this paper constructed a site-directed mutagenesis scheme of multifunctional enzymes, and the aim was to obtain multifunctional and efficient straw lignocellulose-degrading enzymes through protein engineering, which provided a feasible scheme for straw biodegradation.

Million-atom molecular dynamics simulations reveal the interfacial interactions and assembly of plant PSII-LHCII supercomplex

Protein-protein interface interactions dictate efficient excitation energy transfer from light-harvesting antennas to the photosystem II (PSII) core. In this work, we construct a 1.2 million atom-scale model of plant C₂S₂-type PSII-LHCII supercomplex and perform microsecond-scale molecular dynamics (MD) simulations to explore the interactions and assembly mechanisms of the sizeable PSII-LHCII supercomplex. We optimize the nonbonding interactions of the PSII-LHCII cryo-EM structure using microsecond-scale MD simulations. Binding free energy calculations with component decompositions reveal that hydrophobic interactions predominantly drive antenna-core association and the antenna-antenna interactions are relatively weak. Despite the positive electrostatic interaction energies, hydrogen bonds and salt bridges mainly provide directional or anchoring forces for interface binding. Analysis of the roles of small intrinsic subunits of PSII suggests that LHCII and CP26 first interact with small intrinsic subunits and then bind to the core proteins, whereas CP29 adopts a one-step binding process to the PSII core without the assistance of other factors. Our study provides insights into the molecular underpinnings of the self-organization and regulation of plant PSII-LHCII. It lays the framework for deciphering the general assembly principles of photosynthetic supercomplexes and possibly other macromolecular structures. The finding also has implications for repurposing photosynthetic systems to enhance photosynthesis.

Engineered Design of the E-Helix Structure on Ferritin Nanoparticles

Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F₁) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E₁) and flexible linker (L₃) residues (to afford F₁L₃E₁). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F₁L₃E₁ were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F₁L₃E₁. However, they are less thermally stable than F₁L₃E₁. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.

High-Risk Mucosal Human Papillomavirus 16 (HPV16) E6 Protein and Cutaneous HPV5 and HPV8 E6 Proteins Employ Distinct Strategies To Interfere with Interferon Regulatory Factor 3-Mediated Beta Interferon Expression

Persistent infection with some mucosal α-genus human papillomaviruses (HPVs; the most prevalent one being HPV16) can induce cervical carcinoma, anogenital cancers, and a subset of head and neck squamous cell carcinoma (HNSCC). Cutaneous β-genus HPVs (such as HPV5 and HPV8) associate with skin lesions that can progress into squamous cell carcinoma with sun exposure in Epidermodysplasia verruciformis patients and immunosuppressed patients. Here, we analyzed mechanisms used by E6 proteins from the α- and β-genus to inhibit the interferon-β (IFNB1) response. HPV16 E6 mediates this effect by a strong direct interaction with interferon regulatory factor 3 (IRF3). The binding site of E6 was localized within a flexible linker between the DNA-binding domain and the IRF-activation domain of IRF3 containing an LxxLL motif. The crystallographic structure of the complex between HPV16 E6 and the LxxLL motif of IRF3 was solved and compared with the structure of HPV16 E6 interacting with the LxxLL motif of the ubiquitin ligase E6AP. In contrast, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3-binding domain (IBiD) of the CREB-binding protein (CBP), a key transcriptional coactivator in IRF3-mediated IFN-β expression. IMPORTANCE Persistent HPV infections can be associated with the development of several cancers. The ability to persist depends on the ability of the virus to escape the host immune system. The type I interferon (IFN) system is the first-line antiviral defense strategy. HPVs carry early proteins that can block the activation of IFN-I. Among mucosal α-genus HPV types, the HPV16 E6 protein has a remarkable property to strongly interact with the transcription factor IRF3. Instead, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3 cofactor CBP. These results highlight the versatility of E6 proteins to interact with different cellular targets. The interaction between the HPV16 E6 protein and IRF3 might contribute to the higher prevalence of HPV16 than that of other high-risk mucosal HPV types in HPV-associated cancers.

Anticoagulant Oligonucleotide-Peptide Conjugates: Identification of Thrombin Aptamer Conjugates with Improved Characteristics

Oligonucleotide–peptide conjugates (OPCs) are a promising class of biologically active compounds with proven potential for improving nucleic acid therapeutics. OPCs are commonly recognized as an efficient instrument to enhance the cellular delivery of therapeutic nucleic acids. In addition to this application field, OPCs have an as yet unexplored potential for the post-SELEX optimization of DNA aptamers. In this paper, we report the preparation of designer thrombin aptamer OPCs with peptide side chains anchored to a particular thymidine residue of the aptamer. The current conjugation strategy utilizes unmodified short peptides and support-bound protected oligonucleotides with activated carboxyl functionality at the T3 thymine nucleobase. The respective modification of the oligonucleotide strand was implemented using N3-derivatized thymidine phosphoramidite. Aptamer OPCs retained the G-quadruplex architecture of the parent DNA structure and showed minor to moderate stabilization. In a series of five OPCs, conjugates bearing T3–Ser–Phe–Asn (SFN) or T3–Tyr–Trp–Asn (YWN) side chains exhibited considerably improved anticoagulant characteristics. Molecular dynamics studies of the aptamer OPC complexes with thrombin revealed the roles of the amino acid nature and sequence in the peptide subunit in modulating the anticoagulant activity.

NbX: Machine Learning-Guided Re-Ranking of Nanobody-Antigen Binding Poses

Modeling the binding pose of an antibody is a prerequisite to structure-based affinity maturation and design. Without knowing a reliable binding pose, the subsequent structural simulation is largely futile. In this study, we have developed a method of machine learning-guided re-ranking of antigen binding poses of nanobodies, the single-domain antibody which has drawn much interest recently in antibody drug development. We performed a large-scale self-docking experiment of nanobody–antigen complexes. By training a decision tree classifier through mapping a feature set consisting of energy, contact and interface property descriptors to a measure of their docking quality of the refined poses, significant improvement in the median ranking of native-like nanobody poses by was achieved eightfold compared with ClusPro and an established deep 3D CNN classifier of native protein–protein interaction. We further interpreted our model by identifying features that showed relatively important contributions to the prediction performance. This study demonstrated a useful method in improving our current ability in pose prediction of nanobodies.

Rational Design of Constrained Peptides as Protein Interface Inhibitors

The lack of progress in developing targeted therapeutics directed at protein–protein complexes has been due to the absence of well-defined ligand-binding pockets and the extensive intermolecular contacts at the protein–protein interface. Our laboratory has developed approaches to dissect protein–protein complexes focusing on the superfamilies of erbB and tumor necrosis factor (TNF) receptors by the combined use of structural biology and computational biology to facilitate small molecule development. We present a perspective on the development and application of peptide inhibitors as well as immunoadhesins to cell surface receptors performed in our laboratory.

A structural signature motif enlightens the origin and diversification of nuclear receptors

Nuclear receptors are ligand-activated transcription factors that modulate gene regulatory networks from embryonic development to adult physiology and thus represent major targets for clinical interventions in many diseases. Most nuclear receptors function either as homodimers or as heterodimers. The dimerization is crucial for gene regulation by nuclear receptors, by extending the repertoire of binding sites in the promoters or the enhancers of target genes via combinatorial interactions. Here, we focused our attention on an unusual structural variation of the α-helix, called π-turn that is present in helix H7 of the ligand-binding domain of RXR and HNF4. By tracing back the complex evolutionary history of the π-turn, we demonstrate that it was present ancestrally and then independently lost in several nuclear receptor lineages. Importantly, the evolutionary history of the π-turn motif is parallel to the evolutionary diversification of the nuclear receptor dimerization ability from ancestral homodimers to derived heterodimers. We then carried out structural and biophysical analyses, in particular through point mutation studies of key RXR signature residues and showed that this motif plays a critical role in the network of interactions stabilizing homodimers. We further showed that the π-turn was instrumental in allowing a flexible heterodimeric interface of RXR in order to accommodate multiple interfaces with numerous partners and critical for the emergence of high affinity receptors. Altogether, our work allows to identify a functional role for the π-turn in oligomerization of nuclear receptors and reveals how this motif is linked to the emergence of a critical biological function. We conclude that the π-turn can be viewed as a structural exaptation that has contributed to enlarging the functional repertoire of nuclear receptors.

The origin of novelties is a central topic in evolutionary biology. A fundamental question is how organisms constrained by natural selection can divert from existing schemes to set up novel structures or pathways. Among the most important strategies are exaptations, which represent pre-adaptation strategies. Many examples exist in biology, at both morphological and molecular levels, such as the one reported here that focuses on an unusual structural feature called the π-turn. It is found in the structure of the most ancestral nuclear receptors RXR and HNF4. The analyses trace back the complex evolutionary history of the π-turn to more than 500 million years ago, before the Cambrian explosion and show that this feature was essential for the heterodimerization capacity of RXR. Nuclear receptor lineages that emerged later in evolution lost the π-turn. We demonstrate here that this loss in nuclear receptors that heterodimerize with RXR was critical for the emergence of high affinity receptors, such as the vitamin D and the thyroid hormone receptors. On the other hand, the conserved π-turn in RXR allowed it to accommodate multiple heterodimer interfaces with numerous partners. This structural exaptation allowed for the remarkable diversification of nuclear receptors.

Swiss-PO: a new tool to analyze the impact of mutations on protein three-dimensional structures for precision oncology

Swiss-PO is a new web tool to map gene mutations on the 3D structure of corresponding proteins and to intuitively assess the structural implications of protein variants for precision oncology. Swiss-PO is constructed around a manually curated database of 3D structures, variant annotations, and sequence alignments, for a list of 50 genes taken from the Ion AmpliSeqTM Custom Cancer Hotspot Panel. The website was designed to guide users in the choice of the most appropriate structure to analyze regarding the mutated residue, the role of the protein domain it belongs to, or the drug that could be selected to treat the patient. The importance of the mutated residue for the structure and activity of the protein can be assessed based on the molecular interactions exchanged with neighbor residues in 3D within the same protein or between different biomacromolecules, its conservation in orthologs, or the known effect of reported mutations in its 3D or sequence-based vicinity. Swiss-PO is available free of charge or login at https://www.swiss-po.ch .

Binding mechanism of a de novo coiled coil complex elucidated from surface forces measurements

We used the Surface Forces Apparatus to elucidate the interaction mechanism between grafted 5 heptad-long peptides engineered to spontaneously form a heterodimeric coiled-coil complex. The results demonstrated that when intimate contact between peptides is reached, binding occurs first via weakly interacting but more mobile distal heptads, suggesting an induced-fit association process. Precise control of the distance between peptide-coated surfaces allowed to quantitatively monitor the evolution of their biding energy. The binding energy of the coiled-coil complex increased in a stepwise fashion rather than monotonically with the overlapping distance, each step corresponding to the interaction between a quantized number of heptads. Surface forces data were corroborated to surface plasmon resonance measurements and molecular dynamics simulations and allowed the calculation of the energetic contribution of each heptad within the coiled-coil complex.

Biomimetic Design of Peptide Neutralizer of Ebola Virus with Molecular Simulation

Ebola virus (EBOV) belongs to the Filoviridae family, which can cause severe hemorrhagic fever in humans and nonprimates. The neutralization of EBOV by monoclonal antibody (mAb) ADI-15946 was reported recently. In the present study, the molecular interactions between the receptor GPcl of EBOV and ADI-15946 were studied by molecular dynamics (MD) simulation and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis. Hydrophobic interaction was identified as the main driving force for the binding of ADI-15946 on EBOV. Moreover, the contribution of each amino acid residue for the binding was evaluated. Then, an affinity binding model (ABM) was constructed using the residues favorable for the binding, including Y107, F108, D109, W110, and R113. The biomimetic design of neutralizer against EBOV according to the ABM of ADI-15946 was then performed, followed by screening using docking, structural similarity. Two neutralizers YFDWHMR and YFDWRYR were obtained, which were proven to be capable of strong binding on GPcl and then neutralizing GPcl. These results would be helpful for the development of neutralizers for Ebola virus.

Harmicines - harmine and cinnamic acid hybrids as novel antiplasmodial hits

Harmicines constitute novel hybrid compounds that combine two agents with reported antiplasmodial properties, namely β-carboline harmine and a cinnamic acid derivative (CAD). Cu(I) catalyzed azide-alkyne cycloaddition was employed for the preparation of three classes of hybrid molecules: N-harmicines 6a-i, O-harmicines 7a-i and N,O-bis-harmicines 8a-g,i. In vitro antiplasmodial activities of harmicines against the erythrocytic stage of Plasmodium falciparum (chloroquine-sensitive Pf3D7 and chloroquine-resistant PfDd2 strains) and hepatic stage of P. berghei, as well as cytotoxicity against human liver hepatocellular carcinoma cell line (HepG2), were evaluated. Remarkably, most of the compounds exerted significant activities against both stages of the Plasmodium life cycle. The conjugation of various CADs to harmine resulted in the increased antiplasmodial activity relative to harmine. In general, O-harmicines 7 exhibited the highest activity against the erythrocytic stage of both P. falciparum strains, whereas N,O-bis harmicines 8 showed the most pronounced activity against P. berghei hepatic stages. For the latter compound, molecular dynamics simulations confirmed binding within the ATP binding site of PfHsp90, while the weaker binders, namely 6b and harmine, were found to be positioned away from this structural element. In addition, decomposition of the computed binding free energies into contributions from individual residues suggested guidelines for further derivatization of harmine towards more efficient compounds. Cytotoxicity screening revealed N-harmicines 6 as the least, and O-harmicines 7 as the most toxic compounds. Harmicines 6g, 8b and 6d exerted the most selective action towards Plasmodium over human cells, respectively. These results establish harmicines as hits for future optimisation and development of novel antiplasmodial agents.

Self-Associating Peptides for Modular Bifunctional Conjugation of Tetramer Macromolecules in Living Cells

Monitoring the assembly of macromolecules to design entities with novel properties can be achieved either chemically creating covalent bonds or by noncovalent connections using appropriate structural motifs. In this report, two self-associating peptides (named K3 and E3) that originate from p53 tetramerization domain were developed as tools for highly specific and noncovalent heterotetramerization of two biomolecules. The pairing/coupling preferences of K3 and E3 were first evaluated by molecular modeling data and confirmed using circular dichroism spectroscopy, size-exclusion chromatography, and biological assays. Regardless of the moieties fused to K3 and E3, these two peptides self-assembled into dimers of dimers to form bivalent heterotetrameric complexes that proved to be extremely stable inside living cells. The benefits of the multivalency in terms of avidity, specificity, and expanded functional activity were strikingly revealed when the proliferating cell nuclear antigen (PCNA), which is essential for DNA replication, was targeted using a heterotetramer presenting both an antibody fragment against PCNA and a specific PCNA binder peptide. In vitro heterotetramerization of these two known PCNA ligands increased their binding efficiencies to PCNA up to 80-fold compared to the best homotetramer counterpart. In cellulo, the heterotetramers were able to efficiently inhibit DNA replication and to trigger cell death. Altogether, we demonstrate that these two biselective self-assembling peptidic domains offer a versatile noncovalent conjugation method that can be easily implemented for protein engineering.

Modulation of RXR-DNA complex assembly by DNA context

Retinoid X Receptors (RXRs) act as dimer partners for several nuclear receptors including itself, binding to genomic DNA response elements and regulating gene transcription with cell and gene specificity. As homodimers, RXRs bind direct repeats of the half-site (A/G)G(G/T)TCA separated by 1 nucleotide (DR1) and little variability of this consensus site is observed for natural DR1s. However, these variations are responsible of the modulation of RXR receptors function through differential binding affinity and conformational changes. To further our understanding of the molecular mechanisms underlying RXR-DNA interactions, we examined how RXR DBDs bind to different DR1s using thermodynamics, X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy. We show that the half-site sequences modulate the binding cooperativity that results from the protein-protein contacts between the two DBDs. Chemical shifts perturbation NMR experiments revealed that sequence variations in half-sites induce changes that propagate from the protein-DNA interface to the dimerization interface throughout the DBD fold.

Alternative dimerization interfaces in the glucocorticoid receptor-α ligand binding domain

BACKGROUND

Nuclear hormone receptors (NRs) constitute a large family of multi-domain ligand-activated transcription factors. Dimerization is essential for their regulation, and both DNA binding domain (DBD) and ligand binding domain (LBD) are implicated in dimerization. Intriguingly, the glucocorticoid receptor-α (GRα) presents a DBD dimeric architecture similar to that of the homologous estrogen receptor-α (ERα), but an atypical dimeric architecture for the LBD. The physiological relevance of the proposed GRα LBD dimer is a subject of debate.

METHODS

We analyzed all GRα LBD homodimers observed in crystals using an energetic analysis based on the PISA and on the MM/PBSA methods and a sequence conservation analysis, using the ERα LBD dimer as a reference point.

RESULTS

Several dimeric assemblies were observed for GRα LBD. The assembly generally taken to be physiologically relevant showed weak binding free energy and no significant residue conservation at the contact interface, while an alternative homodimer mediated by both helix 9 and C-terminal residues showed significant binding free energy and residue conservation. However, none of the GRα LBD assemblies found in crystals are as stable or conserved as the canonical ERα LBD dimer. GRα C-terminal sequence (F-domain) forms a steric obstacle to the canonical dimer assembly in all available structures.

CONCLUSIONS

Our analysis calls for a re-examination of the currently accepted GRα homodimer structure and experimental investigations of the alternative architectures.

GENERAL SIGNIFICANCE

This work questions the validity of the currently accepted architecture. This has implications for interpreting physiological data and for therapeutic design pertaining to glucocorticoid research.

Artificial antibody created by conformational reconstruction of the complementary-determining region on gold nanoparticles

Mimicking protein-like specific interactions and functions has been a long-pursued goal in nanotechnology. The key challenge is to precisely organize nonfunctional surface groups on nanoparticles into specific 3D conformations to function in a concerted and orchestrated manner. Here, we develop a method to graft the complementary-determining regions of natural antibodies onto nanoparticles and reconstruct their “active” conformation to create nanoparticle-based artificial antibodies that recognize the corresponding antigens. Our work demonstrates that it is possible to create functions on nanoparticles by conformational engineering, namely tuning flexible surface groups into specific conformations. Our straightforward strategy could be used further to create other artificial antibodies for various applications and provides a new tool to understand the structure and folding of natural proteins.

To impart biomedical functions to nanoparticles (NPs), the common approach is to conjugate functional groups onto NPs by dint of the functions of those groups per se. It is still beyond current reach to create protein-like specific interactions and functions on NPs by conformational engineering of nonfunctional groups on NPs. Here, we develop a conformational engineering method to create an NP-based artificial antibody, denoted “Goldbody,” through conformational reconstruction of the complementary-determining regions (CDRs) of natural antibodies on gold NPs (AuNPs). The seemingly insurmountable task of controlling the conformation of the CDR loops, which are flexible and nonfunctional in the free form, was accomplished unexpectedly in a simple way. Upon anchoring both terminals of the free CDR loops on AuNPs, we managed to reconstruct the “active” conformation of the CDR loops by tuning the span between the two terminals and, as a result, the original specificity was successfully reconstructed on the AuNPs. Two Goldbodies have been created by this strategy to specifically bind with hen egg white lysozyme and epidermal growth factor receptor, with apparent affinities several orders of magnitude stronger than that of the original natural antibodies. Our work demonstrates that it is possible to create protein-like functions on NPs in a protein-like way, namely by tuning flexible surface groups to the correct conformation. Given the apparent merits, including good stability, of Goldbodies, we anticipate that a category of Goldbodies could be created to target different antigens and thus used as substitutes for natural antibodies in various applications.

Computational study on the unbinding pathways of B-RAF inhibitors and its implication for the difference of residence time: insight from random acceleration and steered molecular dynamics simulations

Random acceleration and steered molecular dynamics simulations reveal the unbinding pathway of B-RAF inhibitors and the difference in the residence time.

Disruption of protein–protein interactions: hot spot detection, structure-based virtual screening and in vitro testing for the anti-cancer drug target – survivin

Hot spot detection at the protein–protein interaction interface using computational tools helped to identify indinavir as survivin inhibitor.

Solvent accessible surface area-based hot-spot detection methods for protein-protein and protein-nucleic acid interfaces

Due to the importance of hot-spots (HS) detection and the efficiency of computational methodologies, several HS detecting approaches have been developed. The current paper presents new models to predict HS for protein-protein and protein-nucleic acid interactions with better statistics compared with the ones currently reported in literature. These models are based on solvent accessible surface area (SASA) and genetic conservation features subjected to simple Bayes networks (protein-protein systems) and a more complex multi-objective genetic algorithm-support vector machine algorithms (protein-nucleic acid systems). The best models for these interactions have been implemented in two free Web tools.

Reversible ketomethylene-based inhibitors of human neutrophil proteinase 3

Neutrophil serine proteases, proteinase 3 (PR3) and human neutrophil elastase (HNE), are considered as targets for chronic inflammatory diseases. Despite sharing high sequence similarity, the two enzymes have different substrate specificities and functions. While a plethora of HNE inhibitors exist, PR3 specific inhibitors are still in their infancy. We have designed ketomethylene-based inhibitors for PR3 that show low micromolar IC50 values. Their synthesis was made possible by amending a previously reported synthesis of ketomethylene dipeptide isosteres to allow for the preparation of derivatives suitable for solid phase peptide synthesis. The best inhibitor (Abz-VADnV[Ψ](COCH2)ADYQ-EDDnp) was found to be selective for PR3 over HNE and to display a competitive and reversible inhibition mechanism. Molecular dynamics simulations show that the interactions between enzyme and ketomethylene-containing inhibitors are similar to those with the corresponding substrates. We also confirm that N- and C-terminal FRET groups are important for securing high inhibitory potency toward PR3.

Molecular dynamics for computational proteomics of methylated histone H3

BACKGROUND

Post-translational modifications of histones, and in particular of their disordered N-terminal tails, play a major role in epigenetic regulation. The identification of proteins and proteic domains that specifically bind modified histones is therefore of paramount importance to understand the molecular mechanisms of epigenetics.

METHODS

We performed an energetic analysis using the MM/PBSA method in order to study known complexes between methylated histone H3 and effector domains of the PHD family. We then developed a simple molecular dynamics based predictive model based on our analysis.

RESULTS

We present a thorough validation of our procedure, followed by the computational predictions of new PHD domains specific for binding histone H3 methylated on lysine 4 (K4).

CONCLUSIONS

PHD domains recognize methylated K4 on histone H3 in the context of a linear interaction motif (LIM) formed by the first four amino acids of histone H3 as opposed to recognition of a single methylated site. PHD domains with different sequences find chemically equivalent solutions for stabilizing the histone LIM and these can be identified from energetic analysis. This analysis, in turn, allows for the identification of new PHD domains that bind methylated H3K4 using information that cannot be retrieved from sequence comparison alone.

GENERAL SIGNIFICANCE

Molecular dynamics simulations can be used to devise computational proteomics protocols that are both easy to implement and interpret, and that yield reliable predictions that compare favorably to and complement experimental proteomics methods. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Differential modes of peptide binding onto replicative sliding clamps from various bacterial origins

Bacterial sliding clamps are molecular hubs that interact with many proteins involved in DNA metabolism through their binding, via a conserved peptidic sequence, into a universally conserved pocket. This interacting pocket is acknowledged as a potential molecular target for the development of new antibiotics. We previously designed short peptides with an improved affinity for the Escherichia coli binding pocket. Here we show that these peptides differentially interact with other bacterial clamps, despite the fact that all pockets are structurally similar. Thermodynamic and modeling analyses of the interactions differentiate between two categories of clamps: group I clamps interact efficiently with our designed peptides and assemble the Escherichia coli and related orthologs clamps, whereas group II clamps poorly interact with the same peptides and include Bacillus subtilis and other Gram-positive clamps. These studies also suggest that the peptide binding process could occur via different mechanisms, which depend on the type of clamp.

Energetic changes caused by antigenic module insertion in a virus-like particle revealed by experiment and molecular dynamics simulations

The success of recombinant virus-like particles (VLPs) for human papillomavirus and hepatitis B demonstrates the potential of VLPs as safe and efficacious vaccines. With new modular designs emerging, the effects of antigen module insertion on the self-assembly and structural integrity of VLPs should be clarified so as to better enabling improved design. Previous work has revealed insights into the molecular energetics of a VLP subunit, capsomere, comparing energetics within various solution conditions known to drive or inhibit self-assembly. In the present study, molecular dynamics (MD) simulations coupled with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were performed to examine the molecular interactions and energetics in a modular capsomere of a murine polyomavirus (MPV) VLP designed to protect against influenza. Insertion of an influenza antigenic module is found to lower the binding energy within the capsomere, and a more active state is observed in Assembly Buffer as compared with that in Stabilization Buffer, which has been experimentally validated through measurements using differential scanning calorimetry. Further in-depth analysis based on free-energy decomposition indicates that destabilized binding can be attributed to electrostatic interaction induced by the chosen antigen module. These results provide molecular insights into the conformational stability of capsomeres and their abilities to be exploited for antigen presentation, and are expected to be beneficial for the biomolecular engineering of VLP vaccines.

Biomimetic design of affinity peptide ligand for capsomere of virus-like particle

Virus-like particle (VLP) of murine polyomavirus (MPV) is a T = 7d icosahedral capsid that self-assembles from 72 capsomeres (Caps), each of which is a pentamer of major coat protein VP1. VLP has great potential in vaccinology, gene therapy, drug delivery, and materials science. However, its application is hindered by high cost downstream processes, leading to an urgent demand of a highly efficient affinity ligand for the separation and purification of Cap by affinity chromatography. Herein a biomimetic design strategy of an affinity peptide ligand of Cap has been developed on the basis of the binding structure of the C-terminus of minor coat protein (VP2-C) on the inner surface of Cap. The molecular interactions between VP2-C and Cap were first examined using all-atom molecular dynamics (MD) simulations coupled with the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method, where V283, P285, D286, W287, L289, and Y296 of VP2-C were identified as the hot spots. An affinity peptide library (DWXLXLXY, X denotes arbitrary amino acids except cysteine) was then constructed for virtual screening sequently by docking with AUTODOCK VINA, binding structure comparison, and final docking with ROSETTA FlexPepDock. Ten peptide candidates were selected and further confirmed by MD simulations and MM/PBSA, where DWDLRLLY was found to have the highest affinity to Cap. In DWDLRLLY, six residues are favorable for the binding, including W2, L4, L6 and Y8 inheriting from VP2-C, and R5 and L7 selected in the virtual screening. This confirms the high efficiency and accuracy of the biomimetic design strategy. DWDLRLLY was then experimentally validated by a one-step purification of Cap from crude cell lysate using affinity chromatography with the octapeptide immobilized on Sepharose gel. The purified Caps were observed to self-assemble into VLP with consistent structure of authentic MPV.

Biomimetic design of platelet adhesion inhibitors to block integrin α2β1-collagen interactions: I. Construction of an affinity binding model

Platelet adhesion on a collagen surface through integrin α2β1 has been proven to be significant for the formation of arterial thrombus. However, the molecular determinants mediating the integrin-collagen complex remain unclear. In the present study, the dynamics of integrin-collagen binding and molecular interactions were investigated using molecular dynamics (MD) simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis. Hydrophobic interaction is identified as the major driving force for the formation of the integrin-collagen complex. On the basis of the MD simulation and MM-PBSA results, an affinity binding model (ABM) of integrin for collagen is constructed; it is composed of five residues, including Y157, N154, S155, R288, and L220. The ABM has been proven to capture the major binding motif contributing 84.8% of the total binding free energy. On the basis of the ABM, we expect to establish a biomimetic design strategy of platelet adhesion inhibitors, which would be beneficial for the development of potent peptide-based drugs for thrombotic diseases.

The N-terminal sequence of tyrosine hydroxylase is a conformationally versatile motif that binds 14-3-3 proteins and membranes

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the synthesis of catecholamine neurotransmitters, and a reduction in TH activity is associated with several neurological diseases. Human TH is regulated, among other mechanisms, by Ser19-phosphorylation-dependent interaction with 14-3-3 proteins. The N-terminal sequence (residues 1-43), which corresponds to an extension to the TH regulatory domain, also interacts with negatively charged membranes. By using X-ray crystallography together with molecular dynamics simulations and structural bioinformatics analysis, we have probed the conformations of the Ser19-phosphorylated N-terminal peptide [THp-(1-43)] bound to 14-3-3γ, free in solution and bound to a phospholipid bilayer, and of the unphosphorylated peptide TH-(1-43) both free and bilayer bound. As seen in the crystal structure of THp-(1-43) complexed with 14-3-3γ, the region surrounding pSer19 adopts an extended conformation in the bound state, whereas THp-(1-43) adopts a bent conformation when free in solution, with higher content of secondary structure and higher number of internal hydrogen bonds. TH-(1-43) in solution presents the highest mobility and least defined structure of all forms studied, and it shows an energetically more favorable interaction with membranes relative to THp-(1-43). Cationic residues, notably Arg15 and Arg16, which are the recognition sites of the kinases phosphorylating at Ser19, are also contributing to the interaction with the membrane. Our results reveal the structural flexibility of this region of TH, in accordance with the functional versatility and conformational adaptation to different partners. Furthermore, this structural information has potential relevance for the development of therapeutics for neurodegenerative disorders, through modulation of TH-partner interactions.

Structure-Based, Rational Design of T Cell Receptors

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1_157–165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8⁺ T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K_D = ∼1 − 5 μM. Beyond the affinity threshold at K_D < 1 μM we observed an attenuation in cellular function, in line with the “half-life” model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

Hot spots in a network of functional sites

It is of significant interest to understand how proteins interact, which holds the key phenomenon in biological functions. Using dynamic fluctuations in high frequency modes, we show that the Gaussian Network Model (GNM) predicts hot spot residues with success rates ranging between S 8–58%, C 84–95%, P 5–19% and A 81–92% on unbound structures and S 8–51%, C 97–99%, P 14–50%, A 94–97% on complex structures for sensitivity, specificity, precision and accuracy, respectively. High specificity and accuracy rates with a single property on unbound protein structures suggest that hot spots are predefined in the dynamics of unbound structures and forming the binding core of interfaces, whereas the prediction of other functional residues with similar dynamic behavior explains the lower precision values. The latter is demonstrated with the case studies; ubiquitin, hen egg-white lysozyme and M2 proton channel. The dynamic fluctuations suggest a pseudo network of residues with high frequency fluctuations, which could be plausible for the mechanism of biological interactions and allosteric regulation.

Structural basis for hijacking of cellular LxxLL motifs by papillomavirus E6 oncoproteins

E6 viral oncoproteins are key players in epithelial tumors induced by papillomaviruses in vertebrates, including cervical cancer in humans. E6 proteins target many host proteins by specifically interacting with acidic LxxLL motifs. We solved the crystal structures of bovine (BPV1) and human (HPV16) papillomavirus E6 proteins bound to LxxLL peptides from the focal adhesion protein paxillin and the ubiquitin ligase E6AP, respectively. In both E6 proteins, two zinc domains and a linker helix form a basic-hydrophobic pocket, which captures helical LxxLL motifs in a way compatible with other interaction modes. Mutational inactivation of the LxxLL binding pocket disrupts the oncogenic activities of both E6 proteins. This work reveals the structural basis of both the multifunctionality and the oncogenicity of E6 proteins.

Design, synthesis, and functional evaluation of leukocyte function associated antigen-1 antagonists in early and late stages of cancer development

The integrin leukocyte function associated antigen 1 (LFA-1) binds the intercellular adhesion molecule 1 (ICAM-1) by its α(L)-chain inserted domain (I-domain). This interaction plays a key role in cancer and other diseases. We report the structure-based design, small-scale synthesis, and biological activity evaluation of a novel family of LFA-1 antagonists. The design led to the synthesis of a family of highly substituted homochiral pyrrolidines with antiproliferative and antimetastatic activity in a murine model of colon carcinoma, as well as potent antiadhesive properties in several cancer cell lines in the low micromolar range. NMR analysis of their binding to the isolated I-domain shows that they bind to the I-domain allosteric site (IDAS), the binding site of other allosteric LFA-1 inhibitors. These results provide evidence of the potential therapeutic value of a new set of LFA-1 inhibitors, whose further development is facilitated by a synthetic strategy that is versatile and fully stereocontrolled.

Specificity and versatility of substrate binding sites in four catalytic domains of human N-terminal acetyltransferases

Nt-acetylation is among the most common protein modifications in eukaryotes. Although thought for a long time to protect proteins from degradation, the role of Nt-acetylation is still debated. It is catalyzed by enzymes called N-terminal acetyltransferases (NATs). In eukaryotes, several NATs, composed of at least one catalytic domain, target different substrates based on their N-terminal sequences. In order to better understand the substrate specificity of human NATs, we investigated in silico the enzyme-substrate interactions in four catalytic subunits of human NATs (Naa10p, Naa20p, Naa30p and Naa50p). To date hNaa50p is the only human subunit for which X-ray structures are available. We used the structure of the ternary hNaa50p/AcCoA/MLG complex and a structural model of hNaa10p as a starting point for multiple molecular dynamics simulations of hNaa50p/AcCoA/substrate (substrate = MLG, EEE, MKG), hNaa10p/AcCoA/substrate (substrate = MLG, EEE). Nine alanine point-mutants of the hNaa50p/AcCoA/MLG complex were also simulated. Homology models of hNaa20p and hNaa30p were built and compared to hNaa50p and hNaa10p. The simulations of hNaa50p/AcCoA/MLG reproduce the interactions revealed by the X-ray data. We observed strong hydrogen bonds between MLG and tyrosines 31, 138 and 139. Yet the tyrosines interacting with the substrate’s backbone suggest that their role in specificity is limited. This is confirmed by the simulations of hNaa50p/AcCoA/EEE and hNaa10p/AcCoA/MLG, where these hydrogen bonds are still observed. Moreover these tyrosines are all conserved in hNaa20p and hNaa30p. Other amino acids tune the specificity of the S1’ sites that is different for hNaa10p (acidic), hNaa20p (hydrophobic/basic), hNaa30p (basic) and hNaa50p (hydrophobic). We also observe dynamic correlation between the ligand binding site and helix that tightens under substrate binding. Finally, by comparing the four structures we propose maps of the peptide-enzyme interactions that should help rationalizing substrate-specificity and lay the ground for inhibitor design.

Exploring PHD fingers and H3K4me0 interactions with molecular dynamics simulations and binding free energy calculations: AIRE-PHD1, a comparative study

PHD fingers represent one of the largest families of epigenetic readers capable of decoding post-translationally modified or unmodified histone H3 tails. Because of their direct involvement in human pathologies they are increasingly considered as a potential therapeutic target. Several PHD/histone-peptide structures have been determined, however relatively little information is available on their dynamics. Studies aiming to characterize the dynamic and energetic determinants driving histone peptide recognition by epigenetic readers would strongly benefit from computational studies. Herein we focus on the dynamic and energetic characterization of the PHD finger subclass specialized in the recognition of histone H3 peptides unmodified in position K4 (H3K4me0). As a case study we focused on the first PHD finger of autoimmune regulator protein (AIRE-PHD1) in complex with H3K4me0. PCA analysis of the covariance matrix of free AIRE-PHD1 highlights the presence of a "flapping" movement, which is blocked in an open conformation upon binding to H3K4me0. Moreover, binding free energy calculations obtained through Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodology are in good qualitative agreement with experiments and allow dissection of the energetic terms associated with native and alanine mutants of AIRE-PHD1/H3K4me0 complexes. MM/PBSA calculations have also been applied to the energetic analysis of other PHD fingers recognizing H3K4me0. In this case we observe excellent correlation between computed and experimental binding free energies. Overall calculations show that H3K4me0 recognition by PHD fingers relies on compensation of the electrostatic and polar solvation energy terms and is stabilized by non-polar interactions.

Comparing aminoglycoside binding sites in bacterial ribosomal RNA and aminoglycoside modifying enzymes

Aminoglycoside antibiotics are used against severe bacterial infections. They bind to the bacterial ribosomal RNA and interfere with the translation process. However, bacteria produce aminoglycoside modifying enzymes (AME) to resist aminoglycoside actions. AMEs form a variable group and yet they specifically recognize and efficiently bind aminoglycosides, which are also diverse in terms of total net charge and the number of pseudo-sugar rings. Here, we present the results of 25 molecular dynamics simulations of three AME representatives and aminoglycoside ribosomal RNA binding site, unliganded and complexed with an aminoglycoside, kanamycin A. A comparison of the aminoglycoside binding sites in these different receptors revealed that the enzymes efficiently mimic the nucleic acid environment of the ribosomal RNA binding cleft. Although internal dynamics of AMEs and their interaction patterns with aminoglycosides differ, the energetical analysis showed that the most favorable sites are virtually the same in the enzymes and RNA. The most copied interactions were of electrostatic nature, but stacking was also replicated in one AME:kanamycin complex. In addition, we found that some water-mediated interactions were very stable in the simulations of the complexes. We show that our simulations reproduce well findings from NMR or X-ray structural studies, as well as results from directed mutagenesis. The outcomes of our analyses provide new insight into aminoglycoside resistance mechanism that is related to the enzymatic modification of these drugs.

Exploring the inter-molecular interactions in amyloid-β protofibril with molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area free energy calculations

Aggregation of amyloid-β (Aβ) peptides correlates with the pathology of Alzheimer's disease. However, the inter-molecular interactions between Aβ protofibril remain elusive. Herein, molecular mechanics Poisson-Boltzmann surface area analysis based on all-atom molecular dynamics simulations was performed to study the inter-molecular interactions in Aβ(17-42) protofibril. It is found that the nonpolar interactions are the important forces to stabilize the Aβ(17-42) protofibril, while electrostatic interactions play a minor role. Through free energy decomposition, 18 residues of the Aβ(17-42) are identified to provide interaction energy lower than -2.5 kcal/mol. The nonpolar interactions are mainly provided by the main chain of the peptide and the side chains of nine hydrophobic residues (Leu17, Phe19, Phe20, Leu32, Leu34, Met35, Val36, Val40, and Ile41). However, the electrostatic interactions are mainly supplied by the main chains of six hydrophobic residues (Phe19, Phe20, Val24, Met35, Val36, and Val40) and the side chains of the charged residues (Glu22, Asp23, and Lys28). In the electrostatic interactions, the overwhelming majority of hydrogen bonds involve the main chains of Aβ as well as the guanidinium group of the charged side chain of Lys28. The work has thus elucidated the molecular mechanism of the inter-molecular interactions between Aβ monomers in Aβ(17-42) protofibril, and the findings are considered critical for exploring effective agents for the inhibition of Aβ aggregation.

Understanding the molecular mechanism of the broad and potent neutralization of HIV-1 by antibody VRC01 from the perspective of molecular dynamics simulation and binding free energy calculations

VRC01 is one of the most broadly and potently neutralizing HIV-1 antibodies known-it has been shown to neutralize 91 % of the tested primary isolate Env pseudoviruses by recognizing the viral envelope glycoprotein gp120. To explore the mechanism of HIV-1 neutralization by VRC01 and thus obtain valuable information for vaccine design, we performed molecular dynamics simulations and binding free energy calculations for apo-VRC01, apo-gp120, and the gp120-VRC01 complex. For gp120, residue energy decomposition analysis showed that the hotspot residues Asn280, Lys282, Asp368, Ile371, and Asp457 are located in three primary loops, including the CD4-binding loop, loop D, and loop V5. For VRC01, the hotspot residues Trp47, Trp50, Asn58, Arg61, Gln64, Trp100, and Tyr91 mainly come from CDR2 of the heavy chain. By decomposing the binding free energy into different components, intermolecular van der Waals interactions and nonpolar solvation were found to dominate the binding process. Principal component analysis of loops D and V5, which are related to neutralization resistance, indicated that these two areas have a larger conformational space in apo-gp120 compared to bound gp120. A comparison of three representative structures from the cluster analysis of loops D and V5 indicated that changes primarily occur at the tip of loop V5, and are caused by fluctuations in the terminal Glu1 residue of the antibody. This information can be used to guide the design of vaccines and small molecule inhibitors.