Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.

Site-directed mutagenesis and molecular dynamics simulations of the alpha 1B-adrenergic receptor (AR) were combined to explore the potential molecular changes correlated with the transition from R (inactive state) to R (active state). Using molecular dynamics analysis we compared the structural/dynamic features of constitutively active mutants with those of the wild type and of an inactive alpha 1B-AR to build a theoretical model which defines the essential features of R and R. The results of site-directed mutagenesis were in striking agreement with the predictions of the model supporting the following hypothesis. (i) The equilibrium between R and R depends on the equilibrium between the deprotonated and protonated forms, respectively, of D142 of the DRY motif. In fact, replacement of D142 with alanine confers high constitutive activity to the alpha 1B-AR. (ii) The shift of R143 of the DRY sequence out of a conserved 'polar pocket' formed by N63, D91, N344 and Y348 is a feature common to all the active structures, suggesting that the role of R143 is fundamental for mediating receptor activation. Disruption of these intramolecular interactions by replacing N63 with alanine constitutively activates the alpha 1B-AR. Our findings might provide interesting generalities about the activation process of G protein-coupled receptors.

Computational fluid dynamics of the right ventricular outflow tract and of the pulmonary artery: a bench model of flow dynamics.

OBJECTIVES: The reconstruction of the right ventricular outflow tract (RVOT) with valved conduits remains a challenge. The reoperation rate at 5 years can be as high as 25% and depends on age, type of conduit, conduit diameter and principal heart malformation. The aim of this study is to provide a bench model with computer fluid dynamics to analyse the haemodynamics of the RVOT, pulmonary artery, its bifurcation, and left and right pulmonary arteries that in the future may serve as a tool for analysis and prediction of outcome following RVOT reconstruction. METHODS: Pressure, flow and diameter at the RVOT, pulmonary artery, bifurcation of the pulmonary artery, and left and right pulmonary arteries were measured in five normal pigs with a mean weight of 24.6 ± 0.89 kg. Data obtained were used for a 3D computer fluid-dynamics simulation of flow conditions, focusing on the pressure, flow and shear stress profile of the pulmonary trunk to the level of the left and right pulmonary arteries. RESULTS: Three inlet steady flow profiles were obtained at 0.2, 0.29 and 0.36 m/s that correspond to the flow rates of 1.5, 2.0 and 2.5 l/min flow at the RVOT. The flow velocity profile was constant at the RVOT down to the bifurcation and decreased at the left and right pulmonary arteries. In all three inlet velocity profiles, low sheer stress and low-velocity areas were detected along the left wall of the pulmonary artery, at the pulmonary artery bifurcation and at the ostia of both pulmonary arteries. CONCLUSIONS: This computed fluid real-time model provides us with a realistic picture of fluid dynamics in the pulmonary tract area. Deep shear stress areas correspond to a turbulent flow profile that is a predictive factor for the development of vessel wall arteriosclerosis. We believe that this bench model may be a useful tool for further evaluation of RVOT pathology following surgical reconstructions.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.

Mutational analysis of the highly conserved arginine within the Glu/Asp-Arg-Tyr motif of the alpha(1b)-adrenergic receptor: effects on receptor isomerization and activation.

We have suggested previously that both the negatively and positively charged residues of the highly conserved Glu/Asp-Arg-Tyr (E/DRY) motif play an important role in the activation process of the alpha(1b)-adreneric receptor (AR). In this study, R143 of the E/DRY sequence in the alpha(1b)-AR was mutated into several amino acids (Lys, His, Glu, Asp, Ala, Asn, and Ile). The charge-conserving mutation of R143 into lysine not only preserved the maximal agonist-induced response of the alpha(1b)-AR, but it also conferred high degree of constitutive activity to the receptor. Both basal and agonist-induced phosphorylation levels were significantly increased for the R143K mutant compared with those of the wild-type receptor. Other substitutions of R143 resulted in receptor mutants with either a small increase in constitutive activity (R143H and R143D), impairment (R143H, R143D), or complete loss of receptor-mediated response (R143E, R143A, R143N, R143I). The R413E mutant displayed a small, but significant increase in basal phosphorylation despite being severely impaired in receptor-mediated response. Interestingly, all the arginine mutants displayed increased affinity for agonist binding compared with the wild-type alpha(1b)-AR. A correlation was found between the extent of the affinity shift and the intrinsic activity of the agonists. The analysis of the receptor mutants using the allosteric ternary complex model in conjunction with the results of molecular dynamics simulations on the receptor models support the hypothesis that mutations of R143 can drive the isomerization of the alpha(1b)-AR into different states, highlighting the crucial role of this residue in the activation process of the receptor.

