Arginine–phosphate salt bridges in protein–DNA complexes: a Car–Parrinello study

Triple lysine and nucleosome-binding motifs of the viral IE19 protein are required for human cytomegalovirus S-phase infections

Herpesvirus genomes are maintained as extrachromosomal plasmids within the nuclei of infected cells. Some herpesviruses persist within dividing cells, putting the viral genome at risk of being lost to the cytoplasm during mitosis because karyokinesis (nuclear division) requires nuclear envelope breakdown. Oncogenic herpesviruses (and papillomaviruses) avoid genome loss during mitosis by tethering their genomes to cellular chromosomes, thereby ensuring viral genome uptake into newly formed nuclei. These viruses use viral proteins with DNA- and chromatin-binding capabilities to physically link viral and cellular genomes together in a process called tethering. The known viral tethering proteins of human papillomavirus (E2), Epstein-Barr virus (EBNA1), and Kaposi's sarcoma-associated herpesvirus (LANA) each contain two independent domains required for genome tethering, one that binds sequence specifically to the viral genome and another that binds to cellular chromatin. This latter domain is called a chromatin tethering domain (CTD). The human cytomegalovirus UL123 gene encodes a CTD that is required for the virus to productively infect dividing fibroblast cells within the S phase of the cell cycle, presumably by tethering the viral genome to cellular chromosomes during mitosis. The CTD-containing UL123 gene product that supports S-phase infections is the IE19 protein. Here, we define two motifs in IE19 required for S-phase infections: an N-terminal triple lysine motif and a C-terminal nucleosome-binding motif within the CTD.IMPORTANCEThe IE19 protein encoded by human cytomegalovirus (HCMV) is required for S-phase infection of dividing cells, likely because it tethers the viral genome to cellular chromosomes, thereby allowing them to survive mitosis. The mechanism through which IE19 tethers viral genomes to cellular chromosomes is not understood. For human papillomavirus, Epstein-Barr virus, and Kaposi's sarcoma-associated herpesvirus, viral genome tethering is required for persistence (latency) and pathogenesis (oncogenesis). Like these viruses, HCMV also achieves latency, and it modulates the properties of glioblastoma multiforme tumors. Therefore, defining the mechanism through which IE19 tethers viral genomes to cellular chromosomes may help us understand, and ultimately combat or control, HCMV latency and oncomodulation.

Neurobiochemical characteristics of arginine-rich peptides explain their potential therapeutic efficacy in neurodegenerative diseases

Neurodegenerative diseases, including Alzheimer̕ s disease (AD), Parkinson̕ s disease (PD), Huntington̕ s disease (HD), and Amyotrophic Lateral Sclerosis (ALS) require special attention to find new potential treatment methods. This review aims to summarize the current knowledge of the relationship between the biochemical properties of arginine-rich peptides (ARPs) and their neuroprotective effects to deal with the harmful effects of risk factors. It seems that ARPs have portrayed a promising and fantastic landscape for treating neurodegeneration-associated disorders. With multimodal mechanisms of action, ARPs play various unprecedented roles, including as the novel delivery platforms for entering the central nervous system (CNS), the potent antagonists for calcium influx, the invader molecules for targeting mitochondria, and the protein stabilizers. Interestingly, these peptides inhibit the proteolytic enzymes and block protein aggregation to induce pro-survival signaling pathways. ARPs also serve as the scavengers of toxic molecules and the reducers of oxidative stress agents. They also have anti-inflammatory, antimicrobial, and anti-cancer properties. Moreover, by providing an efficient nucleic acid delivery system, ARPs can play an essential role in developing various fields, including gene vaccines, gene therapy, gene editing, and imaging. ARP agents and ARP/cargo therapeutics can be raised as an emergent class of neurotherapeutics for neurodegeneration. Part of the aim of this review is to present recent advances in treating neurodegenerative diseases using ARPs as an emerging and powerful therapeutic tool. The applications and progress of ARPs-based nucleic acid delivery systems have also been discussed to highlight their usefulness as a broad-acting class of drugs.

Formation of versus Recruitment to RNA-Rich Condensates: Controlling Effects Exerted by Peptide Side Chain Identity

