Shape Complementarity Research Articles

A conserved epitope III on hepatitis C virus E2 protein has alternate conformations facilitating cell binding or virus neutralization

Epitope III, a highly conserved amino acid motif of 524APTYSW529 on the hepatitis C virus (HCV) E2 glycoprotein, resides in the critical loop that binds to the host receptor CD81, thus making it one of the most important antibody targets for blocking HCV infections. Here, we have determined the X-ray crystal structure of epitope III at a 2.0-Å resolution when it was captured by a site-specific neutralizing antibody, monoclonal antibody 1H8 (mAb1H8). The snapshot of this complex revealed that epitope III has a relatively rigid structure when confined in the binding grooves of mAb1H8, which confers the residue specificity at both ends of the epitope. Such a high shape complementarity is reminiscent of the "lock and key" mode of action, which is reinforced by the incompatibility of an antibody binding with an epitope bearing specific mutations. By subtly positioning the side chains on the three residues of Tyr527, Ser528, and Trp529 while preserving the spatial rigidity of the rest, epitope III in this cocrystal complex adopts a unique conformation that is different from previously described E2 structures. With further analyses of molecular docking and phage display-based peptide interactions, we recognized that it is the arrangements of two separate sets of residues within epitope III that create these discrete conformations for the epitope to interact selectively with either mAb1H8 or CD81. These observations thus raise the possibility that local epitope III conformational dynamics, in conjunction with sequence variations, may act as a regulatory mechanism to coordinate "mAb1H8-like" antibody-mediated immune defenses with CD81-initiated HCV infections.

Molecular recognition of structurally disordered Pro/Ala-rich sequences (PAS) by antibodies involves an Ala residue at the hot spot of the epitope

Pro/Ala-rich sequences (PAS) are polypeptides that were developed as a biological alternative to poly-ethylene glycol (PEG) to generate biopharmaceuticals with extended plasma half-life. Like PEG, PAS polypeptides are conformationally disordered and show high solubility in water. Devoid of any charged or prominent hydrophobic side chains, these biosynthetic polymers represent an extreme case of intrinsically disordered proteins. Despite lack of immunogenicity of PAS tags in numerous animal studies we now succeeded in generating monoclonal antibodies (MAbs) against three different PAS versions. To this end, mice were immunized with a PAS#1, P/A#1 or APSA 40mer peptide conjugated to keyhole limpet hemocyanin as highly immunogenic carrier protein. In each case, one MAb with high binding activity and specificity towards a particular PAS motif was obtained. The apparent affinity was strongly dependent on the avidity effect and most pronounced for the bivalent MAb when interacting with a long PAS repeat. X-ray structural analysis of four representative anti-PAS Fab fragments in complex with their cognate PAS epitope peptides revealed interactions dominated by hydrogen bond networks involving the peptide backbone as well as multiple Van der Waals contacts arising from intimate shape complementarity. Surprisingly, Ala, the L-amino acid with the smallest side chain, emerged as a crucial feature for epitope recognition, contributing specific contacts at the center of the paratope in several anti-PAS complexes. Apart from these insights into how antibodies can recognize feature-less peptides without secondary structure, the MAbs characterized in this study offer valuable reagents for the preclinical and clinical development of PASylated biologics.

In-Silico Evidence for a Two Receptor Based Strategy of SARS-CoV-2.

We propose a computational investigation on the interaction mechanisms between SARS-CoV-2 spike protein and possible human cell receptors. In particular, we make use of our newly developed numerical method able to determine efficiently and effectively the relationship of complementarity between portions of protein surfaces. This innovative and general procedure, based on the representation of the molecular isoelectronic density surface in terms of 2D Zernike polynomials, allows the rapid and quantitative assessment of the geometrical shape complementarity between interacting proteins, which was unfeasible with previous methods. Our results indicate that SARS-CoV-2 uses a dual strategy: in addition to the known interaction with angiotensin-converting enzyme 2, the viral spike protein can also interact with sialic-acid receptors of the cells in the upper airways.

