Synthesis and comparative molecular field analysis (CoMFA) of symmetric and nonsymmetric cyclic sulfamide HIV-1 protease inhibitors

Preventing lipophilic aggregation in cosolvent molecular dynamics simulations with hydrophobic probes using Plumed Automatic Restraining Tool (PART)

Cosolvent molecular dynamics (MD) simulations are molecular dynamics simulations used to identify preferable locations of small organic fragments on a protein target. Most cosolvent molecular dynamics workflows make use of only water-soluble fragments, as hydrophobic fragments would cause lipophilic aggregation. To date the two approaches that allow usage of hydrophobic cosolvent molecules are to use a low (0.2 M) concentration of hydrophobic probes, with the disadvantage of a lower sampling speed, or to use force field modifications, with the disadvantage of a difficult and inflexible setup procedure. Here we present a third alternative, that does not suffer from low sampling speed nor from cumbersome preparation procedures. We have built an easy-to-use open source command line tool PART (Plumed Automatic Restraining Tool) to generate a PLUMED file handling all intermolecular restraints to prevent lipophilic aggregation. We have compared restrained and unrestrained cosolvent MD simulations, showing that restraints are necessary to prevent lipophilic aggregation at hydrophobic probe concentrations of 0.5 M. Furthermore, we benchmarked PART generated restraints on a test set of four proteins (Factor-Xa, HIV protease, P38 MAP kinase and RNase A), showing that cosolvent MD with PART generated restraints qualitatively reproduces binding features of cocrystallised ligands.

Prediction of protein-ligand binding affinity via deep learning models

Accurately predicting the binding affinity between proteins and ligands is crucial in drug screening and optimization, but it is still a challenge in computer-aided drug design. The recent success of AlphaFold2 in predicting protein structures has brought new hope for deep learning (DL) models to accurately predict protein-ligand binding affinity. However, the current DL models still face limitations due to the low-quality database, inaccurate input representation and inappropriate model architecture. In this work, we review the computational methods, specifically DL-based models, used to predict protein-ligand binding affinity. We start with a brief introduction to protein-ligand binding affinity and the traditional computational methods used to calculate them. We then introduce the basic principles of DL models for predicting protein-ligand binding affinity. Next, we review the commonly used databases, input representations and DL models in this field. Finally, we discuss the potential challenges and future work in accurately predicting protein-ligand binding affinity via DL models.

Novel pseudonucleosides and sulfamoyl-oxazolidinone β-D-glucosamine derivative as anti-COVID-19: design, synthesis, and in silico study

New pseudonucleosides containing cyclic sulfamide moiety and sulfamoyl β-D-glucosamine derivative are described. These pseudonucleosides are synthesized in good yields starting from chlorosulfonyl isocyanate and β-D-glucosamine hydrochloride in five steps; (protection, acetylation, removal of the Boc group, sulfamoylation, and cyclization). Further, novel glycosylated sulfamoyloxazolidin-2-one is prepared in three steps; carbamoylation, sulfamoylation, and intramolecular cyclization. The structures of the synthesized compounds were confirmed by usual spectroscopic and spectrometric methods NMR, IR, MS, and EA. Interesting molecular docking of the prepared pseudonucleosides and (Beclabuvir, Remdesivir) drugs with SARS-CoV-2/Mpro (PDB:5R80) was conducted using the same parameters for a fair comparison. A low binding affinity of the synthesized compounds compared to the Beclabuvir and other analysis showed that pseudonucleosides have the ability to inhibit SARS-CoV-2. After the motivating results of molecular docking study, the complex between the SARS-CoV-2 Mpro and compound 7 was subjected to 100 ns molecular dynamics (MD) simulation using Desmond module of Schrodinger suite, during which the receptor-ligand complex showed substantial stability after 10 ns of MD simulation. Also, we studied the prediction of absorption, distribution, properties of metabolism, excretion, and toxicity (ADMET) of the synthesized compounds.Communicated by Ramaswamy H. Sarma.

Novel N-acylsulfonamides: Synthesis, in silico prediction, molecular docking dynamic simulation, antimicrobial and anti-inflammatory activities