Liquid-liquid phase separation (LLPS), the spontaneous formation of contiguous liquid phases with distinct compositions, has been long known in chemical systems and more recently recognized as a ubiquitous feature of cell biology. We describe a system involving biologically relevant components, synthetic peptides, and total yeast RNA, that has enabled us to explore factors that underlie phase separation. Coulombic complementarity between a cationic peptide and anionic RNA is necessary but not sufficient for formation of a condensed phase in our system. In addition to a net positive charge, the peptide must present the proper type of cationic moiety. Guanidinium groups, as found in the Arg side chain, support phase separation, but ammonium groups, as found in the Lys side chain, or dimethylguanidinium groups, as found in post-translationally modified Arg side chains, do not support phase separation in our system. However, the cationic groups that do not support phase separation via interaction with RNA can nevertheless enable recruitment to a condensed phase, which reveals that the network of forces governing condensed phase formation can differ from the network of forces governing recruitment to such a phase. We introduce a new method for measuring the concentrations of components in condensed phases based on fluorine-containing additives and ¹⁹F NMR.

Infrequent metadynamics study of rare-event electrostatic channeling

The efficiency of cascade reactions, which consist of multiple chemical transformations that occur in a single pot without purification steps, is limited by the transport efficiency of intermediates between adjacent steps. Electrostatic channeling is a proven strategy for intermediate transfer in natural chemical cascades, but implementation into artificial cascades remains a challenge. Here, we combine infrequent metadynamics (InMetaD), umbrella sampling (US), and kinetic Monte Carlo (KMC) models to computationally study the transfer mechanism of glucose-6-phosphate (G6P) on a poly-arginine peptide bridging hexokinase (HK) and glucose-6-dehydrogenase (G6PDH). Transport of G6P by hopping in the presence of poly-arginine peptides is shown to be a rare event, and InMetaD is used to compute the hopping activation energy. US simulations capture the configurational change in the desorption process and enable the determination of the desorption energy. Parameterized by these results, a KMC model is used to estimate transport efficiency for the bridged enzyme complex. Results are compared to a similar complex using a poly-lysine bridge, using kinetic lag time as a metric. Even at a high ionic strength of 120 mM, poly-arginine peptides may be capable of more efficient transport as compared to poly-lysine, with a predicted lag time of 6 seconds for poly-arginine, compared to a previously reported lag time of 59 seconds for poly-lysine. This work indicates that poly-arginine peptides may be an improved bridge structure for electrostatic channeling of anionic intermediates.

Anti-EGFR VHH Antibody under Thermal Stress Is Better Solubilized with a Lysine than with an Arginine SEP Tag

In this study, we assessed the potential of arginine and lysine solubility-enhancing peptide (SEP) tags to control the solubility of a model protein, anti-EGFR V_HH-7D12, in a thermally denatured state at a high temperature. We produced V_HH-7D12 antibodies attached with a C-terminal SEP tag made of either five or nine arginines or lysines (7D12-C5R, 7D12-C9R, 7D12-C5K and 7D12-C9K, respectively). The 5-arginine and 5-lysine SEP tags increased the E. coli expression of V_HH-7D12 by over 80%. Biophysical and biochemical analysis confirmed the native-like secondary and tertiary structural properties and the monomeric nature of all V_HH-7D12 variants. Moreover, all V_HH-7D12 variants retained a full binding activity to the EGFR extracellular domain. Finally, thermal stress with 45-minute incubation at 60 and 75 °C, where V_HH-7D12 variants are unfolded, showed that the untagged V_HH-7D12 formed aggregates in all of the four buffers, and the supernatant protein concentration was reduced by up to 35%. 7D12-C5R and 7D12-C9R did not aggregate in Na-acetate (pH 4.7) and Tris-HCl (pH 8.5) but formed aggregates in phosphate buffer (PB, pH 7.4) and phosphate buffer saline (PBS, pH 7.4). The lysine tags (either C5K or C9K) had the strongest solubilization effect, and both 7D12-C5K and 7D12-C9K remained in the supernatant. Altogether, our results indicate that, under a thermal stress condition, the lysine SEP tags solubilization effect is more potent than that of an arginine SEP tags, and the SEP tags did not affect the structural and functional properties of the protein.

