Crystal Structure Of Protease Research Articles

Design of substituted tetrahydrofuran derivatives for HIV-1 protease inhibitors: synthesis, biological evaluation, and X-ray structural studies.

Substituted tetrahydrofuran derivatives were designed and synthesized to serve as the P2 ligand for a series of potent HIV-1 protease inhibitors. Both enantiomers of the tetrahydrofuran derivatives were synthesized stereoselectivity in optically active forms using lipase-PS catalyzed enzymatic resolution as the key step. These tetrahydrofuran derivatives are designed to promote hydrogen bonding and van der Waals interactions with the backbone atoms in the S2 subsite of the HIV-1 protease active site. Several inhibitors displayed very potent HIV-1 protease inhibitory activity. A high-resolution X-ray crystal structure of an inhibitor-bound HIV-1 protease provided important insight into the ligand binding site interactions in the active site.

Evaluation of therapeutic potentials of some bioactive compounds in selected African plants targeting main protease (Mpro) in SARS-CoV-2: a molecular docking study

BackgroundCoronavirus disease 2019 (COVID-19) is an infectious disease brought on by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a global treat in early 2020. Despite worldwide research proving different medications used to treat COVID-19, the infection still affects the human race; we need to continue researching the virus to protect humanity and reduce the complications that some medications might cause. This study focuses on finding another promising therapeutic compound against SARS-CoV-2. Twenty-four (24) bioactive compounds were selected from the following African plants' Adansonia digitata L, Aframomum melegueta K. Schum, Ageratum conyzoides (L.) L, and Boswellia dalzielii, and Remdesivir was used as the control medication. The PubChem web server acquired the 3D structures of bioactive compounds in the plant and the control medication. The SARS-CoV-2 main protease (Mpro) crystal structure was obtained using the Protein Data Bank (PDB). Using the SwissADME web server, the bioactive compounds' drug-likeness was assessed, and AutoDock was employed for the molecular docking with the Mpro. The Proteins Plus and Protein–Ligand Interaction Profiler web servers were used to analyse the docked complexes. Furthermore, the admetSAR website was utilized to predict the ligands' absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties.ResultsBased on the drug-likeness screening, Rutin violated more than one of the Lipinski rules of five, while Remdesivir violated two. Molecular docking analysis results indicated that Catechin, Epicatechin, Vitexin, Quercetin, Kaempferol, Gamma-Sitosterol, and Kaur-16-ene exhibited a stronger binding affinity with Mpro, with binding scores of − 7.1, − 7.1, − 8.0, − 7.3, − 7.2, − 6.8, and − 6.5 kcal/mol, respectively, compared to Remdesivir's binding score of − 6.3 kcal/mol. Consequently, binding scores of bioactive compounds suggest their potential biological activity against the SARS-CoV-2 main protease. Additionally, these bioactive compounds exhibited favourable ADMET properties. Vitexin also has a plasma protein binding below 90%, a promising medication distribution feature.ConclusionsThis study shows that Catechin, Epicatechin, Vitexin, Quercetin, Kaempferol, Gamma-Sitosterol, and Kaur-16-ene have better binding affinities with Mpro than Remdesivir. Molecular dynamics simulation in vitro and in vivo investigation is required to support this study.

Molecular Docking Study: Identification of Potential Inhibitors in Lavandula Angustifolia Essential Oil Against the Main Protease of Sars-Cov-2

SARS-CoV-2, also known as, Covid-19 is an almost globally endemic disease that has infected hundreds of millions of people and caused ongoing deaths. Researchers are still trying to search for effective drugs against Covid-19. Currently, with the development of science and technology, the process of drug discovery through simulated environments is becoming increasingly accurate and yielding rapid results. The protease crystal structure of Covid-19 plays a crucial role in viral replication and represents promising primary targets for antiviral drug discovery and therapy. Lavender and N3 compound ligands were used as test ligands. A total of 9 lavender ligand compounds were subjected to testing using the molecular docking method. Molecular docking aims to achieve the optimal conformation of both the protein and ligand, as well as the relative orientation between them. The analysis conducted during molecular docking includes the determination of the inhibition constant and the study of interactions between the protein and ligand.Keywords: Covid-19, docking, N3, ligand

Fourteen immunomodulatory alkaloids and two prenylated phenylpropanoids with dual therapeutic approach for COVID-19: molecular docking and dynamics studies

