Computational Analysis of the Interaction Energies between Amino Acid Residues of the Measles Virus Hemagglutinin and Its Receptors (original) (raw)
Related papers
Structure of the measles virus hemagglutinin bound to the CD46 receptor
Nature Structural & Molecular Biology, 2009
The binding site in MV-H for nectin-4 also overlaps extensively with those of the other two receptors. Finally, a hydrophobic pocket centered in the b4-b5 groove is involved in binding to all three identified measles virus receptors, representing a potential target for antiviral drugs.
Protein Science, 2008
The two terminal complement control protein (CCP) modules of the CD46 glycoprotein mediate measles virus binding. Three-dimensional models for these two domains were derived based on the NMR structures of two CCP modules of factor H. Both CD46 CCP modules are about 35 8, long, and form a five-stranded antiparallel P-barrel structure. Monte Carlo simulations, sampling the backbone torsion angles of the linker peptide and selecting possible orientations on the basis of minimal solvent-exposed hydrophobic area, were used to predict the orientation of CCP-I relative to CCP-11. We tested this procedure successfully for factor H. For CD46, three clusters of structures differing in the tilt angle of the two domains were obtained. To test these models, we mutagenized the CCP modules. Four proteins, two without an oligosaccharide chain and two with mutated short amino acid segments, reached the cell surface efficiently. Only the protein without the CCP-I oligosaccharide chain maintained binding to the viral attachment protein hemagglutinin. These results are consistent with one of our models and suggest that the viral hemagglutinin does not bind at the membrane-distal tip of CD46, but near the concave CCP-1-11 interface region.
… and Crystallization …, 2009
The measles virus (MV) hemagglutinin (MV-H) mediates the attachment of MV particles to cell-surface receptors for entry into host cells. MV uses two receptors for attachment to host cells: the complement-control protein CD46 and the signalling lymphocyte activation molecule (SLAM). The MV-H glycoprotein from an Edmonston MV variant and the MV-binding fragment of the CD46 receptor were overproduced in mammalian cells and used to crystallize an MV-H–CD46 complex. Well diffracting crystals containing two complexes in the asymmetric unit were obtained and the structure of the complex was solved by the molecular-replacement method.
Crystal structure of measles virus hemagglutinin provides insight into effective vaccines
Proceedings of the National Academy of Sciences, 2007
Measles still remains a major cause of childhood morbidity and mortality worldwide. Measles virus (MV) vaccines are highly successful, but the mechanism underlying their efficacy has been unclear. Here we report the crystal structure of the MV attachment protein, hemagglutinin, responsible for MV entry. The receptorbinding head domain exhibits a cubic-shaped -propeller structure and forms a homodimer. N-linked sugars appear to mask the broad regions and cause the two molecules forming the dimer to tilt oppositely toward the horizontal plane. Accordingly, residues of the putative receptor-binding site, highly conserved among MV strains, are strategically positioned in the unshielded area of the protein. These conserved residues also serve as epitopes for neutralizing antibodies, ensuring the serological monotype, a basis for effective MV vaccines. Our findings suggest that sugar moieties in the MV hemagglutinin critically modulate virus-receptor interaction as well as antiviral antibody responses, differently from sugars of the HIV gp120, which allow for immune evasion.
International Journal of Molecular Sciences, 2011
Antibodies play an increasingly important role in both basic research and the pharmaceutical industry. Since their efficiency depends, in ultimate analysis, on their atomic interactions with an antigen, studying such interactions is important to understand how they function and, in the long run, to design new molecules with desired properties. Computational docking, the process of predicting the conformation of a complex from its separated components, is emerging as a fast and affordable technique for the structural characterization of antibody-antigen complexes. In this manuscript, we first describe the different computational strategies for the modeling of antibodies and docking of their complexes, and then predict the binding of two antibodies to the stalk region of influenza hemagglutinin, an important pharmaceutical target. The purpose is two-fold: on a general note, we want to illustrate the advantages and pitfalls of computational docking with a practical example, using different approaches and comparing the results to known experimental structures. On a more specific note, we want to assess if docking can be successful in characterizing the binding to the same influenza epitope of other antibodies with unknown structure, which has practical relevance for pharmaceutical and biological research. The paper clearly shows that some of the computational docking predictions can be very accurate, but the algorithm often fails to discriminate them from inaccurate OPEN ACCESS Int. J. Mol. Sci. 2011, 12 227 solutions. It is of paramount importance, therefore, to use rapidly obtained experimental data to validate the computational results.
Molecular Diversity, 2020
The pandemic outbreak of the Corona viral infection has become a critical global health issue. Biophysical and structural evidence shows that spike protein possesses a high binding affinity towards host angiotensin-converting enzyme 2 and viral hemagglutinin-acetylesterase (HE) glycoprotein receptor. We selected HE as a target in this study to identify potential inhibitors using a combination of various computational approaches such as molecular docking, ADMET analysis, dynamics simulations and binding free energy calculations. Virtual screening of NPACT compounds identified 3,4,5-Trihydroxy-1,8-bis[(2R,3R)-3,5,7-trihydroxy-3,4-dihydro-2H-chromen-2-yl]benzo[7]annulen-6-one, Silymarin, Withanolide D, Spirosolane and Oridonin as potential HE inhibitors with better binding energy. Furthermore, molecular dynamics simulations for 100 ns time scale revealed that most of the key HE contacts were retained throughout the simulations trajectories. Binding free energy calculations using MM/PBSA approach ranked the top-five potential NPACT compounds which can act as effective HE inhibitors.
Journal of Molecular Modeling
Molecular dynamics (MD) simulations were carried out to study the behavior of human receptor molecule in the hemagglutinin (HA) of 1918 and 2009 H1N1 influenza viruses respectively. The 2009 HA model was obtained by virtually mutating the 1918 HA crystal structure based on A/Mexico City/MCIG01/2009(H1N1) segment 4 sequence. We found that human receptor molecule has no binding preference between the 2009 HA and the 1918 HA. In addition, among the four sugar moieties in the human receptor molecule, sialic acid contributes the most to the electrostatic and non-polar interaction energy during binding. Furthermore, the hydrogen bonds between sialic acid and the surrounding residues in 1918 HA are preserved in 2009 HA. We also found that the mutated residues contribute to a more favorable binding of hemagglutinin to the human receptor molecule.
Computational studies of H5N1 hemagglutinin binding with SA-a-2, 3-Gal and SA-a-2, 6-Gal
Biochem Biophys Res Commun, 2006
For influenza H5N1 hemagglutinin, a switch from SA-a-2, 3-Gal to SA-a-2, 6-Gal receptor specificity is a critical step leading to the conversion from avian-to-human to human-to-human infection. Therefore, the understanding of the binding modes of SA-a-2, 3-Gal and SA-a-2, 6-Gal to H5N1 hemagglutinin will be very important for the examination of possible mutations needed for going from an avian to a human flu virus. Based on the available H5N1 hemagglutinin crystal structure, the binding profiles between H5N1 hemagglutinin and two saccharide ligands, SA-a-2, 3-Gal and SA-a-2, 6-Gal, were investigated by ab initio quantum mechanics, molecular docking, molecular mechanics, and molecular dynamics simulations. It was found that SA-a-2, 3-Gal has strong multiple hydrophobic and hydrogen bond interactions in its trans conformation with H5N1 hemagglutinin, whereas the SA-a-2, 6-Gal only shows weak interactions in a different conformation (cis type).
Journal of Biophysical Chemistry, 2011
In this present study, we predicted the neutralizing epitope of a modeled H5N1 hemagglutinin 1046T when interacted with a modeled monoclonal antibody variable fragment 8H5Fv using molecular dynamics simulation. Following the production run of the molecular dynamics simulation, we observed the average change of solvent accessible surface of the antigen alongside the formation of hydrogen bonds between the two structures during the simulation. Based on the acquired data, we predicted the neutralizing epitope of the 1046T antigen to be consisted of residues Asp 84, Glu85, Phe86, Ile87, Asn88, Val89, Pro90, Ile132, Ser136, Val147, Pro152, Tyr153, Leu154, Arg161, and Tyr268. By calculating the RMSD of the Cα backbone chain of the complex during the simulation we found the structure to be generally stable suggesting a well maintained steric hindrance, while RMSD calculation of the predicted neutralizing epitope backbone suggests the stability of the neutralizing epitope itself.
Computational Biology and Chemistry, 2023
Intermolecular interaction between key residue N501 of the epitope on SARS-CoV-2 RBD and screening antibody B38 was studied using the QM/MM and QM approach. The QM/MM optimized geometry shows that angle X-H---Y is 165° for O-H---O between mAb light chain S30 and RBD N501. High level MP2 calculations indicated the interaction between RBD N501 and S30 of B38 Fab light chain provide a relatively strong attractive force of − 3.32 kcal/mol, whereas the hydrogen bond between RBD Q498 and S30 was quantified as 0.10 kcal/mol. The decrease in ESP partial charge on hydrogen atom of hydroxyl group on S30 drops from 0.38 a.u. to 0.31 a.u., exhibiting the sharing of 0.07 a.u. from the lone pair electron oxygen of N501 due to hydrogen bond formation. The NBO occupancy of hydrogen atom also decreases from 25.79 % to 22.93 % in the hydroxyl H-O NBO bond of S30. However, the minor change of NBO hybridization of hydroxyl oxygen of S30 from sp3.00 to sp3.05 implies the rigidity of hydrogen bond tetrahedral geometry in the relative dynamic protein complex. The O-H---O angle is 165° which is close but not exactly linear. The structural requirement for sp3 hybridization of oxygen for hydroxyl group on S30 and dimension of protein likely prevent O-H---O from adopting linear geometry. The hydrogen bond strengths were also calculated using a variety of DFT methods, and the result of − 3.33 kcal/mol from the M06L method is the closest to that of the MP2 calculation. Results of this work may aid in the COVID-19 vaccine and drug screening.Open in a separate window