Overlapping Properties of the Short Membrane-Active Peptide BP100 With (i) Polycationic TAT and (ii) α-helical Magainin Family Peptides

BP100 is a short, designer-made membrane-active peptide with multiple functionalities: antimicrobial, cell-penetrating, and fusogenic. Consisting of five lysines and 6 hydrophobic residues, BP100 was shown to bind to lipid bilayers as an amphipathic α-helix, but its mechanism of action remains unclear. With these features, BP100 embodies the characteristics of two distinctly different classes of membrane-active peptides, which have been studied in detail and where the mechanism of action is better understood. On the one hand, its amphiphilic helical structure is similar to the pore forming magainin family of antimicrobial peptides, though BP100 is much too short to span the membrane. On the other hand, its length and high charge density are reminiscent of the HIV-TAT family of cell penetrating peptides, for which inverted micelles have been postulated as translocation intermediates, amongst other mechanisms. Assays were performed to test the antimicrobial and hemolytic activity, the induced leakage and fusion of lipid vesicles, and cell uptake. From these results the functional profiles of BP100, HIV-TAT, and the magainin-like peptides magainin 2, PGLa, MSI-103, and MAP were determined and compared. It is observed that the activity of BP100 resembles most closely the much longer amphipathic α-helical magainin-like peptides, with high antimicrobial activity along with considerable fusogenic and hemolytic effects. In contrast, HIV-TAT shows almost no antimicrobial, fusogenic, or hemolytic effects. We conclude that the amphipathic helix of BP100 has a similar membrane-based activity as magainin-like peptides and may have a similar mechanism of action.

Impact of Arginine-Phosphate Interactions on the Reentrant Condensation of Disordered Proteins

Re-entrant condensation results in the formation of a condensed protein regime between two critical ion concentrations. The process is driven by neutralization and inversion of the protein charge by oppositely charged ions. Re-entrant condensation of cationic proteins by the polyvalent anions, pyrophosphate and tripolyphosphate, has previously been observed, but not for citrate, which has similar charge and size compared to the polyphosphates. Therefore, besides electrostatic interactions, other specific interactions between the polyphosphate ions and proteins must contribute. Here, we show that additional attractive interactions between arginine and tripolyphosphate determine the re-entrant condensation and decondensation boundaries of the cationic, intrinsically disordered saliva protein, histatin 5. Furthermore, we show by small-angle X-ray scattering (SAXS) that polyvalent anions cause compaction of histatin 5, as would be expected based solely on electrostatic interactions. Hence, we conclude that arginine-phosphate-specific interactions not only regulate solution properties but also influence the conformational ensemble of histatin 5, which is shown to vary with the number of arginine residues. Together, the results presented here provide further insight into an organizational mechanism that can be used to tune protein interactions in solution of both naturally occurring and synthetic proteins.

Bivalirudin and sirolimus co-eluting coronary stent: Potential strategy for the prevention of stent thrombosis and restenosis

Localized drug delivery with sustained elution characteristics from nanocarrier coated stents represents a viable therapeutic approach to circumvent concerns related to coronary stent therapy. We fabricated a Sirolimus (SRL) and Bivalirudin (BIV) releasing nanoparticles (NPs) coated stent for concurrent mitigation of vascular restenosis and acute stent thrombosis. SRL NPs were prepared by nanoprecipitation method whereas the BIV vesicles were generated using hydrophobic ion pair approach followed by micellization phenomenon. MTT assay and confocal microscopic analysis indicated superior anti-proliferative activity and higher cellular uptake of SRL NPs into human coronary artery smooth muscle cells, respectively. DSC and ATR-FTIR techniques confirmed the formation of complex between BIV and phosphatidylglycerol via some weak physical interactions. More than 2 fold rise in log P value was obtained for DSPG-BIV at 3:1 M ratio compared with native BIV solution. The SAXS analysis indicated formation of oligolamellar vesicles of DSPG-BIV complex which was preferentially entrapped into lipophilic lamellae of vesicles. APTT, PT, and TT tests revealed that the BIV vesicles caused significant prolongation of clotting time compared to native BIV solution. The SEM analysis showed uniform and defect free stent coating. In vitro release study demonstrated that SRL and BIV were eluted in a sustained manner from coated stents.

