Design, in silico Evaluation, and Determination of Antitumor Activity of Potential Inhibitors Against Protein Kinases: Application to BCR-ABL Tyrosine Kinase.

Despite significant progress made over the past two decades in the treatment of chronic myeloid leukemia (CML), there is still an unmet need for effective and safe agents to treat patients with resistance and intolerance to the drugs used in clinic. In this work, we designed 2-arylaminopyrimidine amides of isoxazole-3-carboxylic acid, assessed in silico their inhibitory potential against Bcr-Abl tyrosine kinase, and determined their antitumor activity in K562 (CML), HL-60 (acute promyelocytic leukemia), and HeLa (cervical cancer) cells. Based on the analysis of computational and experimental data, three compounds with the antitumor activity against K562 and HL-60 cells were identified. The lead compound efficiently suppressed the growth of these cells, as evidenced by the low IC50 values of 2.8 ± 0.8 µM (K562) and 3.5 ± 0.2 µM (HL-60). The obtained compounds represent promising basic structures for the design of novel, effective, and safe anticancer drugs able to inhibit the catalytic activity of Bcr-Abl kinase by blocking the ATP-binding site of the enzyme.

Virtual screening and identification of promising therapeutic compounds against drug-resistant Mycobacterium tuberculosis ß-ketoacyl-acyl carrier protein synthase I (KasA).

The emergence of new Mycobacterium tuberculosis (Mtb) strains resistant to the key drugs currently used in the clinic for tuberculosis treatment can substantially reduce the probability of therapy success, causing the relevance and importance of studies on the development of novel potent antibacterial agents targeting different vulnerable spots of Mtb. In this study, 28,860 compounds from the library of bioactive molecules were screened to identify novel potential inhibitors of ß-ketoacyl-acyl carrier protein synthase I (KasA), one of the key enzymes involved in the biosynthesis of mycolic acids of the Mtb cell wall. In doing so, we used a structure-based virtual screening approach to drug repurposing that included high-throughput docking of the C171Q KasA enzyme with compounds from the library of bioactive molecules including the FDA-approved drugs and investigational drug candidates, assessment of the binding affinity for the docked ligand/C171Q KasA complexes, and molecular dynamics simulations followed by binding free energy calculations. As a result, post-modeling analysis revealed 6 top-ranking compounds exhibiting a strong attachment to the malonyl binding site of the enzyme, as evidenced by the values of binding free energy which are significantly lower than those predicted for the KasA inhibitor TLM5 used in the calculations as a positive control. In light of the data obtained, the identified compounds are suggested to form a good basis for the development of new antitubercular molecules of clinical significance with activity against the KasA enzyme of Mtb.Communicated by Ramaswamy H. Sarma.

Repurposing Navitoclax to block SARS-CoV-2 fusion and entry by targeting heptapeptide repeat sequence 1 in S2 protein.

Along with the long pandemic of COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has come the dilemma of emerging viral variants of concern (VOC), particularly Omicron and its subvariants, able to deftly escape immune surveillance and the otherwise protective effect of current vaccines and antibody drugs. We previously identified a peptide-based pan-CoV fusion inhibitor, termed as EK1, able to bind the HR1 region in viral spike (S) protein S2 subunit. This effectively blocked formation of the six-helix bundle (6-HB) fusion core and, thus, showed efficacy against all human coronaviruses (HCoVs). EK1 is now in phase 3 clinical trials. However, the peptide drug generally lacks oral availability. Therefore, we herein performed a structure-based virtual screening of the libraries of biologically active molecules and identified nine candidate compounds. One is Navitoclax, an orally active anticancer drug by inhibition of Bcl-2. Like EK1 peptide, it could bind HR1 and block 6-HB formation, efficiently inhibiting fusion and infection of all SARS-CoV-2 variants tested, as well as SARS-CoV and MERS-CoV, with IC50 values ranging from 0.5 to 3.7 µM. These findings suggest that Navitoclax is a promising repurposed drug candidate for development as a safe and orally available broad-spectrum antiviral drug to combat the current SARS-CoV-2 and its variants, as well as other HCoVs.

