Origin of Decrease in Potency of Darunavir and Two Related Antiviral Inhibitors against HIV-2 Compared to HIV-1 Protease

Elucidating the conformational dynamics of histo-blood group antigens and their interactions with the rotavirus spike protein through computational lens

In the present study, we investigated the conformational dynamics of histo-blood group antigens (HBGAs) and their interactions with the VP8* domain of four rotavirus genotypes (P[4], P[6], P[19], and P[11]) utilizing all-atom molecular dynamics simulations in explicit water. Our study revealed distinct changes in the dynamic behavior of the same glycan due to linkage variations. We observed that LNFPI HBGA having a terminal β linkage shows two dominant conformations after complexation, whereas only one was obtained for LNFPI with a terminal α linkage. Interestingly, both variants displayed a single dominant structure in the free state. Similarly, LNT and LNnT show a shift in their dihedral linkage profile between their two terminal monosaccharides because of a change in the linkage from β(1-3) to β(1-4). The molecular mechanics generalized Born surface area (MM/GBSA) calculations yielded the highest binding affinity for LNFPI(β)/P[6] (-13.93 kcal/mol) due to the formation of numerous hydrogen bonds between VP8* and HBGAs. LNnT binds more strongly to P[11] (-12.88 kcal/mol) than LNT (-4.41 kcal/mol), suggesting a single change in the glycan linkage might impact its binding profile significantly. We have also identified critical amino acids and monosaccharides (Gal and GlcNAc) that contributed significantly to the protein-ligand binding through the per-residue decomposition of binding free energy. Moreover, we found that the interaction between the same glycan and different protein receptors within the same rotavirus genogroup influenced the micro-level dynamics of the glycan. Overall, our study helps a deeper understanding of the H-type HBGA and rotavirus spike protein interaction.Communicated by Ramaswamy H. Sarma.

Synthetic ERRα/β/γ Agonist Induces an ERRα-Dependent Acute Aerobic Exercise Response and Enhances Exercise Capacity

Repetitive physical exercise induces physiological adaptations in skeletal muscle that improves exercise performance and is effective for the prevention and treatment of several diseases. Genetic evidence indicates that the orphan nuclear receptors estrogen receptor-related receptors (ERRs) play an important role in skeletal muscle exercise capacity. Three ERR subtypes exist (ERRα, β, and γ), and although ERRβ/γ agonists have been designed, there have been significant difficulties in designing compounds with ERRα agonist activity. Additionally, there are limited synthetic agonists that can be used to target ERRs in vivo. Here, we report the identification of a synthetic ERR pan agonist, SLU-PP-332, that targets all three ERRs but has the highest potency for ERRα. Additionally, SLU-PP-332 has sufficient pharmacokinetic properties to be used as an in vivo chemical tool. SLU-PP-332 increases mitochondrial function and cellular respiration in a skeletal muscle cell line. When administered to mice, SLU-PP-332 increased the type IIa oxidative skeletal muscle fibers and enhanced exercise endurance. We also observed that SLU-PP-332 induced an ERRα-specific acute aerobic exercise genetic program, and the ERRα activation was critical for enhancing exercise endurance in mice. These data indicate the feasibility of targeting ERRα for the development of compounds that act as exercise mimetics that may be effective in the treatment of numerous metabolic disorders and to improve muscle function in the aging.

Finding potential inhibitors against RNA-dependent RNA polymerase (RdRp) of bovine ephemeral fever virus (BEFV): an in-silico study

The bovine ephemeral fever virus (BEFV) is an enzootic agent that affects millions of bovines and causes major economic losses. Though the virus is seasonally reported with a very high morbidity rate (80-100%) from African, Australian, and Asiatic continents, it remains a neglected pathogen in many of its endemic areas, with no proper therapeutic drugs or vaccines presently available for treatment. The RNA-dependent RNA polymerase (RdRp) catalyzes the viral RNA synthesis and is an appropriate candidate for antiviral drug developments. We utilized integrated computational tools to build the 3D model of BEFV-RdRp and then predicted its probable active binding sites. The virtual screening and optimization against these active sites, using several small-molecule inhibitors from a different category of Life Chemical database and FDA-approved drugs from the ZINC database, was performed. We found nine molecules that have docking scores varying between -6.84 to -10.43 kcal/mol. Furthermore, these complexes were analyzed for their conformational dynamics and thermodynamic stability using molecular dynamics simulations in conjunction with the molecular mechanics generalized Born surface area (MM-GBSA) scheme. The binding free energy calculations depict that the electrostatic interactions play a dominant role in the RdRp-inhibitor binding. The hot spot residues, such as Arg565, Asp631, Glu633, Asp740, and Glu707, were found to control the RdRp-inhibitor interaction. The ADMET analysis strongly suggests favorable pharmacokinetics of these compounds that may prove useful for treating the BEFV ailment. Overall, we anticipate that these findings would help explore and develop a wide range of anti-BEFV therapy.Communicated by Ramaswamy H. Sarma.

Decoding the Host-Parasite Protein Interactions Involved in Cerebral Malaria Through Glares of Molecular Dynamics Simulations