Artificial Host Molecules to Covalently Capture 8-Nitro-cGMP in Neutral Aqueous Solutions and in Cells

New 1,3-diazaphenoxazine derivatives (nitroG-Grasp-Guanidine, NGG) have been developed to covalently capture 8-nitro-cGMP in neutral aqueous solutions, which furnish a thiol reactive group to displace the 8-nitro group and a guanidine unit for interaction with the cyclic phosphate. The thiol group was introduced to the 1,3-diazaphenoxazine skeleton through a 2-aminobenzylthiol group (NGG-H) and its 4-methyl (NGG-pMe) and 6-methyl (NGG-oMe) substituted derivatives. The covalent adducts were formed between the NGG derivatives and 8-nitro-cGMP in neutral aqueous solutions. Among the NGG derivatives, the one with the 6-methyl group (NGG-oMe) exhibited the most efficient capture reaction. Furthermore, NGG-H showed a cell permeability into HEK-293 and RAW 264.7 cells and reduced the intracellular 8-nitro-cGMP level. The NGG derivatives developed in this study would become a valuable tool to study the intracellular role of 8-nitro-cGMP.

A solid-state NMR tool box for the investigation of ATP-fueled protein engines

Motor proteins are involved in a variety of cellular processes. Their main purpose is to convert the chemical energy released during adenosine triphosphate (ATP) hydrolysis into mechanical work. In this review, solid-state Nuclear Magnetic Resonance (NMR) approaches are discussed allowing studies of structures, conformational events and dynamic features of motor proteins during a variety of enzymatic reactions. Solid-state NMR benefits from straightforward sample preparation based on sedimentation of the proteins directly into the Magic-Angle Spinning (MAS) rotor. Protein resonance assignment is the crucial and often time-limiting step in interpreting the wealth of information encoded in the NMR spectra. Herein, potentials, challenges and limitations in resonance assignment for large motor proteins are presented, focussing on both biochemical and spectroscopic approaches. This work highlights NMR tools available to study the action of the motor domain and its coupling to functional processes, as well as to identify protein-nucleotide interactions during events such as DNA replication. Arrested protein states of reaction coordinates such as ATP hydrolysis can be trapped for NMR studies by using stable, non-hydrolysable ATP analogues that mimic the physiological relevant states as accurately as possible. Recent advances in solid-state NMR techniques ranging from Dynamic Nuclear Polarization (DNP), ³¹P-based heteronuclear correlation experiments, ¹H-detected spectra at fast MAS frequencies >100 kHz to paramagnetic NMR are summarized and their applications to the bacterial DnaB helicase from Helicobacter pylori are discussed.

Nucleotide Binding Modes in a Motor Protein Revealed by 31 P- and 1 H-Detected MAS Solid-State NMR Spectroscopy

Protein-nucleic acid interactions play important roles not only in energy-providing reactions, such as ATP hydrolysis, but also in reading, extending, packaging, or repairing genomes. Although they can often be analyzed in detail with X-ray crystallography, complementary methods are needed to visualize them in complexes, which are not crystalline. Here, we show how solid-state NMR spectroscopy can detect and classify protein-nucleic interactions through site-specific 1 H- and 31 P-detected spectroscopic methods. The sensitivity of 1 H chemical-shift values on noncovalent interactions involved in these molecular recognition processes is exploited allowing us to probe directly the chemical bonding state, an information, which is not directly accessible from an X-ray structure. We show that these methods can characterize interactions in easy-to-prepare sediments of the 708 kDa dodecameric DnaB helicase in complex with ADP:AlF4- :DNA, and this despite the very challenging size of the complex.

The biomaterial polyphosphate blocks stoichiometric binding of the SARS-CoV-2 S-protein to the cellular ACE2 receptor

The polymer polyphosphate, abundant in blood platelets, blocks the binding of the receptor-binding domain (RBD) of the SARS- spike (S)-protein to the angiotensin-converting enzyme 2 (ACE2) at low concentrations.

The conformational changes coupling ATP hydrolysis and translocation in a bacterial DnaB helicase

