The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies

Direct measurement of the dielectric polarization properties of DNA

The electric polarizability of DNA, represented by the dielectric constant, is a key intrinsic property that modulates DNA interaction with effector proteins. Surprisingly, it has so far remained unknown owing to the lack of experimental tools able to access it. Here, we experimentally resolved it by detecting the ultraweak polarization forces of DNA inside single T7 bacteriophages particles using electrostatic force microscopy. In contrast to the common assumption of low-polarizable behavior like proteins (εr ∼ 2-4), we found that the DNA dielectric constant is ∼ 8, considerably higher than the value of ∼ 3 found for capsid proteins. State-of-the-art molecular dynamic simulations confirm the experimental findings, which result in sensibly decreased DNA interaction free energy than normally predicted by Poisson-Boltzmann methods. Our findings reveal a property at the basis of DNA structure and functions that is needed for realistic theoretical descriptions, and illustrate the synergetic power of scanning probe microscopy and theoretical computation techniques.

Large scale characterization of the LC13 TCR and HLA-B8 structural landscape in reaction to 172 altered peptide ligands: a molecular dynamics simulation study

The interplay between T cell receptors (TCRs) and peptides bound by major histocompatibility complexes (MHCs) is one of the most important interactions in the adaptive immune system. Several previous studies have computationally investigated their structural dynamics. On the basis of these simulations several structural and dynamical properties have been proposed as effectors of the immunogenicity. Here we present the results of a large scale Molecular Dynamics simulation study consisting of 100 ns simulations of 172 different complexes. These complexes consisted of all possible point mutations of the Epstein Barr Virus peptide FLRGRAYGL bound by HLA-B*08:01 and presented to the LC13 TCR. We compare the results of these 172 structural simulations with experimental immunogenicity data. We found that simulations with more immunogenic peptides and those with less immunogenic peptides are in fact highly similar and on average only minor differences in the hydrogen binding footprints, interface distances, and the relative orientation between the TCR chains are present. Thus our large scale data analysis shows that many previously suggested dynamical and structural properties of the TCR/peptide/MHC interface are unlikely to be conserved causal factors for peptide immunogenicity.

Structural basis of protein oxidation resistance: a lysozyme study

Accumulation of oxidative damage in proteins correlates with aging since it can cause irreversible and progressive degeneration of almost all cellular functions. Apparently, native protein structures have evolved intrinsic resistance to oxidation since perfectly folded proteins are, by large most robust. Here we explore the structural basis of protein resistance to radiation-induced oxidation using chicken egg white lysozyme in the native and misfolded form. We study the differential resistance to oxidative damage of six different parts of native and misfolded lysozyme by a targeted tandem/mass spectrometry approach of its tryptic fragments. The decay of the amount of each lysozyme fragment with increasing radiation dose is found to be a two steps process, characterized by a double exponential evolution of their amounts: the first one can be largely attributed to oxidation of specific amino acids, while the second one corresponds to further degradation of the protein. By correlating these results to the structural parameters computed from molecular dynamics (MD) simulations, we find the protein parts with increased root-mean-square deviation (RMSD) to be more susceptible to modifications. In addition, involvement of amino acid side-chains in hydrogen bonds has a protective effect against oxidation Increased exposure to solvent of individual amino acid side chains correlates with high susceptibility to oxidative and other modifications like side chain fragmentation. Generally, while none of the structural parameters alone can account for the fate of peptides during radiation, together they provide an insight into the relationship between protein structure and susceptibility to oxidation.

Base flip in DNA studied by molecular dynamics simulationsof differently-oxidized forms of methyl-Cytosine

Distortions in the DNA sequence, such as damage or mispairs, are specifically recognized and processed by DNA repair enzymes. Many repair proteins and, in particular, glycosylases flip the target base out of the DNA helix into the enzyme’s active site. Our molecular dynamics simulations of DNA with intact and damaged (oxidized) methyl-cytosine show that the probability of being flipped is similar for damaged and intact methyl-cytosine. However, the accessibility of the different 5-methyl groups allows direct discrimination of the oxidized forms. Hydrogen-bonded patterns that vary between methyl-cytosine forms carrying a carbonyl oxygen atom are likely to be detected by the repair enzymes and may thus help target site recognition.

