Reduction of conformational mobility and aggregation in W60G β2‐microglobulin: assessment by 15N NMR relaxation

The amyloid pathology associated with long‐term haemodialysis is due to the deposition of β₂‐microglobulin, the non‐polymorphic light chain of class I major histocompatibility complex, that accumulates at bone joints into amyloid fibrils. Several lines of evidence show the relevance of the tryptophan residue at position 60 for the fibrillogenic transition of the protein. A comparative ¹⁵N NMR relaxation analysis is presented for wild‐type human β₂‐microglobulin and W60G β₂‐microglobulin, i.e. the mutant with a glycyne replacing the natural tryptophan residue at position 60. The experimental data, collected at 11.4 T and 310 K, were analyzed by means of the reduced spectral density approach. Molecular dynamics (MD) simulations and corresponding thermodynamic integration, together with hydrodynamic calculations were performed to support data interpretation. The analysis results for the mutant protein are consistent with a reduced aggregation with respect to the wild‐type counterpart, as a consequence of an increased conformational rigidity probed by either NMR relaxation and MD simulations. Although dynamics in solution is other than fibrillar competence, the assessed properties of the mutant protein can be related with its reduced ability of forming fibrils when seeded in 20% trifluoroethanol. Copyright © 2013 John Wiley & Sons, Ltd.     

Fast,metadynamics‐based method for prediction of the stereochemistry‐dependent relative free energies of ligand–receptor interactions

The computational approach applicable for the molecular dynamics (MD)‐based techniques is proposed to predict the ligand–protein binding affinities dependent on the ligand stereochemistry. All possible stereoconfigurations are expressed in terms of one set of force‐field parameters [stereoconfiguration‐independent potential (SIP)], which allows for calculating all relative free energies by only single simulation. SIP can be used for studying diverse, stereoconfiguration‐dependent phenomena by means of various computational techniques of enhanced sampling. The method has been successfully tested on the β₂‐adrenergic receptor (β₂‐AR) binding the four fenoterol stereoisomers by both metadynamics simulations and replica‐exchange MD. Both the methods gave very similar results, fully confirming the presence of stereoselective effects in the fenoterol‐β₂‐AR interactions. However, the metadynamics‐based approach offered much better efficiency of sampling which allows for significant reduction of the unphysical region in SIP. © 2014 Wiley Periodicals, Inc.     

Functional Dynamics of Deuterated β2‐Adrenergic Receptor in Lipid Bilayers Revealed by NMR Spectroscopy

G‐protein‐coupled receptors (GPCRs) exist in conformational equilibrium between active and inactive states, and the former population determines the efficacy of signaling. However, the conformational equilibrium of GPCRs in lipid bilayers is unknown owing to the low sensitivities of their NMR signals. To increase the signal intensities, a deuteration method was developed for GPCRs expressed in an insect cell/baculovirus expression system. The NMR sensitivities of the methionine methyl resonances from the β₂‐adrenergic receptor (β₂AR) in lipid bilayers of reconstituted high‐density lipoprotein (rHDL) increased by approximately 5‐fold upon deuteration. NMR analyses revealed that the exchange rates for the conformational equilibrium of β₂AR in rHDLs were remarkably different from those measured in detergents. The timescales of GPCR signaling, calculated from the exchange rates, are faster than those of receptor tyrosine kinases and thus enable rapid neurotransmission and sensory perception.     

Site‐Specific Investigation of the Steady‐State Kinetics and Dynamics of the Multistep Binding of Bile Acid Molecules to a Lipid Carrier Protein

The investigation of multi‐site ligand–protein binding and multi‐step mechanisms is highly demanding. In this work, advanced NMR methodologies such as 2D ¹H–¹⁵N line‐shape analysis, which allows a reliable investigation of ligand binding occurring on micro‐ to millisecond timescales, have been extended to model a two‐step binding mechanism. The molecular recognition and complex uptake mechanism of two bile salt molecules by lipid carriers is an interesting example that shows that protein dynamics has the potential to modulate the macromolecule–ligand encounter. Kinetic analysis supports a conformational selection model as the initial recognition process in which the dynamics observed in the apo form is essential for ligand uptake, leading to conformations with improved access to the binding cavity. Subsequent multi‐step events could be modelled, for several residues, with a two‐step binding mechanism. The protein in the ligand‐bound state still exhibits a conformational rearrangement that occurs on a very slow timescale, as observed for other proteins of the family. A global mechanism suggesting how bile acids access the macromolecular cavity is thus proposed.     