AI-Driven De Novo Design and Molecular Modeling for Discovery of Small-Molecule Compounds as Potential Drug Candidates Targeting SARS-CoV-2 Main Protease.

Over the past three years, significant progress has been made in the development of novel promising drug candidates against COVID-19. However, SARS-CoV-2 mutations resulting in the emergence of new viral strains that can be resistant to the drugs used currently in the clinic necessitate the development of novel potent and broad therapeutic agents targeting different vulnerable spots of the viral proteins. In this study, two deep learning generative models were developed and used in combination with molecular modeling tools for de novo design of small molecule compounds that can inhibit the catalytic activity of SARS-CoV-2 main protease (Mpro), an enzyme critically important for mediating viral replication and transcription. As a result, the seven best scoring compounds that exhibited low values of binding free energy comparable with those calculated for two potent inhibitors of Mpro, via the same computational protocol, were selected as the most probable inhibitors of the enzyme catalytic site. In light of the data obtained, the identified compounds are assumed to present promising scaffolds for the development of new potent and broad-spectrum drugs inhibiting SARS-CoV-2 Mpro, an attractive therapeutic target for anti-COVID-19 agents.

In silico design and computational evaluation of novel 2-arylaminopyrimidine-based compounds as potential multi-targeted protein kinase inhibitors: application for the native and mutant (T315I) Bcr-Abl tyrosine kinase.

An integrated computational approach to drug discovery was used to identify novel potential inhibitors of the native and mutant (T315I) Bcr-Abl tyrosine kinase, the enzyme playing a key role in the pathogenesis of chronic myeloid leukemia (CML). This approach included i) design of chimeric molecules based on the 2-arylaminopyrimidine fragment, the main pharmacophore of the Abl kinase inhibitors imatinib and nilotinib used in the clinic for the CML treatment, ii) molecular docking of these compounds with the ATP-binding site of the native and mutant Abl kinase, iii) refinement of the ligand-binding poses by the quantum chemical method PM7, iv) molecular dynamics simulations of the ligand/Abl complexes, and v) prediction of the ligand/Abl binding affinity in terms of scoring functions of molecular docking, machine learning, quantum chemistry, and molecular dynamics. As a result, five top-ranking compounds able to effectively block the enzyme catalytic site were identified. According to the data obtained, these compounds exhibit close modes of binding to the Abl kinase active site that are mainly provided by hydrogen bonds and multiple van der Waals contacts. The identified compounds show high binding affinity to the native and mutant Abl kinase comparable with the one calculated for the FDA-approved kinase-targeted inhibitors imatinib, nilotinib, and ponatinib used in the calculations as a positive control. The results obtained testify to the predicted drug candidates against CML may serve as good scaffolds for the design of novel anticancer agents able to target the ATP-binding pocket of the native and mutant Abl kinase.Communicated by Ramaswamy H. Sarma.

Metabolomic profiling of three Araucaria species, and their possible potential role against COVID-19.

The COVID-19 pandemic in Egypt is a part of the worldwide global crisis of coronavirus 2 (SARS-CoV-2). The contagious life-threatening condition causes acute respiratory syndrome. The present study aimed to assess the compounds identified by LC-MS of the methanolic leaves extracts from three conifers trees cultivated in Egypt (Araucaria bidwillii, Araucaria. cunninghamii and Araucaria heterophylla) via docking technique as potential inhibitor of COVID-19 virus on multiple targets; viral main protease (Mpro, 6LU7), non-structural protein-16 which is a methyl transferase (nsp16, 6W4H) and RNA dependent RNA polymerase (nsp12, 7BV2). Among the three targets, nsp16 was the best target recognized by the tested compounds as can be deduced from docking studies. Moreover, the methanolic extract of A. cunninghamii showed the highest radical-scavenging activity using (DPPH test) with 53.7 µg/mL comparable to ascorbic acid with IC50 = 46 µg/mL The anti-inflammatory potential carried using enzyme linked immunoassay showed the highest activity for A. cunninghamii and A. bidwillii followed by A. heterophylla with IC50 = 23.20 ± 1.17 µg/mL, 82.83 ± 3.21 µg/mL and 221.13 ± 6.7 µg/mL, respectively (Celecoxib was used as a standard drug with IC50 = 141.92 ± 4.52 µg/mL). Moreover, a molecular docking study was carried for the LC-MS annotated metabolites to validate their anti-inflammatory inhibitory effect using Celecoxib as a reference compound and showed a high docking score (-7.7 kcal/mol) for Octadecyl (E) P-coumarate and (-7.3 kcal/mol) for secoisolariciresinol rhamnoside.Communicated by Ramaswamy H. Sarma.