Deracemization of Carbohelicenes by a Chiral Perylene Bisimide Cyclophane Template Catalyst.

Deracemization describes the conversion of a racemic mixture of a chiral molecule into an enantioenriched mixture or an enantiopure compound without structural modifications. Herein, we report an inherently chiral perylene bisimide (PBI) cyclophane whose chiral pocket is capable of transforming a racemic mixture of [5]‐helicene into an enantioenriched mixture with an enantiomeric excess of 66 %. UV/Vis and fluorescence titration studies reveal this cyclophane host composed of two helically twisted PBI dyes has high binding affinities for the respective homochiral carbohelicene guests, with outstanding binding constants of up to 3.9×1010 m −1 for [4]‐helicene. 2D NMR studies and single‐crystal X‐ray analysis demonstrate that the observed strong and enantioselective binding of homochiral carbohelicenes and the successful template‐catalyzed deracemization of [5]‐helicene can be explained by the enzyme‐like perfect shape complementarity of the macrocyclic supramolecular host.

Deracemization of Carbohelicenes by a Chiral Perylene Bisimide Cyclophane Template Catalyst

AbstractDeracemization describes the conversion of a racemic mixture of a chiral molecule into an enantioenriched mixture or an enantiopure compound without structural modifications. Herein, we report an inherently chiral perylene bisimide (PBI) cyclophane whose chiral pocket is capable of transforming a racemic mixture of [5]‐helicene into an enantioenriched mixture with an enantiomeric excess of 66 %. UV/Vis and fluorescence titration studies reveal this cyclophane host composed of two helically twisted PBI dyes has high binding affinities for the respective homochiral carbohelicene guests, with outstanding binding constants of up to 3.9×1010 m−1 for [4]‐helicene. 2D NMR studies and single‐crystal X‐ray analysis demonstrate that the observed strong and enantioselective binding of homochiral carbohelicenes and the successful template‐catalyzed deracemization of [5]‐helicene can be explained by the enzyme‐like perfect shape complementarity of the macrocyclic supramolecular host.

In-silico Screening of Synthetic Inhibitors for Human Poly (Adp-ribose) Polymerase 2 Enzyme Using Patch Dock Software for Ovarian Cancer Therapy

Aim: This research deals with the identification of drug synthetic inhibitors docked with PARP enzyme by In silico screening using patch Dock for targeted ovarian cancer. Materials and methods: The sample size (N= 51) preparation is of two-dimensional structures for 50 synthetic compounds and a reference compound (Rucaparib) that were retrieved and collected from the NCBI-PubChem compound database. The ligand molecule was prepared using the LigPrep wizard of the Schrödinger suite. Patch Dock software was used to conduct a molecular docking study of PARP with plant-derived compounds. This program uses an algorithm to generate output complexes based on biomolecule shape complementarity. The best poses were analyzed for non-covalent interaction using PLIP and discovery studio visualizer. Results: Molecular docking analysis revealed (Lyngabyabelline, Sansalvamide, Dolastatin, Beta-carotene, Butulinic) from Plant derived could bind PARP protein with higher affinity, ability to competitively inhibit PARP as a stable drug candidate in comparison with Rucaparib and also show similar residue interaction patterns when compared with other PARP inhibitors that would report an affirmative prognostic factor. Conclusion: The identified inhibitor complexes from plant derived synthetic molecules are expected to bind with PARP protein with better efficiency in comparison with Rucaparib, hence they can be further considered for in vivo and in vitro analysis.

Analysis of Binding Modes of Antigen-Antibody Complexes by Molecular Mechanics Calculation.

