SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The QM/MM MD and free energy simulations show that the dynamics involving a His residue at the P1 site of the substrate may play an important role in substrate-assisted catalysis and specificity for a serine-carboxyl peptidase.
The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC).d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usualcut-of of the long range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other ligands.
BACKGROUND: Defects in PPP1R12A can lead to genitourinary and/or brain malformation syndrome (GUBS). GUBS is primarily characterized by neurological or genitourinary system abnormalities, but a few reported cases are associated with neonatal seizures. Here, we report a case of a female newborn with neonatal seizures caused by a novel variant in PPP1R12A, aiming to enhance the clinical and variant data of genetic factors related to epilepsy in early life. METHODS: Whole-exome and Sanger sequencing were used for familial variant assessment, and bioinformatics was employed to annotate the variant. A structural model of the mutant protein was simulated using molecular dynamics (MD), and the free binding energy between PPP1R12A and PPP1CB was analyzed. A mutant plasmid was constructed, and mutant protein expression was analyzed using western blotting (WB), and the interaction between the mutant and PPP1CB proteins using co-immunoprecipitation (Co-IP) experiments. RESULTS: The patient experienced tonic-clonic seizures on the second day after birth. Genetic testing revealed a heterozygous variant in PPP1R12A, NM_002480.3:c.2533 C > T (p.Arg845Ter). Both parents had the wild-type gene. MD suggested that loss of the C-terminal structure in the mutant protein altered its structural stability and increased the binding energy with PPP1CB, indicating unstable protein-protein interactions. On WB, a low-molecular-weight band was observed, indicating that the protein was truncated. Co-IP indicated that the mutant protein no longer interacted with PPP1CB, indicating an effect on the structural stability of the myosin phase complex. CONCLUSION: The PPP1R12A c.2533 C > T variant may explain the neonatal seizures in the present case. The findings of this study expand the spectrum of PPP1R12A variants and highlight the potential significance of truncated proteins in the pathogenesis of GUBS.
The two main structural differences between calmodulin (CaM) and skeletal troponin C (sTnC) are the absence in CaM of (i) the short N-terminal helix in TnC and (ii) the triplet KGK (residues 91-93; numbering according to chicken sTnC). It was recently shown that deletion of both structural groups from sTnC imparted to the resulting construct the CaM-like ability to activate phosphodiesterase (PDE) and to regulate force development in smooth muscle. To continue probing of the structural basis of the differential behavior of sTnC and CaM, residue Arg-11 in rabbit sTnC was mutated to Ala because the interactions of Arg-11 with distal residues in the N-terminal domain seem to link the N-terminal helix to the rest of the structure. The mutant exhibits CaM-like function in its ability to activate PDE (about 50% of CaM at 5 microM concentration). If, in addition, the KGK triplet is also deleted, PDE activation increases to about 80%. Both constructs retain their TnC function to nearly 100%. To explore the mechanistic basis of this remarkable observation, computational simulations of the molecular dynamics (MD) were carried out for both wild-type 4Ca2+.sTnC and the 4Ca2+.R11A mutant, and the results were compared to those from earlier simulations of 4Ca2+.CaM. Two types of structural changes observed from such simulations of the molecular dynamics of CaM had been considered to have a functional role: (i) a compaction to a more globular form and (ii) a reorientation of the Ca-binding domains around the central tether helix.(ABSTRACT TRUNCATED AT 250 WORDS)
An efficient Monte Carlo (MC) algorithm including concerted rotations is directly compared to molecular dynamics (MD) in all-atom statistical mechanics folding simulations of small polypeptides. The previously reported algorithm "concerted rotations with flexible bond angles" (CRA) has been shown to successfully locate the native state of small polypeptides. In this study, the folding of three small polypeptides (trpzip2/H1/Trp-cage) is investigated using MC and MD, for a combined sampling time of approximately 10(11) MC configurations and 8 micros, respectively. Both methods successfully locate the experimentally determined native states of the three systems, but they do so at different speed, with 2-2.5 times faster folding of the MC runs. The comparison reveals that thermodynamic and dynamic properties can reliably be obtained by both and that results from folding simulations do not depend on the algorithm used. Similar to previous comparisons of MC and MD, it is found that one MD integration step of 2 fs corresponds to one MC scan, revealing the good sampling of MC. The simplicity and efficiency of the MC method will enable its future use in folding studies involving larger systems and the combination with replica exchange algorithms.
This paper presents recent advances in the use of molecular simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy, which enable us to understand solvated ions in solution. We report and discuss the EXAFS spectra and related properties governing solvation processes of different ions in water and methanol. Molecular dynamics (MD) trajectories are coupled to electron scattering simulations to generate the MD-EXAFS spectra, which are found to be in very good agreement with the corresponding experimental measurements. From these simulated spectra, the ion-oxygen distances for the first hydration shell are in agreement with experiment within 0.05-0.1 A. The ionic species studied range from monovalent to divalent, positive and negative: K+, Ca2+, and Cl-. This work demonstrates that the combination of MD-EXAFS and the corresponding experimental measurement provides a powerful tool in the analysis of the solvation structure of aqueous ionic solutions. We also investigate the value of electronic structure analysis of small aqueous clusters as a benchmark to the empirical potentials. In a novel computational approach, we determine the Debye-Waller factors for Ca2+, K+, and Cl- in water by combining the harmonic analysis of data obtained from electronic structure calculations on finite ion-water clusters, providing excellent agreement with the experimental values, and discuss how they compare with results from a harmonic classical statistical mechanical analysis of an empirical potential.
Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.
Here we report molecular dynamics (MD) and free energy perturbation (FEP) simulations applied to hydroxamate-matrix metalloproteinase-2 (MMP-2) complex systems. We have developed some new force field parameters for the hydroxamate functional group that were not included in the AMBER94 force field but were necessary in our simulations. For the representation of the active zinc center, a bonded model was adopted in which restrained electrostatic potential fitting (RESP) charges were used as the electrostatic representation of this model. Using the resulted bonded model, FEP simulations predict the relative binding free energy in good agreement with the experimental value. By analyzing the molecular dynamics (MD) trajectories of the two complex systems, we can provide an explanation of why one of the two inhibitors is favored over the other. The results provide a chemical insight into the interactions between inhibitor and enzyme, and can indicate changes in the inhibitor that would enhance inhibitor-enzyme interactions.
The BTE Work Simulator (Baltimore, Therapeutic Equipment Co., Hanover, Md.) recreates the work environment by providing a rotatable shaft to which a variety of tools can be attached. The patient can simulate the coordinated muscle and joint movements involved in the work situation, exercising either isometrically or isotonically, and thereby increase the potential for return to work. The BTE Work Simulator has four major components: the exercise head to provide variable resistance, tools that attach to the exercise head, a control console containing a microprocessor that calculates work and power, and the Quest software package (Baltimore Therapeutic Equipment Co.). Quest provides four static/isometric and seven dynamic/isotonic tests that can be used to document the extent of a patient's impairment. In addition, it has daily treatment options and acts as an extensive data base capable of retrieving any stored test or treatment results. This software system greatly facilitates the use of the BTE Work Simulator by providing instant calculations of power and coefficients of variation, immediate patient feedback, graphic expansion of the data, and printed evaluation reports.
Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB-lipid bilayer model in order to characterize the molecular aspects of AmB-membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.
This article presents a molecular dynamics (MD) study of the cdk2 enzyme and its two complexes with the inhibitors isopentenyladenine and roscovitine using the Cornell et al. force field from the AMBER software package. The results show that inserting an inhibitor into the enzyme active site does not considerably change enzyme structure but it seemingly changes the distribution of internal motions. The inhibitor causes differences in the domain motions in free cdk2 and in its complexes. It was found out that repulsion of roscovitine N9 substituent causes conformational change on Lys 33 side chain. Isopentenyladenine forms with Lys 33 side chain terminal amino group a hydrogen bond. It implies that the cavity, where N9 substituent of roscovitine is buried, can adopt larger substituent due to Lys 33 side chain flexibility. The composition of electrostatic and van der Waals interactions between the inhibitor and the enzyme were also calculated along both cdk2/inhibitor MD trajectories together with MM-PB/GBSA analysis. These results show that isopentenyladenine-like inhibitors could be more effective after modifications leading to an increase in their van der Waals contact with the enzyme. We suggest that a way leading to better inhibitors occupying isopentenyladenine binding mode could be: to keep N9 and N7 purine positions free, to keep 3,3-dimethylallylamino group at C6 position, and to add, e.g., benzylamino group at C2 position. The results support the idea that the isopentenyladenine binding mode can be used for cdk2 inhibitors design and that all possibilities to improve this binding mode were not uncovered yet.
The Cray MTA-2 (Multithreaded Architecture) is an unusual parallel supercomputer that promises ease of use and high performance. We describe our experience on the MTA-2 with a molecular dynamics code, SIMU-MD, that we are using to simulate the translocation of DNA through a nanopore in a silicon based ultrafast sequencer. Our sequencer is constructed using standard VLSI technology and consists of a nanopore surrounded by Field Effect Transistors (FETs). We propose to use the FETs to sense variations in charge as a DNA molecule translocates through the pore and thus differentiate between the four building block nucleotides of DNA. We were able to port SIMU-MD, a serial C code, to the MTA with only a modest effort and with good performance. Our porting process needed neither a parallelism support platform nor attention to the intimate details of parallel programming and interprocessor communication, as would have been the case with more conventional supercomputers.
Previously published studies on transfers to restore elbow flexion have rarely compared the outcomes among different techniques. We found no comparison between the Steindler procedure and pectoralis major tendon transfer. This study compared functional outcomes after these two types of elbow flexorplasty performed in conjunction with shoulder arthrodesis for patients with brachial plexus injuries. The same examiner assessed 14 patients who underwent shoulder fusion and elbow flexorplasty. Five had a pectoralis major tendon transfer, and nine underwent a modified Steindler procedure. All patients had brachial plexus lesions that were irreparable or had not responded to neurolysis, repair, or grafting. Average follow-up time was 60 months (range, 12 months-10 years). All patients completed a questionnaire (pain, function) and underwent an objective assessment of range of motion, isometric strength, isotonic power, grip, and pinch strength. The BTE Work Simulator (Baltimore Therapeutic Equipment, Hanover, MD) was used for muscle strength testing. Improvement was noted in daily activities including dressing and bilateral activities and was reflected in a significant improvement over preoperative function. No differences (at p = .05 level of significance) were found in objective tests between the pectoralis major transfer and the Steindler groups. In the patient requiring stabilization of the shoulder and flexorplasty, the pectoralis major tendon transfer was at least equivalent to the modified Steindler flexorplasty in terms of range of motion, strength, and subjective measures.
The effect of amino groups attached at positions 2 and 8 of the hypoxanthine moiety in the structure, reactivity and stability of DNA duplexes and triplexes is studied by means of quantum mechanical calculations, as well as extended molecular dynamics (MD) and thermodynamic integration (MD/TI) simulations. Theoretical estimates of the change in stability related to 2'-deoxyguanosine (G) --> 2'-deoxyinosine (I) --> 8-amino-2'-deoxyinosine (8AI) mutations have been experimentally verified, after synthesis of the corresponding compounds. An amino group placed at position 2 stabilizes the duplex, as expected, and surprisingly also the triplex. The presence of an amino group at position 8 of the hypoxanthine moiety stabilizes the triplex but, surprisingly, destabilizes the duplex. The subtle electronic redistribution occurring upon the introduction of an amino group on the purine seems to be responsible for this surprising behavior. Interesting 'universal base' properties are found for 8AI.
Duocarmycins are a potent class of antitumor agents, whose activity arises through their covalent binding to adenine nucleobases of DNA.(1-3) Here, we perform molecular dynamics (MD) and hybrid Car-Parinello QM/MM simulations to investigate aspects of duocarmycin binding to the d(pGpApCpTpApApTpTpGpApC) oligonucleotide. We focus on the derivatives (+)-duocarmycin SA (DSA) and (+)-duocarmycin SI (DSI), for which structural information of the covalent complex with the oligonucleotide is available, as well as on the related, but less reactive, NBOC-duocarmycin SA (NBOC-DSA), interacting with the same oligonucleotide. Comparison is made with adenine alkylation reaction in water performed by the smallest of these compounds (NBOC-DSA). The MD calculations suggest that, in noncovalent complexes, (i) drug binding causes a partial dehydration of the minor groove, without inducing a significant conformational changes, and (ii) DSA and DSI occupy a more favorable position for nucleophilic attack than NBOC-DSA, consistently with the lower reactivity of the latter. The QM/MM calculations, which are used to investigate the first step of the alkylation reaction, turn out to provide strongly underestimated free energy barriers. Within these approximations, our calculations suggest that an important ingredient for the experimentally observed DNA catalytic power is the polarization of the drug by the biomolecular scaffold.
Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.
Vibrational spectra of a 12-residue beta-hairpin peptide, RYVEVBGKKILQ (HBG), stabilized by an Aib-Gly turn sequence (B = Aib) were investigated theoretically using a combination of molecular dynamics (MD) and density functional theory (DFT) calculations. Selected conformations of HBG were extracted from a classical MD trajectory and used for spectral simulations. DFT calculations, based on the Cartesian coordinate spectral property transfer protocol, were carried out for peptide structures in which all residues are replaced with Ala, except for the Aib and Gly residues, but the backbone (phi, psi, omega) structure of the original configuration is retained. The simulations provide a basis for interpretation of the HBG amide I infrared spectra in terms of structural variables such as detailed secondary structure and thermal conformational fluctuation as well as vibrational coupling as indicated by spectra of 13C isotope-labeled variants. The characteristic amide I band shape of such small beta-hairpin peptides appears to arise from the structure of the short antiparallel beta-sheet strands. The role of structural parameter fluctuation in vibrational coupling is evaluated by comparison of DFT-derived amide coupling constants for selected configurations and from transition dipole coupling calculations of coupling parameters between (13)C isotopically labeled residues for a MD-derived ensemble of configurations. Calculated results were compared with the experimentally obtained spectra for several (13)C isotope-labeled peptides of this sequence.
The purpose of this study was to evaluate the relative performance and usefulness of single dose (SD) and multiple dose (MD) regimens for bioequivalence (BE) determination. Drugs such as indomethacin, procainamide, erythromycin, quinidine, nifedipine were tested for BE under SD and MD dose regimens. Drugs characterized by low accumulation indices (AI) showed virtually no change in the 90% confidence interval (CI) of AUC and CMAX upon multiple dosing. On the other hand, drugs with higher AI appeared to have smaller CI at steady-state. For example, the CI range of AUC and CMAX of quinidine (AI of 1.54) decreased from 26 to 12 and from 22 to 12, respectively, upon multiple dosing. A Monte Carlo simulation study of SD and MD bioequivalence trials was performed. The probability of failing the bioequivalence test was evaluated for several situations defined by different levels of variability and correlation in ka constants, presence or absence of inter- and/or intra-individual variability in clearance (CL) and volume of distribution (V), and different degrees of accumulation. All the possible combinations of these factors were tested with SD and MD study designs. All simulations used 1000 data sets with 30 subjects in each data set for a total of 144 unique designs (total of 144,000 simulations of bioequivalence trials). Upon multiple dosing, narrowing of CI ranges was observed for drugs simulated to have high AI high variability and a large difference in absorption constants (ka) between test and reference formulations.(ABSTRACT TRUNCATED AT 250 WORDS)