Application of deep learning and molecular modeling to identify small drug-like compounds as potential HIV-1 entry inhibitors.

A generative adversarial autoencoder for the rational design of potential HIV-1 entry inhibitors able to block CD4-binding site of the viral envelope protein gp120 was developed. To do this, the following studies were carried out: (i) an autoencoder architecture was constructed; (ii) a virtual compound library of potential anti-HIV-1 agents for training the neural network was formed by the concept of click chemistry allowing one to generate a large number of drug candidates by their assembly from small modular units; (iii) molecular docking of all compounds from this library with gp120 was made and calculations of the values of binding free energy were performed; (iv) molecular fingerprints of chemical compounds from the training dataset were generated; (v) training of the developed autoencoder was implemented followed by the validation of this neural network using more than 21 million molecules from the ZINC15 database. As a result, three small drug-like compounds that exhibited the high-affinity binding to gp120 were identified. According to the data from molecular docking, machine learning, quantum chemical calculations, and molecular dynamics simulations, these compounds show the low values of binding free energy in the complexes with gp120 similar to those calculated using the same computational protocols for the HIV-1 entry inhibitors NBD-11021 and NBD-14010, highly potent and broad anti-HIV-1 agents presenting a new generation of the viral CD4 antagonists. The identified CD4-mimetic candidates are suggested to present good scaffolds for the design of novel antiviral drugs inhibiting the early stages of HIV-1 infection.

Computational discovery of small drug-like compounds as potential inhibitors of SARS-CoV-2 main protease.

A computational approach to in silico drug discovery was carried out to identify small drug-like compounds able to show structural and functional mimicry of the high affinity ligand X77, potent non-covalent inhibitor of SARS-COV-2 main protease (MPro). In doing so, the X77-mimetic candidates were predicted based on the crystal X77-MPro structure by a public web-oriented virtual screening platform Pharmit. Models of these candidates bound to SARS-COV-2 MPro were generated by molecular docking, quantum chemical calculations and molecular dynamics simulations. At the final point, analysis of the interaction modes of the identified compounds with MPro and prediction of their binding affinity were carried out. Calculation revealed 5 top-ranking compounds that exhibited a high affinity to the active site of SARS-CoV-2 MPro. Insights into the ligand - MPro models indicate that all identified compounds may effectively block the binding pocket of SARS-CoV-2 MPro, in line with the low values ​​of binding free energy and dissociation constant. Mechanism of binding of these compounds to MPro is mainly provided by van der Waals interactions with the functionally important residues of the enzyme, such as His-41, Met-49, Cys-145, Met-165, and Gln-189 that play a role of the binding hot spots assisting the predicted molecules to effectively interact with the MPro active site. The data obtained show that the identified X77-mimetic candidates may serve as good scaffolds for the design of novel antiviral agents able to target the active site of SARS-CoV-2 MPro.Communicated by Ramaswamy H. Sarma.

In Silico Identification of Novel Aromatic Compounds as Potential HIV-1 Entry Inhibitors Mimicking Cellular Receptor CD4.

Despite recent progress in the development of novel potent HIV-1 entry/fusion inhibitors, there are currently no licensed antiviral drugs based on inhibiting the critical interactions of the HIV-1 envelope gp120 protein with cellular receptor CD4. In this connection, studies on the design of new small-molecule compounds able to block the gp120-CD4 binding are still of great value. In this work, in silico design of drug-like compounds containing the moieties that make the ligand active towards gp120 was performed within the concept of click chemistry. Complexes of the designed molecules bound to gp120 were then generated by molecular docking and optimized using semiempirical quantum chemical method PM7. Finally, the binding affinity analysis of these ligand/gp120 complexes was performed by molecular dynamic simulations and binding free energy calculations. As a result, five top-ranking compounds that mimic the key interactions of CD4 with gp120 and show the high binding affinity were identified as the most promising CD4-mimemic candidates. Taken together, the data obtained suggest that these compounds may serve as promising scaffolds for the development of novel, highly potent and broad anti-HIV-1 therapeutics.

