α-Bn-o-AMPA as a cis peptide bond mimic in somatostatin analogues

Summary In a previous communication we reported the racemic synthesis of the cis peptide bond mimic α-benzyl-o-aminomethylphenylacetic acid and its incorporation in the cyclic somatostatin analoguesc[α(R andS)Bn-o-AMPA-Phe⁷-d-Trp⁸-Lys⁹-Thr¹⁰]. Since the epimeric peptides exhibit different binding affinities, we completed the structure-activity study with an asymmetric synthesis. A model for the solution conformation ofc[α(R andS)Bn-o-AMPA-Phe⁷-d-Trp⁸-Lys⁹-Thr¹⁰] is proposed on the basis of a 2D NMR study in CD₃OH and restrained molecular dynamics.     

Promiscuity in protein-RNA interactions: conformational ensembles facilitate molecular recognition in the spliceosome: conformational diversity in U2AF⁶⁵ facilitates binding to diverse RNA sequences

Here I discuss findings that suggest a universal mechanism for proteins (and RNA) to recognize and interact with various binding partners by selectively binding to different conformations that pre‐exist in the free protein's conformational ensemble. The tandem RNA recognition motif domains of splicing factor U2AF⁶⁵ fluctuate in solution between a predominately closed conformation in which the RNA binding site of one of the domains is blocked, and a lowly populated open conformation in which both RNA binding pockets are accessible. RNA binding to U2AF⁶⁵ may thus occur through the weakly populated open conformation, and the binding interaction stabilizes the open conformation. The conformational diversity observed in U2AF⁶⁵ might also facilitate binding to diverse RNA sequences as found in the polypyrimidine tracts that help define 3′ splice sites. Similar binding pathways in other systems have important consequences in biological regulation, molecular evolution, and information storage.     

Investigations of Takeout proteins’ ligand binding and release mechanism using molecular dynamics simulation

Takeout (To) proteins exist in a diverse range of insect species. They are involved in many important processes of insect physiology and behaviors. As the ligand carriers, To proteins can transport the small molecule to the target tissues. However, ligand release mechanism of To proteins is unclear so far. In this contribution, the process and pathway of the ligand binding and release are revealed by conventional molecular dynamics simulation, steered molecular dynamics simulation and umbrella sampling methods. Our results show that the α4-side of the protein is the unique gate for the ligand binding and release. The structural analysis confirms that the internal cavity of the protein has high rigidity, which is in accordance with the recent experimental results. By using the potential of mean force calculations in combination with residue cross correlation calculation, we concluded that the binding between the ligand and To proteins is a process of conformational selection. Furthermore, the conformational changes of To proteins and the hydrophobic interactions both are the key factors for ligand binding and release.     

Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid–liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Bródka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water–water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.     

Molecular dynamics investigations of structural and functional changes in Bcl-2 induced by the novel antagonist BDA-366

Apoptosis is a fundamental biological phenomenon, in which anti- or proapoptotic proteins of the Bcl-2 family regulate a committed step. Overexpression of Bcl-2, the prototypical antiapoptotic protein in this family, is associated with therapy resistance in various human cancers. Accordingly, Bcl-2 inhibitors intended for cancer therapy have been developed, typically against the BH3 domain. Recent experimental evidences have shown that the antiapoptotic function of Bcl-2 is not immutable, and that BDA-366, a novel antagonist of the BH4 domain, converts Bcl-2 from a survival molecule to an inducer of cell death. In this study, the underlying mechanisms of this functional conversion were investigated by accelerated molecular dynamics simulation. Results revealed that Pro127 and Trp30 in the BH4 domain rotate to stabilize BDA-366 via π-π interactions, and trigger a series of significant conformational changes of the α3 helix. This rearrangement blocks the hydrophobic binding site (HBS) in the BH3 domain and further prevents binding of BH3-only proteins, which consequently allows the BH3-only proteins to activate the proapoptotic proteins. Analysis of binding free energy confirmed that BDA-366 cross-inhibits BH3-only proteins, implying negative cooperative effects across separate binding sites. The newly identified blocked conformation of the HBS along with the open to closed transition pathway revealed by this study advances the understanding of the Bcl-2 transition from antiapoptotic to proapoptotic function, and yielded new structural insights for novel drug design against the BH4 domain.

Communicated by Ramaswamy H. Sarma     

Molecular mechanisms of activation in CDK2

Cyclin-dependent kinases (CDKs) are enzymes involved in crucial cellular processes. Their biological activity is directly linked to their high conformational variability, which involves large protein conformational rearrangements. We present here the application of an enhancing sampling technique to the study of conformational transitions between the open and closed state of CDKs. The analysis of the conformational intermediates supports the idea that the process is regulated by two important protein regions, which sequentially rearrange in order to allow the protein to reach its final conformation. Furthermore, the two paths involve additional (minor) protein rearrangements which are specific to the paths. Our results show that our procedure can provide reasonable transition pathways between the two protein forms at a very reduced computational cost. The robustness and the simplicity of our approach make it of general application to describe virtually any macromolecular conformational transitions.     