g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations

Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages is described. g_mmpbsa was developed as part of the Open Source Drug Discovery (OSDD) consortium. Its aim is to integrate high-throughput molecular dynamics (MD) simulations with binding energy calculations. The tool provides options to select alternative atomic radii and different nonpolar solvation models including models based on the solvent accessible surface area (SASA), solvent accessible volume (SAV), and a model which contains both repulsive (SASA-SAV) and attractive components (described using a Weeks-Chandler-Andersen like integral method). We showcase the effectiveness of the tool by comparing the calculated interaction energy of 37 structurally diverse HIV-1 protease inhibitor complexes with their experimental binding free energies. The effect of varying several combinations of input parameters such as atomic radii, dielectric constant, grid resolution, solute-solvent dielectric boundary definition, and nonpolar models was investigated. g_mmpbsa can also be used to estimate the energy contribution per residue to the binding energy. It has been used to identify those residues in HIV-1 protease that are most critical for binding a range of inhibitors.

Visualising intrinsic disorder and conformational variation in protein ensembles

Intrinsically disordered regions (IDRs) in proteins are still not well understood, but are increasingly recognised as important in key biological functions, as well as in diseases. IDRs often confound experimental structure determination-however, they are present in many of the available 3D structures, where they exhibit a wide range of conformations, from ill-defined and highly flexible to well-defined upon binding to partner molecules, or upon post-translational modifications. Analysing such large conformational variations across ensembles of 3D structures can be complex and difficult; our goal in this paper is to improve this situation by augmenting traditional approaches (molecular graphics and principal components) with methods from human-computer interaction and information visualisation, especially parallel coordinates. We present a new tool integrating these approaches, and demonstrate how it can dissect ensembles to reveal functional insights into conformational variation and intrinsic disorder.

Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins

The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.

Coiled-coil response to mechanical force: global stability and local cracking

Coiled coils are important structural motifs formed by two or more amphipathic α-helices that twist into a supercoil. These motifs are found in a wide range of proteins, including motor proteins and structural proteins, that are known to transmit mechanical loads. We analyze atomically detailed simulations of coiled-coil cracking under load with Milestoning. Milestoning is an approach that captures the main features of the process in a network, quantifying kinetics and thermodynamics. A 112-residue segment of the β-myosin S2 domain was subjected to constant-magnitude (0-200 pN) and constant-direction tensile forces in molecular dynamics simulations. Twenty 20 ns straightforward simulations at several load levels revealed that initial single-residue cracking events (Ψ > 90°) at loads <100 pN were accompanied by rapid refolding without either intra- or interhelix unfolding propagation. Only initial unfolding events at the highest load (200 pN) regularly propagated along and between helices. Analysis of hydrophobic interactions and of interhelix hydrogen bonds did not show significant variation as a function of load. Unfolding events were overwhelmingly located in the vicinity of E929, a charged residue in a hydrophobic position of the heptad repeat. Milestoning network analysis of E929 cracking determined that the mean first-passage time ranges from 20 ns (200 pN) to 80 ns (50 pN), which is ∼20 times the mean first-passage time of an isolated helix with the same sequence.

Scaling behaviour for the water transport in nanoconfined geometries

The transport of water in nanoconfined geometries is different from bulk phase and has tremendous implications in nanotechnology and biotechnology. Here molecular dynamics is used to compute the self-diffusion coefficient D of water within nanopores, around nanoparticles, carbon nanotubes and proteins. For almost 60 different cases, D is found to scale linearly with the sole parameter θ as D(θ)=D_B[1+(D_C/D_B−1)θ], with D_B and D_C the bulk and totally confined diffusion of water, respectively. The parameter θ is primarily influenced by geometry and represents the ratio between the confined and total water volumes. The D(θ) relationship is interpreted within the thermodynamics of supercooled water. As an example, such relationship is shown to accurately predict the relaxometric response of contrast agents for magnetic resonance imaging. The D(θ) relationship can help in interpreting the transport of water molecules under nanoconfined conditions and tailoring nanostructures with precise modulation of water mobility.

A precise control of water transport in different configurations is crucial for many technological applications. Chiavazzo et al. define a dimensionless scaling parameter to fully describe the water self-diffusion coefficient in various nanoconfined geometries using molecular dynamic simulations.

Reduced model captures Mg(2+)-RNA interaction free energy of riboswitches

