- Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] Molecular Dynamics
by Jean-Didier Pierre Marechal 18 Sep '20

18 Sep '20

Dear Edmund, Elaine is completely right about everything. This interface was first developed for our lab for our own needs (training, short MD and Teaching) then we work with the Chimera’s team to make it look the way it looks today (a lot better than the initial shot). It can do a lot of things but in its actual state it remains a simple tool for first insights on the system and not necessarily convenient for all systems (especially if you look for large systems with long MDs). Typical time scale of MDs in today standards are difficult to performe with this interface. The main application we make of it is on small molecules and for teaching purposes it is fantastic. You could eventually look look for short MDs on a complex protein-ligand but I would reckon to go for standard MD drivers. We did not continue on this project because ChimeraX, that already integrates openMM, was about to be relased. Hope this help, JeanDI ***************************** Prof. Dr. Jean-Didier Maréchal (JeanDi) Associate Professor Insilichem Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es personal webpage: http://gent.uab.cat/jdidier insilichem webpage: http://www.insilichem.com

1 0

Molecular Dynamics
by Edmund Marinelli 18 Sep '20

18 Sep '20

Hello: Very new to Chimera and to modeling in general. 1.) Does the Molecular Dynamics Tool allow simulation of a protein ligand complex? 2.) Does a protein ligand docking pose generated by Autodock, for example, qualify as a 'model' that can be read by the MD Tool and be employed in a MD run? Is there a comprehensive tutorial or manual discussing in detail the setting of MD run within Chimera? Any help would be appreciated. Thanks. Edmund R Marinelli, PhD

2 1

measurement of arc-length (geodesic distance) between atoms on a protein surface?
by Wiener, Michael C (mcw2s) 17 Sep '20

17 Sep '20

New subscriber here, with a question about whether Chimera currently can measure arc-length (geodesic distance) between atoms on a protein surface (or, for example, generate a Calpha-Calpha distance matrix of this). I am (very) interested in measuring/obtaining distances between pairs of Calpha carbons (or pairs of other atoms, I suppose) of different amino acids that are on the surface of a protein. While the difference between arc-length (geodesic) distance and cartesian/Euclidean distance will be negligible if, for example, the two solvent-exposed Calpha atoms are on the same solvent-exposed face of an alpha-helix, the difference could be significant for residues on opposite “sides” of the protein. I did see that there are ways in Chimera to get distances between an atom and a molecular surface, or to place markers on the surface & obtain the lengths between these for a distance, but I am VERY interested in whether one can simply select two atoms and get the arc-length (geodesic) distance or generate an arc-length (geodesic distance) matrix of all surface-exposed Calpha pairs. Several computational methods have been developed for this sort of thing, via graph theoretical or other approximation method, and have been published (Hall et al., Biophysical Chemistry 190-191 [2014] 50-55). However, I have not yet (?) found this to be in any graphics/analysis program. I’m thinking that this could perhaps be useful to someone besides me? Any assistance would be greatly appreciated. My lab and I are (much) more on the user side than the developer side, so are not currently well-positioned to code it in-lab. Regards, -MW (btw, if anyone in the “Chimera Group” wants to know more about why I am interested in this, don’t hesitate to send me an email. This type of distance calculation is related to some experimental technology/assay development that we are doing.) ----- Michael C. Wiener. Ph.D., Professor Department of Molecular Physiology and Biological Physics University of Virginia Charlottesville, VA 22908-0886 434-243-2731 (office); 434-243-2730 (lab)

3 2

generate helical symmetry mates around an axis
by sadeem ahmad 17 Sep '20

17 Sep '20

Hi all, I am trying to generate a helical trimer from the monomer pdb around a central double-stranded RNA. The rise and twist are 46 A and 74.0 deg respectively. I opened my monomer pdb in chimera *(Screenshot-0)* and used the following command to generate an axis around the central RNA double helix: *define axis #0:.X,.Y* where X and Y are the chain IDs for the 2 RNA strands This command helped me generate an axis named *a1* along the RNA duplex *(Screenshot-1)*. Now I want to generate a helical trimer along the RNA duplex with rise=46A and twist=74deg. For this I used the command: *sym #0 group H,46,74,3 axis a1* However, it is not placing the symmetric copies around axis a1 but around a distant axis *(Screenshot-2)*. I am not sure how to generate the helical trimer along the a1 axis. Any help in this regard would be appreciated. Thanks. -Sadeem -- :-)

3 5

Re: [Chimera-users] Chimera-users Digest, Vol 209, Issue 27
by Jatin Kashyap 17 Sep '20

17 Sep '20

