Determining functionally important amino acid residues of the E1 protein of Venezuelan equine encephalitis virus

Application of temperature-dependent and steered molecular dynamics simulation to screen anti-dengue compounds against Marburg virus

Marburg virus infections are extremely fatal with a fatality range of 23% to 90%, therefore there is an urgent requirement to design and develop efficient therapeutic molecules. Here, a comprehensive temperature-dependent molecular dynamics (MD) simulation method was implemented to identify the potential molecule from the anti-dengue compound library that can inhibit the function of the VP24 protein of Marburg. Virtual high throughput screening identified five effective binders of VP24 after screening 484 anti-dengue compounds. These compounds were treated in MD simulation at four different temperatures: 300, 340, 380, and 420 K. Higher temperatures showed dissociation of hit compounds from the protein. Further, triplicates of 100 ns MD simulation were conducted which showed that compounds ID = 118717693, and ID = 5361 showed strong stability with the protein molecule. These compounds were further validated using Δ G binding free energies and they showed: -30.38 kcal/mol, and -67.83 kcal/mol binding free energies, respectively. Later, these two compounds were used in steered MD simulation to detect its dissociation. Compound ID = 5361 showed the maximum pulling force of 199.02 kcal/mol/nm to dissociate the protein-ligand complex while ID = 118717693 had a pulling force of 101.11 kcal/mol/nm, respectively. This ligand highest number of hydrogen bonds with varying occupancies at 89.93%, 69.80%, 57.93%, 52.33%, and 50.63%. This study showed that ID = 5361 can bind with the VP24 strongly and has the potential to inhibit its function which can be validated in the in-vitro experiment.Communicated by Ramaswamy H. Sarma.

In silico characterization, molecular docking, and dynamic simulation of a novel fungal cell-death suppressing effector, MoRlpA as potential cathepsin B-like cysteine protease inhibitor during rice blast infection

The blast fungus Magnaporthe oryzae is one of the most notorious pathogens affecting rice production worldwide. The cereal killer employs a special class of small secreted proteins called effectors to manipulate and perturb the host metabolism. In turn, the host plants trigger effector-triggered immunity (ETI) via localized cell death and hypersensitive response (HR). We have identified and characterized a novel secreted effector MoRlpA from M. oryzae by extensive in silico methods. The localization studies suggested that it is exclusively secreted in the host apoplasts. Interestingly, MoRlpA interacts with a protease, cathepsin B from rice with highest affinity. The 3D structural models of both the proteins were generated. Cathepsin B-like cysteine proteases are usually involved in programmed cell death (PCD) and autophagy in plants which lead to generation of HR upon infection. Our results suggest that MoRlpA interacts with rice cathepsin B-like cysteine protease and demolish the host counter-attack by suppressing cell death and HR during an active blast infection. This was further validated by molecular docking and molecular dynamic simulation analyses. The important residues involved in the rice-blast pathogen interactions were deciphered. Overall, this research highlights stable interactions between MoRlpA-OsCathB during rice blast pathogenesis and providing an insight into how this novel RlpA protease inhibitor-cum-effector modulates the host's apoplast to invade the host tissues and establish a successful infection. Thus, this research will help to develop potential fungicide to block the binding region of MoRlpA target so that the cryptic pathogen would be recognized by the host. HIGHLIGHTSFor the first time, a novel secreted effector protein, MoRlpA has been identified and characterised from M. oryzae in silicoMoRlpA contains a rare lipoprotein A-like DPBB domain which is often an enzymatic domain in other systemsMoRlpA as an apoplastic effector interacts with the rice protease OsCathB to suppress the cell death and hypersensitive response during rice blast infectionThe three-dimensional structures of both the MoRlpA and OsCathB proteins were predictedMoRlpA-OsCathB interactions were analysed by molecular docking and molecular dynamic simulation studiesCommunicated by Ramaswamy H. Sarma.

Designing multivalent immunogens for alphavirus vaccine optimization

