Calculating the electrostatic properties of RNA provides new insights into molecular interactions and function.

Solutions to the nonlinear Poisson-Boltzmann equation were used to obtain the electrostatic potentials of RNA molecules that have known three-dimensional structures. The results are described in terms of isopotential contours and surface electrostatic potential maps. Both representations have unexpected features: 'cavities' within isopotential contours and areas of enhanced negative potential on molecular surfaces. Intriguingly, the sites of unusual electrostatic features correspond to functionally important regions, suggesting that electrostatic properties play a key role in RNA recognition and stabilization. These calculations reveal that the electrostatic potentials generated by RNA molecules have a variety of functionally important characteristics that cannot be discerned by simple visual inspection of the molecular structure.     

Protein Surface Dynamics: Interaction with Water and Small Solutes

Previous time resolved measurements had indicated that protons could propagate on the surface of a protein, or a membrane, by a special mechanism that enhances the shuttle of the proton towards a specific site [1]. It was proposed that a proper location of residues on the surface contributes to the proton shuttling function. In the present study, this notion was further investigated using molecular dynamics, with only the mobile charge replaced by Na⁺ and Cl⁻ ions. A molecular dynamics simulation of a small globular protein (the S6 of the bacterial ribosome) was carried out in the presence of explicit water molecules and four pairs of Na⁺ and Cl⁻ ions. A 10 ns simulation indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. Yet it was noted that, close to some domains, the ions extended their duration near the surface, suggesting that the local electrostatic potential prevented them from diffusing to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.     

Profiling charge complementarity and selectivity for binding at the protein surface

A novel analysis and representation of the protein surface in terms of electrostatic binding complementarity and selectivity is presented. The charge optimization methodology is applied in a probe-based approach that simulates the binding process to the target protein. The molecular surface is color coded according to calculated optimal charge or according to charge selectivity, i.e., the binding cost of deviating from the optimal charge. The optimal charge profile depends on both the protein shape and charge distribution whereas the charge selectivity profile depends only on protein shape. High selectivity is concentrated in well-shaped concave pockets, whereas solvent-exposed convex regions are not charge selective. This suggests the synergy of charge and shape selectivity hot spots toward molecular selection and recognition, as well as the asymmetry of charge selectivity at the binding interface of biomolecular systems. The charge complementarity and selectivity profiles map relevant electrostatic properties in a readily interpretable way and encode information that is quite different from that visualized in the standard electrostatic potential map of unbound proteins.     

Electrostatic orientation of the electron-transfer complex between plastocyanin and cytochrome c.

To understand the specificity and efficiency of protein-protein interactions promoting electron transfer, we evaluated the role of electrostatic forces in precollision orientation by the development of two new methods, computer graphics alignment of protein electrostatic fields and a systematic orientational search of intermolecular electrostatic energies for two proteins at present separation distances. We applied these methods to the plastocyanin/cytochrome c interaction, which is faster than random collision, but too slow for study by molecular dynamics techniques. Significant electrostatic potentials were concentrated on one-fourth (969 A2) of the plastocyanin surface, with the greatest negative potential centered on the Tyr-83 hydroxyl within the acidic patch, and on one-eighth (632 A2) of the cytochrome c surface, with the greatest positive potential centered near the exposed heme edge. Coherent electrostatic fields occurred only over these regions, suggesting that local, rather than global, charge complementarity controls productive recognition. The three energetically favored families of pre-collision orientations all directed the positive region surrounding the heme edge of cytochrome c toward the acidic patch of plastocyanin but differed in heme plane orientation. Analysis of electrostatic fields, electrostatic energies of precollision orientations with 12 and 6 A separation distances, and surface topographies suggested that the favored orientations should converge to productive complexes promoting a single electron-transfer pathway from the cytochrome c heme edge to Tyr-83 of plastocyanin. Direct interactions of the exposed Cu ligand in plastocyanin with the cytochrome c heme edge are not unfavorable sterically or electrostatically but should occur no faster than randomly, indicating that this is not the primary pathway for electron transfer.     

The Molecular Determinants of the Increased Reduction Potential of the Rubredoxin Domain of Rubrerythrin Relative to Rubredoxin