Eric, Following are steps I followed: load the structure> del H > run Dock Prep(Add Hydrogen - ticked, Write Mol2 - unticked) > del :10.P@h :8.T@h > minimize (while deselecting the add hydrogen option) I got a warning that I am not sure about. I am pasting the log below[1], can you please confirm that I followed the correct steps and warning can be ignored in the log. [1] Model 0 (7bv2_nsp12-7.pdb) appears to be a protein without secondary structure assignments. Automatically computing assignments using 'ksdssp' and parameter values: energy cutoff -0.5 minimum helix length 3 minimum strand length 3 Use command 'help ksdssp' for more information. Model 0 (7bv2_nsp12-7.pdb) appears to be a protein without secondary structure assignments. Automatically computing assignments using 'ksdssp' and parameter values: energy cutoff -0.5 minimum helix length 3 minimum strand length 3 Use command 'help ksdssp' for more information. Residues with incomplete side chains: LYS 98.A ASP 100.A ASP 218.A ASP 824.A GLU 50.C Replacing each by 'swapaa same (residue atom spec) lib Dunbrack preserve True' LYS 98.A: phi -78.3, psi 128.2 trans Chi angles for LYS 98.A: ASP 100.A: phi -70.1, psi -6.8 trans Chi angles for ASP 100.A: ASP 218.A: phi -127.8, psi 120.9 trans Chi angles for ASP 218.A: ASP 824.A: phi -92.0, psi 1.6 trans Chi angles for ASP 824.A: GLU 50.C: phi -60.9, psi -45.8 trans Chi angles for GLU 50.C: Applying LYS rotamer (chi angles: -62.5 -75.3 99.0 -88.6) to LYS 98.A Applying ASP rotamer (chi angles: -67.2 -26.9) to ASP 824.A Applying ASP rotamer (chi angles: -70.6 -17.9) to ASP 100.A Applying GLU rotamer (chi angles: -170.8 -83.2 -29.7) to GLU 50.C Applying ASP rotamer (chi angles: -175.9 -11.3) to ASP 218.A No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain A; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain C; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain P; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain T; guessing terminii instead Chain-initial residues that are actual N terminii: VAL 31.A, LYS 2.C, G 10.P, U 8.T Chain-initial residues that are not actual N terminii: Chain-final residues that are actual C terminii: U 20.P, C 21.T Chain-final residues that are not actual C terminii: THR 929.A, GLY 64.C 856 hydrogen bonds Removing spurious proton from 'C' of THR 929.A Removing spurious proton from 'C' of GLY 64.C Hydrogens added Charge model: AMBER ff14SB Non-standard atom names: G H (G 10.P H) U H (U 8.T H) Total charge for #0: -30.384 The following residues had non-integral charges: U 20.P -0.6919 C 21.T -0.6919 Correct charges are unknown for 2 non-standard atom names in otherwise standard residues Charges of 0.0 were assigned to the unknown atoms 1 model(s) had non-integral total charge Details in reply log No incomplete side chains Charge model: AMBER ff14SB Total charge for #0: -30.384 The following residues had non-integral charges: U 20.P -0.6919 C 21.T -0.6919 1 model(s) had non-integral total charge Details in reply log Using main parameter file parm10.dat modified by heme-iron.frcmod, frcmod.ff14SB, frcmod.ionsjc_tip3p Initial energy: -45283.236153 kJ/mol Step 12 Potential energy: -58447.379716, Gradient norm: 758.031360 Updated 15035 atoms. RMSD: 0.037342 Finished 10 of 100 steepest descent minimization steps Step 12 Potential energy: -64213.179928, Gradient norm: 362.328469 Updated 15035 atoms. RMSD: 0.018132 Finished 20 of 100 steepest descent minimization steps Step 12 Potential energy: -65930.345155, Gradient norm: 417.614654 Updated 15035 atoms. RMSD: 0.012606 Finished 30 of 100 steepest descent minimization steps Step 12 Potential energy: -67901.748990, Gradient norm: 190.328174 Updated 15035 atoms. RMSD: 0.010435 Finished 40 of 100 steepest descent minimization steps Step 12 Potential energy: -68928.640793, Gradient norm: 79.670023 Updated 15035 atoms. RMSD: 0.007430 Finished 50 of 100 steepest descent minimization steps Step 12 Potential energy: -70504.551913, Gradient norm: 192.974949 Updated 15035 atoms. RMSD: 0.016225 Finished 60 of 100 steepest descent minimization steps Step 12 Potential energy: -71243.383589, Gradient norm: 95.907393 Updated 15035 atoms. RMSD: 0.006255 Finished 70 of 100 steepest descent minimization steps Step 12 Potential energy: -71661.708391, Gradient norm: 193.911209 Updated 15035 atoms. RMSD: 0.006283 Finished 80 of 100 steepest descent minimization steps Step 12 Potential energy: -72255.311758, Gradient norm: 199.037193 Updated 15035 atoms. RMSD: 0.007288 Finished 90 of 100 steepest descent minimization steps Step 12 Potential energy: -72873.927407, Gradient norm: 84.768251 Updated 15035 atoms. RMSD: 0.005517 Finished 100 of 100 steepest descent minimization steps Step 12 Potential energy: -79350.452189, Gradient norm: 206.369396 Updated 15035 atoms. RMSD: 0.114423 Finished 10 of 10 conjugate gradient minimization steps —— Jatin Kashyap > On Sep 16, 2020, at 6:41 PM, chimera-users-request(a)cgl.ucsf.edu wrote: > > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Re: Error "No MMTK name for atom "H" in standard residue ?G"" > (Eric Pettersen) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 16 Sep 2020 15:41:33 -0700 > From: Eric Pettersen <pett(a)cgl.ucsf.edu> > To: Jatin Kashyap <jk435(a)njit.edu> > Cc: Chimera User Help <chimera-users(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Error "No MMTK name for atom "H" in > standard residue ?G"" > Message-ID: <9B74E727-F918-4A68-9E48-83F05E6B329E(a)cgl.ucsf.edu> > Content-Type: text/plain; charset="utf-8" > > Ah I see. Because you also have missing side chains you have to do things a little differently. Namely run Dock Prep (Tools?Structure Editing?Dock Prep) first, keeping "Add hydrogens" on but "Write Mol2" off; then use the delete command to delete the two unwanted protons; then run Minimize. > > --Eric > >> On Sep 16, 2020, at 2:54 PM, Jatin Kashyap <jk435(a)njit.edu> wrote: >> >> I am sorry for the confusion. Actually I did that step by GUI options (Tools>Structure Editing>AddH). But I did it again by issuing ?addh? from command line along with other steps. >> But still I am getting the same error. >> >> ?? >> Jatin Kashyap >> >> >>> On Sep 16, 2020, at 5:49 PM, Eric Pettersen <pett(a)cgl.ucsf.edu <mailto:pett@cgl.ucsf.edu>> wrote: >>> >>> It's "addh" not "add H". >>> >>> --Eric >>> >>>> On Sep 16, 2020, at 2:44 PM, Jatin Kashyap <jk435(a)njit.edu <mailto:jk435@njit.edu>> wrote: >>>> >>>> Hi Eric, >>>> >>>> Just for the clarification I am double checking it here, first I loaded the structure > del H > add H > del :10.P@h :8.T@h > minimize (while deselecting the add hydrogen option). >>>> >>>> But it is throwing the error [1]. >>>> >>>> >>>> [1] >>>> Model 0 (7bv2_nsp12-7.pdb) appears to be a protein without secondary structure assignments. >>>> Automatically computing assignments using 'ksdssp' and parameter values: >>>> energy cutoff -0.5 >>>> minimum helix length 3 >>>> minimum strand length 3 >>>> Use command 'help ksdssp' for more information. >>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain A; guessing terminii instead >>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain C; guessing terminii instead >>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain P; guessing terminii instead >>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain T; guessing terminii instead >>>> Chain-initial residues that are actual N terminii: VAL 31.A, LYS 2.C, G 10.P, U 8.T >>>> Chain-initial residues that are not actual N terminii: >>>> Chain-final residues that are actual C terminii: U 20.P, C 21.T >>>> Chain-final residues that are not actual C terminii: THR 929.A, GLY 64.C >>>> 852 hydrogen bonds >>>> Removing spurious proton from 'C' of THR 929.A >>>> Removing spurious proton from 'C' of GLY 64.C >>>> Hydrogens added >>>> Residues with incomplete side chains: >>>> LYS 98.A >>>> ASP 100.A >>>> ASP 218.A >>>> ASP 824.A >>>> GLU 50.C >>>> Replacing each by 'swapaa same (residue atom spec) lib Dunbrack preserve True' >>>> LYS 98.A: phi -78.3, psi 128.2 trans >>>> Chi angles for LYS 98.A: >>>> ASP 100.A: phi -70.1, psi -6.8 trans >>>> Chi angles for ASP 100.A: >>>> ASP 218.A: phi -127.8, psi 120.9 trans >>>> Chi angles for ASP 218.A: >>>> ASP 824.A: phi -92.0, psi 1.6 trans >>>> Chi angles for ASP 824.A: >>>> GLU 50.C: phi -60.9, psi -45.8 trans >>>> Chi angles for GLU 50.C: >>>> Applying ASP rotamer (chi angles: -67.2 -26.9) to ASP 824.A >>>> Applying LYS rotamer (chi angles: -178.9 71.3 178.8 68.6) to LYS 98.A >>>> Applying ASP rotamer (chi angles: -175.9 -11.3) to ASP 218.A >>>> Applying GLU rotamer (chi angles: -170.8 -83.2 25.9) to GLU 50.C >>>> Applying ASP rotamer (chi angles: -70.6 -17.9) to ASP 100.A >>>> Charge model: AMBER ff14SB >>>> Total charge for #0: -32.172 >>>> The following residues had non-integral charges: >>>> LYS 98.A -0.6068 >>>> ASP 100.A -1.0636 >>>> ASP 218.A -1.0636 >>>> ASP 824.A -1.0636 >>>> GLU 50.C -0.9909 >>>> U 20.P -0.6919 >>>> C 21.T -0.6919 >>>> 1 model(s) had non-integral total charge >>>> Details in reply log >>>> >>>> Residue #0:98.A (LYS/lysine) is missing atoms H_gamma_3, H_epsilon_2, H_zeta_3, H_delta_2, H_epsilon_3, H_alpha, H_beta_3, H_gamma_2, H_delta_3, H_beta_2, H_zeta_1 and H_zeta_2 >>>> >>>> >>>> ?? >>>> Jatin Kashyap >>>> >>>> >>>>> On Sep 16, 2020, at 5:32 PM, Eric Pettersen <pett(a)cgl.ucsf.edu <mailto:pett@cgl.ucsf.edu>> wrote: >>>>> >>>>> Hi Jatin, >>>>> I you don't delete those phosphates, then Chimera will add spurious protons to them, which is what is causing the problem. In order to get this to work, you will have to run Add Hydrogens "by hand" before your minimization, delete the spurious protons, and then run your minimization -- being careful to skip adding hydrogens as part of the minimization process. >>>>> The command (Favorites?Command Line) to delete the unwanted protons in your case is "del :10.P@h :8.T@h". >>>>> >>>>> --Eric >>>>> >>>>> Eric Pettersen >>>>> UCSF Computer Graphics Lab >>>>> >>>>> >>>>>> On Sep 16, 2020, at 2:10 PM, Jatin Kashyap <jk435(a)njit.edu <mailto:jk435@njit.edu>> wrote: >>>>>> >>>>>> Dear Chimera Community Members, >>>>>> >>>>>> I am facing an error while minimizing a protein that have 5? terminal phosphates attached to it. >>>>>> Chimera is showing the warning that AMBER is unable to parametrize the phosphates and if you do not delete it, 0 charges will be assigned to it. >>>>>> I selected the option not to delete it but then it throws error[1] on the following step. I am pasting my log below [2]. >>>>>> >>>>>> Thank you for helping me to understand this. >>>>>> >>>>>> [1] >>>>>> No MMTK name for atom "H" in standard residue ?G" >>>>>> >>>>>> [2] >>>>>> Model 0 (7bv2_nsp12-7.pdb) appears to be a protein without secondary structure assignments. >>>>>> Automatically computing assignments using 'ksdssp' and parameter values: >>>>>> energy cutoff -0.5 >>>>>> minimum helix length 3 >>>>>> minimum strand length 3 >>>>>> Use command 'help ksdssp' for more information. >>>>>> Residues with incomplete side chains: >>>>>> LYS 98.A >>>>>> ASP 100.A >>>>>> ASP 218.A >>>>>> ASP 824.A >>>>>> GLU 50.C >>>>>> Replacing each by 'swapaa same (residue atom spec) lib Dunbrack preserve True' >>>>>> LYS 98.A: phi -78.3, psi 128.2 trans >>>>>> Chi angles for LYS 98.A: >>>>>> ASP 100.A: phi -70.1, psi -6.8 trans >>>>>> Chi angles for ASP 100.A: >>>>>> ASP 218.A: phi -127.8, psi 120.9 trans >>>>>> Chi angles for ASP 218.A: >>>>>> ASP 824.A: phi -92.0, psi 1.6 trans >>>>>> Chi angles for ASP 824.A: >>>>>> GLU 50.C: phi -60.9, psi -45.8 trans >>>>>> Chi angles for GLU 50.C: >>>>>> Applying LYS rotamer (chi angles: -62.5 -75.3 99.0 -88.6) to LYS 98.A >>>>>> Applying ASP rotamer (chi angles: -67.2 -26.9) to ASP 824.A >>>>>> Applying ASP rotamer (chi angles: -175.9 -11.3) to ASP 218.A >>>>>> Applying ASP rotamer (chi angles: -70.6 -17.9) to ASP 100.A >>>>>> Applying GLU rotamer (chi angles: -170.8 -83.2 -29.7) to GLU 50.C >>>>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain A; guessing terminii instead >>>>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain C; guessing terminii instead >>>>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain P; guessing terminii instead >>>>>> No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain T; guessing terminii instead >>>>>> Chain-initial residues that are actual N terminii: VAL 31.A, LYS 2.C, G 10.P, U 8.T >>>>>> Chain-initial residues that are not actual N terminii: >>>>>> Chain-final residues that are actual C terminii: U 20.P, C 21.T >>>>>> Chain-final residues that are not actual C terminii: THR 929.A, GLY 64.C >>>>>> 856 hydrogen bonds >>>>>> Removing spurious proton from 'C' of THR 929.A >>>>>> Removing spurious proton from 'C' of GLY 64.C >>>>>> Hydrogens added >>>>>> Charge model: AMBER ff14SB >>>>>> Non-standard atom names: >>>>>> G H (G 10.P H) >>>>>> U H (U 8.T H) >>>>>> Total charge for #0: -30.384 >>>>>> The following residues had non-integral charges: >>>>>> U 20.P -0.6919 >>>>>> C 21.T -0.6919 >>>>>> Correct charges are unknown for 2 non-standard atom names in otherwise standard residues >>>>>> >>>>>> Charges of 0.0 were assigned to the unknown atoms >>>>>> >>>>>> 1 model(s) had non-integral total charge >>>>>> Details in reply log >>>>>> >>>>>> No MMTK name for atom "H" in standard residue "G" >>>>>> >>>>>> >>>>>> ?? >>>>>> Jatin Kashyap >>>>>> New Jersey Institute of Technology >>>>>> >>>>>> _______________________________________________ >>>>>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>>>>> Manage subscription: https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >>>>> >>>> >>>> _______________________________________________ >>>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>>> Manage subscription: https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >>> >> > >