The stability of RNA tertiary structures depends heavily on Mg(2+). The Mg(2+)-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ2+, the number of excess Mg(2+) associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg(2+) ions interact closely and strongly with the RNA, unlike monovalent ions such as K(+), suggesting that an explicit treatment of Mg(2+) is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg(2+)-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg(2+) through the explicit representation of Mg(2+) ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K(+) and models its competition with condensed Mg(2+). The model contains one fitted parameter, the number of condensed K(+) ions in the absence of Mg(2+). Values of Γ2+ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ2+ to the Mg(2+)-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg(2+) and RNA dynamics.

Atomistic mechanisms of huntingtin N-terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q(N)) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N-terminal 17-amino-acid sequence (htt(NT)) positioned just before the Q(N) region is important for the binding of huntingtin to membranes. Through all-atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt(NT)Q(N) fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps-approach, reorganization, anchoring, and insertion-that are very diverse at the atomic level. On the membrane, the htt(NT) peptide forms a stable α-helix essentially parallel to the membrane with its nonpolar side-chains-mainly Leu-4, Leu-7, Phe-11 and Leu-14-positioned in the hydrophobic core of the membrane. Salt-bridges involving Glu-5, Glu-12, Lys-6, and Lys-15, as well as hydrogen bonds involving Thr-3 and Ser-13 with the phospholipids also stabilize the structure and orientation of the htt(NT) peptide. These observations do not significantly change upon adding the Q(N) region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt(NT) region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q(N) region could also play a key role for the oligomerization of htt(NT)Q(N) on phospholipid membranes.

Molecular dynamics simulations reveal the assembly mechanism of polysaccharides in marine aerosols

This figure represents the configuration of the marine gel after the polysaccharides undergo self-assembly. The calcium cations are blended inside the assembly and they act as ionic bridges.

Structural features of the regulatory ACT domain of phenylalanine hydroxylase

Phenylalanine hydroxylase (PAH) catalyzes the conversion of L-Phe to L-Tyr. Defects in PAH activity, caused by mutations in the human gene, result in the autosomal recessively inherited disease hyperphenylalaninemia. PAH activity is regulated by multiple factors, including phosphorylation and ligand binding. In particular, PAH displays positive cooperativity for L-Phe, which is proposed to bind the enzyme on an allosteric site in the N-terminal regulatory domain (RD), also classified as an ACT domain. This domain is found in several proteins and is able to bind amino acids. We used molecular dynamics simulations to obtain dynamical and structural insights into the isolated RD of PAH. Here we show that the principal motions involve conformational changes leading from an initial open to a final closed domain structure. The global intrinsic motions of the RD are correlated with exposure to solvent of a hydrophobic surface, which corresponds to the ligand binding-site of the ACT domain. Our results strongly suggest a relationship between the Phe-binding function and the overall dynamic behaviour of the enzyme. This relationship may be affected by structure-disturbing mutations. To elucidate the functional implications of the mutations, we investigated the structural effects on the dynamics of the human RD PAH induced by six missense hyperphenylalaninemia-causing mutations, namely p.G46S, p.F39C, p.F39L, p.I65S, p.I65T and p.I65V. These studies showed that the alterations in RD hydrophobic interactions induced by missense mutations could affect the functionality of the whole enzyme.

Membrane driven spatial organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β₂AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β₂AR compared to the highly homologous and structurally similar β₁AR, consonant with experimentally observed oligomerization patterns of β₂AR and β₁AR. The results provide a mechanistic understanding of the structural context of oligomerization.

Molecular dynamics simulation of SDS and CTAB micellization and prediction of partition equilibria with COSMOmic

Molecular dynamics (MD) simulations of the self-assembly of different ionic surfactants have been performed in order to obtain representative micellar structures. Subsequently, these structures were used to predict the partition behavior of various solutes in these micelles with COSMOmic, an extension of COSMO-RS. This paper includes multiple self-assembled micelles of SDS (sodium dodecyl sulfate, anionic surfactant) and CTAB (cetyltrimethylammoniumbromide, cationic surfactant) at different concentrations. Micellar size, density profiles, and shape (eccentricity) have been investigated. However, the size strongly depends on the functional definition of a micelle. For this reason, we present a method based on the free monomer concentration in aqueous solution as an optimization criterion for the micelle definition. The combination of MD with COSMOmic has the benefit of combining detailed atomistic information from MD with fast calculations of COSMOmic. For the first time the influence of micelle structure on pratition equilibria, predicted with COSMOmic, were investigated. In case of SDS more than 4600 and for CTAB more than 800 single micelles have been studied. The predictions of the partition coefficients with COSMOmic are in good agreement with experimental data. Additionally, the most favorable locations of selected molecules in the micelles as well as probable energy barriers are determined even for complex solutes such as toluene, propanolol, ephedrine, acetone, phenol, lidocaine, syringic acid, coumarin, isovanillin, ferulic acid, and vanillic acid. This method can therefore be applied as a potential screening tool for solutes (e.g., drugs) to find the optimal solute-surfactant combination.

