Structural and energetic analysis on the complexes of clinically isolated subtype C HIV-1 proteases and approved inhibitors by molecular dynamics simulation

Secondary Structure and Conformational Stability of the Antigen Residues Making Contact with Antibodies

Antibodies are crucial biomolecules that bring high therapeutic efficacy in medicine and accurate molecular detection in diagnosis. Many studies have been devoted to analyzing the antigen-antibody interaction from the importance of understanding the antibody recognition mechanism. However, most of the previous studies examined the characteristic of the antibody for interaction. It is also informative to clarify the significant antigen residues contributing to the binding. To characterize the molecular interaction of antigens, we computationally analyzed 350 antigen-antibody complex structures by molecular mechanics (MM) calculations and molecular dynamics (MD) simulations. Based on the MM calculations, the antigen residues contributing to the binding were extracted from all the 350 complexes. The extracted residues are located at the antigen-antibody interface and are responsible for making contact with the antibody. The appearances of the charged polar residues, Asp, Glu, Arg, and Lys, were noticeably large. In contrast, the populations of the hydrophobic residues, Leu, Val, and Ala, were relatively low. The appearance frequencies of the other amino acid residues were almost close to the abundance of general proteins of eukaryotes. The binding score indicated that the hydrophilic interaction was dominant at the antigen-antibody contact instead of the hydrophobic one. The positively charged residues, Arg and Lys, remarkably contributed to the binding compared to the negatively charged ones, Asp and Glu. Considerable contributions were also observed for the noncharged polar residues, Asn and Gln. The analysis of the secondary structures of the extracted antigen residues suggested that there was no marked difference in recognition by antibodies among helix, sheet, turn, and coil. A long helix of the antigen sometimes made contact with antibody complementarity-determining regions, and a large sheet also frequently covered the antibody heavy and light chains. The turn structure was the most popularly observed at the contact with antibody among 350 complexes. Three typical complexes were picked up for each of the four secondary structures. MD simulations were performed to examine the stability of the interfacial structures of the antigens for these 12 complex models. The alterations of secondary structures were monitored through the simulations. The structural fluctuations of the contact residues were low compared with the other domains of antigen molecules. No drastic conversion was observed for every model during the 100 ns simulation. The motions of the interfacial antigen residues were small compared to the other residues on the protein surface. Therefore, diverse molecular conformations are possible for antibody recognition as long as the target areas are polar, nonflexible, and protruding on the protein surface.

Analysis of Binding Modes of Antigen-Antibody Complexes by Molecular Mechanics Calculation

Antibodies are one of the most important protein molecules in biopharmaceutics. Due to the recent advance in technology for producing monoclonal antibodies, many structural data are available on the antigen-antibody complexes. To characterize the molecular interaction in antigen-antibody recognition, we computationally analyzed 500 complex structures by molecular mechanics calculations. The presence of Ser and Tyr is markedly large in the complementarity-determining regions (CDRs). Although Ser is abundant in CDRs, its contribution to the binding score is not large. Instead, Tyr, Asp, Glu, and Arg significantly contribute to the molecular interaction from the viewpoint of the binding score. The decomposition of the binding score suggests that the hydrophilic interaction is predominant in all CDRs compared with the hydrophobic one. The contribution of the heavy chain is larger than that of the light chain. In particular, H2 and H3 largely contribute to the binding interaction. Tyr is a main contributing residue both in H2 and H3. The positively charged residue Arg also significantly contributes to the binding score in H3, while the contribution of Lys is small. The appearance of Ser is remarkable in H2, and Asp is abundant in H3. The non-charged polar residues, Thr, Asn, and Gln, appear much in H2, compared to appearing in H3. The negatively charged residues Asp and Glu significantly contribute to the binding score in H3. The contributions of Phe and Trp are not large in spite that the aromatic residues are capable of making the π-π or CH-π interaction. Gly is commonly abundant both in H2 and H3. The average distance of the shortest direct hydrogen bond between the antigen and antibody is longer than that of the hydrogen bonds observed in the complexes between compounds and their target proteins. Therefore, the antigen-antibody interface is not so tight as the compound-target protein interface. The calculation of shape complementarity is consistent with the result of the hydrogen bonds in that the fitness of the antigen-antibody contact is not so high as that of the compound-target protein contact. There exist many water molecules at the antigen-antibody interface. These findings suggest that Tyr, Asp, Glu, and Arg are rich in H3 and work as major contributors for the interaction with the antigen. Ser, Thr, Asn, and Gln are rich in H2 and support the interaction with enhancing molecular fitness. Gly is helpful in increasing flexibility and geometrical diversity. Because the antigen-antibody binding is fundamentally hydrophilic-driven, the non-polar residues are unfavorable for mediating the contact even for the aromatic residues such as Phe and Trp.