There is a pressing need for vaccines against mosquito-borne alphaviruses such as Venezualen and eastern equine encephalitis viruses (VEEV, EEEV). We demonstrate an approach to vaccine development based on physicochemical properties (PCP) of amino acids to design a PCP-consensus sequence of the epitope-rich B domain of the VEEV major antigenic E2 protein. The consensus "spike" domain was incorporated into a live-attenuated VEEV vaccine candidate (ZPC/IRESv1). Mice inoculated with either ZPC/IRESv1 or the same virus containing the consensus E2 protein fragment (VEEVconE2) were protected against lethal challenge with VEEV strains ZPC-738 and 3908, and Mucambo virus (MUCV, related to VEEV), and had comparable neutralizing antibody titers against each virus. Both vaccines induced partial protection against Madariaga virus (MADV), a close relative of EEEV, lowering mortality from 60% to 20%. Thus PCP-consensus sequences can be integrated into a replicating virus that could, with further optimization, provide a broad-spectrum vaccine against encephalitic alphaviruses.

Producing physicochemical property consensus alphavirus protein antigens for broad spectrum vaccine design

There is a pressing need for new vaccines against alphaviruses, which can cause fatal encephalitis (Venezuelan equine encephalitis virus (VEEV) and others) and severe arthralgia (e.g. Chikungunya virus, CHIKV). These positive-strand RNA viruses are diverse and evolve rapidly, meaning that the sequence of any vaccine should cover multiple strains that may be quite different from any previous isolate. Here, consensus proteins were produced to represent the common physicochemical properties (PCPs) of the epitope rich, B domain of the E2 envelope protein. PCP-consensus proteins were based on multiple strains of VEEV (VEEVcon) and CHIKV (CHIKVcon) or the conserved PCPs of 24 different alphaviruses (AllAVcon). The AllAVcon was altered to include binding sites for neutralizing antibodies of both VEEV and CHIKV strains (Mosaikcon). All four designed proteins were produced solubly in E. coli and purified. They formed the β-strand core expected from experimental structures of this region of the wild type E2 proteins as indicated by circular dichroism (CD) spectra. Furthermore, the CHIKVcon protein bound to a structure dependent, CHIKV neutralizing monoclonal antibody. The AllAVcon and Mosaikcon proteins bound to polyclonal antibodies generated during natural infection with either VEEV or CHIKV, indicating they contained epitopes of both serotypes. The Mosaikcon antigen induced antibodies in rabbit sera that recognized both the VEEVcon and CHIKVcon spike proteins. These PCP-consensus antigens are promising starting points for novel, broad-spectrum alphavirus vaccines.

The La-related protein 1-specific domain repurposes HEAT-like repeats to directly bind a 5'TOP sequence

La-related protein 1 (LARP1) regulates the stability of many mRNAs. These include 5'TOPs, mTOR-kinase responsive mRNAs with pyrimidine-rich 5' UTRs, which encode ribosomal proteins and translation factors. We determined that the highly conserved LARP1-specific C-terminal DM15 region of human LARP1 directly binds a 5'TOP sequence. The crystal structure of this DM15 region refined to 1.86 Å resolution has three structurally related and evolutionarily conserved helix-turn-helix modules within each monomer. These motifs resemble HEAT repeats, ubiquitous helical protein-binding structures, but their sequences are inconsistent with consensus sequences of known HEAT modules, suggesting this structure has been repurposed for RNA interactions. A putative mTORC1-recognition sequence sits within a flexible loop C-terminal to these repeats. We also present modelling of pyrimidine-rich single-stranded RNA onto the highly conserved surface of the DM15 region. These studies lay the foundation necessary for proceeding toward a structural mechanism by which LARP1 links mTOR signalling to ribosome biogenesis.

Inhibition of amyloid beta-induced synaptotoxicity by a pentapeptide derived from the glycine zipper region of the neurotoxic peptide