Based on the crystal structures, three possible sequence determinants have been suggested as the cause of a 285 mV increase in reduction potential of the rubredoxin domain of rubrerythrin over rubredoxin by modulating the polar environment around the redox site. Here, electrostatic calculations of crystal structures of rubredoxin and rubrerythrin and molecular dynamics simulations of rubredoxin wild-type and mutants are used to elucidate the contributions to the increased reduction potential. Asn¹⁶⁰ and His¹⁷⁹ in rubrerythrin versus valines in rubredoxins are predicted to be the major contributors, as the polar side chains contribute significantly to the electrostatic potential in the redox site region. The mutant simulations show both side chains rotating on a nanosecond timescale between two conformations with different electrostatic contributions. Reduction also causes a change in the reduction energy that is consistent with a linear response due to the interesting mechanism of shifting the relative populations of the two conformations. In addition to this, a simulation of a triple mutant indicates the side-chain rotations are approximately anticorrelated so whereas one is in the high potential conformation, the other is in the low potential conformation. However, Ala¹⁷⁶ in rubrerythrin versus a leucine in rubredoxin is not predicted to be a large contributor, because the solvent accessibility increases only slightly in mutant simulations and because it is buried in the interface of the rubrerythrin homodimer.     

Molecular dynamics simulation of hepatitis C virus IRES IIId domain: structural behavior,electrostatic and energetic analysis

The dynamic behavior of the HCV IRES IIId domain is analyzed by means of a 2.6-ns molecular dynamics simulation, starting from an NMR structure. The simulation is carried out in explicit water with Na⁺ counterions, and particle-mesh Ewald summation is used for the electrostatic interactions. In this work, we analyze selected patterns of the helix that are crucial for IRES activity and that could be considered as targets for the intervention of inhibitors, such as the hexanucleotide terminal loop (more particularly its three consecutive guanines) and the loop-E motif. The simulation has allowed us to analyze the dynamics of the loop substructure and has revealed a behavior among the guanine bases that might explain the different role of the third guanine of the GGG triplet upon molecular recognition. The accessibility of the loop-E motif and the loop major and minor groove is also examined, as well as the effect of Na⁺ or Mg²⁺ counterion within the simulation. The electrostatic analysis reveals several ion pockets, not discussed in the experimental structure. The positions of these ions are useful for locating specific electrostatic recognition sites for potential inhibitor binding. Figure Superposition of 14 structures representative of the evolution of IRES IIId RNA along 2.6-ns MD simulation     

Effect of the electrostatic potential on the internalization mechanism of cell penetrating peptides derived from TIRAP

In order to develop future therapeutic applications for cell penetrating peptides (CPPs), it is essential to characterize their internalization mechanisms, as they might affect the stability and the accessibility of the carried drug. Several internalization mechanisms have been described in literature, such as endocytosis and transduction. In this work we study the internalization mechanism in HeLa cells of two TIRAP derived peptides: pepTIRAP and pepTIRAP^ALA, where some of the cationic amino acids were replaced with alanines. Detailed analysis of internalization and the peptides electrostatic potential was carried out, to shed light on the internalization mechanism involved. Molecular modeling studies showed that the main difference identified between pepTIRAP and pepTIRAP^ALA is the distribution of their electrostatic potential field. The structure of pepTIRAP displays a predominantly positive potential when compared to pepTIRAP^ALA, which has a more balanced potential distribution. In addition, docking experiments show that interactions between pepTIRAP and negatively charged molecules on the cellular surface such as heparan sulfate are stronger than the ones exhibited by pepTIRAP^ALA. A mathematical model was proposed to quantify the amount of peptide internalized or non-specifically bound to the membrane. The model indicates a stronger interaction of pepTIRAP with the plasma membrane, compared to pepTIRAP^ALA. We propose these discrepancies are related to the differences in the electrostatic potential characteristics of each peptide. In the case of pepTIRAP, these interactions lead to the formation of nucleation zones, which are the first stage of the transduction internalization mechanism. These results should be considered for effective design of a cell penetrating peptide.     

Electrostatic components of drug-receptor recognition. II. The DNA-binding antibiotic actinomycin