2 1

Error "No MMTK name for atom "H" in standard residue “G""
by Jatin Kashyap 16 Sep '20

16 Sep '20

Dear Chimera Community Members, I am facing an error while minimizing a protein that have 5’ terminal phosphates attached to it. Chimera is showing the warning that AMBER is unable to parametrize the phosphates and if you do not delete it, 0 charges will be assigned to it. I selected the option not to delete it but then it throws error[1] on the following step. I am pasting my log below [2]. Thank you for helping me to understand this. [1] No MMTK name for atom "H" in standard residue “G" [2] Model 0 (7bv2_nsp12-7.pdb) appears to be a protein without secondary structure assignments. Automatically computing assignments using 'ksdssp' and parameter values: energy cutoff -0.5 minimum helix length 3 minimum strand length 3 Use command 'help ksdssp' for more information. Residues with incomplete side chains: LYS 98.A ASP 100.A ASP 218.A ASP 824.A GLU 50.C Replacing each by 'swapaa same (residue atom spec) lib Dunbrack preserve True' LYS 98.A: phi -78.3, psi 128.2 trans Chi angles for LYS 98.A: ASP 100.A: phi -70.1, psi -6.8 trans Chi angles for ASP 100.A: ASP 218.A: phi -127.8, psi 120.9 trans Chi angles for ASP 218.A: ASP 824.A: phi -92.0, psi 1.6 trans Chi angles for ASP 824.A: GLU 50.C: phi -60.9, psi -45.8 trans Chi angles for GLU 50.C: Applying LYS rotamer (chi angles: -62.5 -75.3 99.0 -88.6) to LYS 98.A Applying ASP rotamer (chi angles: -67.2 -26.9) to ASP 824.A Applying ASP rotamer (chi angles: -175.9 -11.3) to ASP 218.A Applying ASP rotamer (chi angles: -70.6 -17.9) to ASP 100.A Applying GLU rotamer (chi angles: -170.8 -83.2 -29.7) to GLU 50.C No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain A; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain C; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain P; guessing terminii instead No SEQRES records for 7bv2_nsp12-7.pdb (#0) chain T; guessing terminii instead Chain-initial residues that are actual N terminii: VAL 31.A, LYS 2.C, G 10.P, U 8.T Chain-initial residues that are not actual N terminii: Chain-final residues that are actual C terminii: U 20.P, C 21.T Chain-final residues that are not actual C terminii: THR 929.A, GLY 64.C 856 hydrogen bonds Removing spurious proton from 'C' of THR 929.A Removing spurious proton from 'C' of GLY 64.C Hydrogens added Charge model: AMBER ff14SB Non-standard atom names: G H (G 10.P H) U H (U 8.T H) Total charge for #0: -30.384 The following residues had non-integral charges: U 20.P -0.6919 C 21.T -0.6919 Correct charges are unknown for 2 non-standard atom names in otherwise standard residues Charges of 0.0 were assigned to the unknown atoms 1 model(s) had non-integral total charge Details in reply log No MMTK name for atom "H" in standard residue "G" —— Jatin Kashyap New Jersey Institute of Technology