Conformational Selection in Glycomimetics: Human Galectin‐1 Only Recognizes syn‐Ψ‐Type Conformations of β‐1,3‐Linked Lactose and Its C‐Glycosyl Derivative

The human lectin galectin‐1 (hGal‐1) translates sugar signals, that is, β‐galactosides, into effects on the level of cells, for example, growth regulation, and has become a model for studying binding of biopharmaceutically relevant derivatives. Bound‐state conformations of Galβ‐C‐(1→3)‐Glcβ‐OMe ( 1 ) and its βGal‐(1→3)‐βGlc‐OMe disaccharide parent compound were studied by using NMR spectroscopy (transferred (TR)‐NOESY data), assisted by docking experiments and molecular dynamics (MD) simulations. The molecular recognition process involves a conformational selection event. Although free C‐glycoside access four distinct conformers in solution, hGal‐1 recognizes shape of a local minimum of compound 1 , the syn‐Φ/syn‐Ψ conformer, not the structure at global minimum. MD simulations were run to explain, in structural terms, the observed geometry of the complex.     

On the Use of a Supramolecular Coarse‐Grained Model for the Solvent in Simulations of the Folding Equilibrium of an Octa‐β‐peptide in MeOH and H2O

Molecular dynamics (MD) simulation can give a detailed picture of conformational equilibria of biomolecules, but it is only reliable if the force field used in the simulation is accurate, and the sampling of the conformational space accessible to the biomolecule shows many (un)folding transitions to allow for precise averages of observable quantities. Here, the use of coarse‐grained (CG) solvent MeOH and H₂O models to speed up the sampling of the conformational equilibria of an octa‐β‐peptide is investigated. This peptide is thought to predominantly adopt a 3₁₄‐helical fold when solvated in MeOH, and a hairpin fold when solvated in H₂O on the basis of the NMR data. Various factors such as the chirality of a residue, a force‐field modification for the solute, coarse‐graining of the solvent model, and an extension of the nonbonded interaction cut‐off radius are shown to influence the simulated conformational equilibria and the agreement with the experimental NMR data for the octa‐β‐peptide.     

A rare occurrence of a β‐turn in an amyloid βA4 peptide

The fragment β(25–35) of the amyloid β‐peptide, like its parent βA4, has shown neurotrophic and late neurotoxic activities in cultured cells. The 3D structure of this important peptide was examined by ¹H and ¹³C 2D‐NMR and MD simulations in DMSO‐d₆ and water. The NMR parameters of chemical shift, ³J(N,Hα) coupling constants, temperature coefficients of NH chemical shifts and the pattern of intra and inter‐residue NOEs were used to deduce the structures. In DMSO‐d₆, the peptide was found to take up a type I β‐turn around the C‐terminal residues Ile⁸–Gly⁹–Leu¹⁰–Met¹¹, whereas in water at pH 5.5, it adopts a random coil conformation. This is only the second report of a β‐turn in the β‐amyloid class of peptides. The solution structures generated using restrained molecular dynamics were refined by MARDIGRAS to an R factor of 0.33 in the case of DMSO‐d₆ and to 0.56 for water. Copyright © 2002 John Wiley & Sons, Ltd.     

Structural fluctuation and concerted motions in F1‐ATPase: A molecular dynamics study

F₁‐ATPase is an adenosine tri‐phosphate (ATP)‐driven rotary motor enzyme. We investigated the structural fluctuations and concerted motions of subunits in F₁‐ATPase using molecular dynamics (MD) simulations. An MD simulation for the α₃β₃γ complex was carried out for 30 ns. Although large fluctuations of the N‐terminal domain observed in simulations of the isolated β_E subunit were suppressed in the complex simulation, the magnitude of fluctuations in the C‐terminal domain was clearly different among the three β subunits (β_E, β_TP, and β_DP). Despite fairly similar conformations of the β_TP and β_DP subunits, the β_DP subunit exhibits smaller fluctuations in the C‐terminal domain than the β_TP subunit due to their dissimilar interface configurations. Compared with the β_TP subunit, the β_DP subunit stably interacts with both the adjacent α_DP and α_E subunits. This sandwiched configuration in the β_DP subunit leads to strongly correlated motions between the β_DP and adjacent α subunits. The β_DP subunit exhibits an extensive network of highly correlated motions with bound ATP and the γ subunit, as well as with the adjacent α subunits, suggesting that the structural changes occurring in the catalytically active β_DP subunit can effectively induce movements of the γ subunit. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