The pandemic outbreak of COVID-19 caused by the new severe acute respiratory syndrome coronavirus (SARS-CoV-2) is a global health burden. To date, there is no highly effective antiviral therapy to eradicate the virus; as a result, researchers are racing to introduce new potential therapeutic agents. Alternatively, traditional immunity boosters and symptomatic treatment based on natural bioactive compounds are also an option. The 3-chymotrypsin-like protease (3CLpro) crystal structure, the main proteolytic enzyme of SARS-CoV-2, has been unraveled, allowing the development of effective protease inhibitors via in silico and biological studies. In COVID-19 infected patients, the loss of lung function, and mortality are reported to be linked to several inflammatory mediators and cytokines. In this context, the approach of introducing immunomodulatory agents may be considered a dual lifesaving strategy in combination with antiviral drugs. This study aims to provide immunomodulatory natural products exhibiting potential protease inhibitory activities. Selected groups of alkaloids of different classes and two prenylated phenylpropanoids from the Brazilian green propolis were in silico screened for their ability to inhibit COVID-19 3CLpro protease. Results showed that compounds exhibited binding energy scores with values ranging from −6.96 to −3.70 compared to the reference synthetic protease inhibitor O6K with a binding energy score of −7.57. O6K binding energy was found comparable with lead phytochemicals in our study, while their toxicity and drug-likeness criteria are better than that of O6K. The activities of these molecules are mainly ascribed to their ability to form hydrogen bonding with 3CLpro crucial amino acid residues of the catalytic site. In addition, the molecular dynamics simulations further showed that some of these compounds formed stable complexes as evidenced by the occupancy fraction measurements. The study suggested that the major immunomodulators 3β, 20α-diacetamido-5α-pregnane, (20S)-(benzamido)-3β-(N,N-dimethyamino)-pregnane, and baccharin are 3CLpro inhibitors. Biological screenings of these phytochemicals will be valuable to experimentally validate and consolidate the results of this study before a rigid conclusion is reached, which may pave the way for the development of efficient modulatory bioactive compounds with dual bioactions in COVID-19 intervention. Communicated by Ramaswamy H. Sarma

Virtual and In Vitro Screening of Natural Products Identifies Indole and Benzene Derivatives as Inhibitors of SARS-CoV-2 Main Protease (Mpro).

The rapid spread of the coronavirus disease 2019 (COVID-19) resulted in serious health, social, and economic consequences. While the development of effective vaccines substantially reduced the severity of symptoms and the associated deaths, we still urgently need effective drugs to further reduce the number of casualties associated with SARS-CoV-2 infections. Machine learning methods both improved and sped up all the different stages of the drug discovery processes by performing complex analyses with enormous datasets. Natural products (NPs) have been used for treating diseases and infections for thousands of years and represent a valuable resource for drug discovery when combined with the current computation advancements. Here, a dataset of 406,747 unique NPs was screened against the SARS-CoV-2 main protease (Mpro) crystal structure (6lu7) using a combination of ligand- and structural-based virtual screening. Based on 1) the predicted binding affinities of the NPs to the Mpro, 2) the types and number of interactions with the Mpro amino acids that are critical for its function, and 3) the desirable pharmacokinetic properties of the NPs, we identified the top 20 candidates that could potentially inhibit the Mpro protease function. A total of 7 of the 20 top candidates were subjected to in vitro protease inhibition assay and 4 of them (4/7; 57%), including two beta carbolines, one N-alkyl indole, and one Benzoic acid ester, had significant inhibitory activity against Mpro protease. These four NPs could be developed further for the treatment of COVID-19 symptoms.

In silico screening of potential compounds from begonia genus as 3CL protease (3Cl pro) SARS-CoV-2 inhibitors.

The emergence of Coronavirus disease (COVID-19) has been declared a pandemic and made a medical emergency worldwide. Various attempts have been made, including optimizing effective treatments against the disease or developing a vaccine. Since the SARS-CoV-2 protease crystal structure has been discovered, searching for its inhibitors by in silico technique becomes possible. This study aims to virtually screen the potential of phytoconstituents from the Begonia genus as 3Cl pro-SARS-CoV- 2 inhibitors, based on its crucial role in viral replication, hence making these proteases "promising" for the anti-SARS-CoV-2 target. In silico screening was carried out by molecular docking on the web-based program DockThor and validated by a retrospective method. Predictive binding affinity (Dock Score) was used for scoring the compounds. Further molecular dynamics on Desmond was performed to assess the complex stability. Virtual screening protocol was valid with the area under curve value 0.913. Molecular docking revealed only β-sitosterol- 3-O-β-D-glucopyranoside with a lower docking score of - 9.712 kcal/mol than positive control of indinavir. The molecular dynamic study showed that the compound was stable for the first 30 ns simulations time with Root Mean Square Deviation <3 Å, despite minor fluctuations observed at the end of simulation times. Root Mean Square Fluctuation of catalytic sites HIS41 and CYS145 was 0.756 Å and 0.773 Å, respectively. This result suggests that β-sitosterol-3-O-β-Dglucopyranoside might be a prospective metabolite compound that can be developed as anti-SARS-CoV-2.