Malaria causes millions of deaths every year. The malaria parasite spends a substantial part of its life cycle inside human erythrocytes. Inside erythrocytes, it synthesizes and displays various proteins onto the erythrocyte surface, such as Plasmodium falciparum erythrocytic membrane protein-1 (PfEMP1). This protein contains cysteine-rich interdomain region (CIDR) domains which have many subtypes based on sequence diversity and can cross-talk with host molecules. The CIDRα1.4 subtype can attach host endothelial protein C receptor (EPCR). This interaction facilitates infected erythrocyte adherence to brain endothelium and subsequent development of cerebral malaria. Through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, we explored the mechanism of interaction in the CIDRα1-EPCR complex. We examined the structural behavior of two CIDRα1 molecules (encoded by HB3-isolate var03-gene and IT4-isolate var07-gene) with EPCR unbound and bound (complex) forms. ^HB3var03CIDRα1 in apo and complexed with EPCR was comparatively more stable than ^IT4var07CIDRα1. Both of the complexes adopted two distinct conformational energy states. The hydrophobic residues played a crucial role in the binding of both complexes. For ^HB3var03CIDRα1-EPCR, the dominant energetic components were total polar interactions, while in ^IT4var07CIDRα1-EPCR, the primary interaction was van der Waals and nonpolar solvation energy. The study also revealed details such as correlated conformational motions and secondary structure evolution. Further, it elucidated various hotspot residues involved in protein-protein recognition. Overall, our study provides additional information on the structural behavior of CIDR molecules in unbound and receptor-bound states, which will help to design potent inhibitors.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study

A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.

Identification of Food Compounds as inhibitors of SARS-CoV-2 main protease using molecular docking and molecular dynamics simulations

SARS-CoV-2 has rapidly emerged as a global pandemic with high infection rate. At present, there is no drug available for this deadly disease. Recently, M^pro (Main Protease) enzyme has been identified as essential proteins for the survival of this virus. In the present work, Lipinski's rules and molecular docking have been performed to identify plausible inhibitors of M^pro using food compounds. For virtual screening, a database of food compounds was downloaded and then filtered using Lipinski's rule of five. Then, molecular docking was accomplished to identify hits using M^pro protein as the target enzyme. This led to identification of a Spermidine derivative as a hit. In the next step, Spermidine derivatives were collected from PubMed and screened for their binding with M^pro protein. In addition, molecular dynamic simulations (200 ns) were executed to get additional information. Some of the compounds are found to have strong affinity for M^pro, therefore these hits could be used to develop a therapeutic agent for SARS-CoV-2.

Elucidating biophysical basis of binding of inhibitors to SARS-CoV-2 main protease by using molecular dynamics simulations and free energy calculations

The recent outbreak of novel "coronavirus disease 2019" (COVID-19) has spread rapidly worldwide, causing a global pandemic. In the present work, we have elucidated the mechanism of binding of two inhibitors, namely α-ketoamide and Z31792168, to SARS-CoV-2 main protease (Mpro or 3CLpro) by using all-atom molecular dynamics simulations and free energy calculations. We calculated the total binding free energy (ΔGbind) of both inhibitors and further decomposed ΔGbind into various forces governing the complex formation using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method. Our calculations reveal that α-ketoamide is more potent (ΔGbind= - 9.05 kcal/mol) compared to Z31792168 (ΔGbind= - 3.25 kcal/mol) against COVID-19 3CLpro. The increase in ΔGbind for α-ketoamide relative to Z31792168 arises due to an increase in the favorable electrostatic and van der Waals interactions between the inhibitor and 3CLpro. Further, we have identified important residues controlling the 3CLpro-ligand binding from per-residue based decomposition of the binding free energy. Finally, we have compared ΔGbind of these two inhibitors with the anti-HIV retroviral drugs, such as lopinavir and darunavir. It is observed that α-ketoamide is more potent compared to lopinavir and darunavir. In the case of lopinavir, a decrease in van der Waals interactions is responsible for the lower binding affinity compared to α-ketoamide. On the other hand, in the case of darunavir, a decrease in the favorable intermolecular electrostatic and van der Waals interactions contributes to lower affinity compared to α-ketoamide. Our study might help in designing rational anti-coronaviral drugs targeting the SARS-CoV-2 main protease.Communicated by Ramaswamy H. Sarma.

Structure, Synthesis and Inhibition Mechanism of Nucleoside Analogues as HIV-1 Reverse Transcriptase Inhibitors (NRTIs)

Acquired immunodeficiency syndrome (AIDS) is caused by infection with the human immunodeficiency virus (HIV). Although treatments against HIV infection are available, AIDS remains a serious disease that causes many deaths annually. Although a variety of anti-HIV drugs have been synthesized and marketed to treat HIV-infected patients, nucleoside analogue reverse transcriptase inhibitors (NRTIs), which mimic nucleosides, are used extensively and remain a subject of interest to medicinal chemists. However, HIV has acquired drug resistance against NRTIs, and thus the struggle to find novel therapies continues. In this review, we trace the trajectory of NRTIs, focusing on the synthesis, mechanisms of action and applications of NRTIs that have been developed.

Structural Impacts of Drug-Resistance Mutations Appearing in HIV-2 Protease

The use of antiretroviral drugs is accompanied by the emergence of HIV-2 resistances. Thus, it is important to elucidate the mechanisms of resistance to antiretroviral drugs. Here, we propose a structural analysis of 31 drug-resistant mutants of HIV-2 protease (PR2) that is an important target against HIV-2 infection. First, we modeled the structures of each mutant. We then located structural shifts putatively induced by mutations. Finally, we compared wild-type and mutant inhibitor-binding pockets and interfaces to explore the impacts of these induced structural deformations on these two regions. Our results showed that one mutation could induce large structural rearrangements in side-chain and backbone atoms of mutated residue, in its vicinity or further. Structural deformations observed in side-chain atoms are frequent and of greater magnitude, that confirms that to fight drug resistance, interactions with backbone atoms should be favored. We showed that these observed structural deformations modify the conformation, volume, and hydrophobicity of the binding pocket and the composition and size of the PR2 interface. These results suggest that resistance mutations could alter ligand binding by modifying pocket properties and PR2 stability by impacting its interface. Our results reinforce the understanding of the effects of mutations that occurred in PR2 and the different mechanisms of PR2 resistance.

