Comparing ligand interactions with multiple receptors via serial docking

Standard uses of ligand-receptor docking typically focus on the association of candidate ligands with a single targeted receptor, but actual applications increasingly require comparisons across multiple receptors. This study demonstrates that comparative docking to multiple receptors can help to select homology models for virtual compound screening and to discover ligands that bind to one set of receptors but not to another, potentially similar, set. A serial docking algorithm is furthermore described that reduces the computational costs of such calculations by testing compounds against a series of receptor structures and discarding a compound as soon as it fails to satisfy specified bind/no bind criteria for each receptor. The algorithm also realizes substantial efficiencies by taking advantage of the fact that a ligand typically binds in similar conformations to similar receptors. Thus, once detailed docking has been used to fit a ligand into the first of a series of similar receptors, much less extensive calculations can be used for the remaining structures.     

Enhanced docking with the mining minima optimizer: acceleration and side-chain flexibility

The ligand-protein docking algorithm based on the Mining Minima method has been substantially enhanced. First, the basic algorithm is accelerated by: (1) adaptively determining the extent of each energy well to help avoid previously discovered energy minima; (2) biasing the search away from ligand positions at the surface of the receptor to prevent the ligand from staying at the surface when large sampling regions are used; (3) quickly testing multiple different ligand positions and orientations for each ligand conformation; and (4) tuning the source code to increase computational efficiency. These changes markedly shorten the time needed to discover an accurate result, especially when large sampling regions are used. The algorithm now also allows user-selected receptor sidechains to be treated as mobile during the docking procedure. The energies associated with the mobile side chains are computed as if they belonged to the ligand, except that atoms at the boundary between side chains and the rigid backbone are treated specially. This new capability is tested for several well-known ligand/protein systems, and preliminary application to an enzyme whose substrate is unknown--the recently solved hypothetical protein YecO (HI0319) from Haemophilus influenzae--indicates that side-chains relaxations allow candidate substrates of various sizes to be accommodated.     

A theoretical investigation of adsorbate-induced surface relaxation effects using cluster models: Al on Si(111)

Differently sized cluster calculations are used to investigate theoretically the preferred adsorption site for an Al adatom on the Si(111) surface. By performing partial geometry optimizations at the Hartree-Fock level on AlSi_n subclusters around the site of interest we find significant Al adatom-induced surface relaxation effects distorting the Si atoms from their bulk lattice positions. The largest relaxation effects take place at the T₄ site resulting in Al adsorption at the T₄ site to be 5 kcal/mol more stable than at the H₃ site and considerably more stable than adsorption at the T₁ site. However, we only have confidence in this result after performing for the T₄ site a partial geometry optimization on the AlSi₅ subcluster in the AlSi₂₆H₂₄ cluster and by including appropriate correlation corrections.     

HIV-1 protease inhibitors from inverse design in the substrate envelope exhibit subnanomolar binding to drug-resistant variants

The acquisition of drug-resistant mutations by infectious pathogens remains a pressing health concern, and the development of strategies to combat this threat is a priority. Here we have applied a general strategy, inverse design using the substrate envelope, to develop inhibitors of HIV-1 protease. Structure-based computation was used to design inhibitors predicted to stay within a consensus substrate volume in the binding site. Two rounds of design, synthesis, experimental testing, and structural analysis were carried out, resulting in a total of 51 compounds. Improvements in design methodology led to a roughly 1000-fold affinity enhancement to a wild-type protease for the best binders, from a Ki of 30-50 nM in round one to below 100 pM in round two. Crystal structures of a subset of complexes revealed a binding mode similar to each design that respected the substrate envelope in nearly all cases. All four best binders from round one exhibited broad specificity against a clinically relevant panel of drug-resistant HIV-1 protease variants, losing no more than 6-13-fold affinity relative to wild type. Testing a subset of second-round compounds against the panel of resistant variants revealed three classes of inhibitors: robust binders (maximum affinity loss of 14-16-fold), moderate binders (35-80-fold), and susceptible binders (greater than 100-fold). Although for especially high-affinity inhibitors additional factors may also be important, overall, these results suggest that designing inhibitors using the substrate envelope may be a useful strategy in the development of therapeutics with low susceptibility to resistance.     

Screening drug-like compounds by docking to homology models: a systematic study

In the absence of an experimentally solved structure, a homology model of a protein target can be used instead for virtual screening of drug candidates by docking and scoring. This approach poses a number of questions regarding the choice of the template to use in constructing the model, the accuracy of the screening results, and the importance of allowing for protein flexibility. The present study addresses such questions with compound screening calculations for multiple homology models of five drug targets. A central result is that docking to homology models frequently yields enrichments of known ligands as good as that obtained by docking to a crystal structure of the actual target protein. Interestingly, however, standard measures of the similarity of the template used to build the homology model to the targeted protein show little correlation with the effectiveness of the screening calculations, and docking to the template itself often is as successful as docking to the corresponding homology model. Treating key side chains as mobile produces a modest improvement in the results. The reasons for these sometimes unexpected results, and their implications for future methodologic development, are discussed.     

Study on the cyclization of 6-arylethynylpyrimidine-5-carbaldehydes with tert-butylamine: microwave versus thermal preparation of pyrido[4,3-d]pyrimidines

Thermal and microwave initiated cyclization of 2,4-disubstituted 6-arylethynylpyrimidine-5-carbaldehydes with tert-butylamine has been studied. A novel high-yielding preparation of 2,4-disubstituted 7-arylpyrido[4,3-d]pyrimidines has been developed. The intermediate compounds were isolated and possible mechanism of the reactions is discussed.