An alpha-helical peptidomimetic inhibitor of the HIV-1 Rev-RRE interaction

The interaction between the HIV-1 Rev protein and the Rev-Responsive Element (RRE) RNA is an attractive target for anti-viral therapy. We have designed alpha-helical peptidomimetics of Rev-like peptides using side chain-side chain macrolactam formation between positions i and i+4. One peptidomimetic having an appropriate location, orientation, and length of the macrolactam exhibited both significant helical character and specific RRE binding. This molecule displays 2-fold greater RNA specificity than the wild-type Rev peptide and more than 20-fold greater specificity than an uncyclized control peptide. Thus, specific, high affinity recognition of the RRE is feasible utilizing a small, relatively rigid peptidomimetic scaffold.     

Homology model-guided 3D-QSAR studies of HIV-1 integrase inhibitors

In the present study, we report the exploration of binding modes of potent HIV-1 integrase (IN) inhibitors MK-0518 (raltegravir) and GS-9137 (elvitegravir) as well as chalcone and related amide IN inhibitors we recently synthesized and the development of 3D-QSAR models for integrase inhibition. Homology models of DNA-bound HIV-1 IN were constructed on the basis of the X-ray crystal structure of the foamy virus IN-DNA complex (PDB ID: 3L2T ) and used for docking. The binding modes of raltegravir and elvitegravir in our homology models are in accordance with those in the foamy virus structure revealing interactions important for inhibitor-IN binding. To gain further insights into the structural requirements for IN inhibition, three-dimensional quantitative structure activity relationship (3D-QSAR) studies were conducted using raltegravir, elvitegravir, and their analogs; our synthesized 3-keto salicylic acid IN inhibitor series; as well as other structurally related HIV-1 IN inhibitors. In the first part of the study with 103 compounds, atom-fit alignments, I and II, and docking-based alignment, III, were used to develop 3D-QSAR models 1, 2, and 3, respectively, each comprising comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSARs. This initial analysis indicated that the docking-based (structure-based) model 3 performed better than the atom-fit (ligand-based) models 1 and 2, in terms of statistical significance and robustness. Thus, the docking-based alignment was then subsequently used with an expanded data set of 296 compounds for building a more comprehensive 3D-QSAR, model 4. Model 4 afforded good q2 values of 0.70 and 0.75 for CoMFA and CoMSIA 3D-QSARs, respectively, and showed good predictive performance on an external validation test set of 59 compounds with predictive r2 values up to 0.71. The HIV IN-DNA homology model of biological relevance and the comprehensive 3D-QSAR models developed in the present study provide insights and new predictive tools for structure-based design and optimization of IN inhibitors.     

HIV-1外膜蛋白在毕赤酵母中的构建

为了进一步研究HIV-1外膜蛋白的结构和功能,获得足够量的糖蛋白,对HIV-1的外膜蛋白Env进行了修饰,利用PCR技术克隆目的基因,把Env基因构建到毕赤酵母Pichia pastoris表达载体中,把pPICZαB-Env表达载体转化到酵母中,根据基因重组技术,使Env基因和酵母基因组重组,筛选阳性重组子.     

Synthesis and stability evaluation of novel peptidomimetic Caspase-1 inhibitors for topical application

During our search for topically-active Caspase-1 inhibitors, we identified a novel class of potent inhibitors based on a 1,3,5-trisubstituted uracil motif equipped with an l-aspartate semi-aldehyde derived warhead. In the literature, the majority of Caspase-1 inhibitors possessing the same warhead have been designed and evaluated for oral administration as the ethyl acetal pro-drug form. For our topical program, the pro-drug acetal form was not fully hydrolysed in the skin and was unstable in many of our standard topical excipients, therefore, we were obliged to focus on the actual hemiacetal drug form of the molecule during our drug discovery program. Our work focuses on both the synthesis and achiral and chiral stability of the final drug molecules in topical excipients.     

Theoretical studies on the interactions and interferences of HIV-1 glycoprotein gp120 and its coreceptor CCR5

The interaction between the HIV gp120 protein and coreceptor CCR5 or CXCR4 of the host cell is critical in mediating the HIV entry process. A model for the CCR5-gp120 complex has been developed. In the model, the N-terminus of CCR5 binds to three discontinuous domains of gp120, including the fourth conserved (C4) region, β19/β20 connecting loop, and V3 loop. The second extra-cellular loop (ECL2) of CCR5 also interacts with the crown part of the gp120 V3 loop. The bindings of the three CCR5 antagonists, maraviroc, aplaviroc, and vicriviroc, to the trans-membrane domain of CCR5 have been modeled. The bindings are found to affect the conformation of the ECL2 domain, which in turn drives the N-terminus of CCR5 to an altered state. Aplaviroc is more hydrophilic than maraviroc and vicriviroc, and its binding is more interfered by solvent, resulting in a quite different effect to the structure of CCR5 compared with those of the other two molecules. The above results are in accord with experimental observations and provide a structural basis for further design of CCR5 antagonists.     