The semiempirical molecular orbital method PCILO has been used to study the electronic charge distribution and conformation of the antibiotic actinomycin. Molecular electrostatic potential fields have been constructed in three-dimensional space round the molecule and displayed stereoscopically together with contour maps in the orthogonal planes relative to that of the chromophore. It was found that the drug has a considerable dipole moment and that the electrostatic fields can be separated into two large regions; an area of positive potential surrounds the chromophore and one of negative potential is directed away from the peptide rings. In this and the preceding paper these observations are discussed with respect to drug-receptor recognition phenomena. Pattern-matching of complementary electrostatic fields between the drug and polynucleotide receptor can be discerned.     

Molecular dynamics studies on the NMR structures of rabbit prion protein wild type and mutants: surface electrostatic charge distributions

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species such as sheep and goats, cattle, deer and elk, and humans. But for rabbits, studies have shown that they have a low susceptibility to be infected by prion diseases. This paper does molecular dynamics (MD) studies of rabbit NMR structures (of the wild type and its two mutants of two surface residues), in order to understand the specific mechanism of rabbit prion proteins (RaPrP^C). Protein surface electrostatic charge distributions are specially focused to analyze the MD trajectories. This paper can conclude that surface electrostatic charge distributions indeed contribute to the structural stability of wild-type RaPrP^C; this may be useful for the medicinal treatment of prion diseases.     

Revealing substitution effects on the strength and nature of halogen-hydride interactions: a theoretical study

A quantum chemistry study was carried out to investigate the strength and nature of halogen bond interactions in HXeH···XCCY complexes, where X = Cl, Br and Y = H, F, Cl, Br, CN, NC, C₂H, CH₃, OH, SH, NH₂. Examination of the electrostatic potentials V(r) of the XCCY molecules reveals that the addition of substituents has a significant effect upon the most positive electrostatic potential on the surface of the interacting halogen atom. We found that the magnitude of atomic charges and multipole moments depends upon the halogen atom X and is rather sensitive to the electron-withdrawing/donating power of the remainder of the molecule. An excellent correlation was found between the most positive electrostatic potentials on the halogen atom and the interaction energies. For either HXeH···ClCCY or HXeH···BrCCY complexes, an approximate linear correlation between the interaction energies and halogens multipole moments are established, indicating that the electrostatic and polarization interactions are responsible for the stability of the complexes. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the HXeH···XCCY complexes. Furthermore, the changes in the electrostatic term are mainly responsible for the dependence of interaction energy on the halogen atom.

Crystal structure of Streptococcus pneumoniae Sp1610, a putative tRNA methyltransferase,in complex with S‐adenosyl‐L‐methionine

Streptococcus pneumoniae Sp1610, a Class‐I fold S‐adenosylmethionine (AdoMet)‐dependent methyltransferase, is a member of the COG2384 family in the Clusters of Orthologous Groups database, which catalyzes the methylation of N¹‐adenosine at position 22 of bacterial tRNA. We determined the crystal structure of Sp1610 in the ligand‐free and the AdoMet‐bound forms at resolutions of 2.0 and 3.0 Å, respectively. The protein is organized into two structural domains: the N‐terminal catalytic domain with a Class I AdoMet‐dependent methyltransferase fold, and the C‐terminal substrate recognition domain with a novel fold of four α‐helices. Observations of the electrostatic potential surface revealed that the concave surface located near the AdoMet binding pocket was predominantly positively charged, and thus this was predicted to be an RNA binding area. Based on the results of sequence alignment and structural analysis, the putative catalytic residues responsible for substrate recognition are also proposed.     

Theoretical Study of Antagonists and Inhibitors of Mammalian Adenosine Deaminase: II. Isomeric Aza-deazaanalogues of Adenosine

Isomeric aza-deazaanalogues of adenosine and their N1-protonated forms (except for that of 8-aza-1-deazaadenosine) were studied by computer modeling to find a relationship between their molecular structures and the properties as substrates for the mammalian adenosine deaminase. The atomic charge distribution and maps of electrostatic potential around their van der Waals molecular surface were calculated using the ab initioSTO-3G method. The conformational studies were carried out by the MM⁺ method of molecular mechanics. The previously proposed mechanism of the substrate acceptance in the active site of mammalian adenosine deaminase was refined, and the potential substrate properties were predicted for two previously unstudied adenosine analogues, 5-aza-9-deazaadenosine and 8-aza-3-deazaadenosine.     