Microbial resistance to drugs currently traded in the market is a serious problem in modern medicine. In this field of research, we synthesized a novel N-acylsulfonamides (NAS) derivatives starting from commercially available compounds; morpholine, isocyanate of chlorosulfonyl and alcohols. The in vitro antimicrobial potential of synthesized compounds was screened against 04 Gram-negative bacteria; Escherichia coli, Pseudomonas aeruginosa, Klebsiella pneumoniae, Acinetobacter baumannii, 02 Gram-positive bacteria: Streptococcus sp, Staphylococcus aureus and 07 yeasts and fungi: Candida albicans, Candida spp, Penicillum spp, Aspegillus sp, Aspergillus flavus, Fusarium sp, and Cladosporium spp. The results of inhibition growth were compared with standard antimicrobial drugs with the goal of exploring their potential antimicrobial activity. In addition, the anti-inflammatory activity of the synthesized compounds was determined in-vitro by protein denaturation method. The obtained bioactivity results were further validated by in silico DFT (Density Functional Theory), ADME (Absorption-Distribution-Métabolisation-Excrétion), molecular docking studies and molecular dynamics simulations.Communicated by Ramaswamy H. Sarma.

Investigation on the Anticancer Activity of Symmetric and Unsymmetric Cyclic Sulfamides

The sulfamide functional group has been extensively employed in organic synthesis to discover probes and drugs in various applications such as cancer, human immunodeficiency virus (HIV), virus, and diabetes. Herein, we describe the synthesis of 7-membered symmetric and unsymmetric sulfamide compounds and their biological evaluation through the National Cancer Institute (NCI) panel of 60 human tumor cell lines (NCI-60) and the mechanism of action study. The results of a study from the NCI-60 cell line exhibited that many synthesized cyclic sulfamide compounds inhibited breast cancer (MDA-MB-468). The mechanism of action study of a representative compound 18 showed the inhibition of proliferation and apoptosis in A549 lung cancer cells.

Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson–Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ_RMSE = 1.25 kcal mol⁻¹. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses.

Synthesis of N-Substituted Sulfamate Esters from Sulfamic Acid Salts by Activation with Triphenylphosphine Ditriflate

A general approach to access sulfamate esters through preparation of sulfamic acid salts, subsequent activation with triphenylphosphine ditriflate, and nucleophilic trapping is disclosed. The method proceeds in modest to excellent yields to incorporate nucleophiles derived from aliphatic alcohols and phenols. This approach can be employed to furnish differentially substituted sulfamides.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Pharmacophore modeling using site-identification by ligand competitive saturation (SILCS) with multiple probe molecules

Receptor-based pharmacophore modeling is an efficient computer-aided drug design technique that uses the structure of the target protein to identify novel leads. However, most methods consider protein flexibility and desolvation effects in a very approximate way, which may limit their use in practice. The Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling protocol (SILCS-Pharm) was introduced recently to address these issues, as SILCS naturally takes both protein flexibility and desolvation effects into account by using full molecular dynamics simulations to determine 3D maps of the functional group-affinity patterns on a target receptor. In the present work, the SILCS-Pharm protocol is extended to use a wider range of probe molecules including benzene, propane, methanol, formamide, acetaldehyde, methylammonium, acetate and water. This approach removes the previous ambiguity brought by using water as both the hydrogen-bond donor and acceptor probe molecule. The new SILCS-Pharm protocol is shown to yield improved screening results, as compared to the previous approach based on three target proteins. Further validation of the new protocol using five additional protein targets showed improved screening compared to those using common docking methods, further indicating improvements brought by the explicit inclusion of additional feature types associated with the wider collection of probe molecules in the SILCS simulations. The advantage of using complementary features and volume constraints, based on exclusion maps of the protein defined from the SILCS simulations, is presented. In addition, reranking using SILCS-based ligand grid free energies is shown to enhance the diversity of identified ligands for the majority of targets. These results suggest that the SILCS-Pharm protocol will be of utility in rational drug design.

Inclusion of multiple fragment types in the site identification by ligand competitive saturation (SILCS) approach

