NMR study of non-structural proteins–part III: <Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N backbone and side-chain resonance assignment of macro domain from <Emphasis Type="Italic">Chikungunya</Emphasis> virus (CHIKV)

Macro domains are conserved protein domains found in eukaryotic organisms, bacteria, and archaea as well as in certain viruses. They consist of 130–190 amino acids and can bind ADP-ribose. Although the exact role of these domains is not fully understood, the conserved binding affinity for ADP-ribose indicates that this ligand is important for the function of the domain. Such a macro domain is also present in the non-structural protein 3 (nsP3) of Chikungunya Alphavirus (CHIKV) and consists of 160 amino acids. In this study we describe the high yield expression of the macro domain from CHIKV and its preliminary structural analysis via solution NMR spectroscopy. The macro domain seems to be folded in solution and an almost complete backbone assignment was achieved. In addition, the α/β/α sandwich topology with 4 α-helices and 6 β-strands was predicted by TALOS+.     

Minimal models of multi-site ligand-binding kinetics

Mathematical models based on the current understanding of co-operativity in ligand binding to the (macro) molecule and relating the dose-response (saturation) curve of the (macro) molecule ligation to intrinsic dissociation constants characterizing the affinities of ligand for binding sites of both unliganded and partly liganded (macro) molecule have been developed. The simplified models disregarding the structural properties and considerations concerning conformational changes of the (macro) molecule retain the ability to yield sigmoid curves of ligand binding and reflect the co-operativity. Model 1 contains only three parameters, parameter κ (a multiplier characterising the change in the affinity) reflects also the existence and type of co-operativity of ligand binding: κ<1 corresponds to positive co-operativity, κ>1 to the negative and κ=1 to the absence of any co-operativity. Model 2 contains an extra parameter, ω, equilibrium constant for the T₀↔R₀ transition but fails to produce dose-response, which would suggest negative co-operativity. For any fixed n>1, the deviation of the dose-response (saturation) curve from the Henri hyperbola depends either solely on parameter κ (Model 1) or also on parameter ω (Model 2). The (macro) molecule being a receptor, both models yield a diversity of dose-response curves due to possible variety of efficacies of the (macro) molecule. The models may be considered as extensions of the Henri model: in case the dissociation constants remain unchanged, the proposed models are reduced to the latter.     

Molecular dynamics simulations reveal structural insights into inhibitor binding modes and mechanism of casein kinase II inhibitors

Casein kinase-II, a member of protein kinase family, plays significant role in different cellular processes such as cell growth, differentiation, proliferation, gene expression, and embryogenesis. Being a potent suppressor of apoptosis, it serves as a significant link for its association with various types of malignancies such as colorectal and breast cancer. To overcome its pathological role in various cancerous diseases, CK-II procures great consideration as a therapeutic target. This study aimed at understanding the binding mechanism and structural properties of benzimidazole derivatives by utilizing various computational tools including docking simulation, three-dimensional quantitative structure activity relationships and molecular dynamic simulation. Structure-based 3D-QSAR techniques such as CoMFA and CoMSIA models, were established from the conformations gained by protein–ligand docking approach. The attained models have showed a good extrapolative power for internal as well as external validation. Moreover, MD simulation was carried out to explain the detailed binding mechanism and the comparison of inhibitor’s binding mode with diverse biological activities. A good correlation was observed among docking studies, MD results, and contour map analysis. Interestingly new molecules were designed using detail structural information from MD simulation, showed higher potency of inhibition (pIC₅₀ 7.6–7.7) compare to the most active compound of the series.     

An ecological niche model to predict the geographic distribution of Haemagogus janthinomys,Dyar, 1921 a yellow fever and Mayaro virus vector,in South America

Yellow fever virus (YFV) has a long history of impacting human health in South America. Mayaro virus (MAYV) is an emerging arbovirus of public health concern in the Neotropics and its full impact is yet unknown. Both YFV and MAYV are primarily maintained via a sylvatic transmission cycle but can be opportunistically transmitted to humans by the bites of infected forest dwelling Haemagogus janthinomys Dyar, 1921. To better understand the potential risk of YFV and MAYV transmission to humans, a more detailed understanding of this vector species’ distribution is critical. This study compiled a comprehensive database of 177 unique Hg. janthinomys collection sites retrieved from the published literature, digitized museum specimens and publicly accessible mosquito surveillance data. Covariate analysis was performed to optimize a selection of environmental (topographic and bioclimatic) variables associated with predicting habitat suitability, and species distributions modelled across South America using a maximum entropy (MaxEnt) approach. Our results indicate that suitable habitat for Hg. janthinomys can be found across forested regions of South America including the Atlantic forests and interior Amazon.     