The Effect of Fluoro Substitution upon the β‐Hairpin Fold of a β‐Tetrapeptide in Methanol

The importance of β‐peptides lies in their ability to mimic the conformational behavior of α‐peptides, even with a much shorter chain length, and in their resistance to proteases. To investigate the effect of substitution of β‐peptides on their dominant fold, we have carried out a molecular‐dynamics (MD) simulation study of two tetrapeptides, Ac‐(2R,3S)‐β^2,3hVal(αMe)‐(2S)‐β²hPhe‐(R)‐β³hLys‐(2R,3S)‐β^2,3‐Ala(αMe)‐NH₂, differing in the substitution at the C_α of Phe2 (pepF with F, and pepH with H). Three simulations, unrestrained (UNRES), using ³J‐coupling biasing with local elevation in combination with either instantaneous (INS) or time‐averaging (AVE) NOE distance restraining, were carried out for each peptide. In the unrestrained simulations, we find three (pepF) and two (pepH) NOE distance bound violations of maximally 0.22 nm that involve the terminal residues. The restrained simulations match both the NOE distance bounds and ³J‐values derived from experiment. The fluorinated peptide shows a slightly larger conformational variability than the non‐fluorinated one.     

High‐expression β1 adrenergic receptor/cell membrane chromatography method based on a target receptor to screen active ingredients from traditional Chinese medicines

β‐Adrenergic receptors are important targets for drug discovery. We have developed a new β₁‐adrenergic receptor cell membrane chromatography (β₁AR‐CMC) with offline ultra‐performance LC (UPLC) and MS method for screening active ingredients from traditional Chinese medicines. In this study, Chinese hamster ovary‐S cells with high β₁AR expression levels were established and used to prepare a cell membrane stationary phase in a β₁AR‐CMC model. The retention fractions were separated and identified by the UPLC–MS system. The screening results found that isoimperatorin from Rhizoma et Radix Notopterygii was the targeted component that could act on β₁AR in similar manner of metoprolol as a control drug. In addition, the biological effects of active component were also investigated in order to search for a new type of β₁AR antagonist. It will be a useful method for drug discovery as a leading compound resource.     

Validation of the GROMOS 54A7 Force Field Regarding Mixed α/β‐Peptide Molecules

A molecular‐dynamics (MD) simulation study of two heptapeptides containing α‐ and β‐amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β‐hairpin‐like conformations, while peptide 8 adopts a 14/15‐helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atom? atom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and ³J‐coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β‐peptides in MeOH.     

Packing forces in dichloridobis(3,5‐diphenyl‐1H‐pyrazole‐κN2)cobalt(II) dichloromethane hemisolvate

The title compound, [CoCl₂(C₁₅H₁₂N₂)₂]·0.5CH₂Cl₂, was crystallized from a binary mixture of dichloromethane and hexane and a dimeric supramolecular structure was isolated. The Co^II centre exhibits a distorted tetrahedral geometry, with two independent pyrazole‐based ligands occupying two coordination sites and two chloride ligands occupying the third and fourth coordination sites. The supramolecular structure is supported by complementary hydrogen bonding between the pyrazole NH group and the chloride ligand of an adjacent molecule. This hydrogen‐bonding motif yields a ten‐membered hydrogen‐bonded ring. Density functional theory (DFT) simulations at the PBE/6‐311G level of theory were used to probe the solid‐state structure. These simulations suggest that the chelate undergoes a degree of conformational distortion from the lowest‐energy geometry to allow for optimal hydrogen bonding in the solid state.     

Synthesis and Conformational Analysis of Hybrid α/β‐Dipeptides Incorporating S‐Glycosyl‐β2,2‐Amino Acids