Antibodies are one of the most important protein molecules in biopharmaceutics. Due to the recent advance in technology for producing monoclonal antibodies, many structural data are available on the antigen-antibody complexes. To characterize the molecular interaction in antigen-antibody recognition, we computationally analyzed 500 complex structures by molecular mechanics calculations. The presence of Ser and Tyr is markedly large in the complementarity-determining regions (CDRs). Although Ser is abundant in CDRs, its contribution to the binding score is not large. Instead, Tyr, Asp, Glu, and Arg significantly contribute to the molecular interaction from the viewpoint of the binding score. The decomposition of the binding score suggests that the hydrophilic interaction is predominant in all CDRs compared with the hydrophobic one. The contribution of the heavy chain is larger than that of the light chain. In particular, H2 and H3 largely contribute to the binding interaction. Tyr is a main contributing residue both in H2 and H3. The positively charged residue Arg also significantly contributes to the binding score in H3, while the contribution of Lys is small. The appearance of Ser is remarkable in H2, and Asp is abundant in H3. The non-charged polar residues, Thr, Asn, and Gln, appear much in H2, compared to appearing in H3. The negatively charged residues Asp and Glu significantly contribute to the binding score in H3. The contributions of Phe and Trp are not large in spite that the aromatic residues are capable of making the π-π or CH-π interaction. Gly is commonly abundant both in H2 and H3. The average distance of the shortest direct hydrogen bond between the antigen and antibody is longer than that of the hydrogen bonds observed in the complexes between compounds and their target proteins. Therefore, the antigen-antibody interface is not so tight as the compound-target protein interface. The calculation of shape complementarity is consistent with the result of the hydrogen bonds in that the fitness of the antigen-antibody contact is not so high as that of the compound-target protein contact. There exist many water molecules at the antigen-antibody interface. These findings suggest that Tyr, Asp, Glu, and Arg are rich in H3 and work as major contributors for the interaction with the antigen. Ser, Thr, Asn, and Gln are rich in H2 and support the interaction with enhancing molecular fitness. Gly is helpful in increasing flexibility and geometrical diversity. Because the antigen-antibody binding is fundamentally hydrophilic-driven, the non-polar residues are unfavorable for mediating the contact even for the aromatic residues such as Phe and Trp.

Simple biochemical features underlie transcriptional activation domain diversity and dynamic, fuzzy binding to Mediator

Gene activator proteins comprise distinct DNA-binding and transcriptional activation domains (ADs). Because few ADs have been described, we tested domains tiling all yeast transcription factors for activation in vivo and identified 150 ADs. By mRNA display, we showed that 73% of ADs bound the Med15 subunit of Mediator, and that binding strength was correlated with activation. AD-Mediator interaction in vitro was unaffected by a large excess of free activator protein, pointing to a dynamic mechanism of interaction. Structural modeling showed that ADs interact with Med15 without shape complementarity (“fuzzy” binding). ADs shared no sequence motifs, but mutagenesis revealed biochemical and structural constraints. Finally, a neural network trained on AD sequences accurately predicted ADs in human proteins and in other yeast proteins, including chromosomal proteins and chromatin remodeling complexes. These findings solve the longstanding enigma of AD structure and function and provide a rationale for their role in biology.

Simple biochemical features underlie transcriptional activation domain diversity and dynamic, fuzzy binding to Mediator.

Gene activator proteins comprise distinct DNA-binding and transcriptional activation domains (ADs). Because few ADs have been described, we tested domains tiling all yeast transcription factors for activation in vivo and identified 150 ADs. By mRNA display, we showed that 73% of ADs bound the Med15 subunit of Mediator, and that binding strength was correlated with activation. AD-Mediator interaction in vitro was unaffected by a large excess of free activator protein, pointing to a dynamic mechanism of interaction. Structural modeling showed that ADs interact with Med15 without shape complementarity ('fuzzy' binding). ADs shared no sequence motifs, but mutagenesis revealed biochemical and structural constraints. Finally, a neural network trained on AD sequences accurately predicted ADs in human proteins and in other yeast proteins, including chromosomal proteins and chromatin remodeling complexes. These findings solve the longstanding enigma of AD structure and function and provide a rationale for their role in biology.

