Functional motions of influenza virus hemagglutinin: a structure-based analytical approach

Breathing and Tilting: Mesoscale Simulations Illuminate Influenza Glycoprotein Vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from convalescent human donor, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Biosensing epidemic and pandemic respiratory viruses: Internet of Things with Gaussian noise channel algorithmic model

The current world condition is dire due to epidemics and pandemics as a result of novel viruses, such as influenza and the coronavirus, causing acute respiratory syndrome. To overcome these critical situations, the current research seeks to generate a common surveillance system with the assistance of a controlled Internet of Things operated under a Gaussian noise channel. To create the model system, a study with an analysis of H1N1 influenza virus determination on an interdigitated electrode (IDE) sensor was validated by current-volt measurements. The preliminary data were generated using hemagglutinin as the target against gold-conjugated aptamer/antibody as the probe, with the transmission pattern showing consistency with the Gaussian noise channel algorithm. A good fit with the algorithmic values was found, displaying a similar pattern to that output from the IDE, indicating reliability. This study can be a model for the surveillance of varied pathogens, including the emergence and reemergence of novel strains.

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

COVID-19: Myths and Reality

COVID-19, a new human respiratory disease that has killed nearly 3 million people in a year since the start of the pandemic, is a global public health challenge. Its infectious agent, SARS-CoV-2, differs from other coronaviruses in a number of structural features that make this virus more pathogenic and transmissible. In this review, we discuss some important characteristics of the main SARS-CoV-2 surface antigen, the spike (S) protein, such as (i) ability of the receptor-binding domain (RBD) to switch between the "standing-up" position (open pre-fusion conformation) for receptor binding and the "lying-down" position (closed pre-fusion conformation) for immune system evasion; (ii) advantage of a high binding affinity of the RBD open conformation to the human angiotensin-converting enzyme 2 (ACE2) receptor for efficient cell entry; and (iii) S protein preliminary activation by the intracellular furin-like proteases for facilitation of the virus spreading across different cell types. We describe interactions between the S protein and cellular receptors, co-receptors, and antagonists, as well as a hypothetical mechanism of the homotrimeric spike structure destabilization that triggers the fusion of the viral envelope with the cell membrane at physiological pH and mediates the viral nucleocapsid entry into the cytoplasm. The transition of the S protein pre-fusion conformation to the post-fusion one on the surface of virions after their treatment with some reagents, such as β-propiolactone, is essential, especially in relation to the vaccine production. We also compare the COVID-19 pathogenesis with that of severe outbreaks of "avian" influenza caused by the A/H5 and A/H7 highly pathogenic viruses and discuss the structural similarities between the SARS-CoV-2 S protein and hemagglutinins of those highly pathogenic strains. Finally, we touch on the prospective and currently used COVID-19 antiviral and anti-pathogenetic therapeutics, as well as recently approved conventional and innovative COVID-19 vaccines and their molecular and immunological features.

Normal mode analysis of membrane protein dynamics using the vibrational subsystem analysis

The vibrational subsystem analysis is a useful approach that allows for evaluating the spectrum of modes of a given system by integrating out the degrees of freedom accessible to the environment. The approach could be utilized for exploring the collective dynamics of a membrane protein (system) coupled to the lipid bilayer (environment). However, the application to membrane proteins is limited due to high computational costs of modeling a sufficiently large membrane environment unbiased by end effects, which drastically increases the size of the investigated system. We derived a recursive formula for calculating the reduced Hessian of a membrane protein embedded in a lipid bilayer by decomposing the membrane into concentric cylindrical domains with the protein located at the center. The approach allows for the design of a time- and memory-efficient algorithm and a mathematical understanding of the convergence of the reduced Hessian with respect to increasing membrane sizes. The application to the archaeal aspartate transporter GltPh illustrates its utility and efficiency in capturing the transporter's elevator-like movement during its transition between outward-facing and inward-facing states.

Potent sialic acid inhibitors that target influenza A virus hemagglutinin

Eradicating influenza A virus (IAV) is difficult, due to its genetic drift and reassortment ability. As the infectious cycle is initiated by the influenza glycoprotein, hemagglutinin (HA), which mediates the binding of virions to terminal sialic acids moieties, HA is a tempting target of anti-influenza inhibitors. However, the complexity of the HA structure has prevented delineation of the structural characterization of the HA protein-ligand complex. Our computational strategy efficiently analyzed > 200,000 records of compounds held in the United States National Cancer Institute (NCI) database and identified potential HA inhibitors, by modeling the sialic acid (SA) receptor binding site (RBS) for the HA structure. Our modeling revealed that compound NSC85561 showed significant antiviral activity against the IAV H1N1 strain with EC50 values ranging from 2.31 to 2.53 µM and negligible cytotoxicity (CC50 > 700 µM). Using the NSC85561 compound as the template to generate 12 derivatives, robust bioassay results revealed the strongest antiviral efficacies with NSC47715 and NSC7223. Virtual screening clearly identified three SA receptor binding site inhibitors that were successfully validated in experimental data. Thus, our computational strategy has identified SA receptor binding site inhibitors against HA that show IAV-associated antiviral activity.

High-Yield Expression and Purification of Recombinant Influenza Virus Proteins from Stably-Transfected Mammalian Cell Lines

