Biomedical subjects

Ricardo L Mancera

Publications and source records attributed to Ricardo L Mancera.

13 recordsLinked to original sources

PubMed2006

Novel structural features of CDK inhibition revealed by an ab initio computational method combined with dynamic simulations.

The rational development of specific inhibitors for the approximately 500 protein kinases encoded in the human genome is impeded by a poor understanding of the structural basis for the activity and selectivity of small molecules that compete for ATP binding. Combining classical dynamic simulations with a novel ab initio computational approach linear-scalable to molecular interactions involving thousands of atoms, we have investigated the binding of five distinct inhibitors to the cyclin-dependent kinase CDK2. We report here that polarization and dynamic hydrogen bonding effects, so far undetected by crystallography, affect both their activity and selectivity. The effects arise from the specific solvation patterns of water molecules in the ATP binding pocket or the intermittent formation of hydrogen bonds during the dynamics of CDK/inhibitor interactions and explain the unexpectedly high potency of certain inhibitors such as 3-(3H-imidazol-4-ylmethylene)-5-methoxy-1,3-dihydro-indol-2-one (SU9516). The Lys89 residue in the ATP-binding pocket of CDK2 is observed to form temporary hydrogen bonds with the three most potent inhibitors. This residue is replaced in CDK4 by Thr89, whose shorter side-chain cannot form similar bonds, explaining the relative selectivity of the inhibitors for CDK2. Our results provide a generally applicable computational method for the analysis of biomolecular structures and reveal hitherto unrecognized features of the interaction between protein kinases and their inhibitors.

Adenosine Triphosphate↗

PubMed2006

Comparative analysis of the surface interaction properties of the binding sites of CDK2, CDK4, and ERK2.

Recently developed hydrogen-bonding and hydrophobic analysis algorithms were used to investigate the interaction properties of the ATP binding sites of CDK2, CDK4, and ERK2. We were able to prioritise those hydrogen-bonding groups that are observed to bind the native ATP ligand, as well as to identify other important groups found to bind inhibitors of these enzymes. However, as the hydrogen-bonding groups in the ATP binding sites of these enzymes are fairly well-conserved, we have confirmed that inhibitor selectivity may be predominantly due to differences in either the hydrophobic or steric properties of their binding sites. In particular, the hydrophobic properties of regions outside the specificity surface were observed to provide a rationale for the differences in specificity between various inhibitors to these enzymes. Our method was thus able to identify variations in hydrophobicity. The greater hydrophobicity of certain regions of CDK4 over analogous regions in CDK2 was detectable; likewise, it was possible to distinguish variations in hydrophobicity for regions of CDK2 against those in ERK2, despite the fact that these regions are largely composed of similar residue types.

Adenosine Triphosphate↗

PubMed2005

The effect of a tightly bound water molecule on scaffold diversity in the computer-aided de novo ligand design of CDK2 inhibitors.

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily. [Figure: see text].

Binding Sites↗

PubMed2005

A new method for estimating the importance of hydrophobic groups in the binding site of a protein.

Interactions between the hydrophobic regions of a binding site and those of a complementary ligand are often observed to provide the driving force for binding. We present a new method for the analysis of hydrophobic regions in the binding site of a protein that considers not only atom type but also the nonadditive effects arising from the shape and extent of a nonpolar region. The method has been parametrized using a purpose-built genetic algorithm to optimize its ability to identify those regions that are more likely to form a strong interaction with a nonpolar ligand group. We demonstrate the ability of this method to account for changes in the shape and extent of the exposed nonpolar surface, using both artificial and protein examples. The method is also able to rationalize differences in binding affinity for ligand-protein complexes with largely hydrophobic binding sites.

Algorithms↗

PubMed2004

Assessment of multiple binding modes in ligand-protein docking.

Computational ligand-protein docking is routinely used for binding mode prediction. We have quantified the effect of considering multiple docking solutions on the success rate of obtaining the crystallographic binding mode. By selection of a small set of representatives, the experimentally observed binding mode can be predicted with a higher probability after a ligand-protein docking simulation. The proportion of correctly predicted complexes improved from 69% to 87% when five distinct binding modes were considered.

Crystallography, X-Ray↗

PubMed2004

Ligand-protein docking using a quantum stochastic tunneling optimization method.

A novel hybrid optimization method called quantum stochastic tunneling has been recently introduced. Here, we report its implementation within a new docking program called EasyDock and a validation with the CCDC/Astex data set of ligand-protein complexes using the PLP score to represent the ligand-protein potential energy surface and ScreenScore to score the ligand-protein binding energies. When taking the top energy-ranked ligand binding mode pose, we were able to predict the correct crystallographic ligand binding mode in up to 75% of the cases. By using this novel optimization method run times for typical docking simulations are significantly shortened.

Algorithms↗

PubMed2004

Comparative estimation of vibrational entropy changes in proteins through normal modes analysis.

We compare the vibrational entropy changes of proteins calculated using a full and a number of approximate normal modes analysis methods. The vibrational entropy differences for three conformational changes and three protein binding interactions were computed. In general, the approximate methods yield good estimates of the vibrational entropy change in a fraction of the time required by full normal modes analysis. The absolute entropies are either overestimated or greatly underestimated, but the difference is sufficiently accurate for some methods. This indicates that some of the approximate methods can give reasonable estimates of the associated vibrational entropy changes, making them suitable for inclusion in free energy calculations.