Insight into Inhibitory Mechanism of PDE4D by Dietary Polyphenols Using Molecular Dynamics Simulations and Free Energy Calculations.

Phosphodiesterase 4 (PDE4), mainly present in immune, epithelial, and brain cells, represents a family of key enzymes for the degradation of cyclic adenosine monophosphate (cAMP), which modulates inflammatory response. In recent years, the inhibition of PDE4 has been proven to be an effective therapeutic strategy for the treatment of neurological disorders. PDE4D constitutes a high-interest therapeutic target primarily for the treatment of Alzheimer’s disease, as it is highly involved in neuroinflammation, learning ability, and memory dysfunctions. In the present study, a thorough computational investigation consisting of molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations based on the linear response approximation (LRA) method was performed to study dietary polyphenols as potential PDE4D inhibitors. The obtained results revealed that curcumin, 6-gingerol, capsaicin, and resveratrol represent potential PDE4D inhibitors; however, the predicted binding free energies of 6-gingerol, capsaicin, and resveratrol were less negative than in the case of curcumin, which exhibited the highest inhibitory potency in comparison with a positive control rolipram. Our results also revealed that the electrostatic component through hydrogen bonding represents the main driving force for the binding and inhibitory activity of curcumin, 6-gingerol, and resveratrol, while the van der Waals component through shape complementarity plays the most important role in capsaicin’s inhibitory activity. All investigated compounds form hydrophobic interactions with residues Gln376 and Asn602 as well as hydrogen bonds with nearby residues Asp438, Met439, and Ser440. The binding mode of the studied natural compounds is consequently very similar; however, it significantly differs from the binding of known PDE4 inhibitors. The uncovered molecular inhibitory mechanisms of four investigated natural polyphenols, curcumin, 6-gingerol, capsaicin, and resveratrol, form the basis for the design of novel PDE4D inhibitors for the treatment of Alzheimer’s disease with a potentially wider therapeutic window and fewer adverse side effects.

Integrative proteomics identifies thousands of distinct, multi-epitope, and high-affinity nanobodies.

The antibody immune response is essential for the survival of mammals. However, we still lack a systematic understanding of the antibody repertoire. Here, we developed a proteomic strategy to survey, at an unprecedented scale, the landscape of antigen-engaged, circulating camelid heavy-chain antibodies, whose minimal binding fragments are called VHH antibodies or nanobodies. The sensitivity and robustness of this approach were validated with three antigens spanning orders of magnitude in immune responses; thousands of distinct, high-affinity nanobody families were reliably identified and quantified. Using high-throughput structural modeling, cross-linking mass spectrometry, mutagenesis, and deep learning, we mapped and analyzed the epitopes of >100,000 antigen-nanobody complexes. Our results revealed a surprising diversity of ultrahigh-affinity camelid nanobodies for specific antigen binding on various dominant epitope clusters. Nanobodies utilize both shape and charge complementarity to enable highly selective antigen binding. Interestingly, we found that nanobody-antigen binding can mimic conserved intracellular protein-protein interactions. A record of this paper's Transparent Peer Review process is included in the Supplemental information.

Molecular Mechanisms Behind Anti SARS-CoV-2 Action of Lactoferrin.

Despite the huge effort to contain the infection, the novel SARS-CoV-2 coronavirus has rapidly become pandemic, mainly due to its extremely high human-to-human transmission capability, and a surprisingly high viral charge of symptom-less people. While the seek for a vaccine is still ongoing, promising results have been obtained with antiviral compounds. In particular, lactoferrin is regarded to have beneficial effects both in preventing and soothing the infection. Here, we explore the possible molecular mechanisms with which lactoferrin interferes with SARS-CoV-2 cell invasion, preventing attachment and/or entry of the virus. To this aim, we search for possible interactions lactoferrin may have with virus structural proteins and host receptors. Representing the molecular iso-electron surface of proteins in terms of 2D-Zernike descriptors, we 1) identified putative regions on the lactoferrin surface able to bind sialic acid present on the host cell membrane, sheltering the cell from the virus attachment; 2) showed that no significant shape complementarity is present between lactoferrin and the ACE2 receptor, while 3) two high complementarity regions are found on the N- and C-terminal domains of the SARS-CoV-2 spike protein, hinting at a possible competition between lactoferrin and ACE2 for the binding to the spike protein.