DnaB helicases are motor proteins that couple ATP-hydrolysis to the loading of the protein onto DNA at the replication fork and to translocation along DNA to separate double-stranded DNA into single strands during replication. Using a network of conformational states, arrested by nucleotide mimics, we herein characterize the reaction coordinates for ATP hydrolysis, DNA loading and DNA translocation using solid-state NMR spectroscopy. AMP-PCP is used as pre-hydrolytic, ADP:AlF4- as transition state, and ADP as post-hydrolytic ATP mimic. 31P and 13C NMR spectra reveal conformational and dynamic responses to ATP hydrolysis and the resulting DNA loading and translocation with single amino-acid resolution. This allows us to identify residues guiding the DNA translocation process and to explain the high binding affinities for DNA observed for ADP:AlF4-, which turns out to be optimally preconfigured to bind DNA.

Protein-nucleotide contacts in motor proteins detected by DNP-enhanced solid-state NMR

DNP (dynamic nuclear polarization)-enhanced solid-state NMR is employed to directly detect protein-DNA and protein-ATP interactions and identify the residue type establishing the intermolecular contacts. While conventional solid-state NMR can detect protein-DNA interactions in large oligomeric protein assemblies in favorable cases, it typically suffers from low signal-to-noise ratios. We show here, for the oligomeric DnaB helicase from Helicobacter pylori complexed with ADP and single-stranded DNA, that this limitation can be overcome by using DNP-enhanced spectroscopy. Interactions are established by DNP-enhanced ³¹P-¹³C polarization-transfer experiments followed by the recording of a 2D ¹³C-¹³C correlation experiment. The NMR spectra were obtained in less than 2 days and allowed the identification of residues of the motor protein involved in nucleotide binding.

Arginine-phosphate salt bridges between histones and DNA: intermolecular actuators that control nucleosome architecture

Structural bioinformatics and van der Waals density functional theory are combined to investigate the mechanochemical impact of a major class of histone-DNA interactions, namely, the formation of salt bridges between arginine residues in histones and phosphate groups on the DNA backbone. Principal component analysis reveals that the configurational fluctuations of the sugar-phosphate backbone display sequence-specific directionality and variability, and clustering of nucleosome crystal structures identifies two major salt-bridge configurations: a monodentate form in which the arginine end-group guanidinium only forms one hydrogen bond with the phosphate, and a bidentate form in which it forms two. Density functional theory calculations highlight that the combination of sequence, denticity, and salt-bridge positioning enables the histones to apply a tunable mechanochemical stress to the DNA via precise and specific activation of backbone deformations. The results suggest that selection for specific placements of van der Waals contacts, with high-precision control of the spatial distribution of intermolecular forces, may serve as an underlying evolutionary design principle for the structure and function of nucleosomes, a conjecture that is corroborated by previous experimental studies.

Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos

Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.

The different interactions of lysine and arginine side chains with lipid membranes

The basic amino acids lysine (Lys) and arginine (Arg) play important roles in membrane protein activity, the sensing of membrane voltages, and the actions of antimicrobial, toxin, and cell-penetrating peptides. These roles are thought to stem from the strong interactions and disruptive influences of these amino acids on lipid membranes. In this study, we employ fully atomistic molecular dynamics simulations to observe, quantify, and compare the interactions of Lys and Arg with saturated phosphatidylcholine membranes of different thickness. We make use of both charged (methylammonium and methylguanidinium) and neutral (methylamine and methylguanidine) analogue molecules, as well as Lys and Arg side chains on transmembrane helix models. We find that the free energy barrier experienced by a charged Lys crossing the membrane is strikingly similar to that of a charged Arg (to within 2 kcal/mol), despite the two having different chemistries, H-bonding capability, and hydration free energies that differ by ∼10 kcal/mol. In comparison, the barrier for neutral Arg is higher than that for neutral Lys by around 5 kcal/mol, being more selective than that for the charged species. This can be explained by the different transport mechanisms for charged or neutral amino acid side chains in the membrane, involving membrane deformations or simple dehydration, respectively. As a consequence, we demonstrate that Lys would be deprotonated in the membrane, whereas Arg would maintain its charge. Our simulations also reveal that Arg attracts more phosphate and water in the membrane, and can form extensive H-bonding with its five H-bond donors to stabilize Arg-phosphate clusters. This leads to enhanced interfacial binding and membrane perturbations, including the appearance of a trans-membrane pore in a thinner membrane. These results highlight the special role played by Arg as an amino acid to bind to, disrupt, and permeabilize lipid membranes, as well as to sense voltages for a range of peptide and protein activities in nature and in engineered bionanodevices.