Drug Resistance in the HIV-1 Subtype C Protease Enzyme: A High Throughput Virtual Screening Approach in Search of New Ligands with Activity.

HIV-1 subtype C protease is a strategic target for antiretroviral treatment. However, resistance to protease inhibitors appears after months of treatment. Chromones and 2- biscoumarin derivatives show potential for inhibition of the HIV- subtype C protease. Different heterocyclic structures from the ZINC database were docked against Human Immunodeficiency Virus-1 (HIV) subtype C protease crystal structure 2R5Q and 2R5P. The 5 best molecules were selected to be docked against 62 homology models based on HIV-protease sequences from infants failing antiretroviral protease treatment. This experimentation was performed with two molecular docking programs: Autodock and Autodock Vina. These molecules were modified by substituting protons with different moieties, and the derivatives were docked against the same targets. Ligand-protein interactions, physical/chemical proprieties of the molecules, and dynamics simulations were analyzed. Docking of all of the molecules was performed to find out the binding sites of HIV-1 subtype C proteases. An in-house script was made to substitute protons of molecules with different moieties. According to the Lipinski rule of five, physical and chemical properties were determined. Complexes of certain ligands-protease were compared to the protein alone in molecular dynamics simulations. From the first docking results, the 5 best (lowest energy) ligands (dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[a]acridine) were selected. The binding energy of the modified ligands increased, including the poorest-performing molecules. A correlation between nature, the position, and the resulting binding energy was observed. According to the Lipinski rules, the physico-chemical characteristics of the five best-modified ligands are ideal for oral bioavailability. Molecular dynamics simulations show that some lead-protease complexes were stable. Dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[ a]acridine and their derivatives might be considered as promising HIV-1 subtype C protease inhibitors. This could be confirmed through synthesis and subsequent in vitro assays.

Prediction of SARS-CoV-2 3C-like protease (3CLpro) crystal structure to provide COVID-19 inhibitor design through computational studies

Infectious diseases have lately become pandemic, posing a threat to global public health with the introduction of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously provisionally named 2019 novel coronavirus or 2019-nCoV). Technological advancements have increased the possibility of discovering natural inhibitor candidates capable of preventing and controlling COVID-19 infections. The SARS-CoV-2 3C-like protease (3CLpro) is critical for SARS-CoV-2 replication and is a prospective therapeutic target. This study aims to identify, evaluate, and explore the 3CLpro macromolecular structures from SARS-CoV and SARS-CoV-2, as well as their impact on angiotensin-converting enzyme 2 (ACE-2). The discovery of the two 3CLpro macromolecules revealed structural similarities in several regions. These findings were subsequently confirmed by performing protein-protein docking simulations to observe the interaction of 3CLpro with the active site ACE-2. With an ACE score of 701.41 kJ/mol, SARS-COV-2 3CLpro forms the strongest binding with ACE-2. As a result, the findings of this research can be used to guide the development of potential SARS-CoV-2 3CLpro inhibitors for the treatment of COVID-19 infectious diseases.

Prospecting for Cressa cretica to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2