Structural Insights into JAK2 Inhibition by Ruxolitinib, Fedratinib, and Derivatives Thereof.

The discovery that aberrant activity of Janus kinase 2 (JAK2) is a driver of myeloproliferative neoplasms (MPNs) has led to significant efforts to develop small molecule inhibitors for this patient population. Ruxolitinib and fedratinib have been approved for use in MPN patients, while baricitinib, an achiral analogue of ruxolitinib, has been approved for rheumatoid arthritis. However, structural information on the interaction of these therapeutics with JAK2 remains unknown. Here, we describe a new methodology for the large-scale production of JAK2 from mammalian cells, which enabled us to determine the first crystal structures of JAK2 bound to these drugs and derivatives thereof. Along with biochemical and cellular data, the results provide a comprehensive view of the shape complementarity required for chiral and achiral inhibitors to achieve highest activity, which may facilitate the development of more effective JAK2 inhibitors as therapeutics.

π-π Catalysis in Carbon Flatland-Flipping [8]Annulene on Graphene.

Noncovalent interactions are an integral part of the modern catalysis toolbox. Although stronger noncovalent interactions such as hydrogen bonding are commonly the main driving force of catalysis, π-π interactions typically provide smaller additional stabilizations, for example, to afford selectivity enhancements. Here, it is shown computationally that pristine graphene flakes may efficiently catalyze the skeletal inversions of various benzannulated cyclooctatetraene derivatives, providing an example of a catalytic process driven solely by π-π stacking interactions. Hereby, the catalytic effect results from disproportionate shape complementarity between catalyst and transition structure compared with catalyst and reactant. An energy decomposition analysis reveals electrostatic and, especially with increasing system size, to a larger extent, dispersion interactions as the origin of stabilization.

Shape and Size Complementarity-Induced Formation of Supramolecular Protein Assemblies with Metal-Oxo Clusters

The controlled formation of protein supramolecular assemblies is challenging, but it could provide an important route for the development of hybrid biomaterials. In this work, we demonstrate the formation of well-defined complexes between the eightfold symmetrical designer protein Tako8 and soluble metal-oxo clusters from the family of Anderson–Evans, Keggin, and ZrIV-substituted Wells–Dawson polyoxometalates. A combination of X-ray crystallography and solution studies showed that metal-oxo clusters are able to serve as linker nodes for the bottom-up creation of protein-based supramolecular assemblies. Our findings indicate that clusters with larger size and negative charge are capable of modulating the crystal packing of the protein, highlighting the need for a size and shape complementarity with the protein node for optimal alteration of the crystalline self-assembly.

Crystallization and self-assembly of shape-complementary sequence-defined peptoids

Shape complementarity between polymers is a hallmark of biological systems (e.g. DNA base pairing and protein binding interactions). Here we explore the role of shape complementarity between sequence-defined N-alkyl peptoids in crystal lattices.

Capturing Surface Complementarity in Proteins Using Unsupervised Learning and Robust Curvature Measure

Structure of a protein plays a pivotal role in determining its function. Often, protein surface's shape and curvature dictate its nature of interaction with other proteins and biomolecules. However, marked by corrugations and roughness, a protein's surface poses significant challenges for its curvature-based characterization. In the present study, we employ unsupervised machine learning to segment the protein surface into patches. To measure surface curvature of a patch, we present an algebraic sphere fitting method that is fast, accurate, and robust. Moreover, we use local curvatures to show the existence of “shape complementarity” in protein-protein, antigen-antibody, and protein-ligand interfaces. We believe that the present approach could help understand the relationship between protein structure and its biological function and help finding binding partners of a given protein.