Conformational studies of parathyroid hormone (PTH)/PTH‐related protein (PTHrP) point‐mutated hybrids

The N‐terminal 1–34 segments of both parathyroid hormone (PTH) and parathyroid hormone‐related protein (PTHrP) bind and activate the same membrane receptor in spite of major differences between the two hormones in their amino acid sequence. Recently, it was shown that in (1–34)PTH/PTHrP segmental hybrid peptides, the N‐terminal 1–14 segment of PTHrP is incompatible with the C‐terminal 15–34 region of PTH leading to substantial reduction in potency. The sites of incompatibility were identified as positions 5 in PTH and 19 in PTHrP. In the present paper we describe the synthesis, biological evaluation, and conformational characterization of two point‐mutated PTH/PTHrP 1–34 hybrids in which the arginine residues at positions 19 and 21 of the native sequence of PTHrP have been replaced by valine (hybrid V²¹) and glutamic acid (hybrid E¹⁹), respectively, taken from the PTH sequence. Hybrid V²¹ exhibits both high receptor affinity and biological potency, while hybrid E¹⁹ binds weakly and is poorly active. The conformational properties of the two hybrids were studied in aqueous solution containing dodecylphosphocholine (DPC) micelles and in water/2,2,2‐trifluoroethanol (TFE) mixtures. Upon addition of TFE or DPC micelles to the aqueous solution, both hybrids undergo a coil‐helix transition. The maximum helix content in 1 : 1 water/TFE, obtained by CD data for both hybrids, is ∼ 80%. In the presence of DPC micelles, the maximum helix content is ∼ 40%. The conformational properties of the two hybrids in the micellar system were further investigated by combined 2D‐nmr, distance geometry (DG), and molecular dynamics (MD) calculations. The common structural motif, consisting of two helical segments located at N‐ and C‐termini, was observed in both hybrids. However, the biologically potent hybrid V²¹ exhibits two flexible sites, centered at residues 12 and 19 and connecting helical segments, while the flexibility sites in the weakly active hybrid E¹⁹ are located at position 11 and in the sequence 20–26. Our findings support the hypothesis that the presence and location of flexibility points between helical segments are essential for enabling the active analogs to fold into the bioactive conformation upon interaction with the receptor. © 1999 John Wiley & Sons, Inc. Biopoly 50: 525–535, 1999     

Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.     

Molecular simulation study of CYP2B6 polymorphism with and without psoralen

CYP2B6 is a polymorphic enzyme with a large number of variants which may lead to functional changes in enzyme activity and substrate selectivity. In this study, CYP2B6 and its three variants with and without psoralen, a mechanism-based inactivator, were investigated using molecular simulation method. The obtained docking orientation of psoralen was in agreement with previously identified site of metabolism. Stability analysis showed that the three variants displayed more flexibility than CYP2B6.1, and CYP2B6.34 was the most flexible one without psoralen binding. However, in the presence of psoralen, CYP2B6.34 became more rigidity. Tunnel analysis indicates that the bottleneck change of tunnels may be correlated to the increased or decreased activity of variants. Binding free energy analysis shows that van der Waals interaction dominates the binding of psoralen. CYP2B6.34 has the highest affinity to psoralen with lowest binding free energy. Ile114, Phe115 and heme contribute largely to the binding of psoralen with CYP2B6.6, while Phe206 and Leu363 play important roles for CYP2B6.1 and CYP2B6.4. These computational observations suggest that the increased activity of CYP2B6.4 and reduced activity of CYP2B6.6 may be due to changes in regional structures.     

CDR-H3 loop ensemble in solution – conformational selection upon antibody binding

ABSTRACT

We analyzed pairs of protein-binding, peptide-binding and hapten-binding antibodies crystallized as complex and in the absence of the antigen with and without conformational differences upon binding in the complementarity-determining region (CDR)-H3 loop. Here, we introduce a molecular dynamics-based approach to capture a diverse conformational ensemble of the CDR-H3 loop in solution. The results clearly indicate that the inherently flexible CDR-H3 loop indeed needs to be characterized as a conformational ensemble. The conformational changes of the CDR-H3 loop in all antibodies investigated follow the paradigm of conformation selection, because we observe the experimentally determined binding competent conformation without the presence of the antigen within the ensemble of pre-existing conformational states in solution before binding. We also demonstrate for several examples that the conformation observed in the antibody crystal structure without antigen present is actually selected to bind the carboxyterminal tail region of the antigen-binding fragment (Fab). Thus, special care must be taken when characterizing antibody CDR-H3 loops by Fab X-ray structures, and the possibility that pre-existing conformations are present should always be considered.