A major characteristic of Alzheimer's disease is the presence of amyloid beta (Aβ) oligomers and aggregates in the brain. Aβ oligomers interact with the neuronal membrane inducing perforations, causing an influx of calcium ions and increasing the release of synaptic vesicles that leads to a delayed synaptic failure by vesicle depletion. Here, we identified a neuroprotective pentapeptide anti-Aβ compound having the sequence of the glycine zipper region of the C-terminal of Aβ (G33LMVG37). Docking and Förster resonance energy transfer experiments showed that G33LMVG37 interacts with Aβ at the C-terminal region, which is important for Aβ association and insertion into the lipid membrane. Furthermore, this pentapeptide interfered with Aβ aggregation, association, and perforation of the plasma membrane. The synaptotoxicity induced by Aβ after acute and chronic applications were abolished by G33LMVG37. These results provide a novel rationale for drug development against Alzheimer's disease.

Engineering proteins with enhanced mechanical stability by force-specific sequence motifs

Use of atomic force microscopy (AFM) has recently led to a better understanding of the molecular mechanisms of the unfolding process by mechanical forces; however, the rational design of novel proteins with specific mechanical strength remains challenging. We have approached this problem from a new perspective that generates linear physical-chemical properties (PCP) motifs from a limited AFM data set. Guided by our linear sequence analysis, we designed and analyzed four new mutants of the titin I1 domain with the goal of increasing the domain's mechanical strength. All four mutants could be cloned and expressed as soluble proteins. AFM data indicate that at least two of the mutants have increased molecular mechanical strength. This observation suggests that the PCP method is useful to graft sequences specific for high mechanical stability to weak proteins to increase their mechanical stability, and represents an additional tool in the design of novel proteins besides steered molecular dynamics calculations, coarse grained simulations, and ϕ-value analysis of the transition state.

Mechanical stability and differentially conserved physical-chemical properties of titin Ig-domains

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single-molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I-band of a sarcomere is composed of about 40 Ig-domains and is exposed to force under normal physiological conditions and connects the free-hanging ends of the myosin filaments to the Z-disc. Recent single-molecule force spectroscopy data show a mechanical hierarchy in the I-band domains. Domains near the C-terminus in this region unfold at forces two to three times greater than domains near the beginning of the I-band. Though all of these Ig-domains are thought to share a fold and topology common to members of the Ig-like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I-band Ig domain to specific, conserved physical-chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into "strong" and "weak" families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the beta-sandwich fold, and force sensitive residues are not only confined to the A'-G region.

Automated detection of conformational epitopes using phage display Peptide sequences

Background:

Precise determination of conformational epitopes of neutralizing antibodies represents a key step in the rational design of novel vaccines. A powerful experimental method to gain insights on the physical chemical nature of conformational epitopes is the selection of linear peptides that bind with high affinities to a monoclonal antibody of interest by phage display technology. However, the structural characterization of conformational epitopes from these mimotopes is not straightforward, and in the past the interpretation of peptide sequences from phage display experiments focused on linear sequence analysis to find a consensus sequence or common sequence motifs.

Results:

We present a fully automated search method, EpiSearch that predicts the possible location of conformational epitopes on the surface of an antigen. The algorithm uses peptide sequences from phage display experiments as input, and ranks all surface exposed patches according to the frequency distribution of similar residues in the peptides and in the patch. We have tested the performance of the EpiSearch algorithm for six experimental data sets of phage display experiments, the human epidermal growth factor receptor-2 (HER-2/neu), the antibody mAb Bo2C11 targeting the C₂ domain of FVIII, antibodies mAb 17b and mAb b12 of the HIV envelope protein gp120, mAb 13b5 targeting HIV-1 capsid protein and 80R of the SARS coronavirus spike protein. In all these examples the conformational epitopes as determined by the X-ray crystal structures of the antibody-antigen complexes, were found within the highest scoring patches of EpiSearch, covering in most cases more than 50% residues of experimental observed conformational epitopes. Input options of the program include mapping of a single peptide or a set of peptides on the antigen structure, and the results of the calculation can be visualized on our interactive web server.

Availability:

Users can access the EpiSearch from our web server http://curie.utmb.edu/episearch.html

Structural analysis of linear and conformational epitopes of allergens