Potential HIV-1 fusion inhibitors mimicking gp41-specific broadly neutralizing antibody 10E8: In silico discovery and prediction of antiviral potency.

An integrated computational approach to in silico drug design was used to identify novel HIV-1 fusion inhibitor scaffolds mimicking broadly neutralizing antibody (bNab) 10E8 targeting the membrane proximal external region (MPER) of the HIV-1 gp41 protein. This computer-based approach included (i) generation of pharmacophore models representing 3D-arrangements of chemical functionalities that make bNAb 10E8 active towards the gp41 MPER segment, (ii) shape and pharmacophore-based identification of the 10E8-mimetic candidates by a web-oriented virtual screening platform pepMMsMIMIC, (iii) high-throughput docking of the identified compounds with the gp41 MPER peptide, and (iv) molecular dynamics simulations of the docked structures followed by binding free energy calculations. As a result, eight hits-able to mimic pharmacophore properties of bNAb 10E8 by specific and effective interactions with the MPER region of the HIV-1 protein gp41 were selected as the most probable 10E8-mimetic candidates. Similar to 10E8, the predicted compounds target the critically important residues of a highly conserved hinge region of the MPER peptide that provides a conformational flexibility necessary for its functioning in cell-virus membrane fusion process. In light of the data obtained, the identified small molecules may present promising HIV-1 fusion inhibitor scaffolds for the design of novel potent antiviral drugs.

Computational identification of novel entry inhibitor scaffolds mimicking primary receptor CD4 of HIV-1 gp120.

Virtual screening of novel entry inhibitor scaffolds mimicking primary receptor CD4 of HIV-1 gp120 was carried out in conjunction with evaluation of their potential inhibitory activity by molecular modeling. To do this, pharmacophore models presenting different sets of the hotspots of cellular receptor CD4 for its interaction with gp120 were generated. These models were used as the templates for identification of CD4-mimetic candidates by the pepMMsMIMIC screening platform. Complexes of these candidates with gp120 were built by high-throughput ligand docking and their stability was estimated by molecular dynamics simulations and binding free energy calculations. As a result, five top hits that exhibited strong attachment to the two well-conserved hotspots of the gp120 CD4-binding site were selected for the final analysis. In analogy to CD4, the identified compounds make hydrogen bonds with Asp-368gp120 and multiple van der Waals contacts with the gp120 residues that bind to Phe-43CD4, resulting in destruction of the critical interactions of gp120 with Phe-43CD4 and Arg-59CD4. The complexes of the CD4-mimetic candidates with gp120 show relative conformational stability within the molecular dynamics simulations and expose the high percentage occupancies of intermolecular hydrogen bonds, in line with the data on the binding free energy calculations. In light of these findings, the identified compounds are considered as good scaffolds for the development of new functional antagonists of viral entry with broad HIV-1 neutralization.

Computational discovery of novel HIV-1 entry inhibitors based on potent and broad neutralizing antibody VRC01.

Computational prediction of novel HIV-1 entry inhibitors presenting peptidomimetics of broadly neutralizing antibody (bNAb) VRC01 was carried out based on the analysis of the X-ray complex of this antibody antigen-binding fragment with the HIV envelope gp120 core. Using these empirical data, peptidomimetic candidates of bNAb VRC01 were identified by a public web-oriented virtual screening platform (pepMMsMIMIC) and models of these candidates bound to gp120 were generated by molecular docking. At the final point, the stability of the complexes of these molecules with gp120 was estimated by molecular dynamics and binding free energy calculations. The calculations identified six molecules exhibiting a high affinity to the HIV-1 gp120 protein. These molecules were selected as the most probable peptidomimetics of bNAb VRC01. In a mechanism similar to that of bNAb VRC01, these compounds were predicted to block the functionally conserved regions of gp120 critical for the HIV-1 binding to cellular receptor CD4. The docked structures of the identified molecules with gp120 do not undergo substantial rearrangements during the molecular dynamics simulations, in agreement with the low values of free energy of their formation. Based on these findings, the selected compounds are considered as promising basic structures for the rational design of novel, potent, and broad-spectrum anti-HIV-1 therapeutics.