We synthesized and carried out the conformational analysis of several hybrid dipeptides consisting of an α‐amino acid attached to a quaternary glyco‐β‐amino acid. In particular, we combined a S‐glycosylated β^2,2‐amino acid and two different types of α‐amino acid, namely, aliphatic (alanine) and aromatic (phenylalanine and tryptophan) in the sequence of hybrid α/β‐dipeptides. The key step in the synthesis involved the ring‐opening reaction of a chiral cyclic sulfamidate, inserted in the peptidic sequence, with a sulfur‐containing nucleophile by using 1‐thio‐β‐D ‐glucopyranose derivatives. This reaction of glycosylation occurred with inversion of configuration at the quaternary center. The conformational behavior in aqueous solution of the peptide backbone and the glycosidic linkage for all synthesized hybrid glycopeptides was analyzed by using a protocol that combined NMR experiments and molecular dynamics with time‐averaged restraints (MD‐tar). Interestingly, the presence of the sulfur heteroatom at the quaternary center of the β‐amino acid induced θ torsional angles close to 180° (anti). Notably, this value changed to 60° (gauche) when the peptidic sequence displayed aromatic α‐amino acids due to the presence of CH–π interactions between the phenyl or indole ring and the methyl groups of the β‐amino acid unit.     

Binding structures of tri‐N‐acetyl‐β‐glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well‐studied enzyme that hydrolyzes the β‐(1,4)‐glycosidic linkage of N‐acetyl‐β‐glucosamine (NAG)_n oligomers. The active site of hen egg‐white lysozyme (HEWL) is believed to consist of six subsites, A‐F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all‐atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)₃. MD trajectories are applied to analyze structures and conformation of the complex as well as protein–ligand interactions, including the hydrogen‐bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)₃ are investigated independently based on a fixed‐charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM‐GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root‐mean‐squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed‐charge models. © 2012 Wiley Periodicals, Inc.     

Stereoselectivity in Diphos Addition to Molybdenum Complex with Pyridine‐2‐thionate (pyS) and η3‐Methallyl Coordinated Ligands

The conformational isomers endo‐ and exo‐[Mo{η³‐C₃H₄(CH₃)}(η²‐pyS)(CO)(η²‐diphos)] (diphos: dppm = {bis(diphenylphosphino)methane}, 2 ; dppe = {1,2‐bis(diphenylphosphino)ethane}, 3 ) are prepared by reacting the double‐bridged pyridine‐2‐thionate (pyS) complex [Mo{η³‐C₃H₄(CH₃)}(CO)₂]₂(η¹:η²:μ‐pyS)₂, 1 with diphos in refluxing acetonitrile. Stereoselectivity of the methallyl, C₃H₄(CH₃), ligand improves the formation of the exo‐conformation of 2 and 3 . Orientations and spectroscopy of these complexes are discussed.     

Monitoring Conformational Changes in the NDM‐1 Metallo‐β‐lactamase by 19F NMR Spectroscopy

The New Delhi metallo‐β‐lactamase (NDM‐1) is involved in the emerging antibiotic resistance problem. Development of metallo‐β‐lactamases (MBLs) inhibitors has proven challenging, due to their conformational flexibility. Here we report site‐selective labeling of NDM‐1 with 1,1,1‐trifluoro‐3‐bromo acetone (BFA), and its use to study binding events and conformational changes upon ligand–metal binding using ¹⁹F NMR spectroscopy. The results demonstrate different modes of binding of known NDM‐1 inhibitors, including L ‐ and D ‐captopril by monitoring the changing chemical environment of the active‐site loop of NDM‐1. The method described will be applicable to other MBLs and more generally to monitoring ligand‐induced conformational changes.     

Numerical Simulation of the Pressure Denaturation of a Helical β‐Peptide Heptamer Solvated in Methanol

The effect of elevated pressure on the conformational behavior of a β‐peptide heptamer ( 1 ) in MeOH solution was considered. The response of the peptide to elevated pressure was probed by means of molecular dynamics (MD) simulations, and described in atomic terms. The most‐striking features of the response are that the region of the ‘unfolded’ state of the peptide accessible at elevated pressure is narrow, and that thermal and pressure denaturation produce similar ‘unfolded’ states in the case of 1 .     