Computational and Experimental Investigation of the Antimicrobial Peptide Cecropin XJ and its Ligands as the Impact Factors of Antibacterial Activity

Cecropin XJ, as a heat stable antimicrobial peptide (AMP), displayed broad bacteriostatic activities, effectively inhibited proliferation of cancer cells and induced cell apoptosis in vitro. However, it exhibited little hemolytic activity and very low cytotoxicity to erythrocytes and normal cells. Although exerts multiple remarkable bioactivities, the refined molecular conformation of native Cecropin XJ remains unsolved. The aim of the present study is to comprehensively investigate the physicochemical characteristics and structure-function relationship of this antimicrobial peptide by using a series of bioinformatics and experimental approaches. In this study, we revealed that the mature Cecropin XJ consists of 41 amino acids, containing two α-helical structures from Lys⁷ to Lys²⁵ and from Ala²⁹ to Ile³⁹. The phylogenetic tree indicated that Cecropin XJ belongs to the Class I AMPs of cecropin family. Hydrophobic analysis showed Cecropin XJ is a typical amphiphilic molecule. The surface of Cecropin XJ was found to have a much wide range of electrostatic potential from ?83.243 to +83.243. The amphipathicity and surface potential of Cecropin XJ partially supported the AMP pore-forming hypothesis. Scanning electron microscopy experimentally confirmed the damages of Cecropin XJ to microbial membrane. Four predicted docking sites respectively for magnesium ion (Mg²⁺), adenosine diphosphate (ADP), bacteriopheophytin (BPH), and guanosine triphosphate (GTP) were found on the surface of Cecropin XJ. Thereinto, Mg²⁺ was experimentally proved to suppress the antibacterial activity of Cecropin XJ; both GTP and ADP enhanced the bactericidal activities to varying degrees. The present study provides a foundation for further investigation of molecular evolution, structural modification, and functional mechanisms of Cecropin XJ.     

Theoretical Study of Antagonists and Inhibitors of Mammalian Adenosine Deaminase: I. Adenosine and Its Aza- and Deazaanalogues

Aza- and deazaanalogues of adenosine, including their 1-protonated forms (except for that of 1-deazaadenosine), were studied by computer computation to find a relationship between their molecular structures and substrate properties for the mammalian adenosine deaminase. The atomic charge distribution and maps of the electrostatic potential around their van der Waals molecular surface were calculated for these compounds using the ab initio STO-3G method. The conformational studies are carried out by the MM⁺ method of molecular mechanics. The mechanism that determines the substrate selectivity of mammalian adenosine deaminase is discussed.     

Substrate recognition by ribosome-inactivating protein studied by molecular modeling and molecular electrostatic potentials

A computer model of dianthin 30, a type 1 ribosome-inactivating protein (RIP), is constructed by homology modeling using two known X-ray structures; a type 1 RIP, pokeweed antiviral protein (PAP), and chain A of a type 2 RIP, ricin. The 3D structure is refined by molecular dynamics and its binding site compared with those of PAP and ricin using molecular electrostatic potential mapping. The differences in the maps obtained clearly show how, despite the similarity of the topology of the binding site, differences in electrostatic potential can account for the experimentally observed differences in substrate recognition and binding. This demonstrates the potential of these techniques for guiding further experimental analyses.     

Modelling the three-dimensional structure and electrostatic potential field of the two Cu,Zn superoxide dismutase variants from Xenopus laevis