The site identification by ligand competitive saturation (SILCS) method identifies the location and approximate affinities of small molecular fragments on a target macromolecular surface by performing molecular dynamics (MD) simulations of the target in an aqueous solution of small molecules representative of different chemical functional groups. In this study, we introduce a set of small molecules to map potential interactions made by neutral hydrogen bond donors and acceptors and charged donor and acceptor fragments in addition to nonpolar fragments. The affinity pattern is obtained in the form of discretized probability or, equivalently, free energy maps, called FragMaps, which can be visualized with the target surface. We performed SILCS simulations for four proteins for which structural and thermodynamic data is available for multiple diverse ligands. Good overlap is shown between high affinity regions identified by the FragMaps and the crystallographic positions of ligand functional groups with similar chemical functionality, thus demonstrating the validity of the qualitative information obtained from the simulations. To test the ability of FragMaps in providing quantitative predictions, we calculate the previously introduced ligand grid free energy (LGFE) metric and observe its correspondence with experimentally measured binding affinity. LGFE is computed for different conformational ensembles and improvement in prediction is shown with increasing ligand conformational sampling. Ensemble generation includes a Monte Carlo sampling approach that uses the GFE FragMaps directly as the energy function. The results show that some but not all experimental trends are predicted and warrant improvements in the scoring methodology. In addition, the potential utility of atom-based free energy contributions to the LGFE scores and the use of multiple ligands in SILCS to identify displaceable water molecules during ligand design are discussed.

Microwave-assisted synthesis of small molecules targeting the infectious diseases tuberculosis, HIV/AIDS, malaria and hepatitis C

The unique properties of microwave in situ heating offer unparalleled opportunities for medicinal chemists to speed up lead optimisation processes in early drug discovery. The technology is ideal for small-scale discovery chemistry because it allows full reaction control, short reaction times, high safety and rapid feedback. To illustrate these advantages, we herein describe applications and approaches in the synthesis of small molecules to combat four of the most prevalent infectious diseases; tuberculosis, HIV/AIDS, malaria and hepatitis C, using dedicated microwave instrumentation.

Synthesis of new fused benzothiadiazepines and macrocyclic sulfamides starting from n,n-disubstituted sulfamides and n(boc)-sulfamides

Herein, we describe an efficient one-step synthesis of new fused benzothiadiazepine-1,1-dioxides and macrocyclic sulfamides. The synthesis of these compounds was achieved in moderate yields starting from previously described N,N′-disubstituted symmetric sulfamides and N-tert-butoxycarbonyl, N′-alkyl sulfamide. The chemical structures of all the new compounds reported in this work were confirmed by NMR, IR, and mass spectrometry. These compounds are beneficial building blocks that can be used in deriving new chemical entities that exert a wide spectrum of pharmacological activities.

Multiple receptor conformation docking and dock pose clustering as tool for CoMFA and CoMSIA analysis - a case study on HIV-1 protease inhibitors

Multiple receptors conformation docking (MRCD) and clustering of dock poses allows seamless incorporation of receptor binding conformation of the molecules on wide range of ligands with varied structural scaffold. The accuracy of the approach was tested on a set of 120 cyclic urea molecules having HIV-1 protease inhibitory activity using 12 high resolution X-ray crystal structures and one NMR resolved conformation of HIV-1 protease extracted from protein data bank. A cross validation was performed on 25 non-cyclic urea HIV-1 protease inhibitor having varied structures. The comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were generated using 60 molecules in the training set by applying leave one out cross validation method, r (loo) (2) values of 0.598 and 0.674 for CoMFA and CoMSIA respectively and non-cross validated regression coefficient r(2) values of 0.983 and 0.985 were obtained for CoMFA and CoMSIA respectively. The predictive ability of these models was determined using a test set of 60 cyclic urea molecules that gave predictive correlation (r (pred) (2) ) of 0.684 and 0.64 respectively for CoMFA and CoMSIA indicating good internal predictive ability. Based on this information 25 non-cyclic urea molecules were taken as a test set to check the external predictive ability of these models. This gave remarkable out come with r (pred) (2) of 0.61 and 0.53 for CoMFA and CoMSIA respectively. The results invariably show that this method is useful for performing 3D QSAR analysis on molecules having different structural motifs.