Influenza viruses infect millions of people each year, resulting in significant morbidity and mortality in the human population. Therefore, generation of a universal influenza virus vaccine is an urgent need and would greatly benefit public health. Recombinant protein technology is an established vaccine platform and has resulted in several commercially available vaccines. Herein, we describe the approach for developing stable transfected human cell lines for the expression of recombinant influenza virus hemagglutinin (HA) and recombinant influenza virus neuraminidase (NA) proteins for the purpose of in vitro and in vivo vaccine development. HA and NA are the main surface glycoproteins on influenza virions and the major antibody targets. The benefits for using recombinant proteins for in vitro and in vivo assays include the ease of use, high level of purity and the ability to scale-up production. This work provides guidelines on how to produce and purify recombinant proteins produced in mammalian cell lines through either transient transfection or generation of stable cell lines from plasmid creation through the isolation step via Immobilized Metal Affinity Chromatography (IMAC). Collectively, the establishment of this pipeline has facilitated large-scale production of recombinant HA and NA proteins to high purity and with consistent yields, including glycosylation patterns that are very similar to proteins produced in a human host.

Morin Hydrate Inhibits Influenza Virus entry into Host Cells and Has Anti-inflammatory Effect in Influenza-infected Mice

Influenza virus is the major cause of seasonal and pandemic flu. Currently, oseltamivir, a potent and selective inhibitor of neuraminidase of influenza A and B viruses, is the drug of choice for treating patients with influenza virus infection. However, recent emergence of oseltamivir-resistant influenza viruses has limited its efficacy. Morin hydrate (3,5,7,2′,4′-pentahydroxyflavone) is a flavonoid isolated from Morus alba L. It has antioxidant, anti-inflammatory, neuroprotective, and anticancer effects partly by the inhibition of the NF-кB signaling pathway. However, its effects on influenza virus have not been studied. We evaluated the antiviral activity of morin hydrate against influenza A/Puerto Rico/8/1934 (A/PR/8; H1N1) and oseltamivir-resistant A/PR/8 influenza viruses in vitro. To determine its mode of action, we carried out time course experiments, and time of addition, hemolysis inhibition, and hemagglutination assays. The effects of the co-administration of morin hydrate and oseltamivir were assessed using the murine model of A/PR/8 infection. We found that morin hydrate reduced hemagglutination by A/PR/8 in vitro. It alleviated the symptoms of A/PR/8-infection, and reduced the levels of pro-inflammatory cytokines and chemokines, such as TNF-α and CCL2, in infected mice. Co-administration of morin hydrate and oseltamivir phosphate reduced the virus titers and attenuated pulmonary inflammation. Our results suggest that morin hydrate exhibits antiviral activity by inhibiting the entry of the virus.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Tilapia Lake Virus Does Not Hemagglutinate Avian and Piscine Erythrocytes and NH4Cl Does Not Inhibit Viral Replication In Vitro

Tilapia lake virus (TiLV) is a negative-sense single-stranded RNA (-ssRNA) icosahedral virus classified to be the only member in the family Amnoonviridae. Although TiLV segment-1 shares homology with the influenza C virus PB1 and has four conserved motifs similar to influenza A, B, and C polymerases, it is unknown whether there are other properties shared between TiLV and orthomyxovirus. In the present study, we wanted to determine whether TiLV agglutinated avian and piscine erythrocytes, and whether its replication was inhibited by lysosomotropic agents, such as ammonium chloride (NH₄Cl), as seen for orthomyxoviruses. Our findings showed that influenza virus strain A/Puerto Rico/8 (PR8) was able to hemagglutinate turkey (Meleagris gallopavo), Atlantic salmon (Salmo salar L), and Nile tilapia (Oreochromis niloticus) red blood cells (RBCs), while infectious salmon anemia virus (ISAV) only agglutinated Atlantic salmon, but not turkey or tilapia, RBCs. In contrast to PR8 and ISAV, TiLV did not agglutinate turkey, Atlantic salmon, or tilapia RBCs. qRT-PCR analysis showed that 30 mM NH₄Cl, a basic lysosomotropic agent, neither inhibited nor enhanced TiLV replication in E-11 cells. There was no difference in viral quantities in the infected cells with or without NH₄Cl treatment during virus adsorption or at 1, 2, and 3 h post-infection. Given that hemagglutinin proteins that bind RBCs also serve as ligands that bind host cells during virus entry leading to endocytosis in orthomyxoviruses, the data presented here suggest that TiLV may use mechanisms that are different from orthomyxoviruses for entry and replication in host cells. Therefore, future studies should seek to elucidate the mechanisms used by TiLV for entry into host cells and to determine its mode of replication in infected cells.

REP 2139: Antiviral Mechanisms and Applications in Achieving Functional Control of HBV and HDV Infection

Nucleic acid polymers (NAPs) are broad spectrum antiviral agents whose antiviral activity in hepatitis B virus (HBV) infection is derived from their ability to block the release of the hepatitis B virus surface antigen (HBsAg). This pharmacological activity blocks replenishment of HBsAg in the circulation, allowing host mediated clearance. This effect has important clinical significance as the clearance of circulating HBsAg dramatically potentiates the ability of immunotherapies to restore functional control of HBV infection which persists after antiviral therapy is removed. These effects are reproducible in preclinical evaluations and in several clinical trials that have evaluated the activity of the lead NAP, REP 2139, in monotherapy and in combination with immunotherapy in hepatitis B e antigen (HBeAg) negative and HBeAg positive HBV infection and also in HBeAg negative HBV/hepatitis D virus (HDV) coinfection. These antiviral effects of REP 2139 are achieved in the absence of any direct immunostimulatory effect in the liver and also without any discernible direct interaction with viral components. The search for the host protein interaction with NAPs that drives their antiviral effects is ongoing, and the interaction targeted by REP 2139 within infected cells has not yet been elucidated. This article provides an updated review of available data on the effects of REP 2139 in HBV and HDV infection and the ability of REP 2139-based combination therapy to achieve functional control of HBV and HDV infection in patients.