Synthesis and Conformational Analysis of 1′‐ and 3′‐Substituted 2‐Deoxy‐2‐fluoro‐D‐ribofuranosyl Nucleosides

Convergent syntheses of the 9‐(3‐X‐2,3‐dideoxy‐2‐fluoro‐β‐D ‐ribofuranosyl)adenines 5 (X=N₃) and 7 (X=NH₂), as well as of their respective α‐anomers 6 and 8 , are described, using methyl 2‐azido‐5‐O‐benzoyl‐2,3‐dideoxy‐2‐fluoro‐β‐D ‐ribofuranoside ( 4 ) as glycosylating agent. Methyl 5‐O‐benzoyl‐2,3‐dideoxy‐2,3‐difluoro‐β‐D ‐ribofuranoside ( 12 ) was prepared starting from two precursors, and coupled with silylated N⁶‐benzoyladenine to afford, after deprotection, 2′,3′‐dideoxy‐2′,3′‐difluoroadenosine ( 13 ). Condensation of 1‐O‐acetyl‐3,5‐di‐O‐benzoyl‐2‐deoxy‐2‐fluoro‐β‐D ‐ribofuranose ( 14 ) with silylated N²‐palmitoylguanine gave, after chromatographic separation and deacylation, the N⁷‐β‐anomer 17 as the main product, along with 2′‐deoxy‐2′‐fluoroguanosine ( 15 ) and its N⁹‐α‐anomer 16 in a ratio of ca. 42 : 24 : 10. An in‐depth conformational analysis of a number of 2,3‐dideoxy‐2‐fluoro‐3‐X‐D ‐ribofuranosides (X=F, N₃, NH₂, H) as well as of purine and pyrimidine 2‐deoxy‐2‐fluoro‐D ‐ribofuranosyl nucleosides was performed using the PSEUROT (version 6.3) software in combination with NMR studies.     

Electronic spectra and first‐order hyperpolarizability of 4‐(dicyanomethylene)‐2,6‐bis‐ (2′‐thiophene‐vinyl)‐pyran derivatives

On the basis of the ZINDO program, we have designed a program to calculate the first‐order hyperpolarizability β_ijk and β_μ according to the sum‐over‐states (SOS) expression. The first‐order hyperpolarizability of 4‐(dicyanomethylene)‐2,6‐bis‐(2′‐thiophene‐vinyl)‐pyran derivatives were studied. The calculated results were that the 4‐(dicyanomethylene)‐2,6‐bis‐(2′‐thiophene‐vinyl)‐pyran derivatives exhibit good nonlinearity with their β₀ values, which are slightly less than that of the corresponding 2,6‐bis‐styryl‐4‐(dicyanomethylene)‐pyran derivatives. It does not agree with the auxiliary donor–acceptor effects theory. The 4‐(dicyanomethylene)‐2,6‐bis‐(2′‐thiophene‐vinyl)‐pyran derivatives, having two low‐lying electronic excited states that contribute to the molecular hyperpolarizability in an additive manner, are good candidates as chromophores due to their high nonlinearities and good thermal stability. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 65–72, 2001     

Using open‐tubular capillary electrochromatography with part‐coating column for binding constants determination of β2‐adrenergic receptor with seven drugs

An open‐tubular capillary electrochromatography method has been developed for the determination of binding constants between β₂‐adrenergic receptor (β₂‐AR) and seven drugs. β₂‐AR was oriented immobilized onto one part of inner surface of capillary via microwave‐assisted technical synthesis. According to the linear relationship between coating length and the apparent mobility of analyte, the binding constant (K_b) can be obtained by related theories and equations. The order of K_b values between drugs such as adrenaline hydrochloride, norepinephrine bitartrate, and propranolol hydrochloride with β₂‐AR is well consistent with that reported in the literature. By the method, K_b values between four extracts of Radix Paeoniae Rubra and β₂‐AR were also successfully obtained. Subsequently, computer models were applied to interpret the CEC experiments. And the results proved to be in good agreement with the method. The work, herein, demonstrates the potential of the method in drug‐receptor affinity interactions evaluation and screening of lead compounds from natural sources.