In silico design of novel broad anti-HIV-1 agents based on glycosphingolipid ß-galactosylceramide, a high-affinity receptor for the envelope gp120 V3 loop.

Novel anti-Human immunodeficiency virus (HIV)-1 agents targeting the V3 loop of envelope protein gp120 were designed by computer modeling based on glycosphingolipid ß-galactosylceramide (ß-GalCer), which is an alternative receptor allowing HIV-1 entry into CD4-negative cells of neural and colonic origin. Models of these ß-GalCer analogs bound to the V3 loops from five various HIV-1 variants were generated by molecular docking and their stability was estimated by molecular dynamics (MDs) and binding free energy simulations. Specific binding to the V3 loop was accomplished primarily by non-conventional XH…π interactions between CH/OH sugar groups of the glycolipids and the conserved V3 residues with π-conjugated side chains. The designed compounds were found to block the tip and/or the base of the V3 loop, which form invariant structural motifs that contain residues critical for cell tropism. With the MDs calculations, the docked models of the complexes of the ß-GalCer analogs with V3 are energetically stable in all of the cases of interest and exhibit low values of free energy of their formation. Based on the data obtained, these compounds are considered as promising basic structures for the rational design of novel, potent, and broad-spectrum anti-HIV-1 therapeutics.

Discovery of novel anti-HIV-1 agents based on a broadly neutralizing antibody against the envelope gp120 V3 loop: a computational study.

A computer-aided search for novel anti-HIV-1 agents able to mimic the pharmacophore properties of broadly neutralizing antibody (bNAb) 3074 was carried out based on the analysis of X-ray complexes of this antibody Fab with the MN, UR29, and VI191 peptides from the V3 loop of the HIV envelope protein gp120. Using these empirical data, peptidomimetic candidates of bNAb 3074 were identified by a public, web-oriented virtual screening platform (pepMMsMIMIC) and models of these candidates bound to the above V3 peptides were generated by molecular docking. The docking calculations identified four molecules exhibiting a high affinity to all of the V3 peptides. These molecules were selected as the most probable peptidomimetics of bNAb 3074. Finally, the stability of the complexes of these molecules with the MN, UR29, and VI191 V3 peptides was estimated by molecular dynamics and free energy simulations. Specific binding to the V3 loop was accomplished primarily by π-π interactions between the aromatic rings of the peptidomimetics and the conserved Phe-20 and/or Tyr-21 of the V3 immunogenic crown. In a mechanism similar to that of bNAb 3074, these compounds were found to block the tip of the V3 loop forming its invariant structural motif that contains residues critical for cell tropism. Based on these findings, the compounds selected are considered as promising basic structures for the rational design of novel, potent, and broad-spectrum anti-HIV-1 therapeutics.

Hydrophobic amino acids in the hinge region of the 5A apolipoprotein mimetic peptide are essential for promoting cholesterol efflux by the ABCA1 transporter.

The bihelical apolipoprotein mimetic peptide 5A effluxes cholesterol from cells and reduces inflammation and atherosclerosis in animal models. We investigated how hydrophobic residues in the hinge region between the two helices are important in the structure and function of this peptide. By simulated annealing analysis and molecular dynamics modeling, two hydrophobic amino acids, F-18 and W-21, in the hinge region were predicted to be relatively surface-exposed and to interact with the aqueous solvent. Using a series of 5A peptide analogs in which F-18 or W-21 was changed to either F, W, A, or E, only peptides with hydrophobic amino acids in these two positions were able to readily bind and solubilize phospholipid vesicles. Compared with active peptides containing F or W, peptides containing E in either of these two positions were more than 10-fold less effective in effluxing cholesterol by the ABCA1 transporter. Intravenous injection of 5A in C57BL/6 mice increased plasma-free cholesterol (5A: 89.9 ± 13.6 mg/dl; control: 38.7 ± 4.3 mg/dl (mean ± S.D.); P < 0.05) and triglycerides (5A: 887.0 ± 172.0 mg/dl; control: 108.9 ± 9.9 mg/dl; P < 0.05), whereas the EE peptide containing E in both positions had no effect. Finally, 5A increased cholesterol efflux approximately 2.5-fold in vivo from radiolabeled macrophages, whereas the EE peptide was inactive. These results provide a rationale for future design of therapeutic apolipoprotein mimetic peptides and provide new insights into the interaction of hydrophobic residues on apolipoproteins with phospholipids in the lipid microdomain created by the ABCA1 transporter during the cholesterol efflux process.