Theoretical study on receptor-G protein recognition: new insights into the mechanisms of the alpha1b-adrenergic recptor activation

This work compares the structural/dynamics features of the wild-type alb-adrenergic receptor (AR) with those of the D142A active mutant and the agonist-bound state. The two active receptor forms were compared in their isolated states as well as in their ability to form homodimers and to recognize the G alpha q beta 1 gamma 2 heterotrimer. The analysis of the isolated structures revealed that, although the mutation- and agonist-induced active states of the alpha 1b-AR are different, they, however, share several structural peculiarities including (a) the release of some constraining interactions found in the wild-type receptor and (b) the opening of a cytosolic crevice formed by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. Accordingly, also their tendency to form homodimers shows commonalties and differences. In fact, in both the active receptor forms, helix 6 plays a crucial role in mediating homodimerization. However, the homodimeric models result from different interhelical assemblies. On the same line of evidence, in both of the active receptor forms, the cytosolic opened crevice recognizes similar domains on the G protein. However, the docking solutions are differently populated and the receptor-G protein preorientation models suggest that the final complexes should be characterized by different interaction patterns.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Three-dimensional models of 14α-sterol demethylase (Cyp51A) from Aspergillus lentulus and Aspergillus fumigatus: an insight into differences in voriconazole interaction.

Aspergillus lentulus, an Aspergillus fumigatus sibling species, is increasingly reported in corticosteroid-treated patients. Its clinical significance is unknown, but the fact that A. lentulus shows reduced antifungal susceptibility, mainly to voriconazole, is of serious concern. Heterologous expression of cyp51A from A. fumigatus and A. lentulus was performed in Saccharomyces cerevisiae to assess differences in the interaction of Cyp51A with the azole drugs. The absence of endogenous ERG11 was efficiently complemented in S. cerevisiae by the expression of either Aspergillus cyp51A allele. There was a marked difference between azole minimum inhibitory concentration (MIC) values of the clones expressing each Aspergillus spp. cyp51A. Saccharomyces cerevisiae clones expressing A. lentulus alleles showed higher MICs to all of the azoles tested, supporting the hypothesis that the intrinsic azole resistance of A. lentulus could be associated with Cyp51A. Homology models of A. fumigatus and A. lentulus Cyp51A protein based on the crystal structure of Cyp51p from Mycobacterium tuberculosis in complex with fluconazole were almost identical owing to their mutual high sequence identity. Molecular dynamics (MD) was applied to both three-dimensional protein models to refine the homology modelling and to explore possible differences in the Cyp51A-voriconazole interaction. After 20ns of MD modelling, some critical differences were observed in the putative closed form adopted by the protein upon voriconazole binding. A closer study of the A. fumigatus and A. lentulus voriconazole putative binding site in Cyp51A suggested that some major differences in the protein's BC loop could differentially affect the lock-up of voriconazole, which in turn could correlate with their different azole susceptibility profiles.

Amino acids of the alpha1B-adrenergic receptor involved in agonist binding: differences in docking catecholamines to receptor subtypes.

Site-directed mutagenesis and molecular dynamics analysis of the 3-D model of the alpha1B-adrenergic receptor (AR) were combined to identify the molecular determinants of the receptor involved in catecholamine binding. Our results indicate that the three conserved serines in the fifth transmembrane domain (TMD) of the alpha1B-AR play a distinct role in catecholamine binding versus receptor activation. In addition to the amino acids D125 in TMDIII and S207 in TMDV directly involved in ligand binding, our findings identify a large number of polar residues playing an important role in the activation process of the alpha1B-AR thus providing new insights into the structure/function relationship of G protein-coupled receptors.