Density functional theory (DFT) calculations of conformational energies and interconversion pathways in 1,2,7-thiadiazepane

The molecular structure and conformational analysis of 1,2,7-thiadiazapane conformers were investigated by density functional theory (DFT) calculations at the B3LYP/cc-pVDZ level of theory. Four twist-chair (TC), six twist-boat (TB), two boat (B), two chair (C) and four twist (T) conformers were identified as minima and transition states for 1,2,7-thiadiazepane. The TC1 conformer is the most stable conformer and the twist-chair conformers are predicted to be lower in energy than their corresponding boat and chair conformations. DFT predicts a small barrier to pseudo-rotation and a remarkable activation barrier for the conformational interconversion of the twist-chair conformers to their corresponding boat conformers. The simplest conformational process and the one with the lowest barrier is the degenerate interconversion of the twist-chair 3 (TC3) conformation with itself via the CS symmetric chair (C2) transition state. The calculated strain energy barrier for this process is 2.41 kJ mol-1. The highest conformational interconversion barrier is between TC2 and twistboat 3 (TB3) forms, which was found to be 75.62 kJ mol-1.

Molecular docking and 3D-QSAR studies of HIV-1 protease inhibitors

HIV-1 protease is an obligatory enzyme in the replication process of the HIV virus. The abundance of structural information on HIV-1PR has made the enzyme an attractive target for computer-aided drug design strategies. The daunting ability of the virus to rapidly generate resistant mutants suggests that there is an ongoing need for new HIV-1PR inhibitors with better efficacy profiles and reduced toxicity. In the present investigation, molecular modeling studies were performed on a series of 54 cyclic urea analogs with symmetric P2/P2' substituents. The binding modes of these inhibitors were determined by docking. The docking results also provided a reliable conformational superimposition scheme for the 3D-QSAR studies. To gain insight into the steric, electrostatic, hydrophobic and hydrogen-bonding properties of these molecules and their influence on the inhibitory activity, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed. Two different alignment schemes viz. receptor-based and atom-fit alignment, were used in this study to build the QSAR models. The derived 3D-QSAR models were found to be robust with statistically significant r(2) and r(2)(pred) values and have led to the identification of regions important for steric, hydrophobic and electronic interactions. The predictive ability of the models was assessed on a set of molecules that were not included in the training set. Superimposition of the 3D-contour maps generated from these models onto the active site of enzyme provided additional insight into the structural requirements of these inhibitors. The CoMFA and CoMSIA models were used to design some new inhibitors with improved binding affinity. Pharmacokinetic and toxicity predictions were also carried out for these molecules to gauge their ADME and safety profile. The computational results may open up new avenues for synthesis of potent HIV-1 protease inhibitors.

Assessment of QM/MM scoring functions for molecular docking to HIV-1 protease

We explore the ability of four quantum mechanical (QM)/molecular mechanical (MM) models to accurately identify the native pose of six HIV-1 protease inhibitors and compare them with the AMBER force field and ChemScore and GoldScore scoring functions. Three QM/MM scoring functions treated the ligand at the HF/6-31G*, AM1d, and PM3 levels; the fourth QM/MM function modeled the ligand and active site at the PM3-D level. For the discrimination of native from non-native poses, solvent-corrected HF/6-31G*:AMBER and AMBER functions exhibited the best overall performance. While the electrostatic component of the MM and QM/MM functions appears important for discriminating the native pose of the ligand, the polarization contribution in the QM/MM functions was relatively insensitive to a ligand's binding mode and, for one ligand, actually hindered discrimination. The inclusion of a desolvation penalty, here using a generalized Born solvent model, improved discrimination for the MM and QM/MM methods. There appeared to be no advantage to binding mode prediction by incorporating active site polarization at the PM3-D level. Finally, we found that choice of the protonation state of the aspartyl dyad in the HIV-1 protease active site influenced the ability of scoring methods to determine the native binding pose.

Synthesis and evaluation of sulfamide-type indolizidines as glycosidase inhibitors