Animals↗

PubMed2004

Molecular modelling prediction of ligand binding site flexibility.

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.

Algorithms↗

PubMed2003

WaterScore: a novel method for distinguishing between bound and displaceable water molecules in the crystal structure of the binding site of protein-ligand complexes.

We have performed a multivariate logistic regression analysis to establish a statistical correlation between the structural properties of water molecules in the binding site of a free protein crystal structure, with the probability of observing the water molecules in the same location in the crystal structure of the ligand-complexed form. The temperature B-factor, the solvent-contact surface area, the total hydrogen bond energy and the number of protein-water contacts were found to discriminate between bound and displaceable water molecules in the best regression functions obtained. These functions may be used to identify those bound water molecules that should be included in structure-based drug design and ligand docking algorithms. FIGURE The binding site ( thin sticks) of penicillopepsin (3app) with its crystallographically determined water molecules ( spheres) and superimposed ligand (in thick sticks, from complexed structure 1ppk). Water molecules sterically displaced by the ligand upon complexation are shown in cyan. Bound water molecules are shown in blue. Displaced water molecules are shown in yellow. Water molecules removed from the analysis due to a lack of hydrogen bonds to the protein are shown in white. WaterScore correctly predicted waters in blue as Probability=1 to remain bound and waters in yellow as Probability<1x10(-20) to remain bound.

Algorithms↗

PubMed2003

A new method for estimating the importance of hydrogen-bonding groups in the binding site of a protein.

We introduce a new method to estimate the importance of hydrogen-bonding sitepoints in the binding site of a protein as part of a structure-based design strategy. Our method identifies hydrogen-bonding sitepoints within a binding pocket and ranks them according to both the accessibility of their hydrogen bonding regions to incoming ligands and their hydrogen-bonding strength. The combination of these components produces a prioritised list of sitepoints that are more likely to be involved in hydrogen bonding with an incoming ligand. A dataset of known protein-ligand interactions was used to compare the prioritisation of sitepoints identified by our method with those observed to be engaged in hydrogen bonding in their crystal structures. Our method was able to remove those sitepoints unable to bind the ligand due to a low accessibility or an unfavourable orientation and to award significantly higher hydrogen-bonding ranking values to those sitepoints observed to form hydrogen bonds. Our method can thus be used to identify hydrogen-bonding sitepoints that should be targeted preferentially in a drug design strategy.

Binding Sites↗

PubMed2002

De novo ligand design with explicit water molecules: an application to bacterial neuraminidase.

Most computer-aided drug design methods ignore the presence of crystallographically-determined water molecules in the binding site of a target protein. In this paper, our de novo ligand design methods are applied to the X-ray crystal structure of bacterial neuraminidase in the presence of some selected water molecules. We have found that, for this particular protein, the complete removal of all bound water molecules leads to difficulties in generating any potential ligands if the unsatisfied hydrogen-bonding sitepoints left by removing these water molecules are to be satisfied by a ligand. As more of the crystallographically determined water molecules are allowed in the binding site, it becomes much easier to generate ligands in larger numbers and with wider chemical diversity. This example shows that, in some cases, bound water molecules can be more accessible for hydrogen bonding to an incoming ligand than the actual protein binding sitepoints associated with them. From the point of view of de novo ligand design, water molecules can thus act as versatile amphiprotic hydrogen-bonding sitepoints and reduce the conformational constraints of a particular binding site.

Bacteria↗

PubMedDate not supplied

Including tightly-bound water molecules in de novo drug design. Exemplification through the in silico generation of poly(ADP-ribose)polymerase ligands.

Different strategies for the in silico generation of ligand molecules in the binding site of poly(ADP-ribose)polymerase (PARP) were studied in order to observe the effect of the targeting and displacement of tightly bound water molecules. Several molecular scaffolds were identified as having better interactions in the binding site when targeting one or two tightly bound water molecules in the NAD binding site. Energy calculations were conducted in order to assess the ligand-protein and ligand-water-protein interactions of different functional groups of the generated ligands. These calculations were used to evaluate the energetic consequences of the presence of tightly bound water molecules and to identify those that contribute favorably to the binding of ligands.

Catalytic Domain↗

PubMedDate not supplied

Expanded interaction fingerprint method for analyzing ligand binding modes in docking and structure-based drug design.

An expanded interaction fingerprint method has been developed for analyzing the binding modes of ligands in docking and structure-based design methods. Taking the basic premise of representing a ligand in terms of a binary string that denotes its interactions with a target protein, we have expanded the method to include additional interaction-specific information. By considering the hydrogen-bonding strength and/or accessibility of the hydrogen bonding groups within a binding site as well as their geometric arrangement we aim to provide a better representation of a ligand-protein interaction. These expanded methods have been applied to the postprocessing of binding poses generated in a docking study for 220 different proteins and to the analysis of ligands generated by an automated ligand-generation algorithm for the anthrax oedema factor. In the docking study, the application of the interaction fingerprint method as a postprocessing tool resulted in an increased success rate in identifying the crystallographic binding mode. In the analysis of the ligands generated for the anthrax oedema factor, the incorporation of additional interaction-specific information resulted in a more intuitive and comprehensive analysis of automated ligand-generation output.

Journal Article↗