An efficient alpha helix model and simulation framework for stationary electrostatic interaction force estimation.
Ontology highlight
ABSTRACT:
SUBMITTER: Butcher GG
PROVIDER: S-EPMC8079384 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Aligned protein α-helix dipoles have been implicated in protein function and structure. The recent breakthroughs in high-resolution electron microscopy (EM) of macromolecules makes it possible to explore fundamental aspects of structural biology at the detailed molecular level. The electrostatic potential (ESP) generated by aligned protein α-helix dipole should be observable in high-resolution EM maps despite the fact that the effect may be partially screened by induced electric fields. Here, we show that aligned backbone dipoles in protein α-helices account for long-range features in the protein ESP functions. Our results are consistent with experimental EM maps and density functional theory calculations, including direct Fourier summation for proper calculation of the ESP due to the nonlocal nature of the ESP function from aligned dipoles and other partial atomic charges.
| S-EPMC5563131 | biostudies-literature
Project description:In this paper, we put forward a real-time multiple GPUs (multi-GPU) accelerated virtual-reality interaction simulation framework where the reconstructed objects from camera images interact with virtual deformable objects. Firstly, based on an extended voxel-based visual hull (VbVH) algorithm, we design an image-based 3D reconstruction platform for real objects. Then, an improved hybrid deformation model, which couples the geometry constrained fast lattice shape matching method (FLSM) and total Lagrangian explicit dynamics (TLED) algorithm, is proposed to achieve efficient and stable simulation of the virtual objects' elastic deformations. Finally, one-way virtual-reality interactions including soft tissues' virtual cutting with bleeding effects are successfully simulated. Moreover, with the purpose of significantly improving the computational efficiency of each time step, we propose an entire multi-GPU implementation method of the framework using compute unified device architecture (CUDA). The experiment results demonstrate that our multi-GPU accelerated virtual-reality interaction framework achieves real-time performance under the moderate calculation scale, which is a new effective 3D interaction technique for virtual reality applications.
| S-EPMC6459549 | biostudies-literature
Project description:Contact and non-contact based atomic force microscopy (AFM) approaches have been extensively utilized to explore various nanoscale surface properties. In most AFM-based measurements, a concurrent electrostatic effect between the AFM tip/cantilever and sample surface can occur. This electrostatic effect often hinders accurate measurements. Thus, it is very important to quantify as well as remove the impact of the electrostatic effect on AFM-based measurements. In this study, we examine the impact of the electrostatic effect on the electromechanical (EM) response in piezoresponse force microscopy as a model AFM mode. We quantitatively studied the effects of increasing the external electric field and reducing the spring constant of a cantilever. Further, we explored ways to minimize the electrostatic effect. The results provide broad guidelines for quantitatively analyzing the EM response as well as, eventually, for obtaining the electrostatic-free EM response. The conclusions can be applied to other AFM-based measurements that are subject to a strong electrostatic effect between the AFM tip/cantilever and sample surface, regardless of contact and non-contact modes.
| S-EPMC5282565 | biostudies-literature
Project description:While offering high resolution atomic and electronic structure, scanning probe microscopy techniques have found greater challenges in providing reliable electrostatic characterization on the same scale. In this work, we offer electrostatic discovery atomic force microscopy, a machine learning based method which provides immediate maps of the electrostatic potential directly from atomic force microscopy images with functionalized tips. We apply this to characterize the electrostatic properties of a variety of molecular systems and compare directly to reference simulations, demonstrating good agreement. This approach offers reliable atomic scale electrostatic maps on any system with minimal computational overhead.
| S-EPMC8793147 | biostudies-literature
Project description:One of the fundamental events in protein folding is α-helix formation, which involves sequential development of a series of helical hydrogen bonds between the backbone C═O group of residues i and the -NH group of residues i + 4. While we now know a great deal about α-helix folding dynamics, a key question that remains to be answered is where the productive helical nucleation event occurs. Statistically, a helical nucleus (or the first helical hydrogen-bond) can form anywhere within the peptide sequence in question; however, the one that leads to productive folding may only form at a preferred location. This consideration is based on the fact that the α-helical structure is inherently asymmetric, due to the specific alignment of the helical hydrogen bonds. While this hypothesis is plausible, validating it is challenging because there is not an experimental observable that can be used to directly pinpoint the location of the productive nucleation process. Therefore, in this study we combine several techniques, including peptide cross-linking, laser-induced temperature-jump infrared spectroscopy, and molecular dynamics simulations, to tackle this challenge. Taken together, our experimental and simulation results support an α-helix folding mechanism wherein the productive nucleus is formed at the N-terminus, which propagates toward the C-terminal end of the peptide to yield the folded structure. In addition, our results show that incorporation of a cross-linker can lead to formation of differently folded conformations, underscoring the need for all-atom simulations to quantitatively assess the proposed cross-linking design.
| S-EPMC6497059 | biostudies-literature
Project description:The multilevel summation method (MSM) offers an efficient algorithm utilizing convolution for evaluating long-range forces arising in molecular dynamics simulations. Shifting the balance of computation and communication, MSM provides key advantages over the ubiquitous particle–mesh Ewald (PME) method, offering better scaling on parallel computers and permitting more modeling flexibility, with support for periodic systems as does PME but also for semiperiodic and nonperiodic systems. The version of MSM available in the simulation program NAMD is described, and its performance and accuracy are compared with the PME method. The accuracy feasible for MSM in practical applications reproduces PME results for water property calculations of density, diffusion constant, dielectric constant, surface tension, radial distribution function, and distance-dependent Kirkwood factor, even though the numerical accuracy of PME is higher than that of MSM. Excellent agreement between MSM and PME is found also for interface potentials of air–water and membrane–water interfaces, where long-range Coulombic interactions are crucial. Applications demonstrate also the suitability of MSM for systems with semiperiodic and nonperiodic boundaries. For this purpose, simulations have been performed with periodic boundaries along directions parallel to a membrane surface but not along the surface normal, yielding membrane pore formation induced by an imbalance of charge across the membrane. Using a similar semiperiodic boundary condition, ion conduction through a graphene nanopore driven by an ion gradient has been simulated. Furthermore, proteins have been simulated inside a single spherical water droplet. Finally, parallel scalability results show the ability of MSM to outperform PME when scaling a system of modest size (less than 100 K atoms) to over a thousand processors, demonstrating the suitability of MSM for large-scale parallel simulation.
| S-EPMC4325600 | biostudies-literature
Project description:Electrostatic interactions are important for understanding molecular interactions, since they are long-range interactions and can guide binding partners to their correct binding positions. To investigate the role of electrostatic forces in molecular recognition, we calculated electrostatic forces between binding partners separated at various distances. The investigation was done on a large set of 275 protein complexes using recently developed DelPhiForce tool and in parallel, evaluating the total electrostatic force via electrostatic association energy. To accomplish the goal, we developed a method to find an appropriate direction to move one chain of protein complex away from its bound position and then calculate the corresponding electrostatic force as a function of separation distance. It is demonstrated that at large distances between the partners, the electrostatic force (magnitude and direction) is consistent among the protocols used and the main factors contributing to it are the net charge of the partners and their interfaces. However, at short distances, where partners form specific pair-wise interactions or de-solvation penalty becomes significant, the outcome depends on the precise balance of these factors. Based on the electrostatic force profile (force as a function of distance), we group the cases into four distinctive categories, among which the most intriguing is the case termed "soft landing." In this case, the electrostatic force at large distances is favorable assisting the partners to come together, while at short distance it opposes binding, and thus slows down the approach of the partners toward their physical binding.
| S-EPMC6774301 | biostudies-literature
Project description:Metal ions are crucial for nucleic acid folding. From the free energy landscapes, we investigate the detailed mechanism for ion-induced collapse for a paradigm system: loop-tethered short DNA helices. We find that Na+ and Mg2+ play distinctive roles in helix-helix assembly. High [Na+] (>0.3 M) causes a reduced helix-helix electrostatic repulsion and a subsequent disordered packing of helices. In contrast, Mg2+ of concentration >1 mM is predicted to induce helix-helix attraction and results in a more compact and ordered helix-helix packing. Mg2+ is much more efficient in causing nucleic acid compaction. In addition, the free energy landscape shows that the tethering loops between the helices also play a significant role. A flexible loop, such as a neutral loop or a polynucleotide loop in high salt concentration, enhances the close approach of the helices in order to gain the loop entropy. On the other hand, a rigid loop, such as a polynucleotide loop in low salt concentration, tends to de-compact the helices. Therefore, a polynucleotide loop significantly enhances the sharpness of the ion-induced compaction transition. Moreover, we find that a larger number of helices in the system or a smaller radius of the divalent ions can cause a more abrupt compaction transition and a more compact state at high ion concentration, and the ion size effect becomes more pronounced as the number of helices is increased.
| S-EPMC1751542 | biostudies-literature
Project description:Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38?s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.
| S-EPMC6503141 | biostudies-literature
Project description:RNA pseudoknots are structural elements that participate in a variety of biological processes. At -1 ribosomal frameshifting sites, several types of pseudoknot have been identified which differ in their organisation and functionality. The pseudoknot found in infectious bronchitis virus (IBV) is typical of those that possess a long stem 1 of 11-12 bp and a long loop 2 (30-164 nt). A second group of pseudoknots are distinguishable that contain stems of only 5 to 7 bp and shorter loops. The NMR structure of one such pseudoknot, that of mouse mammary tumor virus (MMTV), has revealed that it is kinked at the stem 1-stem 2 junction, and that this kinked conformation is essential for efficient frameshifting. We recently investigated the effect on frameshifting of modulating stem 1 length and stability in IBV-based pseudoknots, and found that a stem 1 with at least 11 bp was needed for efficient frameshifting. Here, we describe the sequence manipulations that are necessary to bypass the requirement for an 11 bp stem 1 and to convert a short non-functional IBV-derived pseudoknot into a highly efficient, kinked frameshifter pseudoknot. Simple insertion of an adenine residue at the stem 1-stem 2 junction (an essential feature of a kinked pseudoknot) was not sufficient to create a functional pseudoknot. An additional change was needed: efficient frameshifting was recovered only when the last nucleotide of loop 2 was changed from a G to an A. The requirement for an A at the end of loop 2 is consistent with a loop-helix contact similar to those described in other RNA tertiary structures. A mutational analysis of both partners of the proposed interaction, the loop 2 terminal adenine residue and two G.C pairs near the top of stem 1, revealed that the interaction was essential for efficient frameshifting. The specific requirement for a 3'-terminal A residue was lost when loop 2 was increased from 8 to 14 nt, suggesting that the loop-helix contact may be required only in those pseudoknots with a short loop 2.
| S-EPMC7141562 | biostudies-literature