Probing the “fingers” domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking,molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π?π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN?Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50?ns MD simulation, the π?π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.

Communicated by Ramaswamy H. Sarma     

Investigation of intermolecular interactions and stability of verubecestat in the active site of BACE1: Development of first model from QM/MM-based charge density and MD analysis

Alzheimer disease (AD) is a cruel neurodegenerative disorder caused by the deposition of amyloid β (Aβ) peptide inside the brain. The β-secretase (beta amyloid precursor protein (APP) cleaving enzyme 1, BACE1) is one of the enzymes involved in the cleavage of APP that leads to the Aβ formation and it is the primary target for the treatment of AD. Recent report outlines that verubecestat molecule strongly inhibits BACE1; however, its structure, binding mechanism and the stability in the active site of BACE1 are not yet known. The present study aims to determine the structure, binding affinity and the stability of verubecestat molecule in the active site of BACE1 from the molecular docking, quantum mechanics/molecular mechanics (QM/MM)-based charge density analysis and molecular dynamics simulation. Verubecestat molecule was docked at BACE1; it shows high binding affinity towards BACE1. Further, the conformational geometry and the intermolecular interactions of verubecestat in the active site of BACE1 were determined. The molecule forms strong interaction with the neighboring amino acids in the active site of BACE1. The onsite QM/MM-based charge density analysis reveals the nature of charge density distribution and the topological properties of intermolecular interactions of verubecestat molecule in the active site of BACE1. The calculated electrostatic potential (ESP) of verubecestat in the active site of BACE1 displays high negative and positive ESP regions of the molecule. This onsite QM/MM analysis is more relevant to the physiological situation. The molecular dynamics simulation has been performed, which confirms the high stability and compactness of verubecestat in the active site of BACE1. The MM-generalized Born surface area and MM-Poisson Boltzmann surface area free energy calculations of verubecestat–BACE1 also confirm the high binding affinity of verubecestat.

Communicated by Ramaswamy H. Sarma     

Mutualistic Biofilm Communities Develop with Porphyromonas gingivalis and Initial,Early, and Late Colonizers of Enamel