2 5

Re: [Chimera-users] Euler angles and rotations
by Elaine Meng 16 Sep '20

16 Sep '20

Hi Stephen, In Chimera, you can write the current transformations of models to a file with "matrixget" and apply transformations from those matrices by reading the file with "matrixset." <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/matrixset.html> It doesn't have anything to do with Euler angles, though. Not sure if it suffices for what you want to do. To send a question for all staff to see, it is best to use the mailing list CC'd here. For ChimeraX, there is similar stuff in the "view matrix" command and saved "positions" files. However, the positions file doesn't include camera position, so in my opinion is not quite as useful for regenerating the view as matrixget(set) in Chimera. <http://rbvi.ucsf.edu/chimerax/docs/user/commands/view.html#matrix> <http://rbvi.ucsf.edu/chimerax/docs/user/formats/positions.html> Best, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco > On Sep 16, 2020, at 12:04 PM, Chin-Bow, Stephen <stephen.chin-bow(a)yale.edu> wrote: > > Hello Elaine, > > I have a math question which hopefully you or the CGL staff can answer quickly. > > I am using the original Chimera (although I recently upgraded my GPU and I am curious to try ChimerX ). > > I am curious about Euler angles and rotations. I would like to learn more about them, but I suspect it has been too many years since I thought about matrices. > > I want another person to first download a published PDB from the RCSB before rotating it using a command line script. > > The final perspective I have was achieved after ten rotations. > > Is there a way I can apply Euler angles and matrix multiplication to summarize the final position using the fewest number of rotations? > > Yes, I can simply share a session file, but my future goal is more general. I want to encourage the other person to learn to use scripts, because the command line opens up the power of Chimera. > > Thanks a lot for the remedial math help! > > Stephen