A creature with a hundred waggly tails: intrinsically disordered proteins in the ribosome

Intrinsic disorder (i.e., lack of a unique 3-D structure) is a common phenomenon, and many biologically active proteins are disordered as a whole, or contain long disordered regions. These intrinsically disordered proteins/regions constitute a significant part of all proteomes, and their functional repertoire is complementary to functions of ordered proteins. In fact, intrinsic disorder represents an important driving force for many specific functions. An illustrative example of such disorder-centric functional class is RNA-binding proteins. In this study, we present the results of comprehensive bioinformatics analyses of the abundance and roles of intrinsic disorder in 3,411 ribosomal proteins from 32 species. We show that many ribosomal proteins are intrinsically disordered or hybrid proteins that contain ordered and disordered domains. Predicted globular domains of many ribosomal proteins contain noticeable regions of intrinsic disorder. We also show that disorder in ribosomal proteins has different characteristics compared to other proteins that interact with RNA and DNA including overall abundance, evolutionary conservation, and involvement in protein–protein interactions. Furthermore, intrinsic disorder is not only abundant in the ribosomal proteins, but we demonstrate that it is absolutely necessary for their various functions.

Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach

Amyloid β (Aβ) peptides are a primary component of fibrils and oligomers implicated in the etiology of Alzheimer's disease (AD). However, the intrinsic flexibility of these peptides has frustrated efforts to investigate the secondary and tertiary structure of Aβ monomers, whose conformational landscapes directly contribute to the kinetics and thermodynamics of Aβ aggregation. In this work, de novo replica exchange molecular dynamics (REMD) simulations on the microseconds-per-replica timescale are used to characterize the structural ensembles of Aβ42, Aβ40, and M35-oxidized Aβ42, three physiologically relevant isoforms with substantially different aggregation properties. J-coupling data calculated from the REMD trajectories were compared to corresponding NMR-derived values acquired through two different pulse sequences, revealing that all simulations converge on the order of hundreds of nanoseconds-per-replica toward ensembles that yield good agreement with experiment. Though all three Aβ species adopt highly heterogeneous ensembles, these are considerably more structured compared to simulations on shorter timescales. Prominent in the C-terminus are antiparallel β-hairpins between L17-A21, A30-L36, and V39-I41, similar to oligomer and fibril intrapeptide models that expose these hydrophobic side chains to solvent and may serve as hotspots for self-association. Compared to reduced Aβ42, the absence of a second β-hairpin in Aβ40 and the sampling of alternate β topologies by M35-oxidized Aβ42 may explain the reduced aggregation rates of these forms. A persistent V24-K28 bend motif, observed in all three species, is stabilized by buried backbone to side-chain hydrogen bonds with D23 and a cross-region salt bridge between E22 and K28, highlighting the role of the familial AD-linked E22 and D23 residues in Aβ monomer folding. These characterizations help illustrate the conformational landscapes of Aβ monomers at atomic resolution and provide insight into the early stages of Aβ aggregation pathways.

Systematic optimization model and algorithm for binding sequence selection in computational enzyme design

A systematic optimization model for binding sequence selection in computational enzyme design was developed based on the transition state theory of enzyme catalysis and graph-theoretical modeling. The saddle point on the free energy surface of the reaction system was represented by catalytic geometrical constraints, and the binding energy between the active site and transition state was minimized to reduce the activation energy barrier. The resulting hyperscale combinatorial optimization problem was tackled using a novel heuristic global optimization algorithm, which was inspired and tested by the protein core sequence selection problem. The sequence recapitulation tests on native active sites for two enzyme catalyzed hydrolytic reactions were applied to evaluate the predictive power of the design methodology. The results of the calculation show that most of the native binding sites can be successfully identified if the catalytic geometrical constraints and the structural motifs of the substrate are taken into account. Reliably predicting active site sequences may have significant implications for the creation of novel enzymes that are capable of catalyzing targeted chemical reactions.

Confinement-dependent friction in peptide bundles