Structural significance of residues 158-160 in the H3N2 hemagglutnin globular head: A computational study with implications in viral evolution and infection

Influenza A H3N2 has been linked to annual outbreaks within the human population attributable to continuous structural changes. H3N2 HA contains well identified antigenic sites and receptor-binding sites (RBS) that are possibly correlated to viral evolution and infection. However, the structural significance of amino acid residues associated with both viral evolution and infection were not fully demonstrated. Throughout this study, we generated and analyzed H3N2 HA models that represented the clade 3C.2 population (comprised of clades 3C.2, 3C.2a, and 3C.21 from the transitioning 2014-2018 H3N2 strains) and 3C.3a (from the 2016 H3N2 strain). Model quality estimation, structural analyses and superimposition, and network analytics of H3N2 HA1 evolution were performed. We found that the structural properties of residues 158-160 could influence the overall HA backbone. More specifically, amino acid substitutions at residues 159-160 affected the amino acid orientation at residue 158, thereby, causing the overall HA backbone structure to vary. Our results were consistent with 1968-2018 HA1 evolution. Taken together, we propose that our results would highlight the structural significance of residues 158-160 in HA1 for both antigenic drift and RBS.

Genetic and antigenic divergence in the influenza A(H3N2) virus circulating between 2016 and 2017 in Thailand

Influenza virus evolves rapidly due to the accumulated genetic variations on the viral sequence. Unlike in North America and Europe, influenza season in the tropical Southeast Asia spans both the rainy and cool seasons. Thus, influenza epidemiology and viral evolution sometimes differ from other regions, which affect the ever-changing efficacy of the vaccine. To monitor the current circulating influenza viruses in this region, we determined the predominant influenza virus strains circulating in Thailand between January 2016 and June 2017 by screening 7,228 samples from patients with influenza-like illness. During this time, influenza A(H3N2) virus was the predominant influenza virus detected. We then phylogenetically compared the hemagglutinin (HA) gene from a subset of these A(H3N2) strains (n = 62) to the reference sequences and evaluated amino acid changes in the dominant antigenic epitopes on the HA protein structure. The divergence of the circulating A(H3N2) from the A/Hong Kong/4801/2014 vaccine strain formed five genetic groups (designated I to V) within the 3C.2a clade. Our results suggest a marked drift of the current circulating A(H3N2) strains in Thailand, which collectively contributed to the declining predicted vaccine effectiveness (VE) from 74% in 2016 down to 48% in 2017.

Uncoupling of an ammonia channel as a mechanism of allosteric inhibition in anthranilate synthase of Serratia marcescens: dynamic and graph theoretical analysis

Network modeling and molecular dynamic studies reveal the perturbation in communication pathways as a mechanism of allosteric inhibition in anthranilate synthase.

Bridging Enzymatic Structure Function via Mechanics: A Coarse-Grain Approach

Flexibility is a central aspect of protein function, and ligand binding in enzymes involves a wide range of structural changes, ranging from large-scale domain movements to small loop or side-chain rearrangements. In order to understand how the mechanical properties of enzymes, and the mechanical variations that are induced by ligand binding, relate to enzymatic activity, we carried out coarse-grain Brownian dynamics simulations on a set of enzymes whose structures in the unbound and ligand-bound forms are available in the Protein Data Bank. Our results show that enzymes are remarkably heterogeneous objects from a mechanical point of view and that the local rigidity of individual residues is tightly connected to their part in the protein's overall structure and function. The systematic comparison of the rigidity of enzymes in their unbound and bound forms highlights the fact that small conformational changes can induce large mechanical effects, leading to either more or less flexibility depending on the enzyme's architecture and the location of its ligand-biding site. These mechanical variations target a limited number of specific residues that occupy key locations for enzymatic activity, and our approach thus offers a mean to detect perturbation-sensitive sites in enzymes, where the addition or removal of a few interactions will lead to important changes in the proteins internal dynamics.

In Silico Functional and Structural Characterization of H1N1 Influenza A Viruses Hemagglutinin, 2010-2013, Shiraz, Iran

Hemagglutinin (HA) is a major virulence factor of influenza viruses and plays an important role in viral pathogenesis. Analysis of amino acid changes, epitopes' regions, glycosylation and phosphorylation sites have greatly contributed to the development of new generations of vaccine. The hemagglutinins of 10 selected isolates, 8 of 2010 and 2 of 2013 samples were sequenced and analyzed by several bioinformatic softwares and the results were compared with those of 3 vaccine isolates. The study detected several amino acid changes related to altered epitopes' sites, modification sites and physico-chemical properties. The results showed some conserved modification sites in HA structure. This study is the first analytical research on isolates obtained from Shiraz, Iran, and our results can be used to better understand the genetic diversity and antigenic variations in Iranian and Asian H1N1 pathogenic strains.

Order and disorder control the functional rearrangement of influenza hemagglutinin

Influenza hemagglutinin (HA), a homotrimeric glycoprotein crucial for membrane fusion, undergoes a large-scale structural rearrangement during viral invasion. X-ray crystallography has shown that the pre- and postfusion configurations of HA2, the membrane-fusion subunit of HA, have disparate secondary, tertiary, and quaternary structures, where some regions are displaced by more than 100 Å. To explore structural dynamics during the conformational transition, we studied simulations of a minimally frustrated model based on energy landscape theory. The model combines structural information from both the pre- and postfusion crystallographic configurations of HA2. Rather than a downhill drive toward formation of the central coiled-coil, we discovered an order-disorder transition early in the conformational change as the mechanism for the release of the fusion peptides from their burial sites in the prefusion crystal structure. This disorder quickly leads to a metastable intermediate with a broken threefold symmetry. Finally, kinetic competition between the formation of the extended coiled-coil and C-terminal melting results in two routes from this intermediate to the postfusion structure. Our study reiterates the roles that cracking and disorder can play in functional molecular motions, in contrast to the downhill mechanical interpretations of the "spring-loaded" model proposed for the HA2 conformational transition.

