Homology modeling by the ICM method

Five models have been built by the ICM method for the Comparative Modeling section of the Meeting on the Critical Assessment of Techniques for Protein Structure Prediction. The targets have homologous proteins with known three-dimensional structure with sequence identity ranging from 25 to 77%. After alignment of the target sequence with the related three-dimensional structure, the modeling procedure consists of two subproblems: side-chain prediction and loop prediction. The ICM method approaches these problems with the following steps: (1) a starting model is created based on the homologous structure with the conserved portion fixed and the noncon-served portion having standard covalent geometry and free torsion angles; (2) the Biased Probability Monte Carlo (BPMC) procedure is applied to search the subspaces of either all the nonconservative side-chain torsion angles or torsion angles in a loop backbone and surrounding side chains. A special algorithm was designed to generate low-energy loop deformations. The BPMC procedure globally optimizes the energy function consisting of ECEPP/3 and solvation energy terms. Comparison of the predictions with the NMR or crystallographic solutions reveals a high proportion of correctly predicted side chains. The loops were not correctly predicted because imprinted distortions of the backbone increased the energy of the near-native conformation and thus made the solution unrecognizable. Interestingly, the energy terms were found to be reliable and the sampling of conformational space sufficient. The implications of this finding for the strategies of future comparative modeling are discussed. © 1995 Wiley-Liss, Inc.     

Structural analysis of human and macaque mAbs 2909 and 2.5B: implications for the configuration of the quaternary neutralizing epitope of HIV-1 gp120

The quaternary neutralizing epitope (QNE) of HIV-1 gp120 is preferentially expressed on the trimeric envelope spikes of intact HIV virions, and QNE-specific monoclonal antibodies (mAbs) potently neutralize HIV-1. Here, we present the crystal structures of the Fabs of human mAb 2909 and macaque mAb 2.5B. Both mAbs have long beta hairpin CDR H3 regions >20 ? in length that are each situated at the center of their respective antigen-binding sites. Computational analysis showed that the paratopes include the whole CDR H3, while additional CDR residues form shallow binding pockets. Structural modeling suggests a way to understand the configuration of QNEs and the antigen-antibody interaction for QNE mAbs. Our data will be useful in designing immunogens that may elicit potent neutralizing QNE Abs.     

Human anti-V3 HIV-1 monoclonal antibodies encoded by the VH5-51/VL lambda genes define a conserved antigenic structure

Preferential usage of immunoglobulin (Ig) genes that encode antibodies (Abs) against various pathogens is rarely observed and the nature of their dominance is unclear in the context of stochastic recombination of Ig genes. The hypothesis that restricted usage of Ig genes predetermines the antibody specificity was tested in this study of 18 human anti-V3 monoclonal Abs (mAbs) generated from unrelated individuals infected with various subtypes of HIV-1, all of which preferentially used pairing of the VH5-51 and VL lambda genes. Crystallographic analysis of five VH5-51/VL lambda-encoded Fabs complexed with various V3 peptides revealed a common three dimensional (3D) shape of the antigen-binding sites primarily determined by the four complementarity determining regions (CDR) for the heavy (H) and light (L) chains: specifically, the H1, H2, L1 and L2 domains. The CDR H3 domain did not contribute to the shape of the binding pocket, as it had different lengths, sequences and conformations for each mAb. The same shape of the binding site was further confirmed by the identical backbone conformation exhibited by V3 peptides in complex with Fabs which fully adapted to the binding pocket and the same key contact residues, mainly germline-encoded in the heavy and light chains of five Fabs. Finally, the VH5-51 anti-V3 mAbs recognized an epitope with an identical 3D structure which is mimicked by a single mimotope recognized by the majority of VH5-51-derived mAbs but not by other V3 mAbs. These data suggest that the identification of preferentially used Ig genes by neutralizing mAbs may define conserved epitopes in the diverse virus envelopes. This will be useful information for designing vaccine immunogen inducing cross-neutralizing Abs.     

Towards a species-selective acetylcholinesterase inhibitor to control the mosquito vector of malaria, Anopheles gambiae

Anopheles gambiae is the major mosquito vector of malaria in sub-Saharan Africa. At present, insecticide-treated nets (ITNs) impregnated with pyrethroid insecticides are widely used in malaria-endemic regions to reduce infection; however the emergence of pyrethroid-resistant mosquitoes has significantly reduced the effectiveness of the pyrethroid ITNs. An acetylcholinesterase (AChE) inhibitor that is potent for An. gambiae but weakly potent for the human enzyme could potentially be safely deployed on a new class of ITNs. In this paper we provide a preliminary pharmacological characterization of An. gambiae AChE, discuss structural features of An. gambiae and human AChE that could lead to selective inhibition, and describe compounds with 130-fold selectivity for inhibition of An. gambiae AChE relative to human AChE.     

