Vladimir Gelev

On-resin N-methylation of cyclic peptides for discovery of orally bioavailable scaffolds

Backbone N-methylation is common among peptide natural products and has a significant impact on both the physical properties and the conformational states of cyclic peptides. However, the specific impact of N-methylation on passive membrane diffusion in cyclic peptides has not been investigated systematically. Here we report a method for the selective, on-resin N-methylation of cyclic peptides to generate compounds with drug-like membrane permeability and oral bioavailability. The selectivity and degree of N-methylation of the cyclic peptide was determined by backbone stereochemistry, suggesting that conformation dictates the regiochemistry of the N-methylation reaction. The permeabilities of the N-methyl variants were corroborated by computational studies on a 1024-member virtual library of N-methyl cyclic peptides. One of the most permeable compounds, a cyclic hexapeptide (MW = 755) with three N-methyl groups, showed an oral bioavailability of 28% in rat.

Cell-Free Expressed Bacteriorhodopsin in Different Soluble Membrane Mimetics: Biophysical Properties and NMR Accessibility

Selecting a suitable membrane-mimicking environment is of fundamental importance for the investigation of membrane proteins. Nonconventional surfactants, such as amphipathic polymers (amphipols) and lipid bilayer nanodiscs, have been introduced as promising environments that may overcome intrinsic disadvantages of detergent micelle systems. However, structural insights into the effects of different environments on the embedded protein are limited. Here, we present a comparative study of the heptahelical membrane protein bacteriorhodopsin in detergent micelles, amphipols, and nanodiscs. Our results confirm that nonconventional environments can increase stability of functional bacteriorhodopsin, and demonstrate that well-folded heptahelical membrane proteins are, in principle, accessible by solution-NMR methods in amphipols and phospholipid nanodiscs. Our data distinguish regions of bacteriorhodopsin that mediate membrane/solvent contacts in the tested environments, whereas the proteins functional inner core remains almost unperturbed. The presented data allow comparing the investigated membrane mimetics in terms of NMR spectral quality and thermal stability required for structural studies.

DNA breathing dynamics in the presence of a terahertz field

We consider the influence of a terahertz field on the breathing dynamics of double-stranded DNA. We model the spontaneous formation of spatially localized openings of a damped and driven DNA chain, and find that linear instabilities lead to dynamic dimerization, while true local strand separations require a threshold amplitude mechanism. Based on our results we argue that a specific terahertz radiation exposure may significantly affect the natural dynamics of DNA, and thereby influence intricate molecular processes involved in gene expression and DNA replication.

Structure and mechanism of the influenza A M218–60 dimer of dimers

We report a magic angle spinning (MAS) NMR structure of the drug-resistant S31N mutation of M218-60 from Influenza A. The protein was dispersed in diphytanoyl-sn-glycero-3-phosphocholine lipid bilayers, and the spectra and an extensive set of constraints indicate that M218-60 consists of a dimer of dimers. In particular, ∼280 structural constraints were obtained using dipole recoupling experiments that yielded well-resolved (13)C-(15)N, (13)C-(13)C, and (1)H-(15)N 2D, 3D, and 4D MAS spectra, all of which show cross-peak doubling. Interhelical distances were measured using mixed (15)N/(13)C labeling and with deuterated protein, MAS at ωr/2π = 60 kHz, ω0H/2π = 1000 MHz, and (1)H detection of methyl-methyl contacts. The experiments reveal a compact structure consisting of a tetramer composed of four transmembrane helices, in which two opposing helices are displaced and rotated in the direction of the membrane normal relative to a four-fold symmetric arrangement, yielding a two-fold symmetric structure. Side chain conformations of the important gating and pH-sensing residues W41 and H37 are found to differ markedly from four-fold symmetry. The rmsd of the structure is 0.7 Å for backbone heavy atoms and 1.1 Å for all heavy atoms. This two-fold symmetric structure is different from all of the previous structures of M2, many of which were determined in detergent and/or with shorter constructs that are not fully active. The structure has implications for the mechanism of H(+) transport since the distance between His and Trp residues on different helices is found to be short. The structure also exhibits two-fold symmetry in the vicinity of the binding site of adamantyl inhibitors, and steric constraints may explain the mechanism of the drug-resistant S31N mutation.

Molecular Framework for the Activation of RNA-dependent Protein Kinase

The RNA-dependent protein kinase (PKR) plays an integral role in the antiviral response to cellular infection. PKR contains three distinct domains consisting of two conserved N-terminal double-stranded RNA (dsRNA)-binding domains, a C-terminal Ser-Thr kinase domain, and a central 80-residue linker. Despite rich structural and biochemical data, a detailed mechanistic explanation of PKR activation remains unclear. Here we provide a framework for understanding dsRNA-dependent activation of PKR using nuclear magnetic resonance spectroscopy, dynamic light scattering, gel filtration, and autophosphorylation kinetics. In the latent state, PKR exists as an extended monomer, with an increase in self-affinity upon dsRNA association. Subsequent phosphorylation leads to efficient release of dsRNA followed by a greater increase in self-affinity. Activated PKR displays extensive conformational perturbations within the kinase domain. We propose an updated model for PKR activation in which the communication between RNA binding, central linker, and kinase domains is critical in the propagation of the activation signal and for PKR dimerization.

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.