A systematic review and meta-analysis of the potential non-human animal reservoirs and arthropod vectors of the Mayaro virus

Improving our understanding of Mayaro virus (MAYV) ecology is critical to guide surveillance and risk assessment. We conducted a PRISMA-adherent systematic review of the published and grey literature to identify potential arthropod vectors and non-human animal reservoirs of MAYV. We searched PubMed/MEDLINE, Embase, Web of Science, SciELO and grey-literature sources including PAHO databases and dissertation repositories. Studies were included if they assessed MAYV virological/immunological measured occurrence in field-caught, domestic, or sentinel animals or in field-caught arthropods. We conducted an animal seroprevalence meta-analysis using a random effects model. We compiled granular georeferenced maps of non-human MAYV occurrence and graded the quality of the studies using a customized framework. Overall, 57 studies were eligible out of 1523 screened, published between the years 1961 and 2020. Seventeen studies reported MAYV positivity in wild mammals, birds, or reptiles and five studies reported MAYV positivity in domestic animals. MAYV positivity was reported in 12 orders of wild-caught vertebrates, most frequently in the orders Charadriiformes and Primate. Sixteen studies detected MAYV in wild-caught mosquito genera including Haemagogus, Aedes, Culex, Psorophora, Coquillettidia, and Sabethes. Vertebrate animals or arthropods with MAYV were detected in Brazil, Panama, Peru, French Guiana, Colombia, Trinidad, Venezuela, Argentina, and Paraguay. Among non-human vertebrates, the Primate order had the highest pooled seroprevalence at 13.1% (95% CI: 4.3–25.1%). From the three most studied primate genera we found the highest seroprevalence was in Alouatta (32.2%, 95% CI: 0.0–79.2%), followed by Callithrix (17.8%, 95% CI: 8.6–28.5%), and Cebus/Sapajus (3.7%, 95% CI: 0.0–11.1%). We further found that MAYV occurs in a wide range of vectors beyond Haemagogus spp. The quality of evidence behind these findings was variable and prompts calls for standardization of reporting of arbovirus occurrence. These findings support further risk emergence prediction, guide field surveillance efforts, and prompt further in-vivo studies to better define the ecological drivers of MAYV maintenance and potential for emergence.     

Analysis of the structure and dynamics of human serum albumin

A computational approach to designing a peptide-based ligand for the purification of human serum albumin (HSA) was undertaken using molecular docking and molecular dynamics (MD) simulation. A three-step procedure was performed to design a specific ligand for HSA. Based on the candidate pocket structure of HSA (warfarin binding site), a peptide library was built. These peptides were then docked into the pocket of HSA using the GOLD program. The GOLDscore values were used to determine the affinity of peptides for HSA. Consequently, the dipeptide Trp–Trp, which shows a high GOLDscore value, was selected and linked to a spacer arm of Lys[CO(CH₂)₅NH] on the surface of ECH-lysine sepharose 4 gel. For further evaluation, the Autodock Vina program was used to dock the linked compound into the pocket of HSA. The docking simulation was performed to obtain a first guess of the binding structure of the spacer–Trp–Trp–HSA complex and subsequently analyzed by MD simulations to assess the reliability of the docking results. These MD simulations indicated that the ligand–HSA complex remains stable, and water molecules can bridge between the ligand and the protein by hydrogen bonds. Finally, absorption spectroscopic studies were performed to illustrate the appropriateness of the binding affinity of the designed ligand toward HSA. These studies demonstrate that the designed dipeptide can bind preferentially to the warfarin binding site. Graphical Abstract

Three-step computational approach to the design of a dipeptide ligand for human serum albumin purification exploiting structure-based docking and molecular dynamics simulation     

Influence of the lipid membrane environment on structure and activity of the outer membrane protein Ail from Yersinia pestis

The surrounding environment has significant consequences for the structural and functional properties of membrane proteins. While native structure and function can be reconstituted in lipid bilayer membranes, the detergents used for protein solubilization are not always compatible with biological activity and, hence, not always appropriate for direct detection of ligand binding by NMR spectroscopy. Here we describe how the sample environment affects the activity of the outer membrane protein Ail (attachment invasion locus) from Yersinia pestis. Although Ail adopts the correct β-barrel fold in micelles, the high detergent concentrations required for NMR structural studies are not compatible with the ligand binding functionality of the protein. We also describe preparations of Ail embedded in phospholipid bilayer nanodiscs, optimized for NMR studies and ligand binding activity assays. Ail in nanodiscs is capable of binding its human ligand fibronectin and also yields high quality NMR spectra that reflect the proper fold. Binding activity assays, developed to be performed directly with the NMR samples, show that ligand binding involves the extracellular loops of Ail. The data show that even when detergent micelles support the protein fold, detergents can interfere with activity in subtle ways.     