Dynamics of the rhomboid-like protein RHBDD2 expression in mouse retina and involvement of its human ortholog in retinitis pigmentosa

The novel rhomboid-like protein RHBDD2 is distantly related to rhomboid proteins, a group of highly specialized membrane-bound proteases that catalyze regulated intramembrane proteolysis. In retina, RHBDD2 is expressed from embryonic stages to adulthood, and its levels show age-dependent changes. RHBDD2 is distinctly abundant in the perinuclear region of cells, and it localizes to their Golgi. A glycine zipper motif present in one of the transmembrane domains of RHBDD2 is important for its packing into the Golgi membranes. Its deletion causes dislodgment of RHBDD2 from the Golgi. A specific antibody against RHBDD2 recognizes two forms of the protein, one with low (39 kDa; RHBDD2(L)) and the other with high (117 kDa; RHBDD2H) molecular masses in mouse retinal extracts. RHBDD2(L) seems to be ubiquitously expressed in all retinal cells. In contrast, RHBDD2H seems to be present only in the outer segments of cone photoreceptors and may correspond to a homotrimer of RHBDD2(L). This protein consistently co-localizes with S- and M-types of cone opsins. We identified a homozygous mutation in the human RHBDD2 gene, R85H, that co-segregates with disease in affected members of a family with autosomal recessive retinitis pigmentosa. Our findings suggest that the RHBDD2 protein plays important roles in the development and normal function of the retina.

Intracellular fates of cell-penetrating block copolypeptide vesicles

The block copolypeptide poly(l-homoarginine)(60)-b-poly(l-leucine)(20) (R(60)L(20)) was previously found to self-assemble into versatile vesicles with controllable size and encapsulate hydrophilic cargo. These R(60)L(20) vesicles also demonstrated the ability to cross the cell membrane and transport encapsulated cargo into different cell lines. To assess the potential for using the R(60)L(20) vesicles as drug delivery vehicles further, we have investigated their endocytosis and intracellular trafficking behavior. Using drugs that inhibit different endocytosis pathways, we identified macropinocytosis to be a major process by which the R(60)L(20) vesicles enter HeLa cells. Subsequent immunostaining experiments demonstrated that the vesicles entered the early endosomes but not the lysosomes, suggesting that they recycle back to the cell surface. Overall, our studies indicate that the R(60)L(20) vesicles are able to enter cells intact with their cargos, and although some manage to escape from early endosomes, most are trapped within these intracellular compartments.

Car−Parrinello Molecular Dynamics Study of Anharmonic Systems: A Mannich Base in Solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schrödinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.

Strengths of hydrogen bonds involving phosphorylated amino acid side chains

Post-translational phosphorylation plays a key role in regulating protein function. Here, we provide a quantitative assessment of the relative strengths of hydrogen bonds involving phosphorylated amino acid side chains (pSer, pAsp) with several common donors (Arg, Lys, and backbone amide groups). We utilize multiple levels of theory, consisting of explicit solvent molecular dynamics, implicit solvent molecular mechanics, and quantum mechanics with a self-consistent reaction field treatment of solvent. Because the approximately 6 pKa of phosphate suggests that -1 and -2 charged species may coexist at physiological pH, hydrogen bonds involving both protonated and deprotonated phosphates for all donor-acceptor pairs are considered. Multiple bonding geometries for the charged-charged interactions are also considered. Arg is shown to be capable of substantially stronger salt bridges with phosphorylated side chains than Lys. A pSer hydrogen-bond acceptor tends to form more stable interactions than a pAsp acceptor. The effect of phosphate protonation state on the strengths of the hydrogen bonds is remarkably subtle, with a more pronounced effect on pAsp than on pSer.