Pharmacomodulation of a ligand targeting the HBV capsid hydrophobic pocket

Hepatitis B virus (HBV) is a small enveloped retrotranscribing DNA virus and an important human pathogen. Its capsid-forming core protein (Cp) features a hydrophobic pocket proposed to be central notably in capsid envelopment. Indeed, mutations in and around this pocket can profoundly modulate, and even abolish, secretion of enveloped virions. We have recently shown that Triton X-100, a detergent used during Cp purification, binds to the hydrophobic pocket with micromolar affinity. We here performed pharmacomodulation of pocket binders through systematic modifications of the three distinct chemical moieties composing the Triton X-100 molecule. Using NMR and ITC, we found that the flat aromatic moiety is essential for binding, while the number of atoms of the aliphatic chain modulates binding affinity. The hydrophilic tail, in contrast, is highly tolerant to changes in both length and type. Our data provide essential information for designing a new class of HBV antivirals targeting capsid–envelope interactions.

Small-molecule binding to the Hepatitis B virus core protein hydrophobic pocket, a possible strategy for targeting viral particle assembly.     

Design, synthesis, and computational studies of inhibitors of Bcl-XL

One of the primary objectives in the design of protein inhibitors is to shape the three-dimensional structures of small molecules to be complementary to the binding site of a target protein. In the course of our efforts to discover potent inhibitors of Bcl-2 family proteins, we found a unique folded conformation adopted by tethered aromatic groups in the ligand that significantly enhanced binding affinity to Bcl-XL. This finding led us to design compounds that were biased by nonbonding interactions present in a urea tether to adopt this bioactive, folded motif. To characterize the key interactions that induce the desired conformational bias, a series of substituted N,N'-diarylureas were prepared and analyzed using X-ray crystallography and quantum mechanical calculations. Stabilizing pi-stacking interactions and destabilizing steric interactions were predicted to work in concert in two of the substitution patterns to promote the bioactive conformation as a global energy minimum and result in a high target binding affinity. Conversely, intramolecular hydrogen bonding present in the third substitution motif promotes a less active, extended conformer as the energetically favored geometry. These findings were corroborated when the inhibition constant of binding to Bcl-XL was determined for fully elaborated analogues bearing these structural motifs. Finally, we obtained the NMR solution structure of the disubstituted N,N'-diarylurea bound to Bcl-XL demonstrating the folded conformation of the urea motif engaged in extensive pi-interactions with the protein.     

HIV-1逆转录酶抑制剂2-氨基6-苯磺酰基苄腈及其类似物活性OSAR研究

首先利用半经验AM I量子化学方法计算了54个2-氨基-6-苯磺酰基苄腈及其类似物的物理化学、电子结构、指示变量等共28个参数,然后使用偏最小二乘,穷举回归和混沌遗传乘法训练的人工神经网络方法建立了这些参数和其抑制H IV-1逆转录酶活性之间的定量构效关系模型,为设计、合成更高生物活性的该类化合物提供了理论参考。     

Structural investigation of the HIV-1 envelope glycoprotein gp160 cleavage site

The selective proteolytic activation of the HIV-1 envelope glycoprotein gp160 by furin and other precursor convertases (PCs) occurs at the carboxyl side of the sequence Arg508-Glu-Lys-Arg511 (site 1), in spite of the presence of another consensus sequence: Lys500-Ala-Lys-Arg503 (site 2). We report on the solution structural analysis of a 19-residue synthetic peptide, p498, which spans the two gp160-processing sites 1 and 2, and is properly digested by furin at site 1. A molecular model is obtained for p498, by means of molecular dynamics simulations, from NMR data collected in trifluoroethanol/water. The peptide N-terminal side presents a 9-residue helical segment, enclosing the processing site 2; the C-terminal segment can be described as a loop exposing the processing site 1. A hypothesis for the docking of p498 onto the catalytic domain of human furin, modeled by homology and fitting previous site-directed mutagenesis studies, is also presented. p498 site 1 is shown to have easy access to the furin catalytic site, unlike the nonphysiological site 2. Finally, on the basis of available data, we suggest a possible structural motif required for the gp160-PCs recognition.     

Eigen value analysis of HIV-1 integrase inhibitors.