1 0

Zone in pixel or Angstrom?
by Dieter Blaas 14 Sep '20

14 Sep '20

Hi, is the zone around selected atoms cut out with "Volume Viewer > Features > Zone" measured in Angstroms or pixels? Thank you, regards, Dieter ------------------------------------------------------------------------ Dieter Blaas, Max Perutz Laboratories Medical University of Vienna, Inst. Med. Biochem., Vienna Biocenter (VBC), Dr. Bohr Gasse 9/3, A-1030 Vienna, Austria, Tel: 0043 1 4277 61630, Mobile: 0043 699 1942 1659 e-mail: dieter.blaas(a)meduniwien.ac.at ------------------------------------------------------------------------

2 2

issue with heme minimization
by Rio Aquarius 14 Sep '20

14 Sep '20

Hi, I have an issue with heme minimization, Fe atom leaves a heme plane. I saw at least one user had the same issue and got help to solve it ( https://www.cgl.ucsf.edu/pipermail/chimera-users/2019-July/015972.html) I could try the protocol described there, but I don't have the heme.defattr file. I hope you can help me with it too. I have attached a pdb file with this issue.

2 6

Re: [Chimera-users] Is it possible to make a dimer from a monomer protein in chimera?
by Elaine Meng 14 Sep '20

14 Sep '20

On Sep 14, 2020, at 12:07 PM, Tarazkar, Maryam <Maryam.Tarazkar(a)ucsf.edu> wrote: > > Hi Elaine, > I am wondering can you refer me to a tutorial in Chimera where it explains how to make a dimer protein from a monomer? > Thank you, > Maryam Hi Maryam, Yes, of course you can do it, but there is no specific tutorial: how to do it really depends what information is in the input file, what you know about the structure, etc. If it is an entry in the PDB database you could just try fetching the biounit directly and seeing if it is a dimer. E.g. menu File... Fetch by ID, use the "PDB (biounit)" choice. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/fetch.html> Some files have symmetry information in their files so that you can build the multimer with the "sym" command, e.g. open 1fav sym #0 #0 <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/sym.html> If it is not yet a deposited structure, but you already have a related structure that is a dimer, you can just open two copies of your protein and match it onto the two monomers of that other structure, e.g. with the Matchmaker tool. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/superposition.html> I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco

3 10

PDB saving
by Henry Shook 14 Sep '20

14 Sep '20

Chimera Developers, I am working on a research project modeling and experimenting on biomolecules and I have just started using Chimera but whenever I am trying to save a PDB, it says that my permission is denied. Is there a way that this can be fixed? I am currently using a school computer and there are certain restrictions may be put in place so that might be the problem as well. Here is a video of what happens if it may help you: Thanks! Henry Shook -- ***** Please note: The Sender of this email is a Student of Pine Lake Preparatory and is not an employee or agent of Pine Lake Preparatory ***** -- ***************************************************************************** CONFIDENTIALITY NOTICE: This electronic mail message, including any attachment(s), contains information which may be confidential pursuant to the Family Educational Rights & Privacy Act (FERPA), 20 U.S.C. 1232g, and/or N.C.G.S. 115C-319, and/or proprietary non-public information. If you are not the intended recipient or the authorized assistant to the intended recipient, please notify the sender by replying to this message and deleting the e-mail from your system. Use, dissemination, distribution and/or reproduction of this message and/or its attachment(s), if any, by unauthorized recipients is not authorized and may be unlawful. *****************************************************************************

2 2

Question regarding saving protein-ligand and docking in Chimera
by Tarazkar, Maryam 11 Sep '20

11 Sep '20

Dear Sir or Madam, I did molecular docking to find the protein-ligand interactions using Chimera. Now I am trying to save the results of docking ( I took the screenshot here) in a text file and save the structure of protein-ligand in pdb file so I do the analysis using Pymol. When I click on file--> save pdb a page appears where I have to select which models to be saved. ( I took the screenshot of the page here) but when I select ligand, protein and docked options all tigetehr, the program gives me error in saving the files. I appreciate if you help me and let me know how I can resolve it and save the protein-ligand structure in best binding site in pdb so I do further analysis on this. Best Maryam [cid:42d36ba9-f8f9-4ae7-9884-45480afe5dc2] [cid:6467c35e-bbaa-44c4-9a56-8b329bc8f3c1]

2 8

installation issues
by Breves, Alexandra 10 Sep '20

10 Sep '20

Hello, I am having an issue installing chimera on my windows 10 (64-bit) computer. It is saying that I do not have 64x but I do. Can someone please assist me? Thank you, Alexandra Breves

1 0

Saving relaxed system in pdbqt format for Vina
by Jatin Kashyap 09 Sep '20

09 Sep '20

Dear Chimera Community Members, I am new to this code, and tried to navigate through this wonderful piece of code by myself. However, I am not able to find a way to save the relaxed system(relaxation performed in UCSF Chimera) in pdbqt format. Can anybody please help me to understand how I can do it, OR am I doing something completely wrong? Thank you for your time, and efforts. —— Jatin Kashyap Department of Mechanical and Industrial Engineering New Jersey Institute of Technology (NJIT) University Heights Newark, NJ 07102-1982 Phone- (201)889-5783 Email- jk435(a)njit.edu

2 4

Saving marker file on command line
by Rene Frank 08 Sep '20

08 Sep '20

Hi, I would like to save a .cmm marker file on the command line by specifying the model number. Please can you help me with what is the command to do this. Thank you, Rene

2 1

how to set initial threshold when opening a mrc
by Shuo 04 Sep '20

04 Sep '20

Hi, When I open an .mrc map, Chimera will show the map with an initial threshold. But the initial threshold is not so good, so I changed it to high value with the scale bar in Volume Viewer. I want to save the map, so that when I open the map again, the initial threshold will be the high value, not the lower one. How could I do this? Best, Shuo

2 1

Re: [Chimera-users] Chimera autodock vina for ligand docking on RNA
by Reeki Emrizal 03 Sep '20

03 Sep '20

Hi Elaine, Sorry about the subject. Just realised I forgot to add one after I hit send. I’ve tried autodock vina (not through chimera) and it works. However, I would like to automate both the receptor and ligand preparation steps and docking process and it seems to me that chimera is capable of it. Is there a help page that contains the command to turn off those options (Ignore chains of non-standard residues and Ignore all non-standard residues). The help page you provide doesn’t seem to provide a command that I can used in a docking script. Thanks for your help On Thursday, September 3, 2020, Elaine Meng <meng(a)cgl.ucsf.edu> wrote: > Hi Reeki, > I changed the subject line to something useful. > > The "Autodock Vina" tool in Chimera is just an interface that calls > Autodock Vina as a separate program. You should check the "receptor > options" in this tool: there are two options to ignore residues that are > not standard amino acids. I assume you should change those options to > "false". The options are explained in the Help page, which you can see by > clicking the Help button, or view the copy on our website here: > <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/vina/vina.html> > > However, even if you change these options to false and the docking > calculation finishes, I do not know whether Autodock Vina is expected to > give good results for docking to RNA. Also, the tool in Chimera does not > allow very much sampling with Autodock Vina; for more thorough sampling of > space, you may want to use Autodock Vina directly (not from Chimera). > > The Chimera developers are not the developers of the Autodock Vina > program, so I am not an expert on this topic. Maybe you would find more > information about whether Autodock Vina is meant to work for RNA receptors > at the Autodock Vina website: > <http://vina.scripps.edu/> > > I hope this helps, > Elaine > ----- > Elaine C. Meng, Ph.D. > UCSF Chimera(X) team > Department of Pharmaceutical Chemistry > University of California, San Francisco > > > On Sep 2, 2020, at 5:02 PM, Reeki Emrizal <reeki(a)mfrlab.org> wrote: > > > > Dear Chimera team, > > > > I'm attempting to dock small molecule ligand on a RNA molecule using > chimera autodock vina. However, the docking process ended prematurely after > the receptor.pdb (the RNA molecule) file is made. Please find attached a > screenshot of the reply log. I'm suspecting the error comes from the > receptor.pdb file made by chimera. Please find attached the screenshots of > receptor.pdb files for a RNA molecule and a protein molecule. The > receptor.pdb for protein has the column for residue names parsed correctly > (e.g. "THR") but not for the receptor.pdb for RNA (e.g. "G" which I > believed should be " G"). > > > > My questions are: (1) Is chimera allows docking of small molecule ligand > on RNA molecule and (2) Is the a parameter that I've overlooked to allow > this? > > > > Kind Regards > > <receptor.pdb protein.png><reply log.png><receptor.pdb rna.png> >

2 1

flipping model in chimera
by Dmitry Semchonok 03 Sep '20

03 Sep '20

3 4

(no subject)
by Reeki Emrizal 02 Sep '20

02 Sep '20

Dear Chimera team, I'm attempting to dock small molecule ligand on a RNA molecule using chimera autodock vina. However, the docking process ended prematurely after the receptor.pdb (the RNA molecule) file is made. Please find attached a screenshot of the reply log. I'm suspecting the error comes from the receptor.pdb file made by chimera. Please find attached the screenshots of receptor.pdb files for a RNA molecule and a protein molecule. The receptor.pdb for protein has the column for residue names parsed correctly (e.g. "THR") but not for the receptor.pdb for RNA (e.g. "G" which I believed should be " G"). My questions are: (1) Is chimera allows docking of small molecule ligand on RNA molecule and (2) Is the a parameter that I've overlooked to allow this? Kind Regards

2 1

Why is molmap changing the grid spacing and origin of simulated map
by Ouadou, Anes (MU-Student) 01 Sep '20

01 Sep '20

Hello, I am using molmap to simulate cryo-em density maps from protein pdb file. I follow the two steps outlined in Chimera tutorial (http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html#onGrid) I first create an empty map with specific size, grid spacing, and origin. vop new #1 size 116,96,152 gridspacing 1 origin -99.0,-39.0,-58.0 Then simulate the pdb file on the empty map using the option the grid spacing molmap #0 3 onGrid #1 the resulting simulated map has grid spacing changed from 1A to 2A and the origin coordinates are also changed to half the values I provided. This does not happen with all cases. The following commands work perfectly and generate a simulated map with grid spacing of 1A vop new #1 size 78,88,121 gridspacing 1 origin -32.73,-28.64,-63.42 molmap #0 3 onGrid #1 I can re-sample the map to obtain 1A grid spacing, but it adds extra work (I doing this to 100s of examples) and seems to affect the resolution of the simulated map. Could you please adive me to why is this phenomenon is happening and how to solve it. best regards Anes Ouadou CS PhD Student University of Missouri-Columbia aomqc(a)mail.missouri.edu Mobile: (573 )-289-9073

2 1

on mrc map splitting without losing resolution
by Smith Lee 01 Sep '20

01 Sep '20

Dear All, Suppose my protein complex pdb contains A,B,C and D subunits, and its mrc map is available. With Chimera, can we and how to get the mrc map for this protein complex with the part for subunit D cut out, without losing the resolution of the protein complex with subunit A, B and C? Smith

2 1

Chimera, Modeller, and socket modules
by Casey Becker 31 Aug '20

31 Aug '20

Hello, I'm attempting to do some loop refinements using Modeller in Chimera while making use of Modeller's parallel job feature to use multiple CPUs. When I run the scripts directly through python's IDLE interface, it works. However, when I attempt to run a custom script including parallel jobs directly through the Chimera interface, I receive an error stating that I'm missing a socket module. I've contacted the Modeller help email and they referred me to you after suggesting I replace the mod9.24.exe with python.exe in the "Modeller executable," but that didn't work. In fact, it froze my Chimera. I've attached copies of my scripts and the error I get. Any assistance appreciated. Thank you.

2 2

User-defined graphics in chimera: NMR tensor plotting
by Albert Smith-Penzel 31 Aug '20

31 Aug '20

Hello, I have a question about improving user defined graphics in Chimera. For background (maybe not necessary for answering): I’m working with nuclear magnetic resonance and MD simulation, and so I want to plot dynamically averaged spin-spin (H–C) dipole tensors, obtained from MD, on top of the corresponding H–C bonds. I’ve been relatively successful, but while one has a very high-resolution picture of the molecule itself, the tensors themselves are visibly triangulated. How I get this plot: Each tensor is expressed first in spherical coordinate (theta, phi, r). In python, I run Delauney from scipy.spatial on theta and phi to get a list of vertices that can be used for triangulation of the tensor. Then, for each tensor, I convert theta, phi, and r into cartesian coordinates, shift the center so that the tensor sits on top of the correct bond. Finally, I create a .bild file where I write out each triangle using the coordinates and vertices. Tensors have negative and positive values (red and blue in the plot), and so I only take triangles for which r for all vertices is positive, and then in a second step, take only triangles where all vertices are negative. This all gets written to the .bild file, and then opened in chimera on top of the molecule (see picture). The question: The result is certainly acceptable, however, there is clearly a quality difference in the molecule rendering vs. the tensor rendering (obviously, because the tensors are just made of triangles that I’ve created in python). Is there a better way to bring the surface data describing the tensors into chimera that would then allow usage of some sort of smoothing within chimera for a nicer image? (And possibly with a relatively simple data format….) Thanks for your help! Cheers, Albert Smith

4 5

Error installing
by David Butler 28 Aug '20

28 Aug '20

Hi I get the error attached when trying to install the stable version of the software. My OS is Catalina 10.15.6 Thanks Dave David Butler, PhD. VP, Head of Therapeutics Development at Flagship Labs 63 600/700 One Kendall Square, Cambridge, MA 02139 781-866-1936 Disclaimer The information contained in this communication from the sender is confidential. It is intended solely for use by the recipient and others authorized to receive it. If you are not the recipient, you should delete this message. You are hereby notified that any disclosure, copying, distribution or taking action in relation to the contents of this information is strictly prohibited and may be unlawful.

2 2

On installation of Chimera linux version
by Smith Lee 27 Aug '20

27 Aug '20

Dear All, Fromhttps://www.cgl.ucsf.edu/chimera/download.html I have downloaded thechimera-1.14-linux_x86_64.bin, and have tried to install it exactly based onthe accompanied exact "Instructions”. After the installation completed, Igot to the bin folder, and try to launch chimera with ./chimera. However her俄 it says the installed/UCSF-chimera64-1.14/lib/python2.7/lib-tk/Tkinter.py was not configured for mypython. Will you please let me know howto configure TK for my python, which was a downloaded component of the chimerapackage? Best wishes, Smith

2 13