Binding Adaptation of GS-441524 Diversifies Macro Domains and Downregulates SARS-CoV-2 de-MARylation Capacity

Viral infection in cells triggers a cascade of molecular defense mechanisms to maintain host-cell homoeostasis. One of these mechanisms is ADP-ribosylation, a fundamental post-translational modification (PTM) characterized by the addition of ADP-ribose (ADPr) on substrates. Poly(ADP-ribose) polymerases (PARPs) are implicated in this process and they perform ADP-ribosylation on host and pathogen proteins. Some viral families contain structural motifs that can reverse this PTM. These motifs known as macro domains (MDs) are evolutionarily conserved protein domains found in all kingdoms of life. They are divided in different classes with the viral belonging to Macro-D-type class because of their properties to recognize and revert the ADP-ribosylation. Viral MDs are potential pharmaceutical targets, capable to counteract host immune response. Sequence and structural homology between viral and human MDs are an impediment for the development of new active compounds against their function. Remdesivir, is a drug administrated in viral infections inhibiting viral replication through RNA-dependent RNA polymerase (RdRp). Herein, GS-441524, the active metabolite of the remdesivir, is tested as a hydrolase inhibitor for several viral MDs and for its binding to human homologs found in PARPs. This study presents biochemical and biophysical studies, which indicate that GS-441524 selectively modifies SARS-CoV-2 MD de-MARylation activity, while it does not interact with hPARP14 MD2 and hPARP15 MD2. The structural investigation of MD?GS-441524 complexes, using solution NMR and X-ray crystallography, discloses the impact of certain amino acids in ADPr binding cavity suggesting that F360 and its adjacent residues tune the selective binding of the inhibitor to SARS-CoV-2 MD.     

Functional and structural characterization of the integrase from the prototype foamy virus

Establishment of the stable provirus is an essential step in retroviral replication, orchestrated by integrase (IN), a virus-derived enzyme. Until now, available structural information was limited to the INs of human immunodeficiency virus type 1 (HIV-1), avian sarcoma virus (ASV) and their close orthologs from the Lentivirus and Alpharetrovirus genera. Here, we characterized the in vitro activity of the prototype foamy virus (PFV) IN from the Spumavirus genus and determined the three-dimensional structure of its catalytic core domain (CCD). Recombinant PFV IN displayed robust and almost exclusively concerted integration activity in vitro utilizing donor DNA substrates as short as 16 bp, underscoring its significance as a model for detailed structural studies. Comparison of the HIV-1, ASV and PFV CCD structures highlighted both conserved as well as unique structural features such as organization of the active site and the putative host factor binding face. Despite possessing very limited sequence identity to its HIV counterpart, PFV IN was sensitive to HIV IN strand transfer inhibitors, suggesting that this class of inhibitors target the most conserved features of retroviral IN-DNA complexes.     

Crystal Structure of the Mp1p Ligand Binding Domain 2 Reveals Its Function as a Fatty Acid-binding Protein

Penicillium marneffei is a dimorphic, pathogenic fungus in Southeast Asia that mostly afflicts immunocompromised individuals. As the only dimorphic member of the genus, it goes through a phase transition from a mold to yeast form, which is believed to be a requisite for its pathogenicity. Mp1p, a cell wall antigenic mannoprotein existing widely in yeast, hyphae, and conidia of the fungus, plays a vital role in host immune response during infection. To understand the function of Mp1p, we have determined the x-ray crystal structure of its ligand binding domain 2 (LBD2) to 1.3 Å. The structure reveals a dimer between the two molecules. The dimer interface forms a ligand binding cavity, in which electron density was observed for a palmitic acid molecule interacting with LBD2 indirectly through hydrogen bonding networks via two structural water molecules. Isothermal titration calorimetry experiments measured the ligand binding affinity (K_d) of Mp1p at the micromolar level. Mutations of ligand-binding residues, namely S313A and S332A, resulted in a 9-fold suppression of ligand binding affinity. Analytical ultracentrifugation assays demonstrated that both LBD2 and Mp1p are mostly monomeric in vitro, no matter with or without ligand, and our dimeric crystal structure of LBD2 might be the result of crystal packing. Based on the conformation of the ligand-binding pocket in the dimer structure, a model for the closed, monomeric form of LBD2 is proposed. Further structural analysis indicated the biological importance of fatty acid binding of Mp1p for the survival and pathogenicity of the conditional pathogen.     

Exploration of human serum albumin binding sites by docking and molecular dynamics flexible ligand–protein interactions