Friction within globular proteins or between adhering macromolecules crucially determines the kinetics of protein folding, the formation, and the relaxation of self-assembled molecular systems. One fundamental question is how these friction effects depend on the local environment and in particular on the presence of water. In this model study, we use fully atomistic MD simulations with explicit water to obtain friction forces as a single polyglycine peptide chain is pulled out of a bundle of k adhering parallel polyglycine peptide chains. The whole system is periodically replicated along the peptide axes, so a stationary state at prescribed mean sliding velocity V is achieved. The aggregation number is varied between k = 2 (two peptide chains adhering to each other with plenty of water present at the adhesion sites) and k = 7 (one peptide chain pulled out from a close-packed cylindrical array of six neighboring peptide chains with no water inside the bundle). The friction coefficient per hydrogen bond, extrapolated to the viscous limit of vanishing pulling velocity V → 0, exhibits an increase by five orders of magnitude when going from k = 2 to k = 7. This dramatic confinement-induced friction enhancement we argue to be due to a combination of water depletion and increased hydrogen-bond cooperativity.

Protein molecular surface mapped at different geometrical resolutions

Many areas of biochemistry and molecular biology, both fundamental and applications-orientated, require an accurate construction, representation and understanding of the protein molecular surface and its interaction with other, usually small, molecules. There are however many situations when the protein molecular surface gets in physical contact with larger objects, either biological, such as membranes, or artificial, such as nanoparticles. The contribution presents a methodology for describing and quantifying the molecular properties of proteins, by geometrical and physico-chemical mapping of the molecular surfaces, with several analytical relationships being proposed for molecular surface properties. The relevance of the molecular surface-derived properties has been demonstrated through the calculation of the statistical strength of the prediction of protein adsorption. It is expected that the extension of this methodology to other phenomena involving proteins near solid surfaces, in particular the protein interaction with nanoparticles, will result in important benefits in the understanding and design of protein-specific solid surfaces.

Intrinsically disordered regions of p53 family are highly diversified in evolution

Proteins of the p53 family are expressed in vertebrates and in some invertebrate species. The main function of these proteins is to control and regulate cell cycle in response to various cellular signals, and therefore to control the organism's development. The regulatory functions of the p53 family members originate mostly from their highly-conserved and well-structured DNA-binding domains. Many human diseases (including various types of cancer) are related to the missense mutations within this domain. The ordered DNA-binding domains of the p53 family members are surrounded by functionally important intrinsically disordered regions. In this study, substitution rates and propensities in different regions of p53 were analyzed. The analyses revealed that the ordered DNA-binding domain is conserved, whereas disordered regions are characterized by high sequence diversity. This diversity was reflected both in the number of substitutions and in the types of substitutions to which each amino acid was prone. These results support the existence of a positive correlation between protein intrinsic disorder and sequence divergence during the evolutionary process. This higher sequence divergence provides strong support for the existence of disordered regions in p53 in vivo for if they were structured, they would evolve at similar rates as the rest of the protein.

The effect of a G:T mispair on the dynamics of DNA

Distortions in the DNA sequence such as damages or mispairs are specifically recognized and processed by DNA repair enzymes. A particular challenge for the enzymatic specificity is the recognition of a wrongly-placed native nucleotide such as thymine in T:G mispairs. An important step of substrate binding which is observed in many repair proteins is the flipping of the target base out of the DNA helix into the enzyme's active site. In this work we investigate how much the intrinsic dynamics of mispaired DNA is changed compared to canonical DNA. Our molecular dynamics simulations of DNA with and without T:G mispairs show significant differences in the conformation of paired and mispaired DNA. The wobble pair T:G shows local distortions such as twist, shear and stretch which deviate from canonical B form values. Moreover, the T:G mispair is found to be kinetically less stable, exhibiting two states with respect to base opening: a closed state comparable to the canonical base pairs, and a more open state, indicating a proneness for base flip. In addition, we observe that the thymine base in a T:G mispair is significantly more probable to be flipped than thymine in a T:A pair or cytosine in a C:G pair. Such local deformations and in particular the existence of a second, more-open state can be speculated to help the target-site recognition by repair enzymes.

A solvated ligand rotamer approach and its application in computational protein design