The shape selectivity of corannulene dimers based on concave-convex and convex-convex shape complementarity as hosts for C60 and C70.

In the formation of noncovalent complexes, the stacking arrangements of corannulene and fullerene are diverse, most of which are combinations of multiple corannulenes and fullerene. Here, a composition ratio of 2 : 1 was selected for the complex between corannulene and fullerene (C60 and C70) to investigate the effects of different superposition modes, including concave-convex and convex-convex interactions, on the stability and third-order nonlinear optical (NLO) properties of the composite materials. It was found that the concave-convex interaction was stronger and it was reported to stabilize the charge-transfer (CT) complex more effectively than the convex-convex interaction. The dispersion range of the concave-convex interaction was larger than that of the convex-convex interaction, which is consistent with the interaction energy results. The packing design with the double convex-convex interactions exhibited the largest linear optical response and third-order NLO response, which showed that the convex-convex interaction was more likely to be excited and cause intermolecular CT as compared to the concave-convex interaction. This work confirmed that the packing arrangement significantly affected the NLO response and will advance the development of NLO crystals.

Computational optimization of angiotensin-converting enzyme 2 for SARS-CoV-2 Spike molecular recognition

Since the beginning of the Covid19 pandemic, many efforts have been devoted to identifying approaches to neutralize SARS-CoV-2 replication within the host cell. A promising strategy to block the infection consists of using a mutant of the human receptor angiotensin-converting enzyme 2 (ACE2) as a decoy to compete with endogenous ACE2 for the binding to the SARS-CoV-2 Spike protein, which decreases the ability of the virus to enter the host cell. Here, using a computational framework based on the 2D Zernike formalism we investigate details of the molecular binding and evaluate the changes in ACE2-Spike binding compatibility upon mutations occurring in the ACE2 side of the molecular interface. We demonstrate the efficacy of our method by comparing our results with experimental binding affinities changes upon ACE2 mutations, separating ones that increase or decrease binding affinity with an Area Under the ROC curve ranging from 0.66 to 0.93, depending on the magnitude of the effects analyzed. Importantly, the iteration of our approach leads to the identification of a set of ACE2 mutants characterized by an increased shape complementarity with Spike. We investigated the physico-chemical properties of these ACE2 mutants and propose them as bona fide candidates for Spike recognition.

Structure–activity relationship of ipglycermide binding to phosphoglycerate mutases

Catalysis of human phosphoglycerate mutase is dependent on a 2,3-bisphosphoglycerate cofactor (dPGM), whereas the nonhomologous isozyme in many parasitic species is cofactor independent (iPGM). This mechanistic and phylogenetic diversity offers an opportunity for selective pharmacologic targeting of glycolysis in disease-causing organisms. We previously discovered ipglycermide, a potent inhibitor of iPGM, from a large combinatorial cyclic peptide library. To fully delineate the ipglycermide pharmacophore, herein we construct a detailed structure–activity relationship using 280 substituted ipglycermide analogs. Binding affinities of these analogs to immobilized Caenorhabditis elegans iPGM, measured as fold enrichment relative to the index residue by deep sequencing of an mRNA display library, illuminated the significance of each amino acid to the pharmacophore. Using cocrystal structures and binding kinetics, we show that the high affinity of ipglycermide for iPGM orthologs, from Brugia malayi, Onchocerca volvulus, Dirofilaria immitis, and Escherichia coli, is achieved by a codependence between (1) the off-rate mediated by the macrocycle Cys14 thiolate coordination to an active-site Zn2+ in the iPGM phosphatase domain and (2) shape complementarity surrounding the macrocyclic core at the phosphotransferase–phosphatase domain interface. Our results show that the high-affinity binding of ipglycermide to iPGMs freezes these structurally dynamic enzymes into an inactive, stable complex.