The severe acute respiratory syndrome COVID-19 declared as a global pandemic by the World Health Organization has become the present wellbeing worry to the whole world. There is an emergent need to search for possible medications. Cressa cretica is reported to show antitubercular, antibacterial and expectorant property. In this research, we aim to prospect the COVID-19 main protease crystal structure (Mpro; PDB ID: 6LU7) and the active chemical constituents from Cressa cretica in order to understand the structural basis of their interactions. We examined the binding potential of active constituents of Cressa cretica plant to immensely conserved protein Mpro of SARS-CoV-2 followed by exploration of the vast conformational space of protein–ligand complexes by molecular dynamics (MD) simulations. The results suggest the effectiveness of 3,5-Dicaffeoylquinic acid and Quercetin against standard drug Remdesivir. The active chemical constituents exhibited good docking scores, and interacts with binding site residues of Mpro by forming hydrogen bond and hydrophobic interactions. 3,5-Dicaffeoylquinic acid showed the best affinity towards Mpro receptor which is one of the target enzymes required by SARS CoV-2 virus for replication suggesting it to be a novel research molecule. The potential of the active chemical constituents from Cressa cretica against the SARS-CoV-2 virus has best been highlighted through this study. Therefore, these chemical entities can be further scrutinized and provides direction for further consideration for in-vivo and in-vitro validations for the treatment of covid-19. Communicated by Ramaswamy H. Sarma

An automated protocol for modelling peptide substrates to proteases

BackgroundProteases are key drivers in many biological processes, in part due to their specificity towards their substrates. However, depending on the family and molecular function, they can also display substrate promiscuity which can also be essential. Databases compiling specificity matrices derived from experimental assays have provided valuable insights into protease substrate recognition. Despite this, there are still gaps in our knowledge of the structural determinants. Here, we compile a set of protease crystal structures with bound peptide-like ligands to create a protocol for modelling substrates bound to protease structures, and for studying observables associated to the binding recognition.ResultsAs an application, we modelled a subset of protease–peptide complexes for which experimental cleavage data are available to compare with informational entropies obtained from protease–specificity matrices. The modelled complexes were subjected to conformational sampling using the Backrub method in Rosetta, and multiple observables from the simulations were calculated and compared per peptide position. We found that some of the calculated structural observables, such as the relative accessible surface area and the interaction energy, can help characterize a protease’s substrate recognition, giving insights for the potential prediction of novel substrates by combining additional approaches.ConclusionOverall, our approach provides a repository of protease structures with annotated data, and an open source computational protocol to reproduce the modelling and dynamic analysis of the protease–peptide complexes.

Effect of foot-and-mouth disease virus 3C protease B2 β-strand proline mutagenesis on expression and processing of the P1 polypeptide using a plasmid expression vector.

The production of experimental molecular vaccines against foot-and-mouth disease virus utilizes the viral encoded 3C protease for processing of the P1 polyprotein. Expression of wild type 3C protease is detrimental to host cells. The molecular vaccine constructs containing the 3C protease L127P mutant significantly reduce adverse effects associated with protease expression while retaining the ability to process and assemble virus-like particles. In published 3C protease crystal structures, the L127 residue is contained within the B2 β-strand as part of the A2-B2 β-sheet. To provide insight into the mechanism by which the L127P mutant alters the properties of the 3C protease, we performed scanning proline mutagenesis of residues 123-128 of the B2 β-strand and monitored expression and P1 processing. Simultaneously, we utilized random mutagenesis of the full 3C sequence to identify additional mutations presenting a phenotype similar to the L127P mutation. Six of the tested mutants enhanced expression over wild type, and the I22P, T100P and V124P mutations surpassed the L127P mutation in certain cell lines. These data areinterpreted in conjunction with published 3C protease crystal structures to provide insight into the mechanism by which these mutations enhance expression.

Inhibitors of flavivirus by targeting the NS2B‐NS3 protease interface

Flaviviruses require the activity of their viral protease in order to process the single polyprotein encoded by their genome into individual proteins and assemble an infectious virion. This protease consists of a large catalytic subunit (NS3) and a smaller accessory subunit (NS2B) important for function. Using the previously solved flavivirus protease crystal structures, we identified three putative sites of interaction between NS2B and NS3 which we then queried with a virtual screen of NCI library compounds, followed by subsequent in vitro biochemical screening of the resultant candidates. Of more than one hundred initial candidate compounds that passed the virtual screen, thirty‐six demonstrated significant inhibition in vitro. Three most potent compounds, 227186, 73101, and 65828 were further investigated and found to significantly reduce titers of Dengue virus 2 on A549 cells.Support or Funding InformationAI133219, AI134568, AI131669This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Protein backbone flexibility pattern is evolutionarily conserved in the Flaviviridae family: A case of NS3 protease in Flavivirus and Hepacivirus