Mining of Ebola virus genome for the construction of multi-epitope vaccine to combat its infection

Ebola virus is the primary causative agent of viral hemorrhagic fever that is an epidemic disease and responsible for the massive premature deaths in humans. Despite knowing the molecular mechanism of its pathogenesis, to date, no commercial or FDA approved multiepitope vaccine is available against Ebola infection. The current study focuses on designing a multi-epitope subunit vaccine for Ebola using a novel immunoinformatic approach. The best predicted antigenic epitopes of Cytotoxic-T cell (CTL), Helper-T cells (HTL), and B-cell epitopes (BCL) joined by various linkers were selected for the multi-epitope vaccine designing. For the enhanced immune response, two adjuvants were also added to the construct. Further analysis showed the vaccine to be immunogenic and non-allergenic, forming a stable and energetically favorable structure. The stability of the unbound vaccine construct and vaccine/TLR4 was elucidated via atomistic molecular dynamics simulations. The binding free energy analysis (ΔG_Bind = -194.2 ± 0.5 kcal/mol) via the molecular mechanics Poisson-Boltzmann docking scheme revealed a strong association and thus can initiate the maximal immune response. Next, for the optimal expression of the vaccine construct, its gene construct was cloned in the pET28a + vector system. In summary, the Ebola viral proteome was screened to identify the most potential HTLs, CTLs, and BCL epitopes. Along with various linkers and adjuvants, a multi-epitope vaccine is constructed that showed a high binding affinity with the immune receptor, TLR4. Thus, the current study provides a highly immunogenic multi-epitope subunit vaccine construct that may induce humoral and cellular immune responses against the Ebola infection.Communicated by Ramaswamy H. Sarma.

Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations

The binding of an allosteric inhibitor in WNK1 leads to the inactive state.

Computational Investigation of Structural Dynamics of SARS-CoV-2 Methyltransferase-Stimulatory Factor Heterodimer nsp16/nsp10 Bound to the Cofactor SAM

Recently, a highly contagious novel coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has emerged, posing a global threat to public health. Identifying a potential target and developing vaccines or antiviral drugs is an urgent demand in the absence of approved therapeutic agents. The 5'-capping mechanism of eukaryotic mRNA and some viruses such as coronaviruses (CoVs) are essential for maintaining the RNA stability and protein translation in the virus. SARS-CoV-2 encodes S-adenosyl-L-methionine (SAM) dependent methyltransferase (MTase) enzyme characterized by nsp16 (2'-O-MTase) for generating the capped structure. The present study highlights the binding mechanism of nsp16 and nsp10 to identify the role of nsp10 in MTase activity. Furthermore, we investigated the conformational dynamics and energetics behind the binding of SAM to nsp16 and nsp16/nsp10 heterodimer by employing molecular dynamics simulations in conjunction with the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) method. We observed from our simulations that the presence of nsp10 increases the favorable van der Waals and electrostatic interactions between SAM and nsp16. Thus, nsp10 acts as a stimulator for the strong binding of SAM to nsp16. The hydrophobic interactions were predominately identified for the nsp16-nsp10 interactions. Also, the stable hydrogen bonds between Ala83 (nsp16) and Tyr96 (nsp10), and between Gln87 (nsp16) and Leu45 (nsp10) play a vital role in the dimerization of nsp16 and nsp10. Besides, Computational Alanine Scanning (CAS) mutagenesis was performed, which revealed hotspot mutants, namely I40A, V104A, and R86A for the dimer association. Hence, the dimer interface of nsp16/nsp10 could also be a potential target in retarding the 2'-O-MTase activity in SARS-CoV-2. Overall, our study provides a comprehensive understanding of the dynamic and thermodynamic process of binding nsp16 and nsp10 that will contribute to the novel design of peptide inhibitors based on nsp16.

Plant-derived natural polyphenols as potential antiviral drugs against SARS-CoV-2 via RNA-dependent RNA polymerase (RdRp) inhibition: an in-silico analysis

The sudden outburst of Coronavirus disease (COVID-19) caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) poses a massive threat to global public health. Currently, no therapeutic drug or vaccine exists to treat COVID-19. Due to the time taking process of new drug development, drug repurposing might be the only viable solution to tackle COVID-19. RNA-dependent RNA polymerase (RdRp) catalyzes SARS-CoV-2 RNA replication and hence, is an obvious target for antiviral drug design. Interestingly, several plant-derived polyphenols effectively inhibit the RdRp of other RNA viruses. More importantly, polyphenols have been used as dietary supplementations for a long time and played beneficial roles in immune homeostasis. We were curious to study the binding of polyphenols with SARS-CoV-2 RdRp and assess their potential to treat COVID-19. Herein, we made a library of polyphenols that have shown substantial therapeutic effects against various diseases. They were successfully docked in the catalytic pocket of RdRp. The investigation reveals that EGCG, theaflavin (TF1), theaflavin-3'-O-gallate (TF2a), theaflavin-3'-gallate (TF2b), theaflavin 3,3'-digallate (TF3), hesperidin, quercetagetin, and myricetin strongly bind to the active site of RdRp. Further, a 150-ns molecular dynamic simulation revealed that EGCG, TF2a, TF2b, TF3 result in highly stable bound conformations with RdRp. The binding free energy components calculated by the MM-PBSA also confirm the stability of the complexes. We also performed a detailed analysis of ADME prediction, toxicity prediction, and target analysis for their druggability. Overall, our results suggest that EGCG, TF2a, TF2b, TF3 can inhibit RdRp and represent an effective therapy for COVID-19.Communicated by Ramaswamy H. Sarma.