Cleavage site and Ectodomain of HA2 sub-unit sequence of three equine influenza virus isolated in Morocco

Background

The equine influenza (EI) is an infectious and contagious disease of the upper respiratory tract of horses. Two outbreaks were notified in Morocco during 1997 and 2004 respectively in Nador and Essaouira. The aims of the present study concern the amino acids sequences comparison with reference strain A/equine/Miami/1963(H3N8) of the HA2 subunit including the cleavage site of three equine influenza viruses (H3N8) isolated in Morocco: A/equine/Nador/1/1997(H3N8), A/equine/Essaouira/2/2004 (H3N8) and A/equine/Essaouira/3/2004 (H3N8).

Results

The obtained results demonstrated that the substitutions were located at Ectodomain (ED) and transmembrane domain (TD), and they have only one arginine in cleavage site (HA1-PEKQI-R³²⁹-GI-HA2). In the Ectodomain, the mutation N/154²/T deleted the NGT glycosylation site at position 154 for both strains A/equine/Essaouira/2/2004(H3N8) and A/equine/Essaouira/3/2004(H3N8). Except for mutation D/160²/Y of the A/equine/Nador/1/1997(H3N8) strain, the other mutations were involved in non conserved sites. While the transmembrane domain (TM) of the strain A/equine/Essaouira/3/2004(H3N8) exhibits a substitution at residue C/199²/F. For the A/equine/Nador/1/1997(H3N8) strain the HA2 shows a mutation at residue M/207²/L. Three Moroccan strains reveals a common substitution at the residue E/211²/Q located between transmembrane domain TM and the cytoplasmic domain (CD).

Conclusion

The given nature virulence of three Moroccan strains, the identified and reported mutations certainly played a permissive role of infection viral process.

Homology modeling study toward identifying structural properties in the HA2 B-loop that would influence the HA1 receptor-binding site

Influenza hemagglutinin (HA) consists of a fibrous globular stem (HA2) inserted into the viral membrane supporting a globular head (HA1). HA1 receptor-binding has been hypothesized to be structurally correlated to the HA2 B-loop, however, this was never fully understood. Here, we elucidated the structural relationship between the HA2 B-loop and the HA1 receptor-binding site (RBS). Throughout this study, we analyzed 2486 H1N1 HA homology models obtained from human, swine and avian strains during 1976-2012. Quality of all homology models were verified before further analyses. We established that amino acid residue 882 is putatively strain-conserved and differs in the human (K882), swine (H882) and avian (N882) strains. Moreover, we observed that the amino acid at residue 882 and, similarly, its orientation has the potential to influence the HA1 RBS diameter measurements which we hypothesize may consequentially affect influenza H1N1 viral infectivity, immune escape, transmissibility, and evolution.