Five‐nanosecond molecular dynamics (MD) simulations were performed on human serum albumin (HSA) to study the conformational features of its primary ligand binding sites (I and II). Additionally, 11 HSA snapshots were extracted every 0.5 ns to explore the binding affinity (K_d) of 94 known HSA binding drugs using a blind docking procedure. MD simulations indicate that there is considerable flexibility for the protein, including the known sites I and II. Movements at HSA sites I and II were evidenced by structural analyses and docking simulations. The latter enabled the study and analysis of the HSA–ligand interactions of warfarin and ketoprofen (ligands binding to sites I and II, respectively) in greater detail. Our results indicate that the free energy values by docking (K_d observed) depend upon the conformations of both HSA and the ligand. The 94 HSA–ligand binding K_d values, obtained by the docking procedure, were subjected to a quantitative structure‐activity relationship (QSAR) study by multiple regression analysis. The best correlation between the observed and QSAR theoretical (K_d predicted) data was displayed at 2.5 ns. This study provides evidence that HSA binding sites I and II interact specifically with a variety of compounds through conformational adjustments of the protein structure in conjunction with ligand conformational adaptation to these sites. These results serve to explain the high ligand‐promiscuity of HSA. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 161–170, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Solution NMR studies provide structural basis for endotoxin pattern recognition by the innate immune receptor CD14

CD14 functions as a key pattern recognition receptor for a diverse array of Gram-negative and Gram-positive cell-wall components in the host innate immune response by binding to pathogen-associated molecular patterns (PAMPs) at partially overlapping binding site(s). To determine the potential contribution of CD14 residues in this pattern recognition, we have examined using solution NMR spectroscopy, the binding of three different endotoxin ligands, lipopolysaccharide, lipoteichoic acid, and a PGN-derived compound, muramyl dipeptide to a ¹⁵N isotopically labeled 152-residue N-terminal fragment of sCD14 expressed in Pichia pastoris. Mapping of NMR spectral changes upon addition of ligands revealed that the pattern of residues affected by binding of each ligand is partially similar and partially different. This first direct structural observation of the ability of specific residue combinations of CD14 to differentially affect endotoxin binding may help explain the broad specificity of CD14 in ligand recognition and provide a structural basis for pattern recognition. Another interesting finding from the observed spectral changes is that the mode of binding may be dynamically modulated and could provide a mechanism for binding endotoxins with structural diversity through a common binding site.     

Modeling the Epidermal Growth Factor—Epidermal Growth Factor Receptor L2 Domain Interaction: Implications for the Ligand Binding Process

Abstract

Signaling from the epidermal growth factor (EGF) receptor is triggered by the binding of lig-ands such as EGF or transforming growth factor alpha (TGF-α) and subsequent receptor dimerization. An understanding of these processes has been hindered by the lack of structural information about the ligand-bound, dimerized EGF receptor. Using an NMR-derived structure of EGF and a homology model of the major ligand binding domain of the EGF receptor and experimental data, we modeled the binding of EGF to this EGF receptor fragment. In this low resolution model of the complex, EGF sits across the second face of the EGF receptor L2 domain and EGF residues 10–16, 36–37, 40–47 bind to this face. The structural model is largely consistent with previously published NMR data describing the residues of TGF-α which interact strongly with the EGF receptor. Other EGF residues implicated in receptor binding are accounted by our proposal that the ligand binding is a two-step process with the EGF binding to at least one other site of the receptor. This three-dimensional model is expected to be useful in the design of ligand-based antagonists of the receptor.     

Nucleolin interacts with the rabbit hemorrhagic disease virus replicase RdRp,nonstructural proteins p16 and p23, playing a role in virus replication

Rabbit hemorrhagic disease virus (RHDV) is a member of the Caliciviridae family and cannot be propagated in vitro, which has impeded the progress of investigating its replication mechanism. Construction of an RHDV replicon system has recently provided a platform for exploring RHDV replication in host cells. Here, aided by this replicon system and using two-step affinity purification, we purified the RHDV replicase and identified its associated host factors. We identified rabbit nucleolin (NCL) as a physical link, which mediating the interaction between other RNA-dependent RNA polymerase (RdRp)-related host proteins and the viral replicase RdRp. We found that the overexpression or knockdown of NCL significantly increased or severely impaired RHDV replication in RK-13 cells, respectively. NCL was identified to directly interact with RHDV RdRp, p16, and p23. Furthermore, NCL knockdown severely impaired the binding of RdRp to RdRp-related host factors. Collectively, these results indicate that the host protein NCL is essential for RHDV replication and acts as a physical link between viral replicase and host proteins.