Impacts of drug resistance mutations on the structural asymmetry of the HIV-2 protease

Background

Drug resistance is a severe problem in HIV treatment. HIV protease is a common target for the design of new drugs for treating HIV infection. Previous studies have shown that the crystallographic structures of the HIV-2 protease (PR2) in bound and unbound forms exhibit structural asymmetry that is important for ligand recognition and binding. Here, we investigated the effects of resistance mutations on the structural asymmetry of PR2. Due to the lack of structural data on PR2 mutants, the 3D structures of 30 PR2 mutants of interest have been modeled using an in silico protocol. Structural asymmetry analysis was carried out with an in-house structural-alphabet-based approach.

Results

The systematic comparison of the asymmetry of the wild-type structure and a large number of mutants highlighted crucial residues for PR2 structure and function. In addition, our results revealed structural changes induced by PR2 flexibility or resistance mutations. The analysis of the highlighted structural changes showed that some mutations alter protein stability or inhibitor binding.

Conclusions

This work consists of a structural analysis of the impact of a large number of PR2 resistant mutants based on modeled structures. It suggests three possible resistance mechanisms of PR2, in which structural changes induced by resistance mutations lead to modifications in the dimerization interface, ligand recognition or inhibitor binding.

Investigating Conformational Dynamics of Lewis Y Oligosaccharides and Elucidating Blood Group Dependency of Cholera Using Molecular Dynamics

Cholera is caused by Vibrio cholerae and is an example of a blood-group-dependent disease. Recent studies suggest that the receptor-binding B subunit of the cholera toxin (CT) binds histo-blood group antigens at a secondary binding site. Herein, we studied the conformational dynamics of Lewis Y (Le^Y) oligosaccharides, H-tetrasaccharides and A-pentasaccharides, in aqueous solution by conducting accelerated molecular dynamics (aMD) simulations. The flexible nature of both oligosaccharides was displayed in aMD simulations. Furthermore, aMD simulations revealed that for both oligosaccharides in the free form, ⁴C₁ and ¹C₄ puckers were sampled for all but GalNAc monosaccharides, while either the ⁴C₁ (GlcNAc, Gal, GalNAc) or ¹C₄ (Fuc2, Fuc3) pucker was sampled in the CT-bound forms. In aMD, the complete transition from the ⁴C₁ to ¹C₄ pucker was sampled for GlcNAc and Gal in both oligosaccharides. Further, we have observed a transition from the open to closed conformer in the case of A-pentasaccharide, while H-tetrasaccharide remains in the open conformation throughout the simulation. Both oligosaccharides adopted an open conformation in the CT binding site. Moreover, we have investigated the molecular basis of recognition of Le^Y oligosaccharides by the B subunit of the cholera toxin of classical and El Tor biotypes using the molecular mechanics generalized Born surface area (MM/GBSA) scheme. The O blood group determinant, H-tetrasaccharide, exhibits a stronger affinity to both biotypes compared to the A blood group determinant, A-pentasaccharide, which agrees with the experimental data. The difference in binding free energy between O and A blood group determinants mainly arises due to the increased entropic cost and desolvation energy in the case of A-pentasaccharide compared to that of H-tetrasaccharide. Our study also reveals that the terminal Fuc3 contributes most to the binding free energy compared to other carbohydrate residues as it forms multiple hydrogen bonds with CT. Overall, our study might help in designing glycomimetic drugs targeting the cholera toxin.

Exploring the potency of currently used drugs against HIV-1 protease of subtype D variant by using multiscale simulations

Acquired immune deficiency syndrome (AIDS) is caused by the human immunodeficiency virus (HIV), type 1 and 2. Further, the diversity in HIV-1 has given rise to many serotypes and recombinant strains. The currently used protease inhibitors have been developed for subtype B, although non-B subtype strains account for ∼ 90% of the global HIV infections. Subtype D is spreading rapidly and infecting a large population in North Africa and the Middle East. In the current study, molecular dynamics simulations in conjunction with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) scheme was used to investigate the potency of four drugs, namely atazanavir (ATV), darunavir (DRV), lopinavir (LPV) and tipranavir (TPV) against the subtype D variant. Our calculations predicted that the potency of the inhibitors decreased in the order TPV > ATV > DRV > LPV. TPV was found to be the most potent against subtype D due to an increase in van der Waals and electrostatic interactions and reduction in the desolvation energy compared to other inhibitors. This result is further supported by the hydrogen bond interactions between inhibitors and protease. Furthermore, our analyses suggested that the binding of TPV induced a more closed conformation of the flap compared to apo or other complexes. It was observed that TPV/PR^D has a lower cavity volume relative to the other three complexes leading to a tighter binding. The open conformation of the flap was observed for LPV/PR^D. We expect that this study might be useful for designing more potent inhibitors against HIV-1 subtype D. Communicated by Ramaswamy H. Sarma.

Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.

Characterizing the structural variability of HIV-2 protease upon the binding of diverse ligands using a structural alphabet approach

The HIV-2 protease (PR2) is an important target for designing new drugs against the HIV-2 infection. In this study, we explored the structural backbone variability of all available PR2 structures complexed with various inhibitors using a structural alphabet approach. 77% of PR2 positions are structurally variable, meaning they exhibit different local conformations in PR2 structures. This variability was observed all along the structure, particularly in the elbow and flap regions. A part of these backbone changes observed between the 18 PR2 is induced by intrinsic flexibility, and ligand binding putatively induces others occurring in the binding pocket. These latter changes could be important for PR2 adaptation to diverse ligands and are accompanied by changes outside the binding pocket. In addition, the study of the link between structural variability of the pocket and PR2-ligand interactions allowed us to localize pocket regions important for ligand binding and catalytic function, regions important for ligand recognition that adjust their backbone in response to ligand binding and regions important for the pocket opening and closing that have large intrinsic flexibility. Finally, we suggested that differences in ligand effectiveness for PR2 could be partially explained by different backbone deformations induced by these ligands. To conclude, this study is the first characterization of the PR2 structural variability considering ligand diversity. It provides information about the recognition of PR2 to various ligands and its mechanisms to adapt its local conformation to bound ligands that could help understand the resistance of PR2 to its inhibitors, a major antiretroviral class. Communicated by Ramaswamy H. Sarma.