The structure-based design of protein-ligand interfaces with respect to different small molecules is of great significance in the discovery of functional proteins. By statistical analysis of a set of protein-ligand complex structures, it was determined that water-mediated hydrogen bonding at the protein-ligand interface plays a crucial role in governing the binding between the protein and the ligand. Based on the novel statistic results, a solvated ligand rotamer approach was developed to explicitly describe the key water molecules at the protein-ligand interface and a water-mediated hydrogen bonding model was applied in the computational protein design context to complement the continuum solvent model. The solvated ligand rotamer approach produces only one additional solvated rotamer for each rotamer in the ligand rotamer library and does not change the number of side-chain rotamers at each protein design site. This has greatly reduced the total combinatorial number in sequence selection for protein design, and the accuracy of the model was confirmed by two tests. For the water placement test, 61% of the crystal water molecules were predicted correctly in five protein-ligand complex structures. For the sequence recapitulation test, 44.7% of the amino acid identities were recovered using the solvated ligand rotamer approach and the water-mediated hydrogen bonding model, while only 30.4% were recovered when the explicitly bound waters were removed. These results indicated that the developed solvated ligand rotamer approach is promising for functional protein design targeting novel protein-ligand interactions.

Structural features of [NiFeSe] and [NiFe] hydrogenases determining their different properties: a computational approach

Hydrogenases are metalloenzymes that catalyze the reversible reaction H(2)<->2H(+) + 2e(-), being potentially useful in H(2) production or oxidation. [NiFeSe] hydrogenases are a particularly interesting subgroup of the [NiFe] class that exhibit tolerance to O(2) inhibition and produce more H(2) than standard [NiFe] hydrogenases. However, the molecular determinants responsible for these properties remain unknown. Hydrophobic pathways for H(2) diffusion have been identified in [NiFe] hydrogenases, as have proton transfer pathways, but they have never been studied in [NiFeSe] hydrogenases. Our aim was, for the first time, to characterize the H(2) and proton pathways in a [NiFeSe] hydrogenase and compare them with those in a standard [NiFe] hydrogenase. We performed molecular dynamics simulations of H(2) diffusion in the [NiFeSe] hydrogenase from Desulfomicrobium baculatum and extended previous simulations of the [NiFe] hydrogenase from Desulfovibrio gigas (Teixeira et al. in Biophys J 91:2035-2045, 2006). The comparison showed that H(2) density near the active site is much higher in [NiFeSe] hydrogenase, which appears to have an alternative route for the access of H(2) to the active site. We have also determined a possible proton transfer pathway in the [NiFeSe] hydrogenase from D. baculatum using continuum electrostatics and Monte Carlo simulation and compared it with the proton pathway we found in the [NiFe] hydrogenase from D. gigas (Teixeira et al. in Proteins 70:1010-1022, 2008). The residues constituting both proton transfer pathways are considerably different, although in the same region of the protein. These results support the hypothesis that some of the special properties of [NiFeSe] hydrogenases could be related to differences in the H(2) and proton pathways.

Specific and nonspecific effects of glycosylation

Glycosylation regulates vital cellular processes and dramatically influences protein folding and stability. In particular, experiments have demonstrated that asparagine (N)-linked disaccharides drive a "conformational switch" in a model peptide. The present work investigates this conformational switch via extensive atomically detailed replica exchange molecular dynamics simulations in explicit solvent. To distinguish the effects of specific and nonspecific interactions upon the peptide conformational ensemble, these simulations considered model peptides that were N-linked to a disaccharide and to a steric crowder of the same shape. The simulations are remarkably consistent with experiment and provide detailed insight into the peptide structure ensemble. They suggest that steric crowding by N-linked disaccharides excludes extended conformations, but does not significantly impact the tetrahedral structure of the surrounding solvent or otherwise alter the peptide free energy surface. However, the combination of steric crowding with specific hydrogen bonds and hydrophobic stacking interactions more dramatically impacts the peptide ensemble and stabilizes new structures.

Distinct dimerization for various alloforms of the amyloid-beta protein: Aβ(1-40), Aβ(1-42), and Aβ(1-40)(D23N)