Viruses belonging to the Flaviviridae family have been an important health concern for humans, animals and birds alike. No specific treatment is available yet for many of the viral infections caused by the members of this family. Lack of specific drugs against these viruses is mainly due to lack of protein structure information. It has been known that protein backbone fluctuation pattern is highly conserved in protein pairs with similar folds, in spite of the lack of sequence similarity. We hypothesized that this concept should also hold true for proteins (especially enzymes) of viruses included in different genera of the Flaviviridae family, as we know that the sequence similarity between them is low. Using available NS3 protease crystal structures of the Flaviviridae family, our preliminary results have shown that the Cα (i.e. backbone) fluctuation patterns are highly similar between Flaviviruses and a Hepacivirus (i.e. hepatitis C virus, HCV). This has to be validated further experimentally.

In silico development of a novel putative inhibitor of the 3C protease of Coxsackievirus B3 with a benzene sulfonamide skeleton

Availability of X-ray crystal structure of 3C protease of several enteroviruses provided an opportunity for in silico drug design and development approach. Presented study is aimed at designing a novel compound targeting 3C protease of Coxsackievirus (CVB3), which is reported frequently to cause myocarditis in North America and Europe. A pthalimido-sulfonamide derivative (ZINC13799063) was identified through high-throughput virtual screening (HTVS) approach from the top HITs. A small library of phalimido-sulphonamides was enumerated to find a potential LEAD. Compound 17 from the library was found to inhibit CVB3 selectively in cell based antiviral assay at a concentration of EC50=1.0±0.1 µM with a selectivity index of &gt;140. Molecular dynamics study was performed to investigate the selective inhibition of CVB3 over CVB4.

Simple Design of an Enzyme-Inspired Supported Catalyst Based on a Catalytic Triad

Summary Enzyme active sites afford an intricate interplay of functional groups to mediate complex organic and inorganic reactions. Many hydrolytic enzymes use a catalytic triad comprising three different functional residues—(Ser(-OH), Hist(-imidazole), Asp(-CO 2 H))—that catalyze the hydrolysis of numerous unique substrates. Inspired by this design, we have developed a simple one-step synthesis for preparing a new supported catalytic system in which the three reactive groups of the catalytic triad (alcohol, imidazole, and carboxylate) are incorporated into a single functional unit. These artificial active sites can be coupled to a solid-phase support (Merrifield resin) by copper(I)-catalyzed azide-alkyne cycloaddition "click chemistry," and their effectiveness as esterolysis catalysts was demonstrated. Furthermore, tuning the local hydrophobicity of the resin particles with an approach analogous to the native enzyme hydrophobic pocket increased the catalytic efficiency. Quantum mechanics and molecular dynamics computational modeling were used to probe the catalytic effect and suggested a concerted two-step mechanism and hydrophobic nanoenvironment similar to that of hydrolytic enzymes.

The New High Resolution Crystal Structure of NS2B-NS3 Protease of Zika Virus.

Zika virus (ZIKV) is the cause of a significant viral disease affecting humans, which has spread throughout many South American countries and has also become a threat to Southeastern Asia. This commentary discusses the article “Crystal structure of unlinked NS2B-NS3 protease from Zika virus” published recently in the journal Science by Zhang et al. of Nanyang Technological University, Singapore. They resolved a 1.58 Å resolution structure of the NS2B-NS3 protease of ZIKV and demonstrated how peptide and non-peptide inhibitors interact with this structure, along with the different conformational states that were observed. This protease crystal structure offers new opportunities for the design and development of novel antiviral drugs used for the treatment and control of ZIKV.

The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Noncovalent Interactions.

We propose a methodology to derive pairwise-additive noncovalent force fields from monomer electron densities without any empirical input. Energy expressions are based on the symmetry-adapted perturbation theory (SAPT) decomposition of interaction energies. This ensures a physically motivated force field featuring an electrostatic, exchange-repulsion, dispersion, and induction contribution, which contains two types of parameters. First, each contribution depends on several fixed atomic parameters, resulting from a partitioning of the monomer electron density. Second, each of the last three contributions (exchange-repulsion, dispersion, and induction) contains exactly one linear fitting parameter. These three so-called interaction parameters in the model are initially estimated separately using SAPT reference calculations for the S66x8 database of noncovalent dimers. In a second step, the three interaction parameters are further refined simultaneously to reproduce CCSD(T)/CBS interaction energies for the same database. The limited number of parameters that are fitted to dimer interaction energies (only three) avoids ill-conditioned fits that plague conventional parameter optimizations. For the exchange-repulsion and dispersion component, good results are obtained for all dimers in the S66x8 database using one single value for the associated interaction parameters. The values of those parameters can be considered universal and can also be used for dimers not present in the original database used for fitting. For the induction component such an approach is only viable for the dispersion-dominated dimers in the S66x8 database. For other dimers (such as hydrogen-bonded complexes), we show that our methodology remains applicable. However, the interaction parameter needs to be determined on a case-specific basis. As an external validation, the force field predicts interaction energies in good agreement with CCSD(T)/CBS values for dispersion-dominated dimers extracted from an HIV-II protease crystal structure with a bound ligand (indinavir). Furthermore, experimental second virial coefficients of small alkanes and alkenes are well reproduced.

