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.

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.

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.

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.

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.

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.

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.

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.

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.

Determining the invariant structure elements of the HIV-1 variable V3 loops: insight into the HIV-MN and HIV-Haiti isolates.

The computer approaches that combined the 3D protein structure modeling with the mathematical statistics methods were used to compute the NMR-based 3D structures of the HIV-1 gp120 V3 loop for the HIV-MN and HIV-Haiti isolates in water as well as to compare their conformational characteristics with the purpose of determining the structure elements common for the two virus modifications. As a result, the variability of the amino acid sequence was found to stimulate the considerable structural rearrangements of the V3 loop. However, despite this fact, one functionally crucial stretch of V3 and a greater portion of its residues were shown to preserve the conformations in the viral strains of interest. To reveal the structural motifs and individual amino acids giving rise to the close conformations in the HIV-MN and HIV-Haiti V3 loops regardless of the sequence and environment variability, the simulated structures were collated with those deciphered previously in terms of NMR data in a water/trifluoroethanol mixed solvent. The structure elements and single residues of V3 residing in its biologically significant sites and keeping the conformations in all of the cases at question are considered to be the promising targets for anti-AIDS drugs studies. In this context, the structurally inflexible motifs of V3 presenting the weak units in the virus protection system may be utilized as the most convenient landing-places for molecular docking of the V3 loop and ligand structures followed by selecting chemical compounds suitable as a basis for the design of safe and effective antiviral agents.

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.

Dual spatial folds and different local structures of the HIV-1 immunogenic crown in various virus isolates.

Local and global structural properties of the HIV-1 principal neutralizing epitope were studied in terms of NMR spectroscopy data reported in literature for the HIV-Haiti and HIV-RF isolates. To this effect, the NMR-based method comprising a probabilistic model of protein conformation in conjunction with the molecular mechanics and quantum chemical computations was used for determining the ensembles of conformers matching the NMR requirements and energy criteria. As a matter of record, the high resolution 3D structure models were constructed for the HIV-Haiti and HIV-RF immunogenic crowns, and their geometric parameters were collated with the ones of conformers derived previously for describing the conformational features of immunogenic tip of gp120 from Thailand and MN HIV-1 strains. The HIV-1 neutralization site was demonstrated to constitute in water solution highly flexible system sensitive to its environment. This inference is completely valid for the geometric space of dihedral angles where statistically significant differences in local structures of simulated conformers have been found for all virus isolates of interest. In spite of this fact, the stretch analyzed was shown to manifest a certain conservatism in the space of atomic coordinates, building up in four HIV-1 isolates two spatial folds similar to those observed in crystal for the V3 loop peptides bound to different neutralizing Fabs. The results are discussed in the light of literature data on HIV-1 neutralizing epitope structure.

Local structural properties of the V3 loop of Thailand HIV-1 isolate.

The model of locally accurate conformation for the HIV-Thailand principal neutralizing determinant (PND) located within the V3 loop of the virus envelope protein gp120 was built in terms of NMR spectroscopy data. To this end, the NMR-based conformational analysis of synthetic molecule representing the peptide copy of the fragment under study was carried out using the published sequential d connectivity data and values of spin-spin coupling constants. As a result, (i) the local structure for the V3 loop from Thailand isolate was determined, (ii) the conformations of its irregular segments were analyzed, and the secondary structure elements identified, (iii) the ensemble of conformers matching the experimental and theoretical data was derived for the stretch forming the neutralizing epitope of the HIV-Thailand PND, (iv) to estimate the probability of realizing each of these conformers in solution, the results obtained were collated with the X-ray data for corresponding segments in synthetic molecules imitating the central region of the HIV-MN PND as well as for homologous segments 39-44 in Bence-Jonce REI protein (BJRP), 41-46 in immunoglobulin lambda (Ig lambda), and 50-55 in beta-chain of horse hemoglobin (HH), (v) to find the conserved structural motifs inside diverse HIV-1 isolates, the structure determined was compared with the one derived earlier for the HIV-MN PND from NMR spectroscopy data, (vi) on the basis of all data obtained, the 3D structure model describing the set of biologically relevant conformations, which may present different antigenic determinants to the immune system in various HIV-1 isolates, was proposed for the immunogenic crown of the V3 loop. The results obtained are discussed in conjunction with the data on the structure for the HIV-1 PND reported in literature.

3D structure model of the principal neutralizing epitope of Minnesota HIV-1 isolate.

A hierarchical procedure, using a "bottom-up" strategy and combining (i). a probabilistic approach for estimating all possible starting structures, (ii). restrained molecular mechanics algorithms for preliminary selection of all energetically preferred conformers, as well as (iii). quantum chemical computations for refining their geometry, was used to study the structural properties of the HIV-MN neutralizing epitope in terms of NMR spectroscopy data. As a result, only one of initial structures matching the experimental and theoretical data was found to be well-ground for implementing the function of immunoreactive conformation of the virus immunogenic crown. The geometric parameters of this structure in water solution were shown to correspond to a double beta-turn conformation similar to that revealed in crystal for synthetic molecules imitating the central region of the HIV-MN V3 loop. The following conclusion was drawn from the comparative analysis of simulated structure with the one computed previously: the HIV-MN immunogenic tip has some inherent conformational flexibility that manifests at the alterations of hexapeptide environment and leads to the structural transitions changing the local conformation of the stretch of interest but retaining its spatial main chain fold. As a matter of record, the high resolution 3D structure model for the HIV-MN principal neutralization site was constructed, and its geometric parameters were compared with the corresponding characteristics of conformers derived earlier for describing the conformational features of immunogenic tip of gp120 from Thailand HIV-1 isolate.

Structure and polymorphism of the principal neutralization site of Thailand HIV-1 isolate.

Refining the geometric parameters for the ensemble of conformers, derived earlier in terms of NMR-spectroscopy data for the immunogenic tip of Thailand HIV-1 isolate, was carried out by quantum chemical methods. As a result, (i) the energy characteristics of initial structures were significantly improved, (ii) their relative locations on the scale of formation heats were determined, and (iii) the energy barriers between conformers under study were computed. On the basis of all data obtained, the high resolution 3D structure model, describing the set of stable conformers and containing the biologically active conformation, was proposed for neutralizing epitope of Thailand HIV-1 isolate. The following major conclusions were made based on the analysis of simulated conformations: i) the Gly-Pro-Gly-Gln-Val-Phe stretch forming the immunogenic crown of Thailand HIV-1 isolate exhibits the properties characteristic for metastable oligopeptide that constitutes in solution the dominant structure with other conformations admissible; (ii) three structures out of five NMR-based starting models form the cluster of conformers which adequately describes general conformational features of this functionally important site of gp120; (iii) two structures residing in this cluster are found to be well-ground for implementing the function of immunoreactive conformation of the stretch of interest; (iv) in spite of this observation, the "global" structure which gives rise to inverse gamma-turn in the central Gly-Pro-Gly crest of Thailand HIV-1 gp120 is proposed to be the most probable conformation responsible for the formation of viral antigen-antibody complex in particular case under study.