The Amyloid-beta protein is related to Alzheimer's disease, and various experiments have shown that oligomers as small as the dimer are cytotoxic. Two alloforms are mainly produced: Aβ(1-40) and Aβ(1-42). They have very different oligomer distributions, and it was recently suggested, from experimental studies, that this variation may originate from structural differences in their dimer structures. Little structural information is available on the Aβ dimer, however, and to complement experimental observations, we simulated the folding of the wild-type Aβ(1-40) and Aβ(1-42) dimers as well as the mutated Aβ(1-40)(D23N) dimer using an accurate coarse-grained force field coupled to Hamiltonian-temperature replica exchange molecular dynamics. The D23N variant impedes the salt-bridge formation between D23 and K28 seen in the wild-type Aβ, leading to very different fibrillation properties and final amyloid fibrils. Our results show that the Aβ(1-42) dimer has a higher propensity than the Aβ(1-40) dimer to form β-strands at the central hydrophobic core (residues 17-21) and at the C-terminal (residues 30-42), which are two segments crucial to the oligomerization of Aβ. The free energy landscape of the Aβ(1-42) dimer is also broader and more complex than that of the Aβ(1-40) dimer. Interestingly, D23N also impacts the free energy landscape by increasing the population of configurations with higher β-strand propensities when compared against Aβ(40). In addition, while Aβ(1-40)(D23N) displays a higher β-strand propensity at the C-terminal, its solvent accessibility does not change with respect to the wild-type sequence. Overall, our results show the strong impact of the two amino acids Ile41-Ala42 and the salt-bridge D23-K28 on the folding of the Aβ dimer.

Beta-decomposition for the volume and area of the union of three-dimensional balls and their offsets

Given a set of spherical balls, called atoms, in three-dimensional space, its mass properties such as the volume and the boundary area of the union of the atoms are important for many disciplines, particularly for computational chemistry/biology and structural molecular biology. Despite many previous studies, this seemingly easy problem of computing mass properties has not been well-solved. If the mass properties of the union of the offset of the atoms are to be computed as well, the problem gets even harder. In this article, we propose algorithms that compute the mass properties of both the union of atoms and their offsets both correctly and efficiently. The proposed algorithms employ an approach, called the Beta-decomposition, based on the recent theory of the beta-complex. Given the beta-complex of an atom set, these algorithms decompose the target mass property into a set of primitives using the simplexes of the beta-complex. Then, the molecular mass property is computed by appropriately summing up the mass property corresponding to each simplex. The time complexity of the proposed algorithm is O(m) in the worst case where m is the number of simplexes in the beta-complex that can be efficiently computed from the Voronoi diagram of the atoms. It is known in ℝ(3) that m = O(n) on average for biomolecules and m = O(n(2)) in the worst case for general spheres where n is the number of atoms. The theory is first introduced in ℝ(2) and extended to ℝ(3). The proposed algorithms were implemented into the software BetaMass and thoroughly tested using molecular structures available in the Protein Data Bank. BetaMass is freely available at the Voronoi Diagram Research Center web site.

Hexapeptide fragment of carcinoembryonic antigen which acts as an agonist of heterogeneous ribonucleoprotein M

Colorectal cancers with metastatic potential secrete the glycoprotein carcinoembryonic antigen (CEA). CEA has been implicated in colorectal cancer metastasis by inducing Kupffer cells to produce inflammatory cytokines which, in turn, make the hepatic micro-environment ideal for tumor cell implantation. CEA binds to the heterogeneous ribonucleoprotein M (hnRNP M) which acts as a cell surface receptor in Kupffer cells. The amino acid sequence in CEA, which binds the hnRNP M receptor, is Tyr-Pro-Glu-Leu-Pro-Lys. In this study, the structure of Ac-Tyr-Pro-Glu-Leu-Pro-Lys-NH₂ (YPELPK) was investigated using electronic circular dichroism, vibrational circular dichroism, and molecular dynamics simulations. The binding of the peptide to hnRNP M was also investigated using molecular docking calculations. The biological activity of YPELPK was studied using differentiated human THP-1 cells, which express hnRNP M on their surface and secrete IL-6 when stimulated by CEA. YPELPK forms a stable polyproline-II helix and stimulates IL-6 production of THP-1 cells at micromolar concentrations.

On-the-fly Numerical Surface Integration for Finite-Difference Poisson-Boltzmann Methods

Most implicit solvation models require the definition of a molecular surface as the interface that separates the solute in atomic detail from the solvent approximated as a continuous medium. Commonly used surface definitions include the solvent accessible surface (SAS), the solvent excluded surface (SES), and the van der Waals surface. In this study, we present an efficient numerical algorithm to compute the SES and SAS areas to facilitate the applications of finite-difference Poisson-Boltzmann methods in biomolecular simulations. Different from previous numerical approaches, our algorithm is physics-inspired and intimately coupled to the finite-difference Poisson-Boltzmann methods to fully take advantage of its existing data structures. Our analysis shows that the algorithm can achieve very good agreement with the analytical method in the calculation of the SES and SAS areas. Specifically, in our comprehensive test of 1,555 molecules, the average unsigned relative error is 0.27% in the SES area calculations and 1.05% in the SAS area calculations at the grid spacing of 1/2Å. In addition, a systematic correction analysis can be used to improve the accuracy for the coarse-grid SES area calculations, with the average unsigned relative error in the SES areas reduced to 0.13%. These validation studies indicate that the proposed algorithm can be applied to biomolecules over a broad range of sizes and structures. Finally, the numerical algorithm can also be adapted to evaluate the surface integral of either a vector field or a scalar field defined on the molecular surface for additional solvation energetics and force calculations.