Structural analysis of the envelope gp120 V3 loop for some HIV-1 variants circulating in the countries of Eastern Europe.

The V3 loop on gp120 from human immunodeficiency virus type 1 (HIV-1) is a focus of many research groups involved in anti-AIDS drug development because this region of the protein is a principal target for neutralizing antibodies and a major determinant for cell tropism and syncytium formation. In this study, the nucleotide sequences of the env gene region coding the V3 loop were determined by DNA sequencing methods for four novel HIV-1 strains that circulate in the countries of Eastern Europe, such as Russia, Belarus, Ukraine, etc. Based on the empirical data obtained, the 3D structures of the V3 loops associated with these viral modifications were generated by computer modeling and then compared to discover similarities in the spatial arrangement of this functionally important site of gp120. Despite the HIV-1 genetic variety, several regions of the V3 loop that contain residues critical for cell tropism were shown to be structurally invariant, which may explain its exceptional role in a co-receptor usage. These data together with those on the biological activity of the V3 individual residues clearly show that these conserved structural motifs of gp120 represent potential HIV-1 weak points most suitable for therapeutic intervention.

Discovery of novel promising targets for anti-AIDS drug developments by computer modeling: application to the HIV-1 gp120 V3 loop.

The V3 loop on gp120 from HIV-1 is a focus of many research groups involved in anti-AIDS drug studies, because this region of the protein determines the preference of the virus for T-lymphocytes or primary macrophages. Although the V3 loop governs cell tropism and, for this reason, exhibits one of the most attractive targets for anti-HIV-1 drug developments, its high sequence variability is a major complicating factor. Nevertheless, the data on the spatial arrangement of V3 obtained here for different HIV-1 subtypes by computer modeling clearly show that, despite a wide range of 3D folds, this functionally important site of gp120 forms at least three structurally invariant segments, which contain residues critical for cell tropism. It is evident that these conserved V3 segments represent potential HIV-1 vulnerable spots and, therefore, provide a blueprint for the design of novel, potent and broad antiviral agents able to stop the HIV's spread.

Human immunodeficiency virus-1 gp120 V3 loop for anti-acquired immune deficiency syndrome drug discovery: computer-aided approaches to the problem solving.

INTRODUCTION: The V3 loop on gp120 from HIV-1 is a focus of many research groups involved in anti-AIDS drug development because this region of the protein is the principal target for neutralizing antibodies and determines the preference of the virus for T-lymphocytes or primary macrophages. AREAS COVERED: This review summarizes findings related to the 3D structure, conformational mobility, function, antigenicity and immunogenicity of the HIV-1 V3 loop. Particular consideration is given to the V3 loop core sequence Gly-Pro-Gly-Arg/Gln-Ala-Phe, which forms the HIV-1 gp120 immunogenic tip, the role of which has not been completely determined in the virus pathogenesis. New computer-aided approaches for designing potential HIV-1 entry inhibitors are illustrated by a series of examples in which promising basic structures for the V3-based anti-AIDS drug researches have been constructed. Special focus is given to recent studies aimed at defining the structurally conservative V3 sites that may present the HIV-1 weak points most suitable for therapeutic intervention. Finally, the article also discusses how this information can be used to develop novel, potent and broad anti-AIDS agents. EXPERT OPINION: Data on the structure and function of the HIV-1 V3 loop prove convincingly that, in spite of disappointing progress > 20 years, this mysterious site of gp120 (comprising at least three structural motifs recurring in various viral isolates) still remains one of the most attractive targets for the anti-AIDS drug development.