Characterization of HIV-2 Protease Structure by Studying Its Asymmetry at the Different Levels of Protein Description

HIV-2 protease (PR2) is a homodimer, which is an important target in the treatment of the HIV-2 infection. In this study, we developed an in silico protocol to analyze and characterize the asymmetry of the unbound PR2 structure using three levels of protein description by comparing the conformation, accessibility, and flexibility of each residue in the two PR2 chains. Our results showed that 65% of PR2 residues have at least one of the three studied asymmetries (structural, accessibility, or flexibility) with 10 positions presenting the three asymmetries in the same time. In addition, we noted that structural and flexibility asymmetries are linked indicating that the structural asymmetry of some positions result from their large flexibility. By comparing the structural asymmetry of the crystallographic and energetically minimized structures of the unbound PR2, we confirmed that the structural asymmetry of unbound PR2 is an intrinsic property of this protein with an important role for the PR2 deformation upon ligand binding. This analysis also allowed locating asymmetries corresponding to crystallization artefacts. This study provides insight that will help to better understand the structural deformations of PR2 and to identify key positions for ligand binding.

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4-2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4-1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8-2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.

Exploration of the effect of sequence variations located inside the binding pocket of HIV-1 and HIV-2 proteases

HIV-2 protease (PR2) is naturally resistant to most FDA (Food and Drug Administration)-approved HIV-1 protease inhibitors (PIs), a major antiretroviral class. In this study, we compared the PR1 and PR2 binding pockets extracted from structures complexed with 12 ligands. The comparison of PR1 and PR2 pocket properties showed that bound PR2 pockets were more hydrophobic with more oxygen atoms and fewer nitrogen atoms than PR1 pockets. The structural comparison of PR1 and PR2 pockets highlighted structural changes induced by their sequence variations and that were consistent with these property changes. Specifically, substitutions at residues 31, 46, and 82 induced structural changes in their main-chain atoms that could affect PI binding in PR2. In addition, the modelling of PR1 mutant structures containing V32I and L76M substitutions revealed a cooperative mechanism leading to structural deformation of flap-residue 45 that could modify PR2 flexibility. Our results suggest that substitutions in the PR1 and PR2 pockets can modify PI binding and flap flexibility, which could underlie PR2 resistance against PIs. These results provide new insights concerning the structural changes induced by PR1 and PR2 pocket variation changes, improving the understanding of the atomic mechanism of PR2 resistance to PIs.

Analysis of the HIV-2 protease's adaptation to various ligands: characterization of backbone asymmetry using a structural alphabet

The HIV-2 protease (PR2) is a homodimer of 99 residues with asymmetric assembly and binding various ligands. We propose an exhaustive study of the local structural asymmetry between the two monomers of all available PR2 structures complexed with various inhibitors using a structural alphabet approach. On average, PR2 exhibits asymmetry in 31% of its positions-i.e., exhibiting different backbone local conformations in the two monomers. This asymmetry was observed all along its structure, particularly in the elbow and flap regions. We first differentiated structural asymmetry conserved in most PR2 structures from the one specific to some PR2. Then, we explored the origin of the detected asymmetry in PR2. We localized asymmetry that could be induced by PR2's flexibility, allowing transition from the semi-open to closed conformations and the asymmetry potentially induced by ligand binding. This latter could be important for the PR2's adaptation to diverse ligands. Our results highlighted some differences between asymmetry of PR2 bound to darunavir and amprenavir that could explain their differences of affinity. This knowledge is critical for a better description of PR2's recognition and adaptation to various ligands and for a better understanding of the resistance of PR2 to most PR2 inhibitors, a major antiretroviral class.

Trypsin-Ligand binding affinities calculated using an effective interaction entropy method under polarized force field

Molecular dynamics (MD) simulation in the explicit water is performed to study the interaction mechanism of trypsin-ligand binding under the AMBER force field and polarized protein-specific charge (PPC) force field combined the new developed highly efficient interaction entropy (IE) method for calculation of entropy change. And the detailed analysis and comparison of the results of MD simulation for two trypsin-ligand systems show that the root-mean-square deviation (RMSD) of backbone atoms, B-factor, intra-protein and protein-ligand hydrogen bonds are more stable under PPC force field than AMBER force field. Our results demonstrate that the IE method is superior than the traditional normal mode (Nmode) method in the calculation of entropy change and the calculated binding free energy under the PPC force field combined with the IE method is more close to the experimental value than other three combinations (AMBER-Nmode, AMBER-IE and PPC-Nmode). And three critical hydrogen bonds between trypsin and ligand are broken under AMBER force field. However, they are well preserved under PPC force field. Detailed binding interactions of ligands with trypsin are further analyzed. The present work demonstrates that the polarized force field combined the highly efficient IE method is critical in MD simulation and free energy calculation.

Molecular Mechanism and Energy Basis of Conformational Diversity of Antibody SPE7 Revealed by Molecular Dynamics Simulation and Principal Component Analysis