DNA dynamics play a role as a basal transcription factor in the positioning and regulation of gene transcription initiation

We assess the role of DNA breathing dynamics as a determinant of promoter strength and transcription start site (TSS) location. We compare DNA Langevin dynamic profiles of representative gene promoters, calculated with the extended non-linear PBD model of DNA with experimental data on transcription factor binding and transcriptional activity. Our results demonstrate that DNA dynamic activity at the TSS can be suppressed by mutations that do not affect basal transcription factor binding–DNA contacts. We use this effect to establish the separate contributions of transcription factor binding and DNA dynamics to transcriptional activity. Our results argue against a purely ‘transcription factor-centric’ view of transcription initiation, suggesting that both DNA dynamics and transcription factor binding are necessary conditions for transcription initiation.

Toward a Detailed Description of the Thermally Induced Dynamics of the Core Promoter

Establishing the general and promoter-specific mechanistic features of gene transcription initiation requires improved understanding of the sequence-dependent structural/dynamic features of promoter DNA. Experimental data suggest that a spontaneous dsDNA strand separation at the transcriptional start site is likely to be a requirement for transcription initiation in several promoters. Here, we use Langevin molecular dynamic simulations based on the Peyrard-Bishop-Dauxois nonlinear model of DNA (PBD LMD) to analyze the strand separation (bubble) dynamics of 80-bp-long promoter DNA sequences. We derive three dynamic criteria, bubble probability, bubble lifetime, and average strand separation, to characterize bubble formation at the transcriptional start sites of eight mammalian gene promoters. We observe that the most stable dsDNA openings do not necessarily coincide with the most probable openings and the highest average strand displacement, underscoring the advantages of proper molecular dynamic simulations. The dynamic profiles of the tested mammalian promoters differ significantly in overall profile and bubble probability, but the transcriptional start site is often distinguished by large (longer than 10 bp) and long-lived transient openings in the double helix. In support of these results are our experimental transcription data demonstrating that an artificial bubble-containing DNA template is transcribed bidirectionally by human RNA polymerase alone in the absence of any other transcription factors.

DNA breathing dynamics distinguish binding from nonbinding consensus sites for transcription factor YY1 in cells

The genome-wide mapping of the major gene expression regulators, the transcription factors (TFs) and their DNA binding sites, is of great importance for describing cellular behavior and phenotypic diversity. Presently, the methods for prediction of genomic TF binding produce a large number of false positives, most likely due to insufficient description of the physiochemical mechanisms of protein–DNA binding. Growing evidence suggests that, in the cell, the double-stranded DNA (dsDNA) is subject to local transient strands separations (breathing) that contribute to genomic functions. By using site-specific chromatin immunopecipitations, gel shifts, BIOBASE data, and our model that accurately describes the melting behavior and breathing dynamics of dsDNA we report a specific DNA breathing profile found at YY1 binding sites in cells. We find that the genomic flanking sequence variations and SNPs, may exert long-range effects on DNA dynamics and predetermine YY1 binding. The ubiquitous TF YY1 has a fundamental role in essential biological processes by activating, initiating or repressing transcription depending upon the sequence context it binds. We anticipate that consensus binding sequences together with the related DNA dynamics profile may significantly improve the accuracy of genomic TF binding sites and TF binding-related functional SNPs.

The Use of Amphipols for NMR Structural Characterization of 7-TM Proteins

While amphipols have been proven useful for refolding of seven transmembrane helical (7-TM) proteins including G-protein-coupled receptors (GPCRs) and it could be shown that an amphipol environment is in principle suitable for NMR structural studies of the embedded protein, high-resolution NMR insights into amphipol refolded and isotopically labeled GPCRs are still very limited. Here we report on the recent progress toward NMR structural studies of the melanocortin-2 and -4 receptors, two class A GPCRs which so far have not been reported to be incorporated into an amphipol environment. Making use of the established 7-TM protein bacteriorhodopsin (BR) we initially tested and optimized amphipol refolding conditions. Most promising conditions were transferred to the refolding of the two melanocortin receptors. Analytical-scale refolding experiments on the melanocortin-2 receptor show very similar behavior to the results obtained on BR. Using cell-free protein expression we could generate sufficient amounts of isotopically labeled bacteriorhodopsin as well as melanocortin-2 and -4 receptors for an initial NMR analysis. Upscaling of the amphipol refolding protocol to protein amounts needed for NMR structural studies was, however, not straightforward and impeded detailed NMR insights for the two GPCRs. While well-resolved and dispersed NMR spectra could only be obtained for bacteriorhodopsin, a comparison of NMR data recorded on the melanocortin-4 receptor in SDS and in an amphipol environment indicates that amphipol refolding induces larger structural modifications in the receptor.