Computational model of the HIV-1 subtype A V3 loop: study on the conformational mobility for structure-based anti-AIDS drug design.

The V3 loop of the HIV-1gp120 glycoprotein presenting 35-residue-long, frequently glycosylated, highly variable, and disulfide bonded structure plays the central role in the virus biology and forms the principal target for neutralizing antibodies and the major viral determinant for co-receptor binding. Here we present the computer-aided studies on the 3D structure of the HIV-1 subtype A V3 loop (SA-V3 loop) in which its structurally inflexible regions and individual amino acids were identified and the structure-function analysis of V3 aimed at the informational support for anti-AIDS drug researches was put into practice. To this end, the following successive steps were carried out: (i) using the methods of comparative modeling and simulated annealing, the ensemble of the low-energy structures was generated for the consensus amino acid sequence of the SA-V3 loop and its most probable conformation was defined basing on the general criteria widely adopted as a measure of the quality of protein structures in terms of their 3D folds and local geometry; (ii) the elements of secondary V3 structures in the built conformations were characterized and careful analysis of the corresponding data arising from experimental observations for the V3 loops in various HIV-1 strains was made; (iii) to reveal common structural motifs in the HIV-1 V3 loops regardless of their sequence variability and medium inconstancy, the simulated structures were collated with each other as well as with those of V3 deciphered by NMR spectroscopy and X-ray studies for diverse virus isolates in different environments; (iv) with the object of delving into the conformational features of the SA-V3 loop, molecular dynamics trajectory was computed from its static 3D structure followed by determining the structurally rigid V3 segments and comparing the findings obtained with the ones derived hereinbefore; and (v) to evaluate the masking effect that can occur due to interaction of the SA-V3 loop with the two virtual molecules constructed previously by tools of computational modeling and named FKBP and CycA peptides, molecular docking of V3 with these molecules was implemented and inter-atomic contacts appearing in the simulated complexes were analyzed to specify the V3 stretches keeping in touch with the ligands. As a matter of record, V3 segments 3-7, 15-20, and 28-32 containing the highly conserved and biologically meaningful residues of gp120 were shown to retain their 3D main chain shapes in all the cases of interest presenting the forward-looking targets for anti-AIDS drug researches. From the data on molecular docking, synthetic analogs of the CycA and FKBP peptides were suggested being suitable frameworks for making a reality of the V3-based anti-HIV-1 drug projects.

Immunophilins and HIV-1 V3 loop for structure-based anti-AIDS drug design.

The model of the structural complex of cyclophilin A (CycA) belonging to the immunophilins family with the HIV-MN gp120 V3 loop was generated, and the computer-aided design of the immunophilin-derived peptide able to mask the biologically crucial V3 segments was implemented. To this end, the following problems were solved: (i) the NMR-based conformational analysis of the HIV-MN V3 loop was put into effect, and its low energy structure fitting the input experimental observations was determined; (ii) molecular docking of this V3 structure with the X-ray conformation of CycA was carried out, and the energy refining the simulated structural complex was performed; (iii) the matrix of inter-atomic distances for the amino acids of the molecules forming part of the built over-molecular ensemble was computed, the types of interactions responsible for its stabilization were analyzed, and the CycA stretch, which accounts for the binding to V3, was identified; (iv) the most probable 3D structure for this stretch in the unbound state was predicted, and its collation with the X-ray structure for the corresponding site of CycA was performed; (v) the potential energy function and its constituents were studied for the structural complex generated by molecular docking of the V3 loop with the CycA peptide offering the virtual molecule that imitates the CycA segment, making a key contribution to the interactions of the native protein with the HIV-1 principal neutralizing determinant; (vi) as a result of the studies above, the designed molecule was shown to be capable of the efficacious blockading the functionally crucial V3 sites; and (vii) based on the joint analysis of the evidence obtained previously and in the present study, the composition of the peptide cocktail presenting the promising anti-AIDS pharmacological substance was developed. The molecules simulated here by molecular modeling methods may become the first representatives of a new class of the chemical compounds (immunophilin-derived peptides) offering the looking-forward basic structures for the design of efficacious and safe antiviral agents.