More and more researchers are interested in and focused on how a limited repertoire of antibodies can bind and correspondingly protect against an almost limitless diversity of invading antigens. In this work, a series of 200-ns molecular dynamics (MD) simulations followed by principal component (PC) analysis and free energy calculations were performed to probe potential mechanism of conformational diversity of antibody SPE7. The results show that the motion direction of loops H3 and L3 is different relative to each other, implying that a big structural difference exists between these two loops. The calculated energy landscapes suggest that the changes in the backbone angles ψ and φ of H-Y101 and H-Y105 provide significant contributions to the conformational diversity of SPE7. The dihedral angle analyses based on MD trajectories show that the side-chain conformational changes of several key residues H-W33, H-Y105, L-Y34 and L-W93 around binding site of SPE7 play a key role in the conformational diversity of SPE7, which gives a reasonable explanation for potential mechanism of cross-reactivity of single antibody toward multiple antigens.

Anti-Hemagglutinin Antibody Derived Lead Peptides for Inhibitors of Influenza Virus Binding

Antibodies against spike proteins of influenza are used as a tool for characterization of viruses and therapeutic approaches. However, development, production and quality control of antibodies is expensive and time consuming. To circumvent these difficulties, three peptides were derived from complementarity determining regions of an antibody heavy chain against influenza A spike glycoprotein. Their binding properties were studied experimentally, and by molecular dynamics simulations. Two peptide candidates showed binding to influenza A/Aichi/2/68 H3N2. One of them, termed PeB, with the highest affinity prevented binding to and infection of target cells in the micromolar region without any cytotoxic effect. PeB matches best the conserved receptor binding site of hemagglutinin. PeB bound also to other medical relevant influenza strains, such as human-pathogenic A/California/7/2009 H1N1, and avian-pathogenic A/Mute Swan/Rostock/R901/2006 H7N1. Strategies to improve the affinity and to adapt specificity are discussed and exemplified by a double amino acid substituted peptide, obtained by substitutional analysis. The peptides and their derivatives are of great potential for drug development as well as biosensing.

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.

Four Amino Acid Changes in HIV-2 Protease Confer Class-Wide Sensitivity to Protease Inhibitors

Protease is essential for retroviral replication, and protease inhibitors (PI) are important for treating HIV infection. HIV-2 exhibits intrinsic resistance to most FDA-approved HIV-1 PI, retaining clinically useful susceptibility only to lopinavir, darunavir, and saquinavir. The mechanisms for this resistance are unclear; although HIV-1 and HIV-2 proteases share just 38 to 49% sequence identity, all critical structural features of proteases are conserved. Structural studies have implicated four amino acids in the ligand-binding pocket (positions 32, 47, 76, and 82). We constructed HIV-2ROD9 molecular clones encoding the corresponding wild-type HIV-1 amino acids (I32V, V47I, M76L, and I82V) either individually or together (clone PRΔ4) and compared the phenotypic sensitivities (50% effective concentration [EC50]) of mutant and wild-type viruses to nine FDA-approved PI. Single amino acid replacements I32V, V47I, and M76L increased the susceptibility of HIV-2 to multiple PI, but no single change conferred class-wide sensitivity. In contrast, clone PRΔ4 showed PI susceptibility equivalent to or greater than that of HIV-1 for all PI. We also compared crystallographic structures of wild-type HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir to models of the PRΔ4 enzyme. These models suggest that the amprenavir sensitivity of PRΔ4 is attributable to stabilizing enzyme-inhibitor interactions in the P2 and P2' pockets of the protease dimer. Together, our results show that the combination of four amino acid changes in HIV-2 protease confer a pattern of PI susceptibility comparable to that of HIV-1, providing a structural rationale for intrinsic HIV-2 PI resistance and resolving long-standing questions regarding the determinants of differential PI susceptibility in HIV-1 and HIV-2.

IMPORTANCE

Proteases are essential for retroviral replication, and HIV-1 and HIV-2 proteases share a great deal of structural similarity. However, only three of nine FDA-approved HIV-1 protease inhibitors (PI) are active against HIV-2. The underlying reasons for intrinsic PI resistance in HIV-2 are not known. We examined the contributions of four amino acids in the ligand-binding pocket of the enzyme that differ between HIV-1 and HIV-2 by constructing HIV-2 clones encoding the corresponding HIV-1 amino acids and testing the PI susceptibilities of the resulting viruses. We found that the HIV-2 clone containing all four changes (PRΔ4) was as susceptible as HIV-1 to all nine PI. We also modeled the PRΔ4 enzyme structure and compared it to existing crystallographic structures of HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir. Our findings demonstrate that four positions in the ligand-binding cleft of protease are the primary cause of HIV-2 PI resistance.

Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism

Inhibitors of HIV-1 protease (HIV-1-pr) generally only bind to the active site of the protease. However, for some mutants such as V32I and M46L the TMC114 can bind not only to the active cavity but also to the groove of the flexible flaps. Although the second binding site suggests the higher efficiency of the drug against HIV-1-pr, the drug resistance in HIV-1-pr due to mutations cannot be ignored, which prompts us to investigate the molecular mechanisms of drug resistance and behavior of double bound TMC114 (2T) to HIV-1-pr. The conformational dynamics of HIV-1-pr and the binding of TMC114 to the WT, V32I and M46L mutants were investigated with all-atom molecular dynamic (MD) simulation. The 20 ns MD simulation shows many fascinating effects of the inhibitor binding to the WT and mutant proteases. MM-PBSA calculations explain the binding free energies unfavorable for the M46L and V32I mutants as compared to the WT. For the single binding (1T) the less binding affinity can be attributed to the entropic loss for both V32I-1T and M46L-1T. Although the second binding of TMC114 with flap does increase binding energy for the mutants (V32I-2T and M46L-2T), the considerable entropy loss results in the lower binding Gibbs free energies. Thus, binding of TMC114 in the flap region does not help much in the total gain in binding affinity of the system, which was verified from this study and thereby validating experiments.

Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations

Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.