A practical synthesis of reducing sulfamide-derived iminosugar glycomimetics related to the indolizidine glycosidase inhibitor family is reported. The polyhydroxylated bicyclic system was built from readily accessible hexofuranose derivatives through a synthetic scheme that involves 5,6-cyclic sulfamides. Further intramolecular nucleophilic addition of the sulfamide nitrogen atom to the masked aldehyde group of the monosaccharide in the open chain form afforded the target sugar mimics. By starting from d-glucose and d-mannose precursors, 2-aza-3,3-dioxo-3-thiaindolizidine derivatives with hydroxylation profiles that matched those of (+)-castanospermine and 6-epi-(+)-castanospermine were obtained. In vitro screening against a panel of glycosidases evidenced a high selectivity towards alpha-mannosidase.

Exploring structural variability in X-ray crystallographic models using protein local optimization by torsion-angle sampling

Modeling structural variability is critical for understanding protein function and for modeling reliable targets for in silico docking experiments. Because of the time-intensive nature of manual X-ray crystallographic refinement, automated refinement methods that thoroughly explore conformational space are essential for the systematic construction of structurally variable models. Using five proteins spanning resolutions of 1.0-2.8 A, it is demonstrated how torsion-angle sampling of backbone and side-chain libraries with filtering against both the chemical energy, using a modern effective potential, and the electron density, coupled with minimization of a reciprocal-space X-ray target function, can generate multiple structurally variable models which fit the X-ray data well. Torsion-angle sampling as implemented in the Protein Local Optimization Program (PLOP) has been used in this work. Models with the lowest R(free) values are obtained when electrostatic and implicit solvation terms are included in the effective potential. HIV-1 protease, calmodulin and SUMO-conjugating enzyme illustrate how variability in the ensemble of structures captures structural variability that is observed across multiple crystal structures and is linked to functional flexibility at hinge regions and binding interfaces. An ensemble-refinement procedure is proposed to differentiate between variability that is a consequence of physical conformational heterogeneity and that which reflects uncertainty in the atomic coordinates.

Complexity in modeling and understanding protonation states: computational titration of HIV-1-protease-inhibitor complexes

The computational-titration (CT) algorithm based on the 'natural' Hydropathic INTeractions (HINT) force field is described. The HINT software model is an empirical, non-Newtonian force field derived from experimentally measured partition coefficients for solvent transfer between octanol and H(2)O (log P(o/w)). The CT algorithm allows the identification, modeling, and optimization of multiple protonation states of residues and ligand functional groups at the protein-ligand active site. The importance of taking into account pH and ionization states of residues, which strongly affect the process of ligand binding, for correctly predicting binding free energies is discussed. The application of the CT protocol to a set of six cyclic inhibitors in their complexes with HIV-1 protease is presented, and the advance of HINT as a virtual-screening tool is outlined.

Bromoallenes as Allyl Dication Equivalents in the Presence or Absence of Palladium(0): Direct Construction of Bicyclic Sulfamides Containing Five- to Eight-membered Rings by Tandem Cyclization of Bromoallenes

A highly regioselective synthesis of bicyclic sulfamides is described. Based on our recent discovery that bromoallenes can act as allyl dication equivalents in the presence of a palladium catalyst and alcohol, we investigated tandem cyclization of bromoallenes bearing a sulfamide group. It is found that some bromoallenes act as allyl dication equivalents even in the absence of a palladium(0) catalyst to afford cyclosulfamides containing five- or six-membered rings. While the palladium-free cyclization is dependent on the substrate structure affording the bicyclic sulfamides through the first cyclization onto the proximal or central carbon atom of the bromoallenes, the palladium-catalyzed reaction strongly promotes the first cyclization onto the central allenic carbon atom to afford bicyclic sulfamides containing a seven- or eight-membered ring. Formation of two types of bicyclic sulfamides from single bromoallenes by simply changing the reaction conditions is also described.

Unexpected novel binding mode of pyrrolidine-based aspartyl protease inhibitors: design, synthesis and crystal structure in complex with HIV protease