A three-dimensional quantitative structure activity relationship using the eigen value analysis (EVA) paradigm applied to 41 HIV-1 integrase inhibitors that inhibit integrase mediated cleavage (3'-processing step) and integration (3'-strand transfer step) in vitro was performed. The training set consisted of 35 molecules from five structurally diverse classes: salicylhydrazines, lichen acids, coumarins, quinones, and thiazolothiazepines. Models derived using semiempirical (MOPAC AM1 and PM3) calculated normal-mode frequencies were compared. The predictive ability of each resultant model was evaluated using a test set comprised of six molecules belonging to a different structural class: hydrazides. Models derived using AM1 method showed considerable internal as well as external predictivity (r(2)(cv) = 0.806, r(2)(pred) = 0.761 for 3'-processing and r(2)(cv) = 0.677, r(2)(pred) = 0.591 for 3'-strand transfer).     

Atypical protonation states in the active site of HIV-1 protease: a computational study

The HIV protease (HIVP) is a prominent example for successful structure-based drug design. Besides its pharmaceutical impact, it is a well-studied system for which, as experimentally evidenced, protonation changes in the active site occur upon ligand binding. Therefore, it serves as an ideal candidate for a case study using our newly developed partial charge model, which was optimized toward the application of Poisson-Boltzmann based pK(a) calculations. The charge model suggests reliably experimentally determined protonation states in the active site of HIVP. Furthermore, we perform pKa calculations for two HIVP complexes with novel types of inhibitors developed and synthesized in our group. For these complexes, no experimental knowledge about the protonation states is given. For one of the compounds, containing a central pyrrolidine ring, the calculations predict that both catalytic aspartates should be deprotonated upon ligand binding.     

Microcalorimetric studies of interactions between proteins and hydrophobic ligands in hydrophobic interaction chromatography: effects of ligand chain length, density and the amount of bound protein

Using isothermal titration calorimetry (ITC), this investigation directly measured the adsorption enthalpies of proteins on various hydrophobic adsorbents. Various amounts of butyl and octyl groups were attached onto CM-Sepharose to form C4 and C8, two types of hydrophobic adsorbents. The adsorption enthalpies of both trypsinogen and alpha-chymotrypsinogen A were measured at 4.0 M NaCl and pH 10.0, in which most ionic interaction was suppressed. The adsorption isotherms of both proteins on various adsorbents were also measured, thus allowing us to calculate the Gibbs free energy and entropy of adsorption. Experimental results indicated that the adsorption of both proteins on butyl-containing adsorbents was exothermic, while their adsorption on octyl ones was endothermic. In addition, binding of both proteins with the butyl ligand is basically an adsorption process, while binding with the octyl ligand is adsorption and partition processes. Moreover, on both butyl or octyl, the adsorption enthalpy became increasingly positive as the ligand density increased, while the adsorption entropy became more positive as the alkyl chain length or density of the adsorbent increased. In addition, ITC was used to measure protein-protein interaction. The adsorption enthalpy of both proteins increased as the amount of bound protein increased, and the enthalpy increase of trypsinogen appeared to be higher than that of alpha-chymotrypsinogen A. This observation implies that protein-protein repulsion was stronger among trypsinogen molecules in the experiments.     

Photoswitch inhibitors of alpha-chymotrypsin--increased substitution and peptidic character in peptidomimetic boronate esters

A series of peptidomimetic boronate esters containing the photoisomerisable azobenzene group has been synthesised and assayed against the serine protease alpha-chymotrypsin. Compounds with borophenylalanine inhibition groups were found to be inhibitors with micromolar activity, while those with aryl boronate inhibition groups were inactive. Selected compounds were isomerised by UV or visible light to obtain enriched states of the (Z) or (E) isomers, respectively, and assayed. A change in activity on photoisomerisation was observed, however some decomposition of the boronate group on irradiation was also observed, limiting reversibility.     

The dynamics of water at DNA interfaces: computational studies of Hoechst 33258 bound to DNA

Together, spectroscopy combined with computational studies that relate directly to the experimental measurements have the potential to provide unprecedented insight into the dynamics of important biological processes. Recent time-resolved fluorescence experiments have shown that the time scales for collective reorganization at the interface of proteins and DNA with water are more than an order of magnitude slower than in bulk aqueous solution. The molecular interpretation of this change in the collective response is somewhat controversial some attribute the slower reorganization to dramatically retarded water motion, while others describe rapid water dynamics combined with a slower biomolecular response. To connect directly to solvation dynamics experiments of the fluorescent probe Hoechst 33258 (H33258) bound to DNA, we have generated 770 ns of molecular dynamics (MD) simulations and calculated the equilibrium and nonequilibrium solvation response to excitation of the probe. The calculated time scales for the solvation response of H33258 free in solution (0.17 and 1.4 ps) and bound to DNA (1.5 and 20 ps) are highly consistent with experiment (0.2 and 1.2 ps, 1.4 and 19 ps, respectively). Decomposition of the calculated response revealed that water solvating the probe bound to DNA was still relatively mobile, only slowing by a factor of 2-3, while DNA motion was responsible for the long-time component (approximately 20 ps).