g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations

Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages is described. g_mmpbsa was developed as part of the Open Source Drug Discovery (OSDD) consortium. Its aim is to integrate high-throughput molecular dynamics (MD) simulations with binding energy calculations. The tool provides options to select alternative atomic radii and different nonpolar solvation models including models based on the solvent accessible surface area (SASA), solvent accessible volume (SAV), and a model which contains both repulsive (SASA-SAV) and attractive components (described using a Weeks-Chandler-Andersen like integral method). We showcase the effectiveness of the tool by comparing the calculated interaction energy of 37 structurally diverse HIV-1 protease inhibitor complexes with their experimental binding free energies. The effect of varying several combinations of input parameters such as atomic radii, dielectric constant, grid resolution, solute-solvent dielectric boundary definition, and nonpolar models was investigated. g_mmpbsa can also be used to estimate the energy contribution per residue to the binding energy. It has been used to identify those residues in HIV-1 protease that are most critical for binding a range of inhibitors.

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53-MDM2 interaction.

The N-terminal sequence of tyrosine hydroxylase is a conformationally versatile motif that binds 14-3-3 proteins and membranes

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the synthesis of catecholamine neurotransmitters, and a reduction in TH activity is associated with several neurological diseases. Human TH is regulated, among other mechanisms, by Ser19-phosphorylation-dependent interaction with 14-3-3 proteins. The N-terminal sequence (residues 1-43), which corresponds to an extension to the TH regulatory domain, also interacts with negatively charged membranes. By using X-ray crystallography together with molecular dynamics simulations and structural bioinformatics analysis, we have probed the conformations of the Ser19-phosphorylated N-terminal peptide [THp-(1-43)] bound to 14-3-3γ, free in solution and bound to a phospholipid bilayer, and of the unphosphorylated peptide TH-(1-43) both free and bilayer bound. As seen in the crystal structure of THp-(1-43) complexed with 14-3-3γ, the region surrounding pSer19 adopts an extended conformation in the bound state, whereas THp-(1-43) adopts a bent conformation when free in solution, with higher content of secondary structure and higher number of internal hydrogen bonds. TH-(1-43) in solution presents the highest mobility and least defined structure of all forms studied, and it shows an energetically more favorable interaction with membranes relative to THp-(1-43). Cationic residues, notably Arg15 and Arg16, which are the recognition sites of the kinases phosphorylating at Ser19, are also contributing to the interaction with the membrane. Our results reveal the structural flexibility of this region of TH, in accordance with the functional versatility and conformational adaptation to different partners. Furthermore, this structural information has potential relevance for the development of therapeutics for neurodegenerative disorders, through modulation of TH-partner interactions.

High level of HIV-2 false positivity in KwaZulu-Natal province: a region of South Africa with a very high HIV-1 subtype C prevalence

Human immunodeficiency virus 2 (HIV-2) is found predominantly in West Africa. It is not unlikely, however, that HIV-2 may also be found in South Africa, due to the influx of immigrants into this country. It is important to distinguish between HIV-1 and HIV-2 since the clinical courses and treatment responses of these viruses are different. Routine serological methods for diagnosing HIV do not differentiate between HIV-1 and -2 infections, while rapid tests, viral load quantification and PCR are HIV-type--specific. The objective of this study was to describe the seroprevalence and molecular epidemiology of HIV-2 in KwaZulu-Natal, one of the regions with the highest HIV prevalence in the world and home of the two largest harbors in South Africa. HIV-1 positive samples were screened for antibodies against HIV-2, using a rapid test. The confirmation of HIV-2 positive samples was done by PCR. Of the 2,123 samples screened, 319 (15%) were identified as positive by the rapid test. None of these samples were confirmed positive by PCR. To explore this discrepancy in the results, a subset (n = 52) of the rapid HIV-2 positive samples was subjected to Western blotting. Thirty-seven (71%) of these were positive, yielding an overall HIV-2 seroprevalence of 10.6%. Three out of 28 (10.7%) Western blot positive samples were positive by a Pepti-LAV assay. This discrepancy between serological and molecular confirmation may be attributed to non-specific or cross-reacting antibodies. The use of rapid tests and Western blots for HIV-2 diagnosis in South Africa should be interpreted with caution.

Importance of polar solvation and configurational entropy for design of antiretroviral drugs targeting HIV-1 protease

Both KNI-10033 and KNI-10075 are high affinity preclinical HIV-1 protease (PR) inhibitors with affinities in the picomolar range. In this work, the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method has been used to investigate the potency of these two HIV-1 PR inhibitors against the wild-type and mutated proteases assuming that potency correlates with the affinity of the drugs for the target protein. The decomposition of the binding free energy reveals the origin of binding affinities or mutation-induced affinity changes. Our calculations indicate that the mutation I50V causes drug resistance against both inhibitors. On the other hand, we predict that the mutant I84V causes drug resistance against KNI-10075 while KNI-10033 is more potent against the I84V mutant compared to wild-type protease. Drug resistance arises mainly from unfavorable shifts in van der Waals interactions and configurational entropy. The latter indicates that neglecting changes in configurational entropy in the computation of relative binding affinities as often done is not appropriate in general. For the bound complex PR(I50V)-KNI-10075, an increased polar solvation free energy also contributes to the drug resistance. The importance of polar solvation free energies is revealed when interactions governing the binding of KNI-10033 or KNI-10075 to the wild-type protease are compared to the inhibitors darunavir or GRL-06579A. Although the contributions from intermolecular electrostatic and van der Waals interactions as well as the nonpolar component of the solvation free energy are more favorable for PR-KNI-10033 or PR-KNI-10075 compared to PR-DRV or PR-GRL-06579A, both KNI-10033 and KNI-10075 show a similar affinity as darunavir and a lower binding affinity relative to GRL-06579A. This is because of the polar solvation free energy which is less unfavorable for darunavir or GRL-06579A relative to KNI-10033 or KNI-10075. The importance of the polar solvation as revealed here highlights that structural inspection alone is not sufficient for identifying the key contributions to binding affinities and affinity changes for the design of drugs but that solvation effects must be taken into account. A detailed understanding of the molecular forces governing binding and drug resistance might assist in the design of new inhibitors against HIV-1 PR variants that are resistant against current drugs.