At present nine FDA-approved HIV protease inhibitors have been launched to market, however rapid drug resistance arising under antiviral therapy calls upon novel concepts. Possible strategies are the development of ligands with less peptide-like character or the stabilization of a new and unexpected binding-competent conformation of the protein through a novel ligand-binding mode. Our rational design of pyrrolidinedimethylene diamines was inspired by the idea to incorporate key structural elements from classical peptidomimetics with a non-peptidic heterocyclic core comprising an endocyclic amino function to address the catalytic aspartic acid side chains of Asp 25 and 25'. The basic scaffolds were decorated by side chains already optimized for the recognition pockets of HIV protease or cathepsin D. A multistep synthesis has been established to produce the central heterocycle and to give flexible access to side chain decorations. Depending on the substitution pattern of the pyrrolidine moiety, single-digit micromolar inhibition of HIV-1 protease and cathepsin D has been achieved. Successful design is suggested in agreement with our modelling concepts. The subsequently determined crystal structure with HIV protease shows that the pyrrolidine moiety binds as expected to the pivotal position between both aspartic acid side chains. However, even though the inhibitors have been equipped symmetrically by polar acceptor groups to address the flap water molecule, it is repelled from the complex, and only one direct hydrogen bond is formed to the flap. A strong distortion of the flap region is detected, leading to a novel hydrogen bond which cross-links the flap loops. Furthermore, the inhibitor addresses only three of the four available recognition pockets. It achieves only an incomplete desolvation compared with the similarly decorated amprenavir. Taking these considerations into account it is surprising that the produced pyrrolidine derivatives achieve micromolar inhibition and it suggests extraordinary potency of the new compound class. Most likely, the protonated pyrrolidine moiety experiences strong enthalpic interactions with the enzyme through the formation of two salt bridges to the aspartic acid side chains. This might provide challenging opportunities to combat resistance of the rapidly mutating virus.

New uses for the Burgess reagent in chemical synthesis: methods for the facile and stereoselective formation of sulfamidates, glycosylamines, and sulfamides

Although the Burgess reagent (methoxycarbonylsulfamoyltriethylammonium hydroxide, inner salt) has found significant use in chemical synthesis as a dehydrating agent, almost no work has been directed towards its potential in other synthetic applications. As this article will detail, we have found that the Burgess reagent is remarkably effective at accomplishing a number of non-dehydrative synthetic tasks when applied to appropriate substrates, such as the formation of sulfamidates from 1,2-diols or epoxyalcohols, alpha- and beta-glycosylamines from carbohydrates, and cyclic sulfamides from 1,2-aminoalcohols. Beyond delineating the power of these new reaction manifolds, we also describe the construction of a group of alternative Burgess-type reagents that extends the scope of these new reactions even further.

Electrophysiological Characterization and Modeling of the Structure Activity Relationship of the Human Concentrative Nucleoside Transporter 3 (hCNT3)

We characterized the electrophysiology, kinetics, and quantitative structure-activity relationship (QSAR) of the human concentrative nucleoside transporter 3 (hCNT3) expressed in Xenopus laevis oocytes by measuring substrate-induced inward currents using a two-microelectrode voltage-clamp system. At membrane potentials between -30 and -150 mV, sodium activation of gemcitabine transport was sigmoidal, with a K0.5 of 8.5+/-0.3 mM for Na+ and a Hill coefficient of 2.2+/-0.25 independent of membrane potential. We measured the Imax and K0.5 for substrate at -50 mV for the nucleoside analog drugs gemcitabine (638+/-58 nA, 59.7+/-17.5 microM), ribavirin (546+/-37 nA, 61.0+/-13.2 microM), AZT (420+/-4 nA, 310+/-9 microM), and 3-deazauridine (506+/-30 nA, 50.8+/-9.90 microM). K0.5 and Imax for substrate were dependent on membrane potential (both increasing as the membrane became more hyperpolarized) for all four drugs. hCNT3 also exhibited pre-steady-state currents. The quantitative structure-activity relationship (QSAR) was examined using comparative molecular field analysis and comparative molecular similarity indices analysis of the inward currents induced by 27 nucleoside analogs with substitutions at both the ribose and the nucleobase. Two statistically significant QSAR models identified electrostatic interaction as the major force in hCNT3 transport and attributed a critical role to the 3'-hydroxyl position of hCNT3 substrates. Steric hindrance at the 3-position and positive charge at the 5-position of the pyrimidine ring were favorable for transport. Two hCNT3 pharmacophore models revealed the minimal features required for hCNT3 transport as two hydrogen bond acceptors at 3'-OH and 5'-O and the hydrophobic center occupied by the base ring.