In many countries regulatory agencies have adopted safety guidelines, based on bioinformatics rules from the WHO/FAO and EFSA recommendations, to prevent potentially allergenic novel foods or agricultural products from reaching consumers. We created the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to combine data that had previously been available only as flat files on Web pages or in the literature. SDAP was designed to be user friendly, to be of maximum use to regulatory agencies, clinicians, as well as to scientists interested in assessing the potential allergenic risk of a protein. We developed methods, unique to SDAP, to compare the physicochemical properties of discrete areas of allergenic proteins to known IgE epitopes. We developed a new similarity measure, the property distance (PD) value that can be used to detect related segments in allergens with clinical observed cross-reactivity. We have now expanded this work to obtain experimental validation of the PD index as a quantitative predictor of IgE cross-reactivity, by designing peptide variants with predetermined PD scores relative to known IgE epitopes. In complementary work we show how sequence motifs characteristic of allergenic proteins in protein families can be used as fingerprints for allergenicity.

Identification and analysis of conserved sequence motifs in cytochrome P450 family 2. Functional and structural role of a motif 187RFDYKD192 in CYP2B enzymes

Using a multiple alignment of 175 cytochrome P450 (CYP) family 2 sequences, 20 conserved sequence motifs (CSMs) were identified with the program PCPMer. Functional importance of the CSM in CYP2B enzymes was assessed from available data on site-directed mutants and genetic variants. These analyses suggested an important role of the CSM 8, which corresponds to(187)RFDYKD(192) in CYP2B4. Further analysis showed that residues 187, 188, 190, and 192 have a very high rank order of conservation compared with 189 and 191. Therefore, eight mutants (R187A, R187K, F188A, D189A, Y190A, K191A, D192A, and a negative control K186A) were made in an N-terminal truncated and modified form of CYP2B4 with an internal mutation, which is termed 2B4dH/H226Y. Function was examined with the substrates 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC), 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), 7-benzyloxy-4-(trifluoromethyl)coumarin (7-BFC), and testosterone and with the inhibitors 4-(4-chlorophenyl)imidazole (4-CPI) and bifonazole (BIF). Compared with the template and K186A, the mutants R187A, R187K, F188A, Y190A, and D192A showed > or =2-fold altered substrate specificity, k(cat), K(m), and/or k(cat)/K(m) for 7-MFC and 7-EFC and 3- to 6-fold decreases in differential inhibition (IC(50,BIF)/IC(50,4-CPI)). Subsequently, these mutants displayed 5-12 degrees C decreases in thermal stability (T(m)) and 2-8 degrees C decreases in catalytic tolerance to temperature (T(50)) compared with the template and K186A. Furthermore, when R187A and D192A were introduced in CYP2B1dH, the P450 expression and thermal stability were decreased. In addition, R187A showed increased activity with 7-EFC and decreased IC(50,BIF)/IC(50,4-CPI) compared with 2B1dH. Analysis of long range residue-residue interactions in the CYP2B4 crystal structures indicated strong hydrogen bonds involving Glu(149)-Asn(177)-Arg(187)-Tyr(190) and Asp(192)-Val(194), which were significantly-reduced/abolished by the Arg(187)-->Ala and Asp(192)-->Alasubstitutions, respectively.

Dynamic interaction of the measles virus hemagglutinin with its receptor signaling lymphocytic activation molecule (SLAM, CD150)

The interaction of measles virus with its receptor signaling lymphocytic activation molecule (SLAM) controls cell entry and governs tropism. We predicted potential interface areas of the measles virus attachment protein hemagglutinin to begin the investigation. We then assessed the relevance of individual amino acids located in these areas for SLAM-binding and SLAM-dependent membrane fusion, as measured by surface plasmon resonance and receptor-specific fusion assays, respectively. These studies identified one hemagglutinin protein residue, isoleucine 194, which is essential for primary binding. The crystal structure of the hemagglutinin-protein localizes Ile-194 at the interface of propeller blades 5 and 6, and our data indicate that a small aliphatic side chain of residue 194 stabilizes a protein conformation conducive to binding. In contrast, a quartet of residues previously shown to sustain SLAM-dependent fusion is not involved in binding. Instead, our data prove that after binding, this quartet of residues on propeller blade 5 conducts conformational changes that are receptor-specific. Our study sets a structure-based stage for understanding how the SLAM-elicited conformational changes travel through the H-protein ectodomain before triggering fusion protein unfolding and membrane fusion.