Molecular dynamics investigation on a series of HIV protease inhibitors: assessing the performance of MM-PBSA and MM-GBSA approaches

The binding free energies (ΔG(Bind)) obtained from molecular mechanics with Poisson-Boltzmann surface area (MM-PBSA) or molecular mechanics with Generalized Born surface area (MM-GBSA) calculations using molecular dynamics (MD) trajectories are the most popular procedures to measure the strength of interactions between a ligand and its receptor. Several attempts have been made to correlate the ΔG(Bind) and experimental IC(50) values in order to observe the relationship between binding strength of a ligand (with its receptor) and its inhibitory activity. The duration of MD simulations seems very important for getting acceptable correlation. Here, we are presenting a systematic study to estimate the reasonable MD simulation time for acceptable correlation between ΔG(Bind) and experimental IC(50) values. A comparison between MM-PBSA and MM-GBSA approaches is also presented at various time scales. MD simulations (10 ns) for 14 HIV protease inhibitors have been carried out by using the Amber program. MM-PBSA/GBSA based ΔG(Bind) have been calculated and correlated with experimental IC(50) values at different time scales (0-1 to 0-10 ns). This study clearly demonstrates that the MM-PBSA based ΔG(Bind) (ΔG(Bind)-PB) values provide very good correlation with experimental IC(50) values (quantitative and qualitative) when MD simulation is carried out for a longer time; however, MM-GBSA based ΔG(Bind) (ΔG(Bind)-GB) values show acceptable correlation for shorter time of simulation also. The accuracy of ΔG(Bind)-PB increases and ΔG(Bind)-GB remains almost constant with the increasing time of simulation.

Proteochemometric modeling of the bioactivity spectra of HIV-1 protease inhibitors by introducing protein-ligand interaction fingerprint

HIV-1 protease is one of the main therapeutic targets in HIV. However, a major problem in treatment of HIV is the rapid emergence of drug-resistant strains. It should be particularly helpful to clinical therapy of AIDS if one method can be used to predict antivirus capability of compounds for different variants. In our study, proteochemometric (PCM) models were created to study the bioactivity spectra of 92 chemical compounds with 47 unique HIV-1 protease variants. In contrast to other PCM models, which used Multiplication of Ligands and Proteins Descriptors (MLPD) as cross-term, one new cross-term, i.e. Protein-Ligand Interaction Fingerprint (PLIF) was introduced in our modeling. With different combinations of ligand descriptors, protein descriptors and cross-terms, nine PCM models were obtained, and six of them achieved good predictive abilities (Q(2)(test)>0.7). These results showed that the performance of PCM models could be improved when ligand and protein descriptors were complemented by the newly introduced cross-term PLIF. Compared with the conventional cross-term MLPD, the newly introduced PLIF had a better predictive ability. Furthermore, our best model (GD & P & PLIF: Q(2)(test) = 0.8271) could select out those inhibitors which have a broad antiviral activity. As a conclusion, our study indicates that proteochemometric modeling with PLIF as cross-term is a potential useful way to solve the HIV-1 drug-resistant problem.

Molecular dynamics simulation in virus research

Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus–host interactions and evolution, and present future perspectives on this technique.

Mutation-induced loop opening and energetics for binding of tamiflu to influenza N8 neuraminidase

Tamiflu, also known as oseltamivir (OTV), binds to influenza A neuraminidase (H5N1) with very high affinity (0.32 nM). However, this inhibitor binds to other neuraminidases as well. In the present work, a systematic computational study is performed to investigate the mechanism underlying the binding of oseltamivir to N8 neuraminidase (NA) in "open" and "closed" conformations of the 150-loop through molecular dynamics simulations and the popular and well established molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation method. Whereas the closed conformation is stable for wild type N8, it transforms into the open conformation for the mutants Y252H, H274Y, and R292K, indicating that bound to oseltamivir these mutants are preferentially in the open conformation. Our calculations show that the binding of wild type oseltamivir to the closed conformation of N8 neuraminidase is energetically favored compared to the binding to the open conformation. We observe water mediated binding of oseltamivir to the N8 neuraminidase in both conformations which is not seen in the case of binding of the same drug to the H5N1 neuraminidase. The decomposition of the binding free energy reveals the mechanisms underlying the binding and changes in affinity due to mutations. Considering the mutant N8 variants in the open conformation adopted during the simulations, we observe a significant loss in the size of the total binding free energy for the N8(Y252H)-OTV, N8(H274Y)-OTV, and N8(R292K)-OTV complexes compared to N8(WT)-OTV, mainly due to the decrease in the size of the intermolecular electrostatic energy. For R292K, an unfavorable shift in the van der Waals interactions also contributes to the drug resistance. The mutations cause a significant expansion in the active site cavity, increasing its solvent accessible surface compared to the crystal structures of both the open and closed conformations. Our study underscores the need to consider dynamics in rationalizing the structure-function relationships of various antiviral inhibitor-NA complexes.