Identifying ligand binding conformations of the β2-adrenergic receptor by using its agonists as computational probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β(2)-adrenergic receptor (β(2)AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β(2)AR. We show that the select ligands bind preferentially to different predicted conformers of β(2)AR, and identify a role of β(2)AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance.

Potential strategies and biosafety protocols used for dual-use research on highly pathogenic influenza viruses

Influenza A viruses (IAVs), particularly the highly pathogenic avian influenza H5N1, have posed a substantial threat to public health worldwide. Although the laboratory generation of the mutant influenza virus H5N1 with airborne transmissibility among mammals, which has been considered as a dual-use research, may benefit the development of effective vaccines and therapeutics against the emerging infectious agents, it may also pose threats to national biosecurity, laboratory biosafety, and/or public health. This review introduces the classification and characterization of IAVs, pinpoints historic pandemics and epidemics caused by IAVs, emphasizes the significance and necessity of biosafety, summarizes currently established biosafety-related protocols for IAV research, and provides potential strategies to improve biosafety protocols for dual-use research on the highly pathogenic avian influenza viruses and other emerging infectious agents.

Insight into the structure, dynamics and the unfolding property of amylosucrases: implications of rational engineering on thermostability

Amylosucrase (AS) is a kind of glucosyltransferases (E.C. 2.4.1.4) belonging to the Glycoside Hydrolase (GH) Family 13. In the presence of an activator polymer, in vitro, AS is able to catalyze the synthesis of an amylose-like polysaccharide composed of only α-1,4-linkages using sucrose as the only energy source. Unlike AS, other enzymes responsible for the synthesis of such amylose-like polymers require the addition of expensive nucleotide-activated sugars. These properties make AS an interesting enzyme for industrial applications. In this work, the structures and topology of the two AS were thoroughly investigated for the sake of explaining the reason why Deinococcus geothermalis amylosucrase (DgAS) is more stable than Neisseria polysaccharea amylosucrase (NpAS). Based on our results, there are two main factors that contribute to the superior thermostability of DgAS. On the one hand, DgAS holds some good structural features that may make positive contributions to the thermostability. On the other hand, the contacts among residues of DgAS are thought to be topologically more compact than those of NpAS. Furthermore, the dynamics and unfolding properties of the two AS were also explored by the gauss network model (GNM) and the anisotropic network model (ANM). According to the results of GNM and ANM, we have found that the two AS could exhibit a shear-like motion, which is probably associated with their functions. What is more, with the discovery of the unfolding pathway of the two AS, we can focus on the weak regions, and hence designing more appropriate mutations for the sake of thermostability engineering. Taking the results on structure, dynamics and unfolding properties of the two AS into consideration, we have predicted some novel mutants whose thermostability is possibly elevated, and hopefully these discoveries can be used as guides for our future work on rational design.

Identification of motions in membrane proteins by elastic network models and their experimental validation

Identifying the functional motions of membrane proteins is difficult because they range from large-scale collective dynamics to local small atomic fluctuations at different timescales that are difficult to measure experimentally due to the hydrophobic nature of these proteins. Elastic Network Models, and in particular their most widely used implementation, the Anisotropic Network Model (ANM), have proven to be useful computational methods in many recent applications to predict membrane protein dynamics. These models are based on the premise that biomolecules possess intrinsic mechanical characteristics uniquely defined by their particular architectures. In the ANM, interactions between residues in close proximity are represented by harmonic potentials with a uniform spring constant. The slow mode shapes generated by the ANM provide valuable information on the global dynamics of biomolecules that are relevant to their function. In its recent extension in the form of ANM-guided molecular dynamics (MD), this coarse-grained approach is augmented with atomic detail. The results from ANM and its extensions can be used to guide experiments and thus speedup the process of quantifying motions in membrane proteins. Testing the predictions can be accomplished through (a) direct observation of motions through studies of structure and biophysical probes, (b) perturbation of the motions by, e.g., cross-linking or site-directed mutagenesis, and (c) by studying the effects of such perturbations on protein function, typically through ligand binding and activity assays. To illustrate the applicability of the combined computational ANM-experimental testing framework to membrane proteins, we describe-alongside the general protocols-here the application of ANM to rhodopsin, a prototypical member of the pharmacologically relevant G-protein coupled receptor family.

Molecular basis of the structure and function of H1 hemagglutinin of influenza virus

Influenza virus hemagglutinin (HA) contains antigenic sites recognized by the host immune system, cleavage sites cleaved by host proteases, receptor binding sites attaching to sialyl receptors on the target cell, and fusion peptides mediating membrane fusion. Change in an amino acid(s) in these sites may affect the potential of virus infection and spread within and between hosts. Influenza viruses with H1 HA infect birds, pigs and humans and have caused two of the four pandemics in the past 100 years: 1918 pandemic that killed 21-50 million people and 2009 pandemic that caused more than 18,000 deaths. Understanding the relationship between antigenic structure and immune specificity, the receptor binding specificity in virus transmission, how the cleavage site controls pathogenicity, and how the fusion peptide causes membrane fusion for the entry of influenza virus into the host cell should provide information to find more effective ways to prevent and control influenza.

The elastic network model reveals a consistent picture on intrinsic functional dynamics of type II restriction endonucleases

The vibrational dynamics of various type II restriction endonucleases, in complex with cognate/non-cognate DNA and in the apo form, are investigated with the elastic network model in order to reveal common functional mechanisms in this enzyme family. Scissor-like and tong-like motions observed in the slowest modes of all enzymes and their complexes point to common DNA recognition and cleavage mechanisms. Normal mode analysis further points out that the scissor-like motion has an important role in differentiating between cognate and non-cognate sequences at the recognition site, thus implying its catalytic relevance. Flexible regions observed around the DNA-binding site of the enzyme usually concentrate on the highly conserved β-strands, especially after DNA binding. These β-strands may have a structurally stabilizing role in functional dynamics for target site recognition and cleavage. In addition, hot spot residues based on high-frequency modes reveal possible communication pathways between the two distant cleavage sites in the enzyme family. Some of these hot spots also exist on the shortest path between the catalytic sites and are highly conserved.

Structural insight into the role of thrombospondin-1 binding to calreticulin in calreticulin-induced focal adhesion disassembly

Thrombospondin-1 (TSP1) binding to calreticulin (CRT) on the cell surface stimulates association of CRT with LDL receptor-related protein (LRP1) to signal focal adhesion disassembly and engagement of cellular activities. The structural basis for this phenomenon is unknown. We studied the binding thermodynamics of the TSP1-CRT complex and the conformational changes in CRT induced by binding to TSP1 with combined binding free energy analysis, molecular dynamics simulation, and anisotropic network model restrained molecular dynamics simulation. Results showed that mutations of Lys 24 and Lys 32 in TSP1 to Ala and of amino acids 24-26 and 32-34 in CRT to Ala significantly weakened the binding of TSP1 and CRT, which is consistent with experimental results. Upon validation of the calculated binding affinity changes of the TSP1-CRT complex by mutations in key residues in TSP1 and CRT with the experimental results, we performed conformational analyses to understand the role of TSP1 binding to CRT in the induction of conformational changes in CRT. Conformational analyses showed that TSP1 binding to CRT resulted in a more "open" conformation and a significant rotational change for the CRT N-domain with respect to the CRT P-domain, which could expose the potential binding site(s) in CRT for binding to LRP1 to signal focal adhesion disassembly. Results offer structural insight into the role of TSP1 binding to CRT in CRT-induced focal adhesion disassembly.

Vibrational subsystem analysis: A method for probing free energies and correlations in the harmonic limit

A new vibrational subsystem analysis (VSA) method is presented for coupling global motion to a local subsystem while including the inertial effects of the environment. The premise of the VSA method is a partitioning of a system into a smaller region of interest and a usually larger part referred to as environment. This method allows the investigation of local-global coupling, a more accurate estimation of vibrational free energy contribution for parts of a large system, and the elimination of the "tip effect" in elastic network model calculations. Additionally, the VSA method can be used as a probe of specific degrees of freedom that may contribute to free energy differences. The VSA approach can be employed in many ways, but it will likely be most useful for estimating activation free energies in QM/MM reaction path calculations. Four examples are presented to demonstrate the utility of this method.

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.

Protein mechanics: a route from structure to function

In order to better understand the mechanical properties of proteins, we have developed simulation tools which enable these properties to be analysed on a residue-by-residue basis. Although these calculations are relatively expensive with all-atom protein models, good results can be obtained much faster using coarse-grained approaches. The results show that proteins are surprisingly heterogeneous from a mechanical point of view and that functionally important residues often exhibit unusual mechanical behaviour. This finding offers a novel means for detecting functional sites and also potentially provides a route for understanding the links between structure and function in more general terms.

Allosteric transitions in the chaperonin GroEL are captured by a dominant normal mode that is most robust to sequence variations

The Escherichia coli chaperonin GroEL, which helps proteins to fold, consists of two heptameric rings stacked back-to-back. During the reaction cycle GroEL undergoes a series of allosteric transitions triggered by ligand (substrate protein, ATP, and the cochaperonin GroES) binding. Based on an elastic network model of the bullet-shaped double-ring chaperonin GroEL-(ADP)(7)-GroES structure (R''T state), we perform a normal mode analysis to explore the energetically favorable collective motions encoded in the R''T structure. By comparing each normal mode with the observed conformational changes in the R''T --> TR'' transition, a single dominant normal mode provides a simple description of this highly intricate allosteric transition. A detailed analysis of this relatively high-frequency mode describes the structural and dynamic changes that underlie the positive intra-ring and negative inter-ring cooperativity. The dynamics embedded in the dominant mode entails highly concerted structural motions with approximate preservation of sevenfold symmetry within each ring and negatively correlated ones between the two rings. The dominant normal mode (in comparison with the other modes) is robust to parametric perturbations caused by sequence variations, which validates its functional importance. Response of the dominant mode to local changes that mimic mutations using the structural perturbation method technique leads to a wiring diagram that identifies a network of key residues that regulate the allosteric transitions. Many of these residues are located in intersubunit interfaces, and may therefore play a critical role in transmitting allosteric signals between subunits.

Locating the active sites of enzymes using mechanical properties

We have applied the calculation of mechanical properties to a dataset of almost 100 enzymes to determine the extent to which catalytic residues have distinct properties. Specifically, we have calculated force constants describing the ease of moving any given amino acid residue with respect to the other residues in the protein. The results show that catalytic residues are invariably associated with high force constants. Choosing an appropriate cutoff enables the detection of roughly 80% of catalytic residues with only 25% of false positives. It is shown that neither multidomain structures, nor the presence or absence of bound ligands hinder successful detections. It is however noted that active sites near the protein surface are more difficult to detect and that non-catalytic, but structurally key residues may also exhibit high force constants.

Predisposition of the dark state of rhodopsin to functional changes in structure

As the only member of the family of G-protein-coupled receptors for which atomic coordinates are available, rhodopsin is widely studied for insight into the molecular mechanism of G-protein-coupled receptor activation. The currently available structures refer to the inactive, dark state, of rhodopsin, rather than the light-activated metarhodopsin II (Meta II) state. A model for the Meta II state is proposed here by analyzing elastic network normal modes in conjunction with experimental data. Key mechanical features and interactions broken/formed in the proposed model are found to be consistent with the experimental data. The model is further tested by using a set of Meta II fluorescence decay rates measured to empirically characterize the deactivation of rhodopsin mutants. The model is found to correctly predict 93% of the experimentally observed effects in 119 rhodopsin mutants for which the decay rates and misfolding data have been measured, including a systematic analysis of Cys-->Ser replacements reported here. Based on the detailed comparison between model and experiments, a cooperative activation mechanism is deduced that couples retinal isomerization to concerted changes in conformation, facilitated by the intrinsic dynamics of rhodopsin. A global hinge site is identified near the retinal-binding pocket that ensures the efficient propagation of signals from the central transmembrane region to both cytoplasmic and extracellular ends. The predicted activation mechanism opens the transmembrane helices at the critical G-protein binding cytoplasmic domain. This model provides a detailed, mechanistic description of the activation process, extending experimental observations and yielding new insights for further tests.

The Role of histidine residues in low-pH-mediated viral membrane fusion

A central event in the invasion of a host cell by an enveloped virus is the fusion of viral and cell membranes. For many viruses, membrane fusion is driven by specific viral surface proteins that undergo large-scale conformational rearrangements, triggered by exposure to low pH in the endosome upon internalization. Here, we present evidence suggesting that in both class I (helical hairpin proteins) and class II (beta-structure-rich proteins) pH-dependent fusion proteins the protonation of specific histidine residues triggers fusion via an analogous molecular mechanism. These histidines are located in the vicinity of positively charged residues in the prefusion conformation, and they subsequently form salt bridges with negatively charged residues in the postfusion conformation. The molecular surfaces involved in the corresponding structural rearrangements leading to fusion are highly conserved and thus might provide a suitable common target for the design of antivirals, which could be active against a diverse range of pathogenic viruses.

The relevance of salt bridges for the stability of the influenza virus hemagglutinin

Hemagglutinin (HA) of influenza virus undergoes an irreversible conformational change at acidic pH, mediating viral fusion with the host endosomal membrane. To unravel the molecular basis of the pH-dependent stability of HA, we demonstrate by mutagenesis of the prototype HA of virus strain X31 (H3 subtype) that salt bridges, especially a tetrad salt bridge within the monomers, are crucial for folding and stability of the trimeric ectodomain. This complex (tetrad) salt bridge is highly conserved among influenza virus subtypes. Introducing additional sites of electrostatic attraction between monomers in the distal region enhanced the stability of ectodomain at low pH mimicking the natural variant H2 subtype. We propose that distinct salt bridges in the distal domain may contribute to the enhanced stability of HA of natural virus variants. This hypothesis may provide clues to understanding adaptations of virus strains (for example, avian influenza viruses) in order to preserve stability of the protein in the host-specific environment.

Collective dynamics of EcoRI-DNA complex by elastic network model and molecular dynamics simulations

Anisotropic network model (ANM) is used to analyze the collective motions of restriction enzyme EcoRI in free form and in complex with DNA. For comparison, three independent molecular dynamics (MD) simulations, each of 1.5 ns duration, are also performed for the EcoRI-DNA complex in explicit water. Although high mobility (equilibrium fluctuations) of inner and outer loops that surround the DNA is consistent in both methods and experiments, MD runs sample different conformational subspaces from which reliable collective dynamics cannot be extracted. However, ANM employed on different conformations from MD simulations indicates very similar collective motions. The stems of the inner loops are quite immobile even in the free enzyme and form a large, almost fixed, pocket for DNA binding. As a result, the residues that make specific and non-specific interactions with the DNA exhibit very low fluctuations in the free enzyme. The vibrational entropy difference between the EcoRI complex and free protein + unkinked DNA is positive (favorable), which may partially counteract the unfavorable enthalpy difference of DNA kink formation. Dynamic domains in EcoRI complex and cross-correlations between residue fluctuations indicate possible means of communication between the distal active sites.

Anisotropic network model: systematic evaluation and a new web interface

MOTIVATION

The Anisotropic Network Model (ANM) is a simple yet powerful model for normal mode analysis of proteins. Despite its broad use for exploring biomolecular collective motions, ANM has not been systematically evaluated to date. A lack of a convenient interface has been an additional obstacle for easy usage.

RESULTS

ANM has been evaluated on a large set of proteins to establish the optimal model parameters that achieve the highest correlation with experimental data and its limits of accuracy and applicability. Residue fluctuations in globular proteins are shown to be more accurately predicted than those in nonglobular proteins, and core residues are more accurately described than solvent-exposed ones. Significant improvement in agreement with experiments is observed with increase in the resolution of the examined structure. A new server for ANM calculations is presented, which offers flexible options for controlling model parameters and output formats, interactive animation of collective modes and advanced graphical features.

AVAILABILITY

ANM server (http://www.ccbb.pitt.edu/anm)

Loop motions of triosephosphate isomerase observed with elastic networks

The internal dynamics of triosephosphate isomerase have been investigated with elastic networks, with and without a substrate bound. The slowest modes of motion involve large domain motions but also a loop motion that conforms to the changes observed between the crystal structures and . Our computations confirm that the different motions of this loop are important in several of the computed slowest modes. We have shown that elastic network computations on this protein system can combine atoms for the functional parts of the structure with coarse-grained (cg) representations of the remainder of the structure in several different ways. Similar loop motions are seen with elastic network models for atomistic and mixed cg models. The loop motions are reproduced with an overlap of 0.75-0.79 by combining the four slowest modes of motion for the free and complex forms of the enzyme.

Identification of key residues in proteins by using their physical characters

Key residues in proteins are important to their stability, folding, and functions. They usually are highly conserved and can be identified by sequence or structure alignments. However, these methods can only determine the locations of key residues in sequences and structures and give less information about their physical characters. In this paper, we try to identify key residues by analyzing their inter-residue interactions. The model we study is the protein domain from transducin. We show that the usual Gaussian network analysis and distance-based contact analysis have difficulty identifying the key residues in this protein, but the contact energies can do it well. We find that most key residues can be located by the lowest contact energies. This enables us to predict and analyze the key residues in other proteins. Our results suggest that contact energy analysis may provide an alternative approach to investigating the folding and stability of proteins.

Specificity of trypsin and chymotrypsin: loop-motion-controlled dynamic correlation as a determinant

Trypsin and chymotrypsin are both serine proteases with high sequence and structural similarities, but with different substrate specificity. Previous experiments have demonstrated the critical role of the two loops outside the binding pocket in controlling the specificity of the two enzymes. To understand the mechanism of such a control of specificity by distant loops, we have used the Gaussian network model to study the dynamic properties of trypsin and chymotrypsin and the roles played by the two loops. A clustering method was introduced to analyze the correlated motions of residues. We have found that trypsin and chymotrypsin have distinct dynamic signatures in the two loop regions, which are in turn highly correlated with motions of certain residues in the binding pockets. Interestingly, replacing the two loops of trypsin with those of chymotrypsin changes the motion style of trypsin to chymotrypsin-like, whereas the same experimental replacement was shown necessary to make trypsin have chymotrypsin's enzyme specificity and activity. These results suggest that the cooperative motions of the two loops and the substrate-binding sites contribute to the activity and substrate specificity of trypsin and chymotrypsin.

Collective motions of RNA polymerases. Analysis of core enzyme, elongation complex and holoenzyme

The anisotropic network model (ANM), a coarse-grained normal mode analysis, is used to study the vibrational dynamics of RNA polymerases (RNAP) around the native states. The theoretical temperature factors obtained from ANM are in conformity with the experimental values for yeast and bacterial RNAP structures in free and complex forms. In the low-frequency collective modes that are related to biological function, both bacterial and yeast RNAPs with a crab claw shape display an opening/closing of the cleft due to the rigid-body motion of the clamp (bottom pincer), which has been also predicted by experiments, together with the motion of the top pincer. Even though slightly lower fluctuations are observed in the elongation complex of yeast RNAP, similar clamp motion still exists in collective modes, which should be concerted with the flexible switches and the bridge helix in driving the transcription process, pointing at the possibility of a ratchet-like mechanism. Two different bacterial holoenzyme (HE) structures are studied, which may have functional significance at different stages of transcription initiation. In a specific closed conformation of the HE, the clamp and top pincer are highly immobilized due to interactions with the sigma subunit. In contrast, the deformation of the top pincer is not inhibited in a relatively open conformation of another HE, which may help load the DNA into the cleft during transcription initiation, even though the clamp motion is still inhibited.

Probing protein mechanics: residue-level properties and their use in defining domains

It is becoming clear that, in addition to structural properties, the mechanical properties of proteins can play an important role in their biological activity. It nevertheless remains difficult to probe these properties experimentally. Whereas single-molecule experiments give access to overall mechanical behavior, notably the impact of end-to-end stretching, it is currently impossible to directly obtain data on more local properties. We propose a theoretical method for probing the mechanical properties of protein structures at the single-amino acid level. This approach can be applied to both all-atom and simplified protein representations. The probing leads to force constants for local deformations and to deformation vectors indicating the paths of least mechanical resistance. It also reveals the mechanical coupling that exists between residues. Results obtained for a variety of proteins show that the calculated force constants vary over a wide range. An analysis of the induced deformations provides information that is distinct from that obtained with measures of atomic fluctuations and is more easily linked to residue-level properties than normal mode analyses or dynamic trajectories. It is also shown that the mechanical information obtained by residue-level probing opens a new route for defining so-called dynamical domains within protein structures.

Normal-modes-based prediction of protein conformational changes guided by distance constraints

Based on the elastic network model, we develop a novel method that predicts the conformational change of a protein complex given its initial-state crystal structure together with a small set of pairwise distance constraints for the end state. The predicted conformational change, which is a linear combination of multiple low-frequency normal modes that are solved from the elastic network model, is computed as a response displacement induced by a perturbation to the system Hamiltonian that incorporates the given distance constraints. For a list of test cases, we find that the computed response displacement overlaps significantly with the measured conformational changes, when only a handful of pairwise constraints are used (</=10). The performance of this method is also shown to be robust against different choices of pairwise distance constraints and errors in their values. This method, if supplied with the experimentally derived distance constraints (for example, from NMR or other spectroscopic measurements), can be applied to the analysis of protein conformational changes toward transient states.

How does protein architecture facilitate the transduction of ATP chemical-bond energy into mechanical work? The cases of nitrogenase and ATP binding-cassette proteins

Transduction of adenosine triphosphate (ATP) chemical-bond energy into work to drive large-scale conformational changes is common in proteins. Two specific examples of ATP-utilizing proteins are the nitrogenase iron protein and the ATP binding-cassette transporter protein, BtuCD. Nitrogenase catalyzes biological nitrogen fixation whereas BtuCD transports vitamin B(12) across membranes. Both proteins drive their reactions with ATP. To interpret how the mechanical force generated by ATP binding and hydrolysis is propagated in these proteins, a coarse-grained elastic network model is employed. The analysis shows that subunits of the proteins move against each other in a concerted manner. The lowest-frequency modes of the nitrogenase iron protein and of the ATP binding-cassette transporter BtuCD protein are found to link the functionally critical domains, and these modes are suggested to be responsible for (at least the initial stages) large-scale ATP-coupled conformational changes.

In vitro and in vivo inhibitory effects of disodium cromoglycate on influenza virus infection

Disodium cromoglycate (DSCG) is one of the safest drugs for the prevention of bronchial asthma and allergic rhinitis attacks. The effect of DSCG on acute upper respiratory tract viral infection is still controversial. Here we investigated DSCG inhibition of influenza virus infection in vivo and in vitro. In vivo effects of DSCG on viral infection were assessed using a murine model of respiratory tract infection. Intranasal administration of DSCG protected mice from death induced by infection with influenza virus A/PR/8/34. We analyzed DSCG anti-viral effects in vitro by either (i) treating cells prior to viral adsorption, (ii) treating cells concurrently with viral adsorption, or (iii) treating cells after viral adsorption. DSCG treatment of cells during or after, but not before, viral adsorption significantly inhibited influenza viral infection, indicating DSCG acts on events late in viral infection. DSCG exerts anti-influenza effect both in vitro and in vivo at the doses compatible with treatment for asthma. DSCG marginally inhibited influenza viral neuraminidase and membrane fusion functions, suggesting that DSCG inhibition of viral neuraminidase and fusion activities may partially mediate this anti-influenza effect. Our results indicate that treatment of patients including children with DSCG may take advantages for prevention from influenza virus infection.

Global ribosome motions revealed with elastic network model

The motions of large systems such as the ribosome are not fully accessible with conventional molecular simulations. A coarse-grained, less-than-atomic-detail model such as the anisotropic network model (ANM) is a convenient informative tool to study the cooperative motions of the ribosome. The motions of the small 30S subunit, the larger 50S subunit, and the entire 70S assembly of the two subunits have been analyzed using ANM. The lowest frequency collective modes predicted by ANM show that the 50S subunit and 30S subunit are strongly anti-correlated in the motion of the 70S assembly. A ratchet-like motion is observed that corresponds well to the experimentally reported ratchet motion. Other slow modes are also examined because of their potential links to the translocation steps in the ribosome. We identify several modes that may facilitate the E-tRNA exiting from the assembly. The A-site t-RNA and P-site t-RNA are found to be strongly coupled and positively correlated in these slow modes, suggesting that the translocations of these two t-RNAs occur simultaneously, while the motions of the E-site t-RNA are less correlated, and thus less likely to occur simultaneously. Overall the t-RNAs exhibit relatively large deformations. Animations of these slow modes of motion can be viewed at.