The crystallographic structure of bovine superoxide dismutase has been used as a template for the graphic reconstruction of the three-dimensional structures of the two Xenopus laevis variants (Schininà, M.E. et al. Arch. Biochem. Biophys. 272:507-515, 1989). In these models the structure-essential residues maintain their position and their structural role, and the interactions between the subunits and the close packing within the beta-barrel are maintained with conservative substitutions and even increased with "aromatic pairs." Because of the same topological motif and surface location of charges, arising from the model building of the two variants with respect to the bovine enzyme, we have calculated the electrostatic potential fields around the models of the two Xenopus laevis variants by numerically solving the Poisson-Boltzmann equation. We show that conservation of a specific space-relationship of charges maintains the potential field pattern already observed in the bovine enzyme, where a negative potential field surrounds the protein surface and specific positive regions wrap up the copper center active site. This electrostatic potential field distribution supports the idea that electrostatic interactions control, like in the bovine enzyme, the mechanism of enzyme-substrate recognition in the Xenopus laevis Cu,Zn superoxide dismutases, suggesting that coordinated mutation of charged residues has occurred in the evolution of this enzyme.     

Molecular similarity analysis between insect juvenile hormone and N, N-diethyl-m-toluamide (DEET) analogs may aid design of novel insect repellents

Molecular similarity analysis of stereoelectronic properties between natural insect juvenile hormone (JH), -a synthetic insect juvenile hormone mimic (JH-mimic, undecen-2-yl carbamate), and N, N-diethyl-m-toluamide (DEET) and its analogs reveals similarities that may aid the design of more efficacious insect repellents and give a better insight into the mechanism of repellent action. The study involves quantum chemical calculations using the AM1 semi-empirical computational method enabling a conformational search for the lowest and most abundant energy conformers of JH, JH-mimic, and 15 DEET compounds, followed by complete geometry optimization of the conformers. Similarity analyses of stereoelectronic properties such as structural parameters, atomic charges, dipole moments, molecular electrostatic potentials, and highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were performed on JH, JH-mimic and the DEET compounds. The similarity of stereoelectronic attributes of the amide/ester moiety, the negative electrostatic potential regions beyond the van der Waals surface, and the large distribution of hydrophobic regions in the compounds appear to be the three important factors leading to a similar interaction with the JH receptor. The similarity of electrostatic profiles beyond the van der Waals surface is likely to play a crucial role in molecular recognition interaction with the JH receptor from a distance. This also suggests electrostatic bioisosterism of the amide group of the DEET compounds and JH-mimic and, thus, a model for molecular recognition at the JH receptor. The insect repellent property of the DEET analogs may thus be attributed to a conflict of complementarity for the JH receptor binding sites.     

Hydration effects on the electrostatic potential around tuftsin

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr–Lys–Pro–Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson–Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed. © 1999 John Wiley & Sons, Inc. Biopoly 50: 133–143, 1999     

Interrogating the mechanism of a tight binding inhibitor of AIR carboxylase

The enzyme aminoimidazole ribonucleotide (AIR) carboxylase catalyzes the synthesis of the purine intermediate, 4-carboxy-5-aminoimidazole ribonucleotide (CAIR). Previously, we have shown that the compound 4-nitro-5-aminoimidazole ribonucleotide (NAIR) is a slow, tight binding inhibitor of the enzyme with a Ki of 0.34 nM. The structural attributes and the slow, tight binding characteristics of NAIR implicated this compound as a transition state or reactive intermediate analog. However, it is unclear what molecular features of NAIR contribute to the mimetic properties for either of the two proposed mechanisms of AIR carboxylase. In order to gain additional information regarding the mechanism for the potent inhibition of AIR carboxylase by NAIR, a series of heterocyclic analogs were prepared and evaluated. We find that all compounds are weaker inhibitors than NAIR and that CAIR analogs are not alternative substrates for the enzyme. Surprisingly, rather subtle changes in the structure of NAIR can lead to profound changes in binding affinity. Computational investigations of enzyme intermediates and these inhibitors reveal that NAIR displays an electrostatic potential surface similar to a proposed reaction intermediate. The result indicates that AIR carboxylase is likely sensitive to the electrostatic surface of reaction intermediates and thus compounds which mimic these surfaces should possess tight binding characteristics. Given the evolutionary relationship between AIR carboxylase and N⁵-CAIR mutase, we believe that this concept extends to the mutase enzyme as well. The implications of this hypothesis for the design of selective inhibitors of the N⁵-CAIR mutase are discussed.