Stable conformations of a single stranded deprotonated DNA i-motif

We present molecular dynamics simulations of a single stranded deprotonated DNA i-motif in explicit solvent. Our results indicate that hairpin structures are stable equilibrium conformations at 300 K. The entropic preference of these configurations is explained by strong water ordering effects due to the present number of hydrogen bonds. We observe a full unfolding at higher temperatures in good agreement with experimental results.

Comparison of the structural basis for thermal stability between archaeal and bacterial proteins

In this study, the structural basis for thermal stability in archaeal and bacterial proteins was investigated. There were many common factors that confer resistance to high temperature in both archaeal and bacterial proteins. These factors include increases in the Lys content, the bends and blanks of secondary structure, the Glu content of salt bridge; decreases in the number of main-side chain hydrogen bond and exposed surface area, and changes in the bends and blanks of amino acids. Certainly, the utilization of charged amino acids to form salt bridges is a primary factor. In both heat-resistant archaeal and bacterial proteins, most Glu and Asp participate in the formation of salt bridges. Other factors may influence either archaeal or bacterial protein thermostability, which includes the more frequent occurrence of shorter 3(10)-helices and increased hydrophobicity in heat-resistant archaeal proteins. However, there were increases in average helix length, the Glu content in salt bridges, temperature factors and decreases in the number of main-side chain hydrogen bonds, uncharged-uncharged hydrogen bonds, hydrophobicity, and buried and exposed polar surface area in heat-resistant bacterial proteins. Evidently, there are few similarities and many disparities between the heat-resistant mechanisms of archaeal and bacterial proteins.

Binding affinity of a small molecule to an amorphous polymer in a solvent. Part 2: preferential binding to local sites on a surface

Crystallization, a separation and purification process, is commonly used to produce a wide range of materials in various industries, and it usually begins with heterogeneous nucleation on a foreign surface in industrial practice and most other circumstances. Recent studies show that amorphous polymeric substrates are useful in controlling crystallization and selectively producing pharmaceutical polymorphs. In our previous publication, we investigated the possible correlation of the binding affinity of one molecule to key binding sites (local binding), and the possibility of using this binding affinity to guide the selection of polymers promoting heterogeneous nucleation. The studied systems were aspirin binding to four nonporous cross-linked polymers in ethanol-water 38 v% mixture. Cross-linked with divinylbenzene (DVB), these polymers were poly(4-acryloylmorpholine) (PAM), poly(2-carboxyethyl acrylate) (PCEA), poly(4-hydroxylbutyl acrylate) (PHBA), and polystyrene (PS). We discovered that the trend of the magnitudes of the average free energies of binding to the best sites is very similar to that of heterogeneous nucleation activities. This Article aims to investigate whether or not local binding to key sites is the important variable to describe heterogeneous nucleation as opposed to the overall/average binding affinity of molecules to a surface, and to investigate the possibility of using the overall binding affinity to guide the selection of polymers. We used the polymer surfaces generated from our previous study to calculate the overall binding affinity of aspirin molecules to the surface as measured by the preferential interaction coefficients of aspirin (1 m) to these polymers. We discovered that the trend of the average preferential interaction coefficients does not correlate as well to that of heterogeneous nucleation activities as the free energies of binding to the best sites. We also computed the average numbers of aspirin molecules associated with the areas of the surfaces' best binding sites and found that they correlate better to heterogeneous nucleation activities than the average preferential interaction coefficients. These results further support that local binding is indicative of heterogeneous nucleation. Moreover, we found a weak trend of the distance order parameters of the aspirin molecules to be similar that of heterogeneous nucleation activities. Our results from the two-part study suggest the importance of local binding to heterogeneous nucleation as well as the possibility of using the binding affinity to the local area (the free energy of binding to the best site and the number of nucleating molecules associated with the area of the best binding site) and the distance order parameters to guide the selection of polymers.