Co-localization of glyceraldehyde-3-phosphate dehydrogenase with ferredoxin-NADP reductase in pea leaf chloroplasts

In immunogold double-labeling of pea leaf thin sections with antibodies raised against ferredoxin-NADP reductase (EC 1.18.1.2, FNR) and antibodies directed against the A or B subunits of the NADP-linked glyceraldehyde-3-P dehydrogenase (GAPD) (EC 1.2.1.13), many small and large gold particles were found together over the chloroplasts. Nearest neighbor analysis of the distribution of the gold particles indicates that FNR and the NADP-linked GAPD are co-localized, in situ. This suggests that FNR might carry FADH2 or NADPH from the thylakoid membrane to GAPD, or that ferredoxin might carry electrons to FNR co-localized with GAPD in the stroma. Crystal structures of the spinach enzymes are available. When they are docked computationally, the proteins appear, as modeled, to be able to form at least two different complexes. One involves a single GAPD monomer and an FNR monomer (or dimer). The amino acid residues located at the putative interface are highly conserved on the chloroplastic forms of both enzymes. The other potential complex involves the GAPD A2B2 tetramer and an FNR monomer (or dimer). The interface residues are conserved in this model as well. Ferredoxin is able to interact with FNR in either complex.

Statistical analysis of physical-chemical properties and prediction of protein-protein interfaces

We have developed a fully automated method, InterProSurf, to predict interacting amino acid residues on protein surfaces of monomeric 3D structures. Potential interacting residues are predicted based on solvent accessible surface areas, a new scale for interface propensities, and a cluster algorithm to locate surface exposed areas with high interface propensities. Previous studies have shown the importance of hydrophobic residues and specific charge distribution as characteristics for interfaces. Here we show differences in interface and surface regions of all physical chemical properties of residues as represented by five quantitative descriptors. In the current study a set of 72 protein complexes with known 3D structures were analyzed to obtain interface propensities of residues, and to find differences in the distribution of five quantitative descriptors for amino acid residues. We also investigated spatial pair correlations of solvent accessible residues in interface and surface areas, and compared log-odds ratios for interface and surface areas. A new scoring method to predict potential functional sites on the protein surface was developed and tested for a new dataset of 21 protein complexes, which were not included in the original training dataset. Empirically we found that the algorithm achieves a good balance in the accuracy of precision and sensitivity by selecting the top eight highest scoring clusters as interface regions. The performance of the method is illustrated for a dimeric ATPase of the hyperthermophile, Methanococcus jannaschii, and the capsid protein of Human Hepatitis B virus. An automated version of the method can be accessed from our web server at http://curie.utmb.edu/prosurf.html.

InterProSurf: a web server for predicting interacting sites on protein surfaces

A new web server, InterProSurf, predicts interacting amino acid residues in proteins that are most likely to interact with other proteins, given the 3D structures of subunits of a protein complex. The prediction method is based on solvent accessible surface area of residues in the isolated subunits, a propensity scale for interface residues and a clustering algorithm to identify surface regions with residues of high interface propensities. Here we illustrate the application of InterProSurf to determine which areas of Bacillus anthracis toxins and measles virus hemagglutinin protein interact with their respective cell surface receptors. The computationally predicted regions overlap with those regions previously identified as interface regions by sequence analysis and mutagenesis experiments.

AVAILABILITY

The InterProSurf web server is available at http://curie.utmb.edu/

Flavitrack: an annotated database of flavivirus sequences

MOTIVATION

Properly annotated sequence data for flaviviruses, which cause diseases, such as tick-borne encephalitis (TBE), dengue fever (DF), West Nile (WN) and yellow fever (YF), can aid in the design of antiviral drugs and vaccines to prevent their spread. Flavitrack was designed to help identify conserved sequence motifs, interpret mutational and structural data and track evolution of phenotypic properties.

SUMMARY

Flavitrack contains over 590 complete flavivirus genome/protein sequences and information on known mutations and literature references. Each sequence has been manually annotated according to its date and place of isolation, phenotype and lethality. Internal tools are provided to rapidly determine relationships between viruses in Flavitrack and sequences provided by the user.