Evaluation of a neural networks QSAR method based on ligand representation using substituent descriptors. Application to HIV-1 protease inhibitors

We present here a neural networks method designed to predict biological activity based on a local representation of the ligand. The compounds of the series are represented by a vector mapping for each of four substituent properties: volume, log P, dipole moment and a simple 'steric' parameter relating to its shape. This ligand representation was tested using neural networks on a set of 42 cyclic-urea derivatives, inhibiting HIV-1 protease. The leave-one-out cross-validation using all descriptors in the input gave a correlation factor between prediction and experiment of 0.76 for the overall set and 0.88 when three outliers were left out. To rank the significance of the four descriptors, we further tested all combinations of two and three parameters for each substituent, using two disjunctive testing sets of five inhibitors. In these sets, vectors with extreme descriptor values were used either in the training or the testing set (sets A and B, respectively). The method is a very good interpolator (set A, 95+/-2% accuracy) but a less effective extrapolator (set B, 85+/-2% accuracy). Generally, the combinations including the 'steric' parameter predict better than average, while those containing the volume are less effective. The best prediction, 98.8+/-1.2%, was obtained when log P, the dipole and the steric parameter were used on set A. At the opposite end, the lowest ranked descriptor set was obtained when replacing log P with the volume, giving 92.3+/-6.7% accuracy over the set A.

N,N'-SUBStituted 1,2,5 thiadiazolidine 1,1-dioxides: synthesis, selected chemical and spectral proprieties and antimicrobial evaluation

The sulfamide functional group is increasingly relevant in both medicinal and bioorganic chemistry. We report here practical access to a series of N2,N5-substituted five-membered cyclosulfamides. The five-membered heterocyclic motif was prepared starting from proteogenic amino acids and chlorosulfonyl isocyanate via the Mitsunobu reaction. Selected chemical and spectral proprieties and the antimicrobial evaluation of these compounds are detailed.

Computational approaches to modeling drug transporters

Computational modeling has advanced our understanding of drug absorption, tissue distribution, excretion and toxicity profiles by providing both direct and indirect knowledge of drug-transporter interactions that would otherwise be unavailable using experimental methods. Currently, two complementary approaches are available in modeling transporters: substrate-based and transporter-based methods. The transporter-based approach directly predicts the transporter's three-dimensional structure to assist in understanding the drug transport process, whereas substrate-based models infer such information by studying a group of substrates or inhibitors with measured activities. In this review, the available strategies in both transporter-based and substrate-based approaches are explained and illustrated with applications and case studies. With increasing computational power and continuously improving modeling algorithms, computational techniques can assist in further understanding transporter-substrate interactions as well as, the optimization of transporter-directed drug design.

Alternative method of Boc-removal from sulfamide using silica-phenyl sulfonic acid in conjunction with microwave heating

A general method was used to prepare an array of unsymmetric sulfamides. This was accomplished by the stepwise addition of CSI to tert-butanol followed by the addition of amines. To increase diversity, nitrogen group of Boc-sulfamides was alkylated with alcohols using Mitsunobu reaction and Boc-group was removed using Si-TsOH. Microwave heating was used in all the steps. The final sulfamides were released from Si-TsOH using NH3 in MeOH.

Cyclic sulfamide HIV-1 protease inhibitors, with sidechains spanning from P2/P2' to P1/P1'

Previous studies of HIV protease inhibitors have shown that it is possible to elongate the P1/P1' sidechains to reach the S3/S3' binding sites. By analogy, we expected that it would be possible to design inhibitors reaching between the S1/S1' and S2/S2' binding sites. Molecular modeling suggested that this could be achieved with appropriate ortho-substitution of the P2/P2' benzyl groups in our cyclic sulfamide inhibitors. Four different spacer groups were investigated. The compounds were smoothly prepared from tartaric acid in five steps and exhibit low to moderate activity, the most potent inhibitor possessing a Ki value of 0.53 microM.