Flexible ligand docking to multiple receptor conformations: a practical alternative

State of the art docking algorithms predict an incorrect binding pose for about 50-70% of all ligands when only a single fixed receptor conformation is considered. In many more cases, lack of receptor flexibility results in meaningless ligand binding scores, even when the correct pose is obtained. Incorporating conformational rearrangements of the receptor binding pocket into predictions of both ligand binding pose and binding score is crucial for improving structure-based drug design and virtual ligand screening methodologies. However, direct modeling of protein binding site flexibility remains challenging because of the large conformational space that must be sampled, and difficulties remain in constructing a suitably accurate energy function. Here we show that using multiple fixed receptor conformations, either experimentally determined by crystallography or NMR, or computationally generated, is a practical shortcut that may improve docking calculations. In several cases, such an approach has led to experimentally validated predictions.     

High-throughput docking for lead generation

Recent improvements in flexible docking technology may lead to a bigger role for computational methods in lead discovery. Although fast and accurate computational prediction of binding affinities for an arbitrary molecule is still beyond the limits of current methods, the docking and screening procedures can select small sets of likely lead candidates from large libraries of either commercially or synthetically available compounds.     

Estimated Secondary Structure Propensities within V1/V2 Region of HIV gp120 Are an Important Global Antibody Neutralization Sensitivity Determinant

Background

Neutralization sensitivity of HIV-1 virus to antibodies and anti-sera varies greatly between the isolates. Significant role of V1/V2 domain as a global neutralization sensitivity regulator has been suggested. Recent X-ray structures revealed presence of well-defined tertiary structure within this domain but also demonstrated partial disorder and conformational heterogeneity.

Methods

Correlations of neutralization sensitivity with the conformational propensities for beta-strand and alpha-helix formation over the entire folded V1/V2 domain as well as within sliding 5-residue window were investigated. Analysis was based on a set of neutralization data for 106 HIV isolates for which consistent neutralization sensitivity measurements against multiple pools of human immune sera have been previously reported.

Results

Significant correlation between beta-sheet formation propensity of the folded segments of V1/V2 domain and neutralization sensitivity was observed. Strongest correlation peaks localized to the beta-strands B and C. Correlation persisted when subsets of HIV isolates belonging to clades B, C and circulating recombinant form BC where analyzed individually or in combinations.

Conclusions

Observed correlations suggest that stability of the beta-sheet structure and/or degree of structural disorder in the V1/V2 domain is an important determinant of the global neutralization sensitivity of HIV-1 virus. While specific mechanism is to yet to be investigated, plausible hypothesis is that less ordered V1/V2s may have stronger masking effect on various neutralizing epitopes, perhaps effectively occupying larger volume and thereby occluding antibody access.     

Re-engineering aryl methylcarbamates to confer high selectivity for inhibition of Anopheles gambiae versus human acetylcholinesterase

To identify potential human-safe insecticides against the malaria mosquito we undertook an investigation of the structure-activity relationship of aryl methylcarbamates inhibitors of acetylcholinesterase (AChE). Compounds bearing a β-branched 2-alkoxy or 2-thioalkyl group were found to possess good selectivity for inhibition of Anopheles gambiae AChE over human AChE; up to 530-fold selectivity was achieved with carbamate 11d. A 3D QSAR model is presented that is reasonably consistent with log inhibition selectivity of 34 carbamates. Toxicity of these compounds to live Anopheles gambiae was demonstrated using both tarsal contact (filter paper) and topical application protocols.     

ICM-DISCO docking by global energy optimization with fully flexible side-chains

The ICM-DISCO (Docking and Interface Side-Chain Optimization) protein-protein-docking method is a direct stochastic global energy optimization from multiple starting positions of the ligand. The first step is performed by docking of a rigid all-atom ligand molecule to a set of soft receptor potentials precalculated on a 0.5 A grid from realistic solvent-corrected force-field energies. This step finds the correct solution as the lowest energy conformation in almost 100% of the cases in which interfaces do not change on binding. The second step is needed to deal with the induced changes and includes the global optimization of the interface side-chains of up to 400 best solutions. The CAPRI predictions were performed fully automatically with this method. Available experimental information was included as a filtering step to favor expected docking surfaces. In three of the seven proposed targets, the ICM-DISCO method found a good solution (>50% of correct contacts) within the five submitted models. The procedure is global and fully automated. We demonstrate that the algorithm handles the induced changes of surface side-chains but is less successful if the backbone undergoes large-scale rearrangements.