55 Questions 110 Answers 0 Followers
Questions related from Anthony Nash
Hi all, I am hoping for some assistance. I (the researcher) would like to know whether it is possible and how to go about building a Jupyter Notebook that a collaborator could load in their web...
27 January 2021 1,221 1 View
I've generated a set of maps for my receptor using ADT (AutoDockTools). It has resulted in the files: .A.map .C.map .HD.map .N.map .NA.map .OA.map .SA.map .d.map .e.map But when I try loading the...
21 December 2020 7,919 2 View
Hi all, There are two big questions here. I am extremely grateful for even the smallest contribution that yields a result. This is concerning Raccoon(1) original, not Raccoon2. The original...
18 December 2020 4,522 2 View
Hi all, Please imagine this to be a brand new machine with the latest installation of: Autodock Vina Autodock AutodockTools (MGLTools). I'm following a Nature Protocol for the prediction of...
15 December 2020 2,430 3 View
Hi all, I guess there are probably several ways of doing this. I'm looking to get access to every Pubchem CID. All compounds are published for download, but the files are quite large in size and...
02 November 2020 8,845 3 View
Hi, I'm up against a horrible deadline so I'm hoping for some help figuring out what I'm doing wrong. I have a list of custom atom names in Gromacs of which I want to perform hydrogen bond...
11 January 2020 7,598 6 View
I'm waiting on email conformation for the GMX user group, but this is in reply to Reviewers so I'm pressed for time and asking here. I'm building a small collagen-based system. The crystal...
14 November 2019 5,786 10 View
I'm picking up quantum chemistry after a three year break. I've forgot a lot and I'm now using ORCA (version 4.0.1.2). My system is a non-enzymatic chemical reaction between lysine, arginine...
31 August 2019 4,112 2 View
Hello all, I have a huge dataset of longitudinal medical records. Each record contains indications for drug and disease. I have been able to calculate the Cure Rate in the population given a...
05 October 2018 2,597 3 View
I have a longitudinal retrospective data set of human medical records. They feature CONDITION and DRUG. There is no way of saying why a drug was prescribed other than observing the...
23 September 2018 4,877 3 View
I have a medical longitudinal retrospective dataset, records between the observation period of 2000 and end 2016. For many reasons not every medical record spans that entire time-frame, e.g. the...
31 May 2018 6,182 7 View
I am performing a literature survey and I am looking for as many computational chemistry studies of Advanced Glycation End-products (AGE) as possible. I am interested in any and all AGEs...
14 April 2018 9,225 2 View
I would like some thoughts on a way of quantifying a simple measure of lipid chain order parameter from a gromacs file. The g_order will give me an average and an error term for each of the...
04 February 2018 5,959 8 View
I would like to try, for my first time, modelling the interactions between protein and DNA in an aqueous environment using Molecular Dynamics. I have the latest Gromacs distribution available....
27 January 2018 8,493 2 View
I have two all atom molecule dynamics systems of a protein. One is a wildtype and the other is a mutant. They differ by amino acid sequences. I am running my simulations in Gromacs. Each system...
27 December 2017 4,163 4 View
A rather basic question just to make sure I understand some of the theory. I have an amino acid backbone, with the amine and carboxyl groups positively and negatively charged, respectively...
18 May 2017 8,190 14 View
If you perform a restart on any Gaussian job it will write over the contents of the earlier .log/.out file. Typically, an incomplete log./.out file can be loaded in GView showing the results thus...
17 February 2017 4,063 3 View
I have derived force constants for MD simulations from vibrational frequency calculations of QM Gaussian calculations. The QM to MD force constant conversion is from my own software. Now that...
22 January 2017 6,230 5 View
A simple question, perhaps badly phrased in the title, what is the largest singular biomolecule i.e., not stacks of carbohydrates one after the other, to be involved in post-translation modified...
20 January 2017 3,548 3 View
I always think of a heme cofactor with bound iron as a flat molecule. However, after reparameterizing heme without iron for the purpose of an MD simulation, the cofactor demonstrates quite a...
30 September 2016 5,691 3 View
Out of interest, I would like to know who has reported and for what reason a set of quantitative results using the total energy of an MD study? That is, the total energy reported at the end of the...
29 July 2016 1,554 6 View
What is the latest in free/open software to predict small molecule/ligand binding to proteins? I am running a Mac for GUI and a cluster of unix machines for computational calculations. A very...
27 July 2016 1,276 6 View
I want a good measure for water content around a wild type side chain, and then a measure around the glycated form of that side chain i.e., a comparison of two individual systems. I think that...
23 July 2016 449 6 View
I have two setups for a cluster of rod shaped proteins as a surface with water on the outside (it is a fibril of collagen protein, periodically stacked in the x-y dimension and a long water box in...
29 May 2016 5,794 4 View
I am trying to remove an -OH group from a ring structure using DFT/ab initio after having closed a carbon ring structure (Re->TS1->Int1, below). -OH is extremely difficult to remove as it...
04 May 2016 7,666 13 View
I have a very long pathway to map using DFT. There are approximately 15 intermediates. I am using Gaussian to capture the reactants, intermediates, transitions states and the products. I have so...
28 January 2016 8,952 3 View
Hi all, For a while now I have been building a small application to generate AMBER forcefield compatible force constants (stretch and angle). I use the full Force Constant matrix for a Gaussian...
22 November 2015 7,304 15 View
I have an alpha-helix embedded in a bilayer. I'm using gromacs 5+ and plumed 2.2. I want to encourage the protein to sample a 'rotation', so I can then put it up against another helix to find the...
20 November 2015 4,660 4 View
I'm using plumed 2.2 with Gromacs 5. I am currently struggling to cover a sufficient range of my three collective variables, two torsional (dihedral) and one distance. My system essentially...
15 November 2015 3,082 6 View
This hadn't really occurred to me until an experimentalist with no experience in computational chemistry asked me "to how many significant figures is Gaussian accurate when measuring distance?" I...
29 October 2015 8,140 5 View
I am reconstructing the force constants from the weighted hessian matrix, as calculated using HF/6-311+G(2d,2p) in Gaussian G09. The mathematics simple used a submatrix of the atoms of interest...
26 October 2015 891 7 View
Hi all, I’m a little stuck on deriving an angular force constant for the amber forcefield. Consider in amber, that the potential energy from the angle ABC is: U_{ABC} = K(theta - theta_eq)^2 The...
06 October 2015 1,191 3 View
Hi all, I am trying to work out the angular force constant from the frequency calculation of an ab initio Hessian for a mechanical force field (amber). The force constant of bond stretching is...
05 October 2015 997 3 View
How do I extract the whole Hessian from my chk/output/fchk Gaussian file? Note: The data in the .out file from a FREQ run is displacement and frequencies, NOT hessian data. The only think I can...
28 September 2015 5,014 18 View
I have a frcmod.unknown after parameterising a structure using R.E.D. for an AMBER simulation. Thankfully all of the atom types are known (organic, no metals), however, the bond connectivity is...
26 September 2015 9,805 11 View
Hi all, I am trying to provide the force constants and distance/angles in frcmod.user for a SECOND (follow-up) submission to R.E.D server so I have a complete amber force field. Unfortunately I...
24 September 2015 3,277 3 View
This seems like it should be pretty obvious, but I have missed it completely! What are the units for the height and width (sigma) of the gaussian bias potentials used during a METAD (metadynamics)...
05 August 2015 5,664 2 View
I have performed a 200 ns MD simulation (Gromacs - Amber forcefield) using PLUMED meta-dynamics with a distance based reaction coordinate of protein A entering protein B. I do see some noticeable...
01 August 2015 4,556 7 View
Hi all, Typically, restricted to electron structure of organic molecules and Molecular Dynamics, I've only just started to look into transition metal complexes. My question, to accurately model an...
16 April 2015 380 4 View
Hi all, I have a compound with a particular functional group surrounded by water molecules. I am using Gaussian G09 to optimise, first with HF/6-31+g** and then I will progressively move higher...
24 March 2015 648 12 View
Hi all, I'm having trouble coming up with a design for a 'stable' dihedral measurement along the principle orientation of a helical protein in a GROMACS topology. I will use this measurement as a...
04 March 2015 4,578 1 View
When you set up the initial geometry of a TS in Gaussian, should I be modelling the charge and spin multiplicity of the TS itself or either the product/reactant? For example, in the case of...
02 March 2015 8,602 4 View
I have been approached by a colleague who is having difficulty executing OPT->FREQ->Pop(NBOs) over a number of structures. So I performed the OPT and FREQ separate from the population...
25 February 2015 334 8 View
Dear RG experts, I'm playing with Meta dynamics for the first time and I'm struggling to come up with a way of rotating my ligand in the catalytic site of an enzyme. The ligand is a long narrow...
20 December 2014 4,922 3 View
Hi all, I have successfully derived partial charges using RED for two separate zinc binding domains in a protein. All the settings were carefully selected for the Amber sb99 FF and the structure...
28 November 2014 8,618 1 View
Hi all, I have done this countless time and I've "googled" this a lot. Unfortunately I have a simulation that simply won't behave. I am trying to visualise (and analyse) a trajectory in VMD after...
11 November 2014 6,433 19 View
Hi, I have a peptide fragment with a Cu^2+/Ag^+ accompanying ion. I want to optimise their geometry before I derive partial charges. I thought I would try a couple of different approaches. The...
26 October 2014 5,334 10 View
Hi all, I've decided to try some bench marking with respects to performance and given that I have plenty of memory to play with I thought I would start by preserving all of the integrals using...
19 September 2014 1,684 1 View
I have derived the partial charges for a zinc ion in a catalytic site of an enzyme using R.E.S.P for the Amberffsb force field. For two of the three transition metals (the third belong to the...
22 August 2014 4,529 5 View
I have successfully modelled the reaction pathway for the nucleophilic addition of an electronegative nitrogen attacking a slightly electropositive carbon. What was a double bond of an oxygen to...
06 July 2014 1,398 2 View
I am able to run both ab initio calculations and MD simulation. However, I am relatively new to parameterising non-bonded interactions (transition metal in an active site supported by non-bonded...
17 May 2014 6,489 2 View
I have an enzyme and a protein ligand bound. In the active site of the enzyme is a ZN2+ ion. I am using the latest version of gromacs with the Amber99SBFF. After carefully settling my solvated...
17 May 2014 2,194 3 View
I'm fairly new to the field of cheminformatics so I would be very grateful for experienced suggestions. I have three question (at the bottom), followed by two further questions that I thought of...
01 January 1970 286 3 View
I'm having a go at calculating the free energy of solvation of two small proteins, as separate systems. I am using the Bennett Acceptance Ratio approach by coupling the coloumbic and van der Waals...
01 January 1970 2,706 9 View
Dear all, I'm hoping for a little bit of guidance. I'm using RDKit and sklearn to implement a random forest model trained using our drug screening (1600 drugs) results to predict new drug...
01 January 1970 9,130 2 View