Localization of breakage points in knotted strings

It is a common macroscopic observation that knotted ropes or fishing lines under tension easily break at the knot. However, a more precise localization of the breakage point in knotted macroscopic strings is a difficult task. In the present work, the tightening of knots was numerically simulated, a comparison of strength of different knots was experimentally performed and a high velocity camera was used to precisely localize the site where knotted macroscopic strings break. In the case of knotted spaghetti, the breakage occurs at the position with high curvature at the entry to the knot. This localization results from joint contributions of loading, bending and friction forces into the complex process of knot breakage. The present simulations and experiments are in agreement with recent molecular dynamics simulations of a knotted polymer chain and with experiments performed on actin and DNA filaments. The strength of the knotted string is greatly reduced (down to 50%) by the presence of a knot, therefore reducing the resistance to tension of all materials containing chains of any sort. The present work with macroscopic strings revels some important aspects, which are not accessible by experiments with microscopic chains.

The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling.

The resistance of mosquitoes to chemical insecticides is threatening vector control programmes worldwide. Cytochrome P450 monooxygenases (CYPs) are known to play a major role in insecticide resistance, allowing resistant insects to metabolize insecticides at a higher rate. Among them, members of the mosquito CYP6Z subfamily, like Aedes aegypti CYP6Z8 and its Anopheles gambiae orthologue CYP6Z2, have been frequently associated with pyrethroid resistance. However, their role in the pyrethroid degradation pathway remains unclear. In the present study, we created a genetically modified yeast strain overexpressing Ae. aegypti cytochrome P450 reductase and CYP6Z8, thereby producing the first mosquito P450-CPR (NADPH-cytochrome P450-reductase) complex in a yeast recombinant system. The results of the present study show that: (i) CYP6Z8 metabolizes PBAlc (3-phenoxybenzoic alcohol) and PBAld (3-phenoxybenzaldehyde), common pyrethroid metabolites produced by carboxylesterases, producing PBA (3-phenoxybenzoic acid); (ii) CYP6Z8 transcription is induced by PBAlc, PBAld and PBA; (iii) An. gambiae CYP6Z2 metabolizes PBAlc and PBAld in the same way; (iv) PBA is the major metabolite produced in vivo and is excreted without further modification; and (v) in silico modelling of substrate-enzyme interactions supports a similar role of other mosquito CYP6Zs in pyrethroid degradation. By playing a pivotal role in the degradation of pyrethroid insecticides, mosquito CYP6Zs thus represent good targets for mosquito-resistance management strategies.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Comparison of T1 relaxation times of the neurochemical profile in rat brain at 9.4 tesla and 14.1 tesla.

Knowledge of T(1) relaxation times can be important for accurate relative and absolute quantification of brain metabolites, for sensitivity optimizations, for characterizing molecular dynamics, and for studying changes induced by various pathological conditions. (1)H T(1) relaxation times of a series of brain metabolites, including J-coupled ones, were determined using a progressive saturation (PS) technique that was validated with an adiabatic inversion-recovery (IR) method. The (1)H T(1) relaxation times of 16 functional groups of the neurochemical profile were measured at 14.1T and 9.4T. Overall, the T(1) relaxation times found at 14.1T were, within the experimental error, identical to those at 9.4T. The T(1)s of some coupled spin resonances of the neurochemical profile were measured for the first time (e.g., those of gamma-aminobutyrate [GABA], aspartate [Asp], alanine [Ala], phosphoethanolamine [PE], glutathione [GSH], N-acetylaspartylglutamate [NAAG], and glutamine [Gln]). Our results suggest that T(1) does not increase substantially beyond 9.4T. Furthermore, the similarity of T(1) among the metabolites (approximately 1.5 s) suggests that T(1) relaxation time corrections for metabolite quantification are likely to be similar when using rapid pulsing conditions. We therefore conclude that the putative T(1) increase of metabolites has a minimal impact on sensitivity when increasing B(0) beyond 9.4T.