Molecular recognition mechanism and motion of HCV NS3/4A protease with Faldaprevir analogue

Faldaprevir analogue molecule (FAM) has been reported to effectively inhibit the catalytic activity of HCV NS3/4A protease, making it a potential lead compound against HCV. A series of HCV NS3/4A protease crystal structures were analyzed by bioinformatics methods, and the FAM-HCV NS3/4A protease crystal structure was chosen for this study. A 20.4 ns molecular dynamics simulation of the complex consists of HCV NS3/4A protease and FAM was conducted. The key amino acid residues for interaction and the binding driving force for the molecular recognition between the protease and FAM were identified from the hydrogen bonds and binding free energy analyses. With the driving force of hydrogen bonds and van der Waals, FAM specifically bind to the active pocket of HCV NS3/4A protease, including V130-S137, F152-D166, D77-D79 and V55, which agreed with the experimental data. The effect of R155K, D168E/V and V170T site-directed mutagenesis on FAM molecular recognition was analyzed for their effect on drug resistance, which provided the possible molecular explanation of FAM resistance. Finally, the system conformational change was explored by using free energy landscape and conformational cluster. The result showed four kinds of dominant conformation, which provides theoretical basis for subsequent design of Faldaprevir analogue inhibitors based on the structure of HCV NS3/4A protease.

A combined cheminformatic and bioinformatic approach to address the proteolytic stability challenge in peptide-based drug discovery.

We have created models to predict cleavage sites for several human proteases including caspase-1, caspase-3, caspase-6, caspase-7, cathepsin B, cathepsin D, cathepsin G, cathepsin K, cathepsin L, elastase-2, granzyme A, granzyme B, matrix metallopeptidase-2 (MMP2), MMP7, MMP9, thrombin, and trypsin-1. Rather than representing the sequence pattern around the potential cleavage site through a series of flags with each flag representing one of the 20 standard amino acids, we first represent each amino acid by its calculated properties. For these calculated properties, we use validated cheminformatic descriptors, such as molecular weight, logP, and polar surface area, of the individual amino acids. Finally, the cleavage site-specific descriptors are calculated through various combinations of the individual amino acid descriptors for the residues surrounding the cleavage site. Some of these combinations do not take into account the location of the residue, as long as it is in a prescribed neighborhood of the potential cleavage site, whereas others are sensitive to the precise order of the residues in the sequence. The key advantage of this approach is that it allows one to perform meaningful calculations with nonstandard amino acids for which little or no data exists. Finally, using both docking and molecular dynamics simulations, we examine the potential for and limitations of protease crystal structures to impact the design of proteolytically stable peptides.

Structure-based design of potent HIV-1 protease inhibitors with modified P1-biphenyl ligands: synthesis, biological evaluation, and enzyme-inhibitor X-ray structural studies.

We report the design, synthesis, X-ray structural studies, and biological evaluation of a novel series of HIV-1 protease inhibitors. We designed a variety of functionalized biphenyl derivatives to make enhanced van der Waals interactions in the S1 subsite of HIV-1 protease. These biphenyl derivatives were conveniently synthesized using a Suzuki-Miyaura cross-coupling reaction as the key step. We examined the potential of these functionalized biphenyl-derived P1 ligands in combination with 3-(S)-tetrahydrofuranyl urethane and bis-tetrahydrofuranyl urethane as the P2 ligands. Inhibitor 21e, with a 2-methoxy-1,1'-biphenyl derivative as P1 ligand and bis-THF as the P2 ligand, displayed the most potent enzyme inhibitory and antiviral activity. This inhibitor also exhibited potent activity against a panel of multidrug-resistant HIV-1 variants. A high resolution X-ray crystal structure of related Boc-derivative 17a-bound HIV-1 protease provided important molecular insight into the ligand-binding site interactions of the biphenyl core in the S1 subsite of HIV-1 protease.