The mechanistic overview of SARS-CoV-2 using angiotensin-converting enzyme 2 to enter the cell for replication: possible treatment options related to the renin-angiotensin system

The SARS-CoV-2 pandemic is a healthcare crisis caused by insufficient knowledge applicable to effectively combat the virus. Therefore, different scientific discovery strategies need to be connected, to generate a rational treatment which can be made available as rapidly as possible. This relies on a solid theoretical understanding of the mechanisms of SARS-CoV-2 infection and host responses, which is coupled to the practical experience of clinicians that are treating patients. Because SARS-CoV-2 enters the cell by binding to angiotensin-converting enzyme 2 (ACE2), targeting ACE2 to prevent such binding seems an obvious strategy to combat infection. However, ACE2 performs its functions outside the cell and was found to enter the cell only by angiotensin II type 1 receptor (AT1R)-induced endocytosis, after which ACE2 is destroyed. This means that preventing uptake of ACE2 into the cell by blocking AT1R would be a more logical approach to limit entry of SARS-CoV-2 into the cell. Since ACE2 plays an important protective role in maintaining key biological processes, treatments should not disrupt the functional capacity of ACE2, to counterbalance the negative effects of the infection. Based on known mechanisms and knowledge of the characteristics of SARS-CoV we propose the hypothesis that the immune system facilitates SARS-CoV-2 replication which disrupts immune regulatory mechanisms. The proposed mechanism by which SARS-CoV-2 causes disease immediately suggests a possible treatment, since the AT1R is a key player in this whole process. AT1R antagonists appear to be the ideal candidate for the treatment of SARS-CoV-2 infection. AT1R antagonists counterbalance the negative consequences of angiotesnin II and, in addition, they might even be involved in preventing the cellular uptake of the virus without interfering with ACE2 function. AT1R antagonists are widely available, cheap, and safe. Therefore, we propose to consider using AT1R antagonists in the treatment of SARS-CoV-2.

New insights into the in silico prediction of HIV protease resistance to nelfinavir

The Human Immunodeficiency Virus type 1 protease enzyme (HIV-1 PR) is one of the most important targets of antiretroviral therapy used in the treatment of AIDS patients. The success of protease-inhibitors (PIs), however, is often limited by the emergence of protease mutations that can confer resistance to a specific drug, or even to multiple PIs. In the present study, we used bioinformatics tools to evaluate the impact of the unusual mutations D30V and V32E over the dynamics of the PR-Nelfinavir complex, considering that codons involved in these mutations were previously related to major drug resistance to Nelfinavir. Both studied mutations presented structural features that indicate resistance to Nelfinavir, each one with a different impact over the interaction with the drug. The D30V mutation triggered a subtle change in the PR structure, which was also observed for the well-known Nelfinavir resistance mutation D30N, while the V32E exchange presented a much more dramatic impact over the PR flap dynamics. Moreover, our in silico approach was also able to describe different binding modes of the drug when bound to different proteases, identifying specific features of HIV-1 subtype B and subtype C proteases.

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.

Identifying the molecular mechanics and binding dynamics characteristics of potent inhibitors to HIV-1 protease

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the primary inhibition targets for chemotherapy of AIDS because of its critical role in the replication cycle of the HIV. In this work, a combinatory coarse-grained and atomistic simulation method was developed for dissecting molecular mechanisms and binding process of inhibitors to the active site of HIV-1 PR, in which 35 typical inhibitors were trialed. We found that the molecular size and stiffness of the inhibitors and the binding energy between the inhibitors and PR play important roles in regulating the binding process. Comparatively, the smaller and more flexible inhibitors have larger binding energy and higher binding rates; they even bind into PR without opening the flaps. In contrast, the larger and stiffer inhibitors have lower binding energy and lower binding rate, and their binding is subjected to the opening and gating of the PR flaps. Furthermore, the components of binding free energy were quantified and analyzed by their dependence on the molecular size, structures, and hydrogen bond networks of inhibitors. Our results also deduce significant dynamics descriptors for determining the quantitative structure and property relationship in potent drug ligands for HIV-1 PR inhibition.