Porphyromonas gingivalis is present in dental plaque as early as 4 h after tooth cleaning, but it is also associated with periodontal disease, a late-developing event in the microbial successions that characterize daily plaque development. We report here that P. gingivalis ATCC 33277 is remarkable in its ability to interact with a variety of initial, early, middle, and late colonizers growing solely on saliva. Integration of P. gingivalis into multispecies communities was investigated by using two in vitro biofilm models. In flow cells, bacterial growth was quantified using fluorescently conjugated antibodies against each species, and static biofilm growth on saliva-submerged polystyrene pegs was analyzed by quantitative real-time PCR using species-specific primers. P. gingivalis could not grow as a single species or together with initial colonizer Streptococcus oralis but showed mutualistic growth when paired with two other initial colonizers, Streptococcus gordonii and Actinomyces oris, as well as with Veillonella sp. (early colonizer), Fusobacterium nucleatum (middle colonizer), and Aggregatibacter actinomycetemcomitans (late colonizer). In three-species flow cells, P. gingivalis grew with Veillonella sp. and A. actinomycetemcomitans but not with S. oralis and A. actinomycetemcomitans. Also, it grew with Veillonella sp. and F. nucleatum but not with S. oralis and F. nucleatum, indicating that P. gingivalis and S. oralis are not compatible. However, P. gingivalis grew in combination with S. gordonii and S. oralis, demonstrating its ability to overcome the incompatibility when cultured with a second initially colonizing species. Collectively, these data help explain the observed presence of P. gingivalis at all stages of dental plaque development.Removal of dental plaque by routine oral hygiene procedures is followed by a repetition of a species succession that starts with initially colonizing streptococci and actinomyces (5, 16). Other species follow as early, middle, and late colonizers, which establishes the following developmental process: successive attachment of saliva-suspended species to already attached bacteria and formation of multispecies communities.Attachment is a critical event essential to preventing the bacteria from being swallowed by salivary flow. Initial colonizers bind to host-derived receptors in the salivary pellicle coating of the tooth enamel. The remainder of typical plaque development occurs by accretion of saliva-suspended species and growth of attached bacteria, thereby increasing the microbial diversity. Adherence of suspended single cells to attached cells is called coadhesion (1). Some suspended cells are already coaggregated and adhere to attached cells as coaggregates; coaggregation is defined as the specific cell-to-cell recognition and adherence of genetically distinct cell types (8). All human oral bacterial species exhibit coaggregation. For example, Streptococcus oralis coaggregates with Streptococcus gordonii (intrageneric coaggregation). Both species pair with Actinomyces oris (intergeneric coaggregation), and all of them coaggregate with Fusobacterium nucleatum (multigeneric coaggregation). Multispecies communities composed of coaggregating species characterize dental plaque biofilms in vivo (3, 17, 18).To increase our understanding of interactions among species, we have employed two in vitro model systems and are testing numerous combinations of seven species for their ability to grow on saliva as their sole nutritional source (20, 21). First, we reported that F. nucleatum (middle colonizer) failed to grow when paired with S. oralis but grew well when A. oris was included in the three-species biofilm (20), indicating specificity by F. nucleatum for the presence of a particular initial colonizer. Recently, we showed that Aggregatibacter actinomycetemcomitans (late colonizer and periodontopathogen) exhibited mutualistic relationships with F. nucleatum and Veillonella sp. (early colonizer and commensal organism), illustrating the ability of commensals and pathogens to grow together (21).Porphyromonas gingivalis, another periodontopathogen, forms three-species communities with F. nucleatum and S. gordonii (11). Proteomics of P. gingivalis in this three-species community revealed a broad increase in proteins involved in protein synthesis, suggesting that a multispecies relationship is advantageous for the porphyromonad (11). This research group had previously reported the presence of differentially regulated porphyromonad genes when P. gingivalis and S. gordonii were together in biofilms (22). Thus, P. gingivalis responds to the presence of other oral species.P. gingivalis is detected in dental plaque samples within 6 h after professional tooth cleaning (5, 13), and its numbers increase in periodontally diseased sites (15). It forms biofilms with S. gordonii but not with Streptococcus mutans (12) or Streptococcus cristatus (23). P. gingivalis required a preformed streptococcal substratum for its incorporation into a biofilm (12). Partner specificity was also noted among four fresh isolates of P. gingivalis, which showed no coaggregation with a variety of oral actinomyces, aggregatibacteria, capnocytophagae, and streptococci (9) but coaggregated with F. nucleatum (7, 10). We show here that P. gingivalis exhibits widespread mutualism with initial, early, middle, and late colonizers but also shows specificity with initially colonizing streptococci, which could help explain its early appearance in the development of dental plaque biofilms. The relationship of porphyromonads with initial, early, middle, and later colonizers during biofilm growth on saliva as a sole nutritional source has not been explored previously. We hypothesize that the ability of P. gingivalis to coaggregate with S. gordonii and A. oris (initial colonizers), Veillonella sp. (early colonizer), F. nucleatum (middle colonizer), and A. actinomycetemcomitans (late colonizer) allows these bacteria to form multispecies biofilm communities.     

A novel 9-bp insertion detected in steroid 21-hydroxylase gene (CYP21A2): prediction of its structural and functional implications by computational methods

Background

Steroid 21-hydroxylase deficiency is the most common cause of congenital adrenal hyperplasia (CAH). Detection of underlying mutations in CYP21A2 gene encoding steroid 21-hydroxylase enzyme is helpful both for confirmation of diagnosis and management of CAH patients. Here we report a novel 9-bp insertion in CYP21A2 gene and its structural and functional consequences on P450c21 protein by molecular modeling and molecular dynamics simulations methods.

Methods