Voltage gating and anions, especially phosphate: a model system

The voltage sensor of voltage gated sodium and potassium channels consists of four sets of transmembrane segments, of which one, called S4, contains at least four arginines; these are presumed to each carry positive charges. The channel opening is usually attributed to the outward (i.e., toward the extracellular side of the membrane) motion of S4. The evidence for this motion is based on certain experiments that appear to show differential access to parts of S4 from the intracellular and extracellular sides of the membrane in the open and closed states. A newly available structure [S.B. Long, E.B. Campbell and R. MacKinnon, Crystal structure of a mammalian voltage-dependent Shaker family K(+) channel. Science 309 (2005) 897-903; S.B. Long, E.B. Campbell, R. MacKinnon, Voltage sensor of Kv1.2: structural basis of electromechanical coupling. Science 309 (2005) 903-908][1,2] has now been used to argue for a large scale motion, although, as a static structure, it is not conclusive. In this paper, we consider the effect of anions in the surrounding medium. Phosphate is present in the intracellular as well as the extracellular fluid, apparently at hundreds of micromolar concentration, or more. There is evidence in the literature suggesting that phosphate-arginine complexes are rather strong. In a recent calculation one of us [M.E. Green, A possible role for phosphate in complexing the arginines of S4 in voltage gated channels J. Theor. Biol. 233 (2005) 337-341][3] has shown that a model peptide with a 2:1 arg:phosphate complex should have a favorable geometry. Here, we present NMR evidence of the existence of phosphate complexes of a model peptide with two arginines separated by two hydrophobic residues, the same spacing as in S4 segments. The complexes (there are different complexes for HPO(4)(2-) and for H(2)PO(4)(-) [3]) form with concentrations of peptide in the range of hundreds of micromolar, making it significant in the biological context. NMR spectra provide changes in chemical shift as functions of both phosphate concentration and pH. The resulting curves show titration of the phosphate, with its standard pK. Possible implications for other anion-S4 interactions, including ion pairs rather than complexes, as with Cl(-), are also discussed.

A possible role for phosphate in complexing the arginines of S4 in voltage gated channels

Phosphate ions are known to complex guanidinium groups, which are the side chains of arginine. Voltage gated channels that allow passage of ions through cell membranes, producing, for example the nerve impulse, are in many cases composed of four domains, each with six transmembrane segments. The S4 transmembrane segments of these channels have arginines placed in such a way that they would be expected to complex phosphate. Known phosphate-arginine complexes are reasonably strong. Here, an ab initio calculation reinforces the expectation that a strong complex could form. As a consequence, if the S4 moved, it would carry either no charge, or at most half of what is expected from fully charged arginines. This suggests that it may be necessary to rethink voltage gating models in which the gating current is produced by physical motion of the S4 transmembrane segments.

The transmembrane domain of the acetylcholine receptor: insights from simulations on synthetic peptide models

We have studied the structure and properties of a bundle of alpha-helical peptides embedded in a 1,2-dimyristoyl-3-phosphatidylcholine phospholipid bilayer by molecular dynamics simulations. The bundle of five transmembrane deltaM2 segments constitutes the model for the pore region of the nicotinic acetylcholine receptor, which is the neurotransmitter-gated ion-channel responsible for the fast propagation of electrical signals between cells at the nerve-muscle synapse. The deltaM2 segments were shown to oligomerize in biomembranes resulting in ion-channel activity with characteristics similar to the native protein, and the structure of the isolated peptides was studied in 1,2-dimyristoyl-3-phosphatidylcholine bilayers and micelles by NMR experiments (Opella, S. J., et al. 1999. Nat. Struct. Biol. 6:374-379). Our analyses indicate that the structure, helix tilt, and the overall shape of the channel are in good agreement with the NMR experiments and the proposed model for the channel, which we show is formed by rings of functional residues. The studied geometry resulted in a closed pore state, where the channel is partially dehydrated at the hydrophobic extracellular half and the extracellular mouth of the channel blocked by the hydrocarbon chains of Arg+ residues. The arginine amino acids form intermolecular salt-bridges with the C-terminus, which contribute as well to the bundle stabilization.