A 30-day-old child was referred to our laboratory for molecular diagnosis of CAH. Sequencing of the entire CYP21A2 gene revealed a novel insertion (duplication) of 9-bp in exon 2 of one allele and a well-known mutation I172N in exon 4 of other allele. Molecular modeling and simulation studies were carried out to understand the plausible structural and functional implications caused by the novel mutation.

Results

Insertion of the nine bases in exon 2 resulted in addition of three valine residues at codon 71 of the P450c21 protein. Molecular dynamics simulations revealed that the mutant exhibits a faster unfolding kinetics and an overall destabilization of the structure due to the triple valine insertion was also observed.

Conclusion

The novel 9-bp insertion in exon 2 of CYP21A2 genesignificantly lowers the structural stability of P450c21 thereby leading to the probable loss of its function.     

A Missense Mutation in the Arabidopsis COPII Coat Protein Sec24A Induces the Formation of Clusters of the Endoplasmic Reticulum and Golgi Apparatus

How the endoplasmic reticulum (ER) and the Golgi apparatus maintain their morphological and functional identity while working in concert to ensure the production of biomolecules necessary for the cell''s survival is a fundamental question in plant biology. Here, we isolated and characterized an Arabidopsis thaliana mutant that partially accumulates Golgi membrane markers and a soluble secretory marker in globular structures composed of a mass of convoluted ER tubules that maintain a connection with the bulk ER. We established that the aberrant phenotype was due to a missense recessive mutation in sec24A, one of the three Arabidopsis isoforms encoding the coat protomer complex II (COPII) protein Sec24, and that the mutation affects the distribution of this critical component at ER export sites. By contrast, total loss of sec24A function was lethal, suggesting that Arabidopsis sec24A is an essential gene. These results produce important insights into the functional diversification of plant COPII coat components and the role of these proteins in maintaining the dynamic identity of organelles of the early plant secretory pathway.     

Micelle-bound structures and dynamics of the hinge deleted analog of melittin and its diastereomer: Implications in cell selective lysis by d-amino acid containing antimicrobial peptides

Melittin, the major component of the honey bee venom, is a 26-residue hemolytic and membrane active peptide. Structures of melittin determined either in lipid environments by NMR or by use of X-ray demonstrated two helical regions at the N- and C-termini connected by a hinge or a bend at the middle. Here, we show that deletion of the hinge residues along with two C-terminal terminal Gln residues (Q25 and Q26), yielding a peptide analog of 19-residue or Mel-H, did not affect antibacterial activity but resulted in a somewhat reduced hemolytic activity. A diastereomer of Mel-H or Mel-^dH containing d-amino acids [^dV5, ^dV8, ^dL11 and ^dK16] showed further reduction in hemolytic activity without lowering antibacterial activity. We have carried out NMR structures, dynamics (H-D exchange and proton relaxation), membrane localization by spin labeled lipids, pulse-field-gradient (PFG) NMR and isothermal titration calorimetry (ITC) in dodecylphosphocholine (DPC) micelles, as a mimic to eukaryotic membrane, to gain insights into cell selectivity of these melittin analogs. PFG-NMR showed Mel-H and Mel-^dH both were similarly partitioned into DPC micelles. ITC demonstrated that Mel-H and Mel-^dH interact with DPC with similar affinity. The micelle-bound structure of Mel-H delineated a straight helical conformation, whereas Mel-^dH showed multiple β-turns at the N-terminus and a short helix at the C-terminus. The backbone amide-proton exchange with solvent D₂O demonstrated a large difference in dynamics between Mel-H and Mel-^dH, whereby almost all backbone protons of Mel-^dH showed a much faster rate of exchange as compared to Mel-H. Proton T₁ relaxation had suggested a mobile backbone of Mel-^dH peptide in DPC micelles. Resonance perturbation by paramagnetic lipids indicated that Mel-H inserted deeper into DPC micelles, whereas Mel-^dH is largely located at the surface of the micelle. Taken together, results presented in this study demonstrated that the poor hemolytic activity of the d-amino acid containing analogs of antimicrobial peptides